文章标题 FWT （快速沃尔什变换）模板

最新推荐文章于 2023-03-25 11:19:43 发布

原创最新推荐文章于 2023-03-25 11:19:43 发布 · 373 阅读

0 ·

CC 4.0 BY-SA版权

模板同时被 2 个专栏收录

27 篇文章

订阅专栏

FWT

2 篇文章

订阅专栏

本文详细介绍了快速沃尔特变换(FWT)及逆变换(UFWT)的实现原理与过程，并通过具体的代码示例展示了如何利用FWT进行数组操作，包括乘法运算。此变换对于处理布尔函数和组合数学问题具有重要意义。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

void FWT(int a[],int n)
{
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;j++)
            {
                int x=a[i+j],y=a[i+j+d];
                a[i+j]=(x+y)%mod,a[i+j+d]=(x-y+mod)%mod;
                //xor:a[i+j]=x+y,a[i+j+d]=(x-y+mod)%mod;
                //and:a[i+j]=x+y;
                //or:a[i+j+d]=x+y;
            }
}

void UFWT(int a[],int n)
{
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;j++)
            {
                int x=a[i+j],y=a[i+j+d];
                a[i+j]=1LL*(x+y)*rev%mod,a[i+j+d]=(1LL*(x-y)*rev%mod+mod)%mod;
                //xor:a[i+j]=(x+y)/2,a[i+j+d]=(x-y)/2;
                //and:a[i+j]=x-y;
                //or:a[i+j+d]=y-x;
            }
}
void solve(int a[],int b[],int n)
{
    FWT(a,n);
    FWT(b,n);
    for(int i=0;i<n;i++) a[i]=1LL*a[i]*b[i]%mod;
    UFWT(a,n);
}