FWT快速沃尔什变换例题

最新推荐文章于 2023-03-25 11:19:43 发布
最新推荐文章于 2023-03-25 11:19:43 发布
阅读量172 收藏
点赞数
原文链接：http://www.cnblogs.com/PaperCloud/p/10027071.html
版权

模板题


传送门 


#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define swap(a,b) (a^=b^=a^=b)
inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}

#define mod 998244353
#define inv2 499122177

ll a[1<<17],b[1<<17],c[1<<17],N;

inline void FWT_or(ll *a,int f)
{
    register int i,j,p,k;
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k)
        a[i+j+k]=(a[i+j+k]+a[j+k]*f+mod)%mod;
}

inline void FWT_and(ll *a,int f)
{
    register int i,j,p,k;
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k)
        a[j+k]=(a[j+k]+a[i+j+k]*f+mod)%mod;
}

inline void FWT_xor(ll *a,bool f)
{
    register int i,j,p,k;
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k)
    {
        int X=a[j+k],Y=a[i+j+k];
        a[j+k]=(X+Y)%mod;a[i+j+k]=(X+mod-Y)%mod;
        if(f==0)a[j+k]=1ll*a[j+k]*inv2%mod,a[i+j+k]=1ll*a[i+j+k]*inv2%mod;
    }
}

int main()
{
    N=1<<read();register int i;
    
    for(i=0;i<N;++i) a[i]=1ll*read()%mod;
    for(i=0;i<N;++i) b[i]=1ll*read()%mod;
    
    FWT_or(a,1);FWT_or(b,1);
    for(i=0;i<N;++i) c[i]=1ll*a[i]*b[i]%mod;
    FWT_or(a,-1);FWT_or(b,-1);FWT_or(c,-1);
    for(i=0;i<N;++i) printf("%lld ",c[i]);puts("");
    
    FWT_and(a,1);FWT_and(b,1);
    for(i=0;i<N;++i) c[i]=1ll*a[i]*b[i]%mod;
    FWT_and(a,-1);FWT_and(b,-1);FWT_and(c,-1);
    for(i=0;i<N;++i) printf("%lld ",c[i]);puts("");
    
    FWT_xor(a,1);FWT_xor(b,1);
    for(i=0;i<N;++i) c[i]=1ll*a[i]*b[i]%mod;
    FWT_xor(a,0);FWT_xor(b,0);FWT_xor(c,0);
    for(i=0;i<N;++i) printf("%lld ",c[i]);puts("");
    
    return 0;
}




（bzoj 4589）Hard Nim


传送门

题意：n堆石子,每堆石子的数量都是不大于m的质数的Nim游戏，求先手必败的局面数量模1000000007


Solution

  • 用线筛求出不大于m的质数
  • \(是质数a[i]=[i是质数]\)
  • 直接对FWT(a)进行快速幂
  • 答案就是a[0]



#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define MN 65540
int prime[MN],tot;
bool mark[MN];
#define mod 1000000007
int a[MN],inv2,N;
inline int fpow(int x,int m)
{
    int res=1;
    while(m)
    {
        if(m&1) res=(int)(1ll*res*x%mod);
        x=(int)(1ll*x*x%mod);m>>=1;
    }
    return res;
}
inline void init()
{
    inv2=fpow(2,mod-2);
    mark[1]=1;register int i,j;
    for(i=2;i<=MN-5;i++)
    {
        if(!mark[i]){prime[++tot]=i;}
        for(j=1;j<=tot&&prime[j]*i<=MN-5;++j)
        {
            mark[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void FWT_xor(int *a,bool f)
{
    register int i,j,p,k;
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k)
    {
        int X=a[j+k],Y=a[i+j+k];
        a[j+k]=(X+Y)%mod;a[i+j+k]=(X+mod-Y)%mod;
        if(f==0)a[j+k]=1ll*a[j+k]*inv2%mod,a[i+j+k]=1ll*a[i+j+k]*inv2%mod;
    }
}
int main()
{
    register int n,m,i,j;
    init();
    while(~scanf("%d%d",&n,&m))
    {
        for(N=1;N<=m;N<<=1);
        memset(a,0,sizeof a);
        for(i=1;i<=m;++i) a[i]=(!mark[i]);
        FWT_xor(a,1);
        for(i=0;i<N;++i) a[i]=fpow(a[i],n);
        FWT_xor(a,0);
        printf("%d\n",a[0]);
    }
    return 0;
}



（hdu 5909）Tree Cutting


传送门

题意：给出一棵树，求异或和为[0, m-1]的非空连通子图的个数


Solution

  • 树形dp
  • f[i][j]表示以i为根的子树内，异或和为j的联通块个数
  • 用FWT优化转移



#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define swap(a,b) (a^=b^=a^=b)
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
#define MN 2005
#define mod 1000000007
#define inv2 500000004
int val[MN],en,hr[MN];
struct edge{int to,nex;}e[MN<<1];
int N,f[MN][1<<11],tmp[1<<11],ans[1<<11];
inline void ins(int f,int t)
{
    e[++en]=(edge){t,hr[f]};hr[f]=en;
    e[++en]=(edge){f,hr[t]};hr[t]=en;
}
inline void FWT_xor(int *a,bool f)
{
    register int i,j,p,k;
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k)
    {
        int X=a[j+k],Y=a[j+i+k];
        a[j+k]=(X+Y)%mod,a[j+i+k]=(X+mod-Y)%mod;
        if(f==0) a[j+k]=1ll*a[j+k]*inv2%mod,a[j+i+k]=1ll*a[j+i+k]*inv2%mod;
    }
}
inline void rw(int *a,int *b)
{
    FWT_xor(a,1);FWT_xor(b,1);register int i;
    for(i=0;i<N;++i) a[i]=1ll*a[i]*b[i]%mod;
    FWT_xor(a,0);
}
inline void dfs(int x,int fa)
{
    f[x][val[x]]=1;
    register int i,j;
    for(i=hr[x];i;i=e[i].nex)if(e[i].to^fa)
    {
        dfs(e[i].to,x);
        for(j=0;j<N;++j) tmp[j]=f[x][j];
        rw(f[x],f[e[i].to]);
        for(j=0;j<N;++j) (f[x][j]+=tmp[j])%=mod;
    }
    for(i=0;i<N;++i) (ans[i]+=f[x][i])%=mod;
}
int main()
{
    register int Case=read(),n,m,i;
    while(Case--)
    {
        n=read(),m=read();N=m;
    //  for(N=1;N<=m;N<<=1);
        en=0;
        memset(hr,0,sizeof hr);
        memset(ans,0,sizeof ans);
        memset(f,0,sizeof f);
        for(i=1;i<=n;++i) val[i]=read();
        for(i=1;i<n;++i) ins(read(),read());
        dfs(1,0);
        for(i=0;i<m-1;i++) printf("%lld ",ans[i]);
        printf("%lld\n",ans[m-1]);
    }
    return 0;
}

（bzoj 4036）[HAOI2015]按位或

传送门

Description

刚开始你有一个数字0，每一秒钟你会随机选择一个[0,2^n-1]的数字，与你手上的数字进行或（c++,c的|,pascal的or）操作。

选择数字i的概率是p[i]。保证0<=p[i]<=1，Σp[i]=1问期望多少秒后，你手上的数字变成2^n-1。

Solution

相当于给你一个集合求最后一个元素出现的时间

可以用\(minmax\) 容斥一波，这样就是求每个子集中第一个元素出现的时间\(min(T)\)啦

我们设\(P(T)\)表示取到\(T\)的子集的概率
\[ P(min(T)==k)=P(S-T)^{k-1}(1-P(S-T)) \]
然后因为：
\[ 若P(x==k)=(1-p)^{k-1}p(k \in N^{+})，则E(x)=\frac{1}{p} \]
所以：
\[ E(min(T))=\frac{1}{1-P(S-T)} \]
问题在于，如何求\(P(T)\)呢

显然
\[ P(T)=\sum_{x⊆T}p(x)，这里我们把一个数都当作一个集合，p(x)就是题目给出的得到x的概率 \]
求子集和？可以用像「PKUWC2018」随机游走 一样的子集和dp，当然，也可以直接\(FWT\)变换一下。


Code 

#include<bits/stdc++.h>
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
int n,N,num[1<<20];
double P[1<<20],Ans=0.;
int main()
{
    scanf("%d",&n);N=1<<n;
    register int i,j,p,k;
    for(i=0;i<N;++i)scanf("%lf",&P[i]);
    for(i=1;i<N;i<<=1)for(p=i<<1,j=0;j<N;j+=p)for(k=0;k<i;++k) P[i+j+k]+=P[j+k];
    for(i=1;i<N;++i)
    {
        num[i]=num[i>>1]+(i&1);
        if(1-P[(N-1)^i]<1e-9) return 0*puts("INF");
        Ans+=((num[i]&1)?1.:-1.)*(1./(1.-P[(N-1)^i]));
    }
    printf("%.9lf\n",Ans);
    return 0;
}




Blog来自PaperCloud，未经允许，请勿转载，TKS！

转载于:https://www.cnblogs.com/PaperCloud/p/10027071.html

weixin_33811539
关注
解题报告（一）快速沃尔什变换FWT（ACM / OI）超高质量题解
繁凡さん的博客
02-23 1506
解题报告（一）快速沃尔什变换FWT（ACM / OI）超高质量题解
快速沃尔什变换（FWT）介绍
HolmiumJiang的博客
07-24 4793
快速沃尔什变换（FWT）介绍 众所周知，FFT可以有效地解决如下的问题 A=(a1,…,an)B=(b1,...,bn)C=(c1,...,cn)ck=∑i+j=kai∗bj A=(a_1,\dots,a_n)\\ B=(b_1,...,b_n)\\ C=(c_1,...,c_n)\\ c_k=\sum_{i+j=k}a_i*b_j\\ A=(a1​,…,an​)B=(b1​,...,bn​)C=(c1​,...,cn​)ck​=i+j=k∑​ai​∗bj​ 但事实上，把求和号里的条件里的“+”换成与、或、
解题报告：CodeForces - 662C:Binary Table FWT(快速沃尔什变换)
Hee的博客~
08-03 806
题目链接 题意： 给定一个至多20*100000的01矩阵，可以选择任意行或列，选择行或列的01值改变，问经过操作能到最少的含1的数量。 思路： 因为行最大为20，考虑将每一列压缩成一个20位的整数，选择的行也最多只有20，同样压缩到一个20位的整数。 改变行上的01值，对应的操作为异或，那么就能写出对应的卷积式： 根据异或运算的性质，可以交换j,k的位置，得到： C[k
FWT(快速沃尔什变换) 试题总结
海边拾贝，沧海一粟
08-21 2927
FWT FWT学习 参考博客 FWT 题目总计 bzoj4589: Hard Nim FWT 模板题 参考代码 E Touring cities 根据线性基的概念推出最多去掉log(sumxor) 个，进行一次FWT然后相乘，求IFWT判断sumxor 是否可以推出 参考代码 3. ...
CSU 1911 Card Game 快速沃尔什变换（FWT）模板题
PrimeXuan
04-29 1346
题意：给你两个二进制数的集合，给出q次询问，输出两个集合之间元素或的值等于查询值的种数。   #include #include #include using namespace std; #define LL long long const int maxn = 1 LL a[maxn],b[maxn]; int n,m; char s[20]; void FWT(LL a
FWT(快速沃尔什变换)零基础详解qaq（ACM/OI）
qq_37136305的博客
08-12 3717
1.前言（废话） 记得一年半之前做SRM518 Nim的时候还不知道FWT，当时自己用分治完美的水过去了。然后昨天的牛客有一道题，是说nim博弈中有n堆石子，请问最多取出多少堆石子可以让先手必败。当时竟然没思路QAQ？？？想了想使劲往字典树靠边靠不上去QAQ，然后就没想出来！！！ 想当年自己手推FWT或运算，真的感叹岁月是把杀猪刀！于是怒写这篇博客QAQ把这个算法总结一下QAQ （其实我觉得...
To_Heart—总结——FWT(快速沃尔什变换)
最新发布
xf2056188203的博客
03-25 448
这个比FFT简单了很多呢，，大概是我可以学懂的水平！好像是叫 快速沃尔什变换？
利用快速沃尔什变换算法对数字图像或矩阵作快速变换
08-01
快速沃尔什变换（Fast Walsh Transform, FWT）是数字信号处理和图像处理领域中的一种重要算法，它在处理离散数据时提供了高效、简洁的方法。本文将深入探讨快速沃尔什变换的基本原理、应用以及实现过程。 沃尔什...
快速沃尔什变换 FWT 学习笔记【多项式】
weixin_30487701的博客
04-29 260
〇、前言 之前看到异或就担心是 FWT，然后才开始想别的。 这次学了 FWT 以后，以后判断应该就很快了吧？ 参考资料 FWT 详解 知识点 by neither_nor 集训队论文 2015 集合幂级数的性质与应用及其快速算法 by 吕凯风 一、FWT 是什么 FWT 是快速沃尔什变换。它和快速傅里叶变换一样，原本都用于物理中的频谱分析。 但是由于它可分治的特点，在算法竞赛中常被用来计算位运...
模板：快速莫比乌斯变换（FMT）+快速沃尔什变换（FWT）（多项式）
BUG_Creater_jie的博客
01-04 1400
下标进行位运算而不是简单加法的卷积时的处理方法
【NOI模拟赛】废物题（快速沃尔什变换，分治NTT）
DD(XYX)的硫氰化氘博客
05-02 398
【图片无法显示】
牛客练习赛41-简单的数学题 快速沃尔什变换
Monster_ixx的博客
03-02 355
题目链接:https://ac.nowcoder.com/acm/contest/373/F 乍一看这个题貌似是个数学题,其实仔细一看,就..嘿嘿 其实题目中的函数都可以进行相应的转换f(x)函数在大于等于2时可以理解为x的所有质因子的乘积,因为题目中说明了最大质因子不超过71,也就是第20个质数。之后我们观察发现g函数和异或函数有些类似，因此我们可以想到可以把每个原数组中的数转化为一个不超过...
【FWT快速沃尔什变换】讲解 + 例题 BZO-J4589 Hard Nim（FWT、博弈论、快速幂）
Floraqiu的博客
01-24 719
参考：FWT快速沃尔什变换学习笔记(证明等都非常详细，虽然我没看证明 ) 一、FWT的作用 FFT 、NTT：计算多项式卷积 Ck=∑i+j=kAi∗BjC_k=\sum_{i+j=k}A_i*B_jCk​=i+j=k∑​Ai​∗Bj​ FWT ：计算多项式的位运算卷积 Ck=∑i∣j=kAi∗BjC_k=\sum_{i|j=k}A_i*B_jCk​=∑i∣j=k​Ai​∗Bj​ Ck=...
FWT快速沃尔什变换及其应用
平凡人Kalzn的博客
08-31 290
有要求我们对于两个序列，求解下列特殊卷积: 1、按位或卷积 Ci=∑j∣k=iAjBk(1)C_i=\sum_{j|k=i}A_jB_k\tag1Ci​=j∣k=i∑​Aj​Bk​(1) 2、按位与 Ci=∑j&k=iAjBk(2)C_i=\sum_{j\&k=i}A_jB_k\tag2Ci​=j&k=i∑​Aj​Bk​(2) 3、异或卷积 Ci=∑j⊕k=iAjBk(3)C_i=\sum_{j\oplus k=i}A_jB_k\tag3Ci​=j⊕k=i∑​Aj​Bk​(3) 三者
[FWT] 快速沃尔什变换学习笔记
HeRaNO@HSAHRBNU
05-05 1917
mdzz终于补坑了…… FWT是一种类似于FFT的变换。我们都知道对于数组aa与数组bb的卷积cc可以表示为： cj=∑i=1jaibj−i\begin{equation} c_j=\sum\limits_{i=1}^ja_ib_{j-i} \end{equation} 那么FWT只是把计算下标的加减运算改为了逻辑运算（与、或、异或等）。 比如与运算下的FWT表示为： ck=∑i&j=ka
快速沃尔什变换FWT
miaom的博客
10-17 962
大概在做CC月赛的时候看到了这样一道题：https://www.codechef.com/OCT17/problems/XORTREEH 题意是要你做一个%330301441的类似FWT的东西。定义数组A，B，C。定义一个操作A⊕B=CA\oplus B=C C[i]=∑u∑vA[u]B[v]  |  u,v在k进制下不进位的和为i  其中k=2,3,4,...,9,10 C[i]=\sum_u
快速沃尔什变换（FWT）
Never give in.
10-24 1165
Preface 我们知道，对于下面的式子 f(n)=∑i⊕j=ng(i)h(j)f(n)=\sum\limits_{i⊕j=n}g(i)h(j)f(n)=i⊕j=n∑​g(i)h(j) 当⊕⊕⊕为+运算时，可以用FFT优化 当⊕⊕⊕为位运算(or,and,xor)时是否也有一种优化呢 答案是肯定的。 FWT 首先，对于多项式，我们定义A∗BA*BA∗B为某一种位运算卷积 显然它满足交换律、结合律 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值