首页 > 代码库 > bzoj 1912 : [Apio2010]patrol 巡逻 树的直径
bzoj 1912 : [Apio2010]patrol 巡逻 树的直径
题目链接
如果k==1, 显然就是直径。
k==2的时候, 把直径的边权变为-1, 然后在求一次直径。 变为-1是因为如果在走一次这条边, 答案会增加1.
学到了新的求直径的方法...
#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 secondtypedef 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;}
bzoj 1912 : [Apio2010]patrol 巡逻 树的直径
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。