把博客园图标替换成自己的图标
把博客园图标替换成自己的图标end

CF1616H Keep XOR Low

题面传送门
终究还是太菜了,赛时没冲出来不说了,到最后搭了个拍才发现自己错在哪里。
看上去这种异或题肯定是扔到Trie上去。然后看上去就是个先确定前多少位相同,然后一位\(1\)变成\(0\)然后随便算?
\(f_{u}\)为Trie上\(u\)这个点的子树内的答案,然后发现好像不太好算?
因为如果\(x\)这一位为\(1\)那么继续\(1\)下去要分割成两个子树各选至少一个点的答案,就不太好做。
那就设\(f_{u,v}\)为这两个点各至少选一个的答案。
如果这一位为\(0\)那么显然\(f_{u,v}=f_{u_0,v_0}+f_{u_1,v_1}\)
如果这一位为\(1\),那么这一位为\(0\)的任意选,如果这一位为\(1\)那么就是\(f_{u_0,v_1}\)\(f_{u_1,v_0}\)
但是这样不对,因为可以剩下的位还可以选,加上就好了。
时间复杂度因为每个节点只会遍历一次,所以显然是\(O(n\log w)\)
code:

#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define RI re int
#define ll long long
#define db double
#define lb long db
#define N (150000+5)
#define M (1<<20)
#define mod 998244353
#define Mod (mod-1)
#define eps (1e-9)
#define U unsigned int
#define it iterator
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) (n*(x-1)+(y))
#define R(n) (rand()*rand()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define PB push_back
using namespace std;
int n,x,A[N+5],Ro;ll Ans,Po[N+5];
namespace Tree{
	int cnt,L[N*40+5],R[N*40+5],F[N*40+5];I void Ins(int x,int &now,int d){!now&&(now=++cnt);F[now]++;if(d==-1) return;Ins(x,x>>d&1?R[now]:L[now],d-1);}
	I ll calc(int N1,int N2,int d){
		if(!N1||!N2) return 0;if(d==-1)return (Po[F[N1]]-1)*(N1^N2?Po[F[N2]]-1:1)%mod;
		if(x>>d&1){
			if(N1==N2) return /*printf("%lld\n",x=calc(L[N1],R[N1],d-1)),*/(calc(L[N1],R[N1],d-1)+Po[F[L[N1]]]-1+Po[F[R[N1]]]-1)%mod;
			else{
				ll F1=calc(L[N1],R[N2],d-1),F2=calc(R[N1],L[N2],d-1);
				return ((F1+1)*(F2+1)-1+F1*(Po[F[R[N1]]]-1+Po[F[L[N2]]]-1)+F2*(Po[F[L[N1]]]-1+Po[F[R[N2]]]-1)+(Po[F[L[N1]]]-1)*(Po[F[L[N2]]]-1)+(Po[F[R[N1]]]-1)*(Po[F[R[N2]]]-1))%mod;
			} //return ((calc(L[N1],R[N2],d-1)+1)*(calc(R[N1],L[N2],d-1)+1)-1+(Po[F[L[N1]]]-1)*(Po[F[L[N2]]]-1)+(Po[F[R[N1]]]-1)*(Po[F[R[N2]]]-1))%mod;
		} 
		else return (calc(L[N1],L[N2],d-1)+calc(R[N1],R[N2],d-1))%mod;
	}
}
int main(){
//	freopen("1.in","r",stdin);//freopen("2.out","w",stdout);
	RI i;scanf("%d%d",&n,&x);for(i=1;i<=n;i++) scanf("%d",&A[i]),Tree::Ins(A[i],Ro,30);
	for(Po[0]=i=1;i<=n;i++) Po[i]=Po[i-1]*2%mod;printf("%lld\n",Tree::calc(Ro,Ro,30)%mod);
}
posted @ 2022-03-16 18:02  275307894a  阅读(56)  评论(0编辑  收藏  举报
浏览器标题切换
浏览器标题切换end