NTT FWT(xor or and) 模板

void nnt(int a[],int len,int on)
{
    for(int i=;i<len;i++)
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=;i<len;i<<=) {
        int wn=mod_pow(,(mod-)/(i<<));
        for(int j=;j<len;j+=(i<<)) {
            int w=;
            for(int k=;k<i;k++,w=1ll*w*wn%mod) {
                int u=a[j+k], v=1ll*w*a[j+k+i]%mod;
                a[j+k]=(u+v)%mod, a[j+k+i]=(u-v+mod)%mod;
            }
        }
    }
    if(on==-) {
        reverse(a+,a+len);
        int inv=mod_pow(len,mod-);
        for(int i=;i<len;i++)
            a[i]=1ll*a[i]*inv%mod;
    }
}

NTT

FWT讲解：

https://blog.csdn.net/no_name233/article/details/52824587

https://blog.csdn.net/qq_34454069/article/details/79524001

void fwtXor(int a[],int len,int on)
{
    for(int i=;i<len;i++)
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=;i<len;i<<=)
        for(int j=;j<len;j+=(i<<))
            for(int k=;k<i;k++){
                int u=a[j+k], v=a[j+k+i];
                if(on)
                    a[j+k]=(u+v)%mod, a[j+k+i]=(u-v+mod)%mod;
                else
                    a[j+k]=(u+v)*inv%mod, a[j+k+i]=(u-v+mod)*inv%mod;
            }
}

fwtXor

void fwtAnd(int a[],int len,int on)
{
    for(int i=;i<n;++i)
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=;i<len;i<<=)
        for(int j=;j<len;j+=(i<<))
            for(int k=;k<i;k++) {
                int u=a[j+k], v=a[j+k+i];
                if(on) a[j+k]=(u+v)%mod, a[j+k+i]=v;
                else a[j+k]=(u-v+mod)%mod, a[j+k+i]=v;
            }
}

fwtAnd

void fwtOr(int a[],int len,int on)
{
    for(int i=;i<n;++i)
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=;i<len;i<<=)
        for(int j=;j<len;j+=(i<<))
            for(int k=;k<i;k++) {
                int u=a[j+k], v=a[j+k+i];
                if(on) a[j+k]=u, a[j+k+i]=(v+u)%mod;
                else a[j+k]=u, a[j+k+i]=(v-u+mod)%mod;
            }
}

fwtOr