CF1413F Roads and Ramen

XVII.CF1413F Roads and Ramen

首先,注意到本题等价于求路径上所有边权的异或和为\(0\)的路径长度的最大值。

然后,我们要猜/证明出一个结论,即任意一条极长合法路径,其必有一个端点是直径端点。

证明:

我们设有一条直径\((S,T)\)。我们再设\(col_i\)表示从\(i\)节点到根节点(不妨设为\(1\)号节点)上所有边权的异或和。

  1. \(col_S=col_T\)

    此时显然有\((S,T)\)合法,又有\((S,T)\)是直径,故\((S,T)\)即为答案。

  2. \(col_S\neq col_T\)

    假设当前我们有一条合法路径\((s,t)\);则显然,必有\(s,t\)的颜色同\(S,T\)中某一个点的颜色相等,假设是\(S\)。则\((s,S),(t,S)\)两条路径中较长的一条的长度肯定大于等于\((s,t)\)的长度,这是直径的性质。

于是我们就直接找到直径\((S,T)\),对二者各维护一棵线段树,查询所有到它距离为偶数的点中深度最大值即可。

时间复杂度\(O(n\log n)\)

代码:

#include<bits/stdc++.h>
using namespace std;
int n,m,head[500100],cnt;
struct node{
	int to,next;
	bool tp;
}edge[1001000];
struct SegTree{
	int dfn[500100],sz[500100],rev[500100],dep[500100],tot;
	bool ini[500100];
	int mx[2000100][2];
	bool tag[2001000];
	#define lson x<<1
	#define rson x<<1|1
	#define mid ((l+r)>>1)
	#define REV(x) tag[x]^=1,swap(mx[x][0],mx[x][1])
	void pushdown(int x){if(tag[x])REV(lson),REV(rson),tag[x]=false;}
	void pushup(int x){mx[x][0]=max(mx[lson][0],mx[rson][0]),mx[x][1]=max(mx[lson][1],mx[rson][1]);}
	void dfs(int x,int fa){
		dfn[x]=++tot,sz[x]=1,rev[tot]=x;
		for(int i=head[x];i!=-1;i=edge[i].next)if(edge[i].to!=fa)dep[edge[i].to]=dep[x]+1,ini[edge[i].to]=ini[x]^edge[i].tp,dfs(edge[i].to,x),sz[x]+=sz[edge[i].to];
	}
	void build(int x,int l,int r){
		if(l==r){mx[x][ini[rev[l]]]=dep[rev[l]];return;}
		build(lson,l,mid),build(rson,mid+1,r),pushup(x);
	}
	void modify(int x,int l,int r,int L,int R){
		if(l>R||r<L)return;
		if(L<=l&&r<=R){REV(x);return;}
		pushdown(x),modify(lson,l,mid,L,R),modify(rson,mid+1,r,L,R),pushup(x);
	}
}seg[2];
void ae(int u,int v,bool w){
	edge[cnt].next=head[u],edge[cnt].to=v,edge[cnt].tp=w,head[u]=cnt++;
	edge[cnt].next=head[v],edge[cnt].to=u,edge[cnt].tp=w,head[v]=cnt++;
}
pair<int,int>dfs(int x,int fa){
	pair<int,int>ret=make_pair(-1,x);
	for(int i=head[x];i!=-1;i=edge[i].next)if(edge[i].to!=fa)ret=max(ret,dfs(edge[i].to,x));
	ret.first++;
	return ret;
}
int main(){
	scanf("%d",&n),memset(head,-1,sizeof(head));
	for(int i=1,x,y,z;i<n;i++)scanf("%d%d%d",&x,&y,&z),ae(x,y,z);
	int S=dfs(1,0).second,T=dfs(S,0).second;
	seg[0].dfs(S,0),seg[1].dfs(T,0);
	seg[0].build(1,1,n),seg[1].build(1,1,n);
	scanf("%d",&m);
	for(int i=1,x;i<=m;i++){
		scanf("%d",&x),x--;
		int u=edge[x<<1].to,v=edge[x<<1|1].to;
		if(seg[0].dep[u]<seg[0].dep[v])swap(u,v);
		seg[0].modify(1,1,n,seg[0].dfn[u],seg[0].dfn[u]+seg[0].sz[u]-1);
		if(seg[1].dep[u]<seg[1].dep[v])swap(u,v);
		seg[1].modify(1,1,n,seg[1].dfn[u],seg[1].dfn[u]+seg[1].sz[u]-1);
		printf("%d\n",max(seg[0].mx[1][0],seg[1].mx[1][0]));
	}
	return 0;
} 

原文地址:https://www.cnblogs.com/Troverld/p/14611360.html