洛谷P4717 【模板】快速沃尔什变换（FWT）

传送门

这玩意儿太骚了……

参考了yyb巨佬的

 //minamoto
 #include<iostream>
 #include<cstdio>
 #define ll long long
 #define add(x,y) ((x+=y)>=mod?x-=mod:x)
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 char sr[<<],z[];int C=-,Z;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(int x){
     if(C><<)Ot();if(x<)sr[++C]=,x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]=' ';
 }
 const int N=(<<)+,mod=,inv=;
 int n,len,A[N],B[N],R[N];
 void FWT_or(int *A,int type){
     for(int i=;i<len;i<<=)
     for(int p=i<<,j=;j<len;j+=p)
     for(int k=;k<i;++k)
     add(A[i+j+k],type==?A[j+k]:mod-A[j+k]);
 }
 void FWT_and(int *A,int type){
     for(int i=;i<len;i<<=)
     for(int p=i<<,j=;j<len;j+=p)
     for(int k=;k<i;++k)
     add(A[j+k],type==?A[i+j+k]:mod-A[i+j+k]);
 }
 void FWT_xor(int *A,int type){
     for(int i=;i<len;i<<=)
     for(int p=i<<,j=;j<len;j+=p)
     for(int k=;k<i;++k){
         int x=A[j+k],y=A[i+j+k];
         A[j+k]=(x+y)%mod,A[i+j+k]=(x+mod-y)%mod;
         if(type==-) A[j+k]=1ll*A[j+k]*inv%mod,A[i+j+k]=1ll*A[i+j+k]*inv%mod;
     }
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     n=read(),len=<<n;
     for(int i=;i<len;++i) A[i]=read();
     for(int i=;i<len;++i) B[i]=read();

     FWT_or(A,),FWT_or(B,);
     for(int i=;i<len;++i) R[i]=1ll*A[i]*B[i]%mod;
     FWT_or(A,-),FWT_or(B,-),FWT_or(R,-);
     for(int i=;i<len;++i) print(R[i]);sr[++C]='\n';

     FWT_and(A,),FWT_and(B,);
     for(int i=;i<len;++i) R[i]=1ll*A[i]*B[i]%mod;
     FWT_and(A,-),FWT_and(B,-),FWT_and(R,-);
     for(int i=;i<len;++i) print(R[i]);sr[++C]='\n';

     FWT_xor(A,),FWT_xor(B,);
     for(int i=;i<len;++i) R[i]=1ll*A[i]*B[i]%mod;
     FWT_xor(R,-);
     for(int i=;i<len;++i) print(R[i]);sr[++C]='\n';

     Ot();
     return ;
 }