BZOJ4811 [Ynoi2017]由乃的OJ

显然用LCT维护一下经过一条链之后每一位原来是1会变成什么原来是0会变成什么即可

这样的话LCT复杂度是n log n * k的,链剖是n log^2 n * k的,都过不了

考虑把这个k给去掉,我们直接开两个unsigned long long,v0和v1,v0的每一位代表这一位原来是0的话经过这条链会变成什么,v1的每一位代表这一位原来是1的话经过这条链会变成什么

那么考虑信息的合并,假如你现在知道了v0[x],v1[x],v0[y],v1[y],先经过x后经过y,那么合并之后新的v0=(v0[x]&v1[y])|(~v0[x]&v0[y]),v1=(v1[x]&v1[y])|(~v1[x]&v0[y])

比较好理解吧……就是比如v0,原来的0变成1那些位就会变成v1[y]中对应的位,变成0的那些位会变成v0[y]中对应的位,那么与一下,然后合起来就好了

链剖和LCT我都写了,结果LCT TLE了,链剖3s过了

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
#define MAXN 100010
#define MAXM 1010
#define ull unsigned long long
#define eps 1e-8
#define MOD 1000000007
#define INF 1000000000
char xB[1<<15],*xS=xB,*xTT=xB;
#define getc() (xS==xTT&&(xTT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xTT)?0:*xS++)
#define isd(c) (c>='0'&&c<='9')
ull read(){
	char xchh;
	ull xaa;
    while(xchh=getc(),!isd(xchh));(xaa=xchh-'0');
    while(xchh=getc(),isd(xchh))xaa=xaa*10+xchh-'0';return xaa;
}
struct vec{
	int to;
	int fro;
};
struct data{
	ull v0;
	ull v1;
	data(){
		
	}
	data(int o,ull x){
		if(o==1){
			v0=0&x;
			v1=(~0)&x;
		}
		if(o==2){
			v0=0|x;
			v1=(~0)|x;
		}
		if(o==3){
			v0=0^x;
			v1=(~0)^x;
		}
	}
	friend data operator +(data x,data y){
		data z;
		z.v0=(x.v0&y.v1)|((~x.v0)&y.v0);
		z.v1=(x.v1&y.v1)|((~x.v1)&y.v0);
		return z;
	}
};
vec mp[MAXN*2];
int tai[MAXN],cnt;
int dfn[MAXN],fa[MAXN],son[MAXN],siz[MAXN],dep[MAXN],ndf[MAXN],tp[MAXN],tim;
int op[MAXN];
ull a[MAXN];
data vl[MAXN<<2],vr[MAXN<<2];
int n,m,k;
inline void be(int x,int y){
	mp[++cnt].to=y;
	mp[cnt].fro=tai[x];
	tai[x]=cnt;
}
inline void bde(int x,int y){
	be(x,y);
	be(y,x);
}
inline void ud(int x){
	vl[x]=vl[x<<1]+vl[x<<1|1];
	vr[x]=vr[x<<1|1]+vr[x<<1];
}
void dfs(int x){
	int i,y;
	siz[x]=1;
	dep[x]=dep[fa[x]]+1;
	for(i=tai[x];i;i=mp[i].fro){
		y=mp[i].to;
		if(!siz[y]){
			fa[y]=x;
			dfs(y);
			siz[x]+=siz[y];
			if(siz[y]>siz[son[x]]){
				son[x]=y;
			}
		}
	}
}
void dfs2(int x,int z){
	int i,y;
	ndf[dfn[x]=++tim]=x;
	tp[x]=z;
	if(son[x]){
		dfs2(son[x],z);
		for(i=tai[x];i;i=mp[i].fro){
			y=mp[i].to;
			if(!dfn[y]){
				dfs2(y,y);
			}
		}
	}
}
void build(int x,int y,int z){
	if(y==z){
		vl[x]=vr[x]=data(op[ndf[y]],a[ndf[y]]);
		return ;
	}
	int mid=y+z>>1;
	build(x<<1,y,mid);
	build(x<<1|1,mid+1,z);
	ud(x);
}
void change(int x,int y,int z,int p){
	if(y==z){
		vl[x]=vr[x]=data(op[ndf[y]],a[ndf[y]]);
		return ;
	}
	int mid=y+z>>1;
	if(p<=mid){
		change(x<<1,y,mid,p);
	}else{
		change(x<<1|1,mid+1,z,p);
	}
	ud(x);
}
data askl(int x,int y,int z,int l,int r){
	if(y==l&&z==r){
		return vl[x];
	}
	int mid=y+z>>1;
	if(r<=mid){
		return askl(x<<1,y,mid,l,r);
	}else if(l>mid){
		return askl(x<<1|1,mid+1,z,l,r);
	}else{
		return askl(x<<1,y,mid,l,mid)+askl(x<<1|1,mid+1,z,mid+1,r);
	}
}
data askr(int x,int y,int z,int l,int r){
	if(y==l&&z==r){
		return vr[x];
	}
	int mid=y+z>>1;
	if(r<=mid){
		return askr(x<<1,y,mid,l,r);
	}else if(l>mid){
		return askr(x<<1|1,mid+1,z,l,r);
	}else{
		return askr(x<<1|1,mid+1,z,mid+1,r)+askr(x<<1,y,mid,l,mid);
	}
}
data toask(int x,int y){
	data vx(3,0),vy(3,0);
	while(tp[x]!=tp[y]){
		if(dep[tp[x]]>=dep[tp[y]]){
			vx=vx+askr(1,1,n,dfn[tp[x]],dfn[x]);
			x=fa[tp[x]];
		}else{
			vy=askl(1,1,n,dfn[tp[y]],dfn[y])+vy;
			y=fa[tp[y]];
		}
	}
	if(dep[x]>dep[y]){
		return vx+askr(1,1,n,dfn[y],dfn[x])+vy;
	}else{
		return vx+askl(1,1,n,dfn[x],dfn[y])+vy;
	}
}
int main(){
	int i,o,x,y;
	ull z;
	n=read();
	m=read();
	k=read();
	for(i=1;i<=n;i++){
		op[i]=read();
		a[i]=read();
	}
	for(i=1;i


你可能感兴趣的:(BZOJ,树链剖分)