NNT FWT(xor or and) 模板

Posted 2020-12-29 zquzjx

 

void nnt(int a[],int len,int on)
{
    for(int i=0;i<len;i++)
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=1;i<len;i<<=1) {
        int wn=mod_pow(3,(mod-1)/(i<<1));
        for(int j=0;j<len;j+=(i<<1)) {
            int w=1;
            for(int k=0;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==-1) {
        reverse(a+1,a+len);
        int inv=mod_pow(len,mod-2);
        for(int i=0;i<len;i++)
            a[i]=1ll*a[i]*inv%mod;
    }
}

NNT

 

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=0;i<n;++i) 
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=1;i<len;i<<=1)
        for(int j=0;j<len;j+=(i<<1))
            for(int k=0;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)%mod;
            }
}

fwtXor

void fwtAnd(int a[],int len,int on)
{
    for(int i=0;i<n;++i) 
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=1;i<len;i<<=1)
        for(int j=0;j<len;j+=(i<<1))
            for(int k=0;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=0;i<n;++i) 
        if(i<r[i]) swap(a[i],a[r[i]]);
    for(int i=1;i<len;i<<=1)
        for(int j=0;j<len;j+=(i<<1))
            for(int k=0;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