zoukankan      html  css  js  c++  java
  • 【UOJ #310】【UNR #2】—黎明前的巧克力(FWT)

    传送门


    考虑说2个人选的集合不相交而且分别异或起来相等

    其实就是有多少个集合划分异或和为0
    对于每个巧克力构建形式幂级数(1+2xai)(1+2x^{a_i})

    显然全部FwtFwt之后乘起来再IFwtIFwt回去后就是a01a_0-1
    但是直接FwtFwt肯定要爆炸

    考虑分析一下这个东西
    对于SSFwtFwt化成的点值实际上是T2U(1)STaTsum_{Tin 2^{U}}(-1)^{|Sigcap T|}a_T

    11对点值的贡献一定是11
    那么每一项要么是33要么是1-1

    又由于FwtFwt的和等于和的FwtFwt
    由于每个集合的贡献是加起来的,所以也很显然

    于是把形式幂级数加起来,做一个FwtFwt
    考虑有多少个1-1和多少个33
    然后就相当于解一个模意义下的一次方程
    k+3(nk)=as-k+3(n-k)=a_s
    解出来个数就完了

    #include<bits/stdc++.h>
    using namespace std;
    #define gc getchar
    inline int read(){
    	char ch=gc();
    	int res=0,f=1;
    	while(!isdigit(ch))f^=ch=='-',ch=gc();
    	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    	return f?res:-res;
    }
    #define re register
    #define pb push_back
    #define cs const
    #define pii pair<int,int>
    #define fi first
    #define se second
    #define ll long long
    #define poly vector<int>
    #define bg begin
    cs int mod=998244353,G=3;
    inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
    inline void Add(int &a,int b){(a+=b)>=mod?(a-=mod):0;}
    inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
    inline void Dec(int &a,int b){(a-=b)<0?(a+=mod):0;}
    inline int mul(int a,int b){return 1ll*a*b>=mod?1ll*a*b%mod:a*b;}
    inline void Mul(int &a,int b){a=mul(a,b);}
    inline int ksm(int a,int b,int res=1){
    	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));return res;
    }
    inline void chemx(int &a,int b){a<b?a=b:0;}
    inline void chemn(int &a,int b){a>b?a=b:0;}
    cs int N=20,M=(1<<N)|5,inv4=ksm(4,mod-2);
    int n,sta,f[M],a[M],p3[M];
    inline void Fwt(int *f,int lim,int kd){
    	for(int mid=1,a0,a1;mid<lim;mid<<=1)
    	for(int i=0;i<lim;i+=(mid<<1))
    	for(int j=0;j<mid;j++)
    	a0=f[i+j],a1=f[i+j+mid],f[i+j]=add(a0,a1),f[i+j+mid]=dec(a0,a1);
    	if(kd==-1)
    	for(int i=0,inv=ksm(lim,mod-2);i<lim;i++)Mul(f[i],inv);
    }
    int main(){
    	n=read(),sta=1<<N;
    	p3[0]=1;for(int i=1;i<M;i++)p3[i]=mul(p3[i-1],3);
    	for(int i=1;i<=n;i++)a[i]=read(),f[0]++,f[a[i]]+=2;
    	Fwt(f,sta,1);
    	for(int i=0;i<sta;i++){
    		int x=f[i],c1=mul(dec(3*n,x),inv4),c3=dec(n,c1);
    		f[i]=(c1&1)?mod-p3[c3]:p3[c3];	
    	}
    	Fwt(f,sta,-1);
    	cout<<dec(f[0],1);
    }
    
  • 相关阅读:
    前端异常上报
    前端异常解析:Source Map
    前端操作剪切板不完全指北
    多系统之间模块相互引用的引发的深思
    浅谈vue原理(四)
    浅谈vue原理(三)
    浅谈vue原理(二)
    浅谈vue原理(一)
    vue中路由嵌套的作用
    常用的学习网站和小工具
  • 原文地址:https://www.cnblogs.com/stargazer-cyk/p/12328664.html
Copyright © 2011-2022 走看看