FWT求解的是一类问题:\( a[i] = \sum\limits_{j\bigoplus k=i}^{} b[j]*c[k] \)

其中,\( \bigoplus \) 可以是 or,and,xor

三种问题的解决思路都是对多项式 \( a \) 构造一个 \( a' \),令 \( a' = b' * c' \);

那么只需要把 \( b \) 变换成 \( b' \),\( c \) 变换成 \( c' \),然后乘出 \( a' \),再逆变换得到 \( a \);

下面问题就变成如何快速(logn)求 \( b \) 到 \( b' \) 的变换,这个变换就是 FWT;

始终要记住进行位运算的是位置(角标)而不是值;

一、or

构造 \( a'[i] = \sum\limits_{j|i=i}^{} a[j] \)

1.正变换

考虑把 \( a \) 分成前后两个部分 \( a0 \) 和 \( a1 \),先分别递归下去做好,得到 \( a0' \) 和 \( a1' \);

可以发现,\( a0' \) 和 \( a1' \) 的位置(角标)数字上唯一不同就是最高位是0或1;

但递归下去做的时候,\( a0' \) 和 \( a1' \) 的位置数字相当与去掉了最高位(因为折半了);

所以合并的时候,关键要考虑到最高位的0和1的不同:

(1) 对于 \( a' \) 的一个位置 \( i \) ,如果它在前半部分,那么它可以直接继承 \( a0'[i] \);

而 \( a1'[i]\) 由于实际上 \( i \) 还应该加上最高位的1,or 运算使它能贡献的位置最高位也是1,但 \( i \) 的最高位是0,所以不贡献给 \( a'[i] \) ;

(2)对于后半部分的 \( i \) ,\( a0'[i] \) 和 \( a1'[i] \) 都会对它产生贡献,因为两部分的位置数字都是 \( i \) 的子集;

所以可以得到:\( a' = \left ( a0' , a0'+a1' \right ) \)

递归的底层,只有一个元素的时候,\( a = a' \) ,于是我们可以递归做出正变换了;

当然,仿照 FFT 的写法即可,并不需要真的写递归函数,而且也不用蝴蝶变换;

2.逆变换

同样先考虑两个部分 \( a'0 \) 和 \( a'1 \) ,表示 \( a' \) 的前后部分;

已经做了 \( a' = \left ( a0' , a0'+a1' \right ) \)

现在要从 \( a' \) 拆出 \( a0' \) 和 \( a1' \)

那么 \( a0' = a'0 \)

而且 \( a1' = a'1 - a'0 \)

知道了 \( a0' \) 和 \( a0' \) ,就可以继续递归求解 \( a0 \) 和 \( a1 \),二者合起来就可以得到 \( a \)

递归的底层,只有一个元素的时候,\( a' = a \) ,于是我们可以递归作出逆变换了;

void fwt1(int *a,int tp)//a'=(a0',a0'+a1')  //a=(a0',a1'-a0')
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
a[j+mid+k]=upt(a[j+mid+k]+tp*a[j+k]);
}

or

二、and

构造 \( a' = \sum\limits_{j \& i=i}^{} a[j] \)

1.正变换

和 or 同理,考虑最高位01的不同,后面继承本身,而前面要加上后面的贡献;

得到 \( a' = \left ( a0'+a1' , a1' \right ) \)

2.逆变换

同理,得到

\( a0' = a'0 - a'1 \)

\( a1' = a'1 \)

void fwt2(int *a,int tp)//a'=(a0'+a1',a1')  //a=(a0'-a1',a1')
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
a[j+k]=upt(a[j+k]+tp*a[j+mid+k]);
}

and

三、xor

设 \( d(i,j) \) 表示 \( i\&j \) 二进制表示中1的个数;

构造 \( a' = \sum\limits_{d(i,j)\%2==0}^{} a[j] - \sum\limits_{d(i,j)\%2==1}^{} a[j] \)

1.正变换

让我们三步走:

(1) \( a' = \left ( a0' + a1' , a0' - a1' \right ) \)

首先明确,\( a' \) 是 \( d(i\&j) \) 为偶数的 \( a[j] \) 求和,减去 \( d(i\&j) \) 为奇数的 \( a[j] \) 求和;

<1> 对于整体的一个位置 \( i \),它在前半部分

对于前半部分(折半)的相同位置 \( i' \),在前半部分的 \( j \) 中,\( d(i'\&j) \) 的奇偶性和 \( d(i\&j) \) 一样,所以继承答案;

对于后半部分(折半)的相同位置 \( i' \),在后半部分的 \( j \) 中,计算 \( i'\&j \) 时是没有考虑最高位的,所以它们的最高位上都是0,

而因为 \( i \) 的最高位是0,\( i\&j \) 的最高位同样是0,所以正好符合,答案可以加上;

也就是,\( a' = \left ( a0' + a1' , ... \right ) \)

<2> 对于整体的一个位置 \( i \),它在后半部分

对于前半部分(折半)的相同位置 \( i' \),在前半部分的 \( j \) 中,\( d(i'\&j) \) 的最高位都是0,

而因为 \( j \) 的最高位是0,\( i\&j \) 的最高位同样是0,所以正好符合,答案可以加上;

对于后半部分(折半)的相同位置 \( i' \),在后半部分的 \( j \) 中,计算 \( i'\&j \) 时是没有考虑最高位的,所以它们的最高位上都是0,

但 \( i\&j \) 的最高位是1,所以奇偶性都反了,答案加上的是负的;

这样,就得到 \( a' = \left ( a0' + a1' , a0' - a1' \right ) \)

(2) \( d(i\&k) \otimes d(j\&k) = d( (i \otimes j)\&k ) \)

因为是 \( \& \) ,我们就看 \( k \) 是1的那些位;

如果 \( d(i\&k) \) 是偶数,说明 \( i\&k \) 有偶数个1和 \( k \) 重合,奇数同理,\( j \) 同理;

<1> 当 \( d(i\&k) \) 和 \( d(j\&k) \) 奇偶性相同时

\( d(i\&k) + d(j\&k) \) 是偶数;

而 \( i \otimes j \) 同时消去 \( i \) 和 \( j \) 相同位置的1,不是 \( k \) 的1就算了,是 \( k \) 的1,消去的也是偶数;

所以 \( d( (i \otimes j)\&k ) \) 是偶数;

<2> 当 \( d(i\&k) \) 和 \( d(j\&k) \) 奇偶性不同时

\( d(i\&k) + d(j\&k) \) 是奇数;

而 \( i \otimes j \) 同时消去 \( i \) 和 \( j \) 相同位置的1,不是 \( k \) 的1就算了,是 \( k \) 的1,消去的是偶数;

所以 \( d( (i \otimes j)\&k ) \) 是奇数;

这样我们就证明了 \( d(i\&k) \otimes d(j\&k) = d( (i \otimes j)\&k ) \)

(3) 若 \( c[i] = \sum_{j \otimes k=i}^{} a[j]*b[k] \) ,有 \( c' = a' * b' \)

因为 \( c[i] = \sum_{j \otimes k=i}^{} a[j]*b[k] \)

又 \( c' = \sum_{d(i,j)\%2==0}^{} c[j] - \sum_{d(i,j)\%2==1}^{} c[j] \)

代入,得到 \( c'[i] = \sum_{d((j \otimes k)\&i)\%2==0}^{} a[j]*b[k] - \sum_{d((j \otimes k)\&i)\%2==1}^{} a[j]*b[k] \)

而 \( a'[i] * b'[i] = ( \sum_{d(i,j)\%2==0}^{} a[j] - \sum_{d(i,j)\%2==1}^{} a[j]) * ( \sum_{d(i,j)\%2==0}^{} b[j] - \sum_{d(i,j)\%2==1}^{} b[j]) \)

拆开再组合,并使用(2)得到的结论,就得到 \( a'[i] * b'[i] = ( \sum_{d((j \otimes k)\&i)\%2==0}^{} a[j]*b[k] - \sum_{d((j \otimes k)\&i))\%2==1}^{} a[j]*b[k]) \)

所以 \( c' = a' * b' \)

综上,我们仍然可以递归求 xor 的正变换,\( a' = \left ( a0' + a1' , a0' - a1' \right ) \)

2.逆变换

根据正变换就可以知道咯:

\( a0' = (a'0 + a'1) / 2 \)

\( a1' = (a'0 - a'1) / 2 \)

void fwt3(int *a,int tp)//a'=(a0'+a1',a0'-a1')  //a=((a0'+a1')/2,(a0'-a1')/2)
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
{
int x=a[j+k],y=a[j+mid+k];
a[j+k]=upt(x+y); a[j+mid+k]=upt(x-y);
if(tp==-)a[j+k]=(ll)a[j+k]*inv%mod,a[j+mid+k]=(ll)a[j+mid+k]*inv%mod;
}
}

xor

看例题:https://www.luogu.org/problemnew/show/P4717

代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=(<<),mod=;
int n,a[xn],b[xn],f[xn],g[xn],lim,inv;
int rd()
{
int ret=,f=; char ch=getchar();
while(ch<''||ch>''){if(ch=='')f=; ch=getchar();}
while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
return f?ret:-ret;
}
ll pw(ll a,int b)
{
ll ret=;
for(;b;b>>=,a=(a*a)%mod)
if(b&)ret=(ret*a)%mod;
return ret;
}
int upt(int x){while(x>=mod)x-=mod; while(x<)x+=mod; return x;}
void fwt1(int *a,int tp)//a'=(a0',a0'+a1') //a=(a0',a1'-a0')
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
a[j+mid+k]=upt(a[j+mid+k]+tp*a[j+k]);
}
void fwt2(int *a,int tp)//a'=(a0'+a1',a1') //a=(a0'-a1',a1')
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
a[j+k]=upt(a[j+k]+tp*a[j+mid+k]);
}
void fwt3(int *a,int tp)//a'=(a0'+a1',a0'-a1') //a=((a0'+a1')/2,(a0'-a1')/2)
{
for(int mid=;mid<lim;mid<<=)
for(int j=,len=(mid<<);j<lim;j+=len)
for(int k=;k<mid;k++)
{
int x=a[j+k],y=a[j+mid+k];
a[j+k]=upt(x+y); a[j+mid+k]=upt(x-y);
if(tp==-)a[j+k]=(ll)a[j+k]*inv%mod,a[j+mid+k]=(ll)a[j+mid+k]*inv%mod;
}
}
int main()
{
n=rd(); lim=(<<n); inv=pw(,mod-);
for(int i=;i<lim;i++)a[i]=f[i]=rd();
for(int i=;i<lim;i++)b[i]=g[i]=rd();
fwt1(f,); fwt1(g,);
for(int i=;i<lim;i++)f[i]=(ll)f[i]*g[i]%mod;
fwt1(f,-);
for(int i=;i<lim;i++)printf("%d ",f[i]); puts(""); for(int i=;i<lim;i++)f[i]=a[i],g[i]=b[i];
fwt2(f,); fwt2(g,);
for(int i=;i<lim;i++)f[i]=(ll)f[i]*g[i]%mod;
fwt2(f,-);
for(int i=;i<lim;i++)printf("%d ",f[i]); puts(""); for(int i=;i<lim;i++)f[i]=a[i],g[i]=b[i];
fwt3(f,); fwt3(g,);
for(int i=;i<lim;i++)f[i]=(ll)f[i]*g[i]%mod;
fwt3(f,-);
for(int i=;i<lim;i++)printf("%d ",f[i]); puts("");
return ;
}

参考博客:https://www.cnblogs.com/ACMLCZH/p/8022502.html

https://blog.csdn.net/neither_nor/article/details/60335099

快速沃尔什变换(FWT)学习笔记 + 洛谷P4717 [模板]的更多相关文章

  1. 一个数学不好的菜鸡的快速沃尔什变换(FWT)学习笔记

    一个数学不好的菜鸡的快速沃尔什变换(FWT)学习笔记 曾经某个下午我以为我会了FWT,结果现在一丁点也想不起来了--看来"学"完新东西不经常做题不写博客,就白学了 = = 我没啥智 ...

  2. 快速沃尔什变换 (FWT)学习笔记

    证明均来自xht37 的洛谷博客 作用 在 \(OI\) 中,\(FWT\) 是用于解决对下标进行位运算卷积问题的方法. \(c_{i}=\sum_{i=j \oplus k} a_{j} b_{k} ...

  3. 快速沃尔什变换 FWT 学习笔记【多项式】

    〇.前言 之前看到异或就担心是 FWT,然后才开始想别的. 这次学了 FWT 以后,以后判断应该就很快了吧? 参考资料 FWT 详解 知识点 by neither_nor 集训队论文 2015 集合幂 ...

  4. 斜率优化dp学习笔记 洛谷P3915[HNOI2008]玩具装箱toy

    本文为原创??? 作者写这篇文章的时候刚刚初一毕业…… 如有错误请各位大佬指正 从例题入手 洛谷P3915[HNOI2008]玩具装箱toy Step0:读题 Q:暴力? 如果您学习过dp 不难推出d ...

  5. 边带权并查集 学习笔记 & 洛谷P1196 [NOI2002] 银河英雄传说 题解

    花了2h总算把边带权并查集整明白了qaq 1.边带权并查集的用途 众所周知,并查集擅长维护与可传递关系有关的信息.然而我们有时会发现并查集所维护的信息不够用,这时"边带权并查集"就 ...

  6. STL Stack(栈)学习笔记 + 洛谷 P1449 后缀表达式

    稍微看了看刘汝佳的白皮书,“实用主义”的STL实在是香到我了,而且在实验室大佬的推荐下我开始了stl的学习. 每篇附带一个题目方便理解,那行,直接开始. 毕竟是实用主义,所以就按照给的题目的例子来理解 ...

  7. STL Queue(队列)学习笔记 + 洛谷 P1540 机器翻译

    队(Queue) 队简单来说就是一个先进先出的“栈”,但是不同于标准“栈”的先进后出. 基本操作: push(x) 将x压入队列的末端 pop() 弹出队列的第一个元素(队顶元素),注意此函数并不返回 ...

  8. 【学习笔鸡】快速沃尔什变换FWT

    [学习笔鸡]快速沃尔什变换FWT OR的FWT 快速解决: \[ C[i]=\sum_{j|k=i} A[j]B[k] \] FWT使得我们 \[ FWT(C)=FWT(A)*FWT(B) \] 其中 ...

  9. FMT/FWT学习笔记

    目录 FMT/FWT学习笔记 FMT 快速莫比乌斯变换 OR卷积 AND卷积 快速沃尔什变换(FWT/XOR卷积) FMT/FWT学习笔记 FMT/FWT是算法竞赛中求or/and/xor卷积的算法, ...

随机推荐

  1. JobConf

    /**  * A map/reduce job configuration. * 翻译:一个map/reduce作业配置 * <p><code>JobConf</code ...

  2. Name和:Name

    http://www.cnblogs.com/tianyutingxy/archive/2011/02/23/1962078.html http://www.tuicool.com/articles/ ...

  3. VMware虚拟机下如何安装一个64位的win7系统

    原文地址:http://www.xitongcheng.com/jiaocheng/win7_article_21001.html VMware虚拟机软件可以在一台电脑上运行多个操作系统,一些网友想在 ...

  4. Ubuntu下安装libsvm

    在安装LibSVM前需要先装 python 和 gnuplot linux 一般都自带了python2.7,所以python的安装不再赘述 在 ubuntu 下安装 gnuplot 不能直接 sudo ...

  5. Perl图书的一些体会

    近期,由于项目须要.又又一次将Perl学习起来. Perl老实说.让我又爱又恨. 爱它.是由于自己写代码的确非常爽. 是代码最少.速度最快的语言. 恨是由于看别人的代码实在太累了. 但,整体体会,在文 ...

  6. Django-缓存的配置

    缓存的介绍 在动态网站中,用户所有的请求,服务器都会去数据库中进行相应的增,删,查,改,渲染模板,执行业务逻辑,最后生成用户看到的页面. 当一个网站的用户访问量很大的时候,每一次的的后台操作,都会消耗 ...

  7. smartUpload注意过程

    操作的过程中一定要注意的几个方面:       1.将smartUpload.jar拷贝到tomcat/lib以及项目的lib下面,最好是只多不少!       2.因为上传的文件一般都很大,所以应该 ...

  8. CSS入门学习

    一.What? CSS的全称是CascadingStyle Sheet,汉语意思是"级联样式表".通常又称为"风格样式表(StyleSheet)".它是用来进行 ...

  9. Python中十六进制和字符串的转换(转载)

    调用Python内置int()函数把该字串转为数字.以下为在Python解释器编程环境下的操作示范: 把十六进制的字串转为十进制数字:Python代码>>> print int('f ...

  10. Open Source Streaming Server--EasyDarwin

    Welcome to EasyDarwin Streaming Server, which is an open source Streaming Server Based On Appple's D ...