uoj207 共价大爷游长沙(lct+Hash)

考虑我们给S中的每条路径一个权值,然后对这段路径异或上这个值。
那么最后询问的那条边如果权值等于S集合的异或和那么就说它被所有路径经过了。这样的正确性是很高的。
然后考虑新加入一条边,就会出现一个环,然后肯定是在这个环上删去一条边,假设是x,y那么在新的树上x->y的路径上异或上x,y这条边的异或值就可以得到新树的正确异或了。(画画图理解一下就是对存在性取反了)
因此维护一棵边权lct,链修改,单点询问即可。
复杂度 O((n+m)logn) O ( ( n + m ) l o g n )

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
#define mod 1000000007
#define pa pair<int,int>
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,fa[N<<1],c[N<<1][2],v[N<<1],lazy[N<<1],rt,sum,owo=0,qq[N<<1];
bool rev[N<<1];map<pa,int>mp;
struct Icefox{
    int x,y,val;
}S[N*3];
inline bool isroot(int x){return x!=c[fa[x]][0]&&x!=c[fa[x]][1];}
inline void rotate(int x){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(!isroot(y)) c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;
}
inline void pushdown(int p){
    int l=c[p][0],r=c[p][1];
    if(rev[p]) rev[l]^=1,rev[r]^=1,rev[p]=0,swap(c[p][0],c[p][1]);
    if(lazy[p]) lazy[l]^=lazy[p],lazy[r]^=lazy[p],v[l]^=lazy[p],v[r]^=lazy[p],lazy[p]=0;
}
inline void splay(int x){
    int top=0;qq[++top]=x;
    for(int xx=x;!isroot(xx);xx=fa[xx]) qq[++top]=fa[xx];
    while(top) pushdown(qq[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(x==c[y][1]^y==c[z][1]) rotate(x);
            else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;while(x){splay(x);c[x][1]=y;y=x;x=fa[x];}
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void link(int x,int y){
    makeroot(x);fa[x]=y;
}
inline void cut(int x,int y){
    makeroot(x);access(y);splay(y);c[y][0]=fa[x]=0;
}
inline void add(int x,int y,int val){
    makeroot(x);access(y);splay(y);lazy[y]^=val;v[y]^=val;
}
inline int ask(int p){
    splay(p);return v[p];
}
int main(){
//  freopen("a.in","r",stdin);
//  freopen("a.out","w",stdout);
    read();n=read();m=read();srand(20000711);
    for(int i=1;i<n;++i){
        int x=read(),y=read();if(x>y) swap(x,y);mp[make_pair(x,y)]=i+n;link(x,i+n);link(i+n,y);
    }while(m--){
        int ty=read();
        if(ty==1){
            int x=read(),y=read(),u=read(),vv=read();if(x>y) swap(x,y);if(u>vv) swap(u,vv);
            int id=mp[make_pair(x,y)];mp.erase(make_pair(x,y));mp[make_pair(u,vv)]=id;
            cut(x,id);cut(y,id);link(u,id);link(vv,id);int tmp=v[id];v[id]=0;add(x,y,tmp);continue;
        }if(ty==2){
            S[++owo].x=read();S[owo].y=read();S[owo].val=rand()%mod;add(S[owo].x,S[owo].y,S[owo].val);sum^=S[owo].val;continue;
        }if(ty==3){int x=read();sum^=S[x].val;add(S[x].x,S[x].y,S[x].val);continue;}
        int x=read(),y=read();if(x>y) swap(x,y);puts(ask(mp[make_pair(x,y)])==sum?"YES":"NO");
    }return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值