【模板】快速沃尔什变换FWT

最新推荐文章于 2024-01-21 21:18:35 发布

ccosi

最新推荐文章于 2024-01-21 21:18:35 发布

阅读量237

点赞数 1

分类专栏： -------MATH------- FWT

本文链接：https://blog.csdn.net/corsica6/article/details/89258497

版权

-------MATH------- 同时被 2 个专栏收录

80 篇文章 1 订阅

订阅专栏

FWT

1 篇文章 0 订阅

订阅专栏

本文介绍了快速沃尔什变换（FWT）的基础知识，包括其在子集运算中的应用，如加法、逆变换、与运算和异或运算。FWT在处理这些问题时展现出高效性，与快速傅里叶变换（FFT）有相似之处。同时，文章提到逆变换的操作，如减去真子集或真包含于的集合，以及异或操作的处理方式。

摘要由CSDN通过智能技术生成

传送门:洛谷【模板】快速沃尔什变换

存板子

或：所有子集(加) 逆变换：减去真子集
与：所有它包含于的集合(加) 逆变换：减去所有真包含与的集合
异或：类似于FFT 逆变换：同正变换，不过值每次要/2

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=(1<<18)+100,mod=998244353,ivt=499122177;

int K,n,a[N],b[N],f[N],g[N],c[N];

inline int ad(int x,int y){x+=y;return x>=mod?x-mod:x;}
inline int dc(int x,int y){x-=y;return x<0?x+mod:x;}

inline void fwt_or(int *a,int pr)
{
    register int i,j,k;
    for(i=1;i<n;i<<=1)
     for(j=0;j<n;j+=(i<<1))
      for(k=0;k<i;++k)
      	 pr?a[i+k+j]=ad(a[i+k+j],a[j+k]):
            a[i+k+j]=dc(a[i+k+j],a[j+k]);
}

inline void fwt_and(int *a,int pr)
{
    register int i,j,k;
    for(i=1;i<n;i<<=1)
     for(j=0;j<n;j+=(i<<1))
      for(k=0;k<i;++k)
      	 pr?a[j+k]=ad(a[j+k],a[i+j+k]):
            a[j+k]=dc(a[j+k],a[i+j+k]);
}

inline void fwt_xor(int *a,int pr)
{
    register int i,j,k,x,y;
    for(i=1;i<n;i<<=1)
     for(j=0;j<n;j+=(i<<1))
      for(k=0;k<i;++k){
          x=a[j+k];y=a[i+j+k];
          pr?(a[j+k]=ad(x,y),a[i+j+k]=dc(x,y)):
          	 (a[j+k]=(ll)ad(x,y)*ivt%mod,a[i+j+k]=(ll)dc(x,y)*ivt%mod);
      }
}

inline int fp(int x,int y)
{
    int re=1;
    for(;y;y>>=1,x=(ll)x*x%mod)
     if(y&1) re=(ll)re*x%mod;
    return re;
}

int main(){
    int i,j;
    scanf("%d",&K);n=1<<K;
    for(i=0;i<n;++i) scanf("%d",&a[i]);
    for(i=0;i<n;++i) scanf("%d",&b[i]);
    
    for(i=0;i<n;++i) f[i]=a[i],g[i]=b[i];
    fwt_or(f,1);fwt_or(g,1);
    for(i=0;i<n;++i) c[i]=(ll)f[i]*g[i]%mod;
    fwt_or(c,0);
    for(i=0;i<n;++i) {printf("%d",c[i]);putchar(i+1<n?' ':'\n');}
    
    for(i=0;i<n;++i) f[i]=a[i],g[i]=b[i];
    fwt_and(f,1);fwt_and(g,1);
    for(i=0;i<n;++i) c[i]=(ll)f[i]*g[i]%mod;
    fwt_and(c,0);
    for(i=0;i<n;++i) {printf("%d",c[i]);putchar(i+1<n?' ':'\n');}
    
    for(i=0;i<n;++i) f[i]=a[i],g[i]=b[i];
    fwt_xor(f,1);fwt_xor(g,1);
    for(i=0;i<n;++i) c[i]=(ll)f[i]*g[i]%mod;
    fwt_xor(c,0);
    for(i=0;i<n;++i) {printf("%d",c[i]);putchar(i+1<n?' ':'\n');}
    
    return 0;
}