bzoj 1912 : [Apio2010]patrol 巡逻 树的直径

#include <bits/stdc++.h>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, n, a) for(int i = a; i<n; i++)
#define fi first
#define se second
typedef complex <double> cmx;
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
const int maxn = 1e5+5;
int head[maxn], num, mx1[maxn], mx2[maxn], n, maxx, p;
struct node
{
int to, nextt, val;
}e[maxn*2];
void add(int u, int v, int val) {
e[num].to = v, e[num].nextt = head[u], e[num].val = val, head[u] = num++;
}
void init() {
num = 0;
mem1(head);
}
int dfs(int u, int fa)
{
int maxx1 = 0, maxx2 = 0;
for(int i = head[u]; ~i; i = e[i].nextt) {
int v = e[i].to;
if(v == fa)
continue;
int tmp = e[i].val + dfs(v, u);
if(tmp > maxx1) {
maxx2 = maxx1;
maxx1 = tmp;
mx2[u] = mx1[u];
mx1[u] = i;
} else if(tmp > maxx2) {
maxx2 = tmp;
mx2[u] = i;
}
}
if(maxx < maxx1 + maxx2) {
maxx = maxx1+maxx2;
p = u;
}
return maxx1;
}
int main()
{
int k, u, v;
cin>>n>>k;
init();
for(int i = 0; i < n - 1; i ++) {
scanf("%d%d", &u, &v);
add(u, v, 1);
add(v, u, 1);
}
dfs(1, 0);
int ans = 2*(n-1);
ans -= (maxx - 1);
if(k == 2) {
for(int i = mx1[p]; i; i = mx1[e[i].to]) e[i].val = e[i^1].val = -1;
for(int i = mx2[p]; i; i = mx1[e[i].to]) e[i].val = e[i^1].val = -1;
maxx = 0;
dfs(1, 0);
ans -= (maxx-1);
}
cout<<ans<<endl;
return 0;
}