2020牛客暑期多校训练营（第二场） E、Exclusive OR（FWT）

题目链接

题面：

题意：
给定一个序列A，对于 i 从1–n，输出每次选 i 个数的时候的异或最大值，其中每个数可重复选。

官方题解：

题解：

先了解一下fwt是干啥的。
因为可以重复选数，我们可以构造一个a数组和一个b数组，初始时在令 b [ A [ i ] ] = a [ A [ i ] ] = 1。
初始时，a 数组是 i = 1的答案
然后进行a=a^b，此时a数组是 i = 2 的答案。
再进行一次a=a^b，那么a数组就是 i == 3的答案。
因为有可能在取模时，导致不为0的变成0了，所以每次可以给a数组不为0的重新赋值1即可。

因为Ai<2¹⁸ ，那么如果能满秩的话，在 i <= 18就会满秩。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#include<bitset>
#include<map>
#include<unordered_map>
#include<set>
#define ui unsigned int
#define ll long long
#define llu unsigned ll
#define ld long double
#define pr make_pair
#define pb push_back
#define lc (cnt<<1)
#define rc (cnt<<1|1)
#define len(x)  (t[(x)].r-t[(x)].l+1)
#define tmid ((l+r)>>1)
using namespace std;

const int inf=0x3f3f3f3f;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const double dnf=1e18;
const int mod=1e9+7;
const double eps=1e-8;
const double pi=acos(-1.0);
const int hp=13331;
const int maxn=400100;
const int maxm=100100;
const int up=100000;

int a[maxn],b[maxn],c[maxn],ans[maxn];
int inv2=500000004;

void fwt_or_mod(int *x,int len,int opt)
{
    for(int i=1;i<len;i<<=1)
    {
        for(int p=i<<1,j=0;j<len;j+=p)
        {
            for(int k=0;k<i;k++)
            {
                if(opt==1) x[i+j+k]=(x[j+k]+x[i+j+k])%mod;
                else x[i+j+k]=(x[i+j+k]+mod-x[j+k])%mod;
            }
        }
    }
}

void fwt_and_mod(int *x,int len,int opt)
{
    for(int i=1;i<len;i<<=1)
    {
        for(int p=i<<1,j=0;j<len;j+=p)
        {
            for(int k=0;k<i;k++)
            {
                if(opt==1) x[j+k]=(x[j+k]+x[i+j+k])%mod;
                else x[j+k]=(x[j+k]+mod-x[i+j+k])%mod;
            }
        }
    }
}

void fwt_xor_mod(int *x,int len,int opt)
{
    for(int i=1;i<len;i<<=1)
    {
        for(int p=i<<1,j=0;j<len;j+=p)
        {
            for(int k=0;k<i;k++)
            {
                int xx=x[j+k],yy=x[i+j+k];
                x[j+k]=(xx+yy)%mod;
                x[i+j+k]=(xx+mod-yy)%mod;
                if(opt==-1)
                    x[j+k]=1ll*x[j+k]*inv2%mod,x[i+j+k]=1ll*x[i+j+k]*inv2%mod;
            }
        }
    }
}

int main(void)
{
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&c[i]);
        a[c[i]]=b[c[i]]=1;
    }
    int len=(1<<18);
    fwt_xor_mod(b,len,1);
    for(int i=1;i<=19;i++)
    {
        for(int j=0;j<len;j++)
        {
            if(a[j])
                ans[i]=j,a[j]=1;
        }
        fwt_xor_mod(a,len,1);
        for(int j=0;j<len;j++)
            a[j]=1ll*a[j]*b[j]%mod;
        fwt_xor_mod(a,len,-1);
    }
    for(int i=20;i<=n;i++)
        ans[i]=ans[i-2];
    for(int i=1;i<=n;i++)
        printf("%d ",ans[i]);
    putchar('\n');
    return 0;
}