其他分享
首页 > 其他分享> > P3629 [APIO2010]巡逻 题解

P3629 [APIO2010]巡逻 题解

作者:互联网

题面

\(k=1\) 的时候就是边数 \(\times 2\) 再减去直径长度。考虑 \(k=2\) 的时候,减去的相当于是选出的两条路径的对称差。考虑一个贪心,第一次取某一条直径一定不会更劣,这个感性理解就好。然后把这个直径的边权都设成 \(-1\) 再跑一遍直径就行了。

点击查看代码
const int N=1e5+13;
struct Edge{int v,nxt;}e[N<<1];
int n,k,h[N],etot=1;
inline void add_edge(int u,int v){e[++etot]=(Edge){v,h[u]};h[u]=etot;}

namespace sub1{
int ans,dep[N];
void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].v;if(v==fa) continue;
		dfs(v,u);
	}
}
inline void solve(){
	dfs(1,0);
	int rt=1;
	for(int i=2;i<=n;++i)
		if(dep[i]>dep[rt]) rt=i;
	dfs(rt,0);
	for(int i=1;i<=n;++i)
		if(dep[i]-1>ans) ans=dep[i]-1;
	println(2*(n-1)-ans+1);
}
}

namespace sub2{
int ans,dep[N],mx1[N],mx2[N];
bool vis[N<<1];
void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].v;if(v==fa) continue;
		dfs(v,u);
	}
}
bool dfs1(int u,int fa,int ed){
	if(u==ed) return 1;
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].v;if(v==fa) continue;
		if(dfs1(v,u,ed)) return vis[i]=vis[i^1]=1;
	}
	return 0;
}
void dfs2(int u,int fa){
	mx1[u]=0,mx2[u]=-INF;
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].v;if(v==fa) continue;
		int w=(vis[i]?-1:1);
		dfs2(v,u);
		if(mx1[v]+w>mx1[u]) mx2[u]=mx1[u],mx1[u]=mx1[v]+w;
		else if(mx1[v]+w>mx2[u]) mx2[u]=mx1[v]+w;
	}
	ans=max(ans,mx1[u]+mx2[u]);
}
inline void solve(){
	dfs(1,0);
	int x=1;
	for(int i=2;i<=n;++i)
		if(dep[i]>dep[x]) x=i;
	dfs(x,0);
	int y=1;
	for(int i=2;i<=n;++i)
		if(dep[i]>dep[y]) y=i;
	dfs1(x,0,y);
	dfs2(1,0);
	println(2*(n-1)-(ans+dep[y]-1)+2);
}
}

int main(){
	read(n),read(k);
	for(int i=1;i<n;++i){int u,v;read(u),read(v);add_edge(u,v),add_edge(v,u);}
	if(k==1) sub1::solve();
	else sub2::solve();
	return 0;
}

标签:rt,int,题解,APIO2010,mx2,dfs,mx1,P3629,ans
来源: https://www.cnblogs.com/winterfrost/p/p3629-solution.html