FWT模板代码

opt是1表示正变换，-1表示逆变换。原理证明参考https://www.cnblogs.com/cjyyb/p/9065615.html
实际这个板子也是从那里抄的。但是知道原理没啥用，过一段就忘了。
题目链接：https://ac.nowcoder.com/acm/contest/5667/E

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int N;
const int NN=2e5+10;
const int MM=(1<<19)+10;
long long a[MM],b[MM];
void FWT_or(long long *a,int opt)
{
    for(int i=1;i<N;i<<=1)
        for(int p=i<<1,j=0;j<N;j+=p)
            for(int k=0;k<i;++k)
                if(opt==1)a[i+j+k]=(a[j+k]+a[i+j+k]);
                else a[i+j+k]=(a[i+j+k]-a[j+k]);
}
void FWT_and(long long *a,int opt)
{
    for(int i=1;i<N;i<<=1)
        for(int p=i<<1,j=0;j<N;j+=p)
            for(int k=0;k<i;++k)
                if(opt==1)a[j+k]=(a[j+k]+a[i+j+k]);
                else a[j+k]=(a[j+k]-a[i+j+k]);
}
void FWT_xor(long long *a,int opt)
{
    for(int i=1;i<N;i<<=1)
        for(int p=i<<1,j=0;j<N;j+=p)
            for(int k=0;k<i;++k)
            {
                long long X=a[j+k],Y=a[i+j+k];
                a[j+k]=(X+Y);a[i+j+k]=(X-Y);
                if(opt==-1)a[j+k]=1ll*a[j+k]/2,a[i+j+k]=1ll*a[i+j+k]/2;
            }
}
int ans[NN];
int main() {
    int n;scanf("%d",&n);
    int maxn=0;
    for(int i=1;i<=n;i++){
        int x;scanf("%d",&x);
        maxn=max(maxn,x);
        b[x]++;
    }
    for(N=1;N<maxn;N<<=1);N<<=1;
    a[0]=1;
    FWT_xor(b,1);
    for(int i=1;i<=19;i++){
        FWT_xor(a,1);
        for(int j=0;j<N;j++)a[j]=a[j]*b[j];
        FWT_xor(a,-1);
        for(int j=0;j<N;j++){
            if(a[j]>0){a[j]=1;ans[i]=j;}
        }
    }
    for(int i=1;i<=n;i++){
        if(i>19){
            if((i-19)%2){
                printf("%d ",ans[18]);
            }
            else{
                printf("%d ",ans[19]);
            }
        }
        else{
            printf("%d ",ans[i]);
        }
    }
    return 0;
}
/*
23
1 2 4 8 16 32 64 128 256 512 1024 2048 4096 8192 16384 32768 65536 131072 0 0 0 0 0
*/