zoukankan      html  css  js  c++  java
  • NTT FWT(xor or and) 模板

    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;
        }
    }
    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=0;i<len;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)*inv%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
  • 相关阅读:
    linux命令大全
    IP协议
    TCP、IP、ARP协议之间的工作关系
    程序出现问题后
    HTTP(超文本传输协议)
    多线程
    syslog.conf文件
    logger命令
    gdb
    二、数据的存储结构
  • 原文地址:https://www.cnblogs.com/zquzjx/p/9473670.html
Copyright © 2011-2022 走看看