学习笔记第三十九节:快速沃尔什变换(FWT)

正题

      我们学会了快速傅立叶变换,狄利克雷卷积,现在我们来解决一下快速沃尔什变换。

      它主要是解决这样的问题的:C(k)=\sum_{i\oplus j=k} A(i)*B(j)

      也就是说C=A\times B

      其中那个你看不懂的符号指的是三则运算的其中一种:or,and,xor。

      怎么做?

      我们来模仿一下解决多项式乘法的过程,我们先是把原数组做一次傅立叶变换,然后再相乘,最后再做一次逆变换。

      是否可以把这种方法运用到这里呢?

      显然是可以的:我们设一个数组tf(A),它表示A数组FWT变换之后的结果。

      我们要使得tf(C)=tf(A)*tf(B),其中乘法是按位相乘。

      然后我们做一次tf(C)逆变换即可。

      那么这个tf(C)应该怎么找呢?

      三种运算就有三种tf

  or

       我们设tf(A)=(tf(A_0),tf(A_0)+tf(A_1))

       其中+表示按位相加,二元组表示前后相接,也就是说A=(A_0,A_1)其中0,1表示二进制下标的第一位恰好分为两部分。

       发现,当我们使得tf这样定义的时候,tf(C)=tf(A)*tf(B)

       为什么?

       首先来看CC=A\times B,明显的当我们把A,B拆开时,贡献到C_0的只有A_0*B_0,而贡献到C_1的就是A_0*B_1+A_1*B_0+A_1*B_1。也就是说C=A\times B=(A_0,A_1)\times (B_0,B_1)=(A_0*B_0,A_0*B_1+A_1*B_0+A_1*B_1)。(因为是或运算啊

       接着我们就可以知道

      \\tf(C)=(tf(C_0),tf(C_0+C_1)) \\=(tf(A_0*B_0),tf(A_0 * B_0+A_0* B_1+A_1* B_0+A_1* B_1))) \\=(tf(A_0*B_0),tf(A_0* B_0)+tf(A_0*B_1)+tf(A_1*B_0)+tf(A_1*B_1))

      为什么可以拆?因为考虑每个多个数组 加起来和分开 并不影响 每个位置上面的数对后面的影响,具体可以试一下手算。
      再看tf(A)*tf(B),这个按位相乘又等于什么呢?相当于前面和前面的相乘,后面的和后面的相乘。

      \small \\tf(A) * tf(B)=(tf(A_0),tf(A_0)+tf(A_1))*(tf(B_0),tf(B_0)+tf(B_1))) \\=(tf(A_0)*tf(B_0),tf(A_0)*tf(B_0)+tf(A_0)*tf(B_1)+tf(A_1)*tf(B_0+tf(A_1)*tf(B_1)))

      我们知道当只有一个数的时候,tf(C)=tf(A)*tf(B)

      然后我们就可以知道,只要我们能满足当

       \\tf(A_0\times B_0)=tf(A_0)*tf(B_0) \\tf(A_0\times B_1)=tf(A_0)*tf(B_1) \\tf(A_1\times B_0)=tf(A_1)*tf(B_0) \\tf(A_1\times B_1)=tf(A_1)*tf(B_1)

       时,那么我们就可以证明tf(A\times B)=tf(A)*tf(B)

       把上面的四条式子换到上面tf(C)的式子中,可以发现上面所化简的tf(C=A\times B)tf(A)*tf(B)是相等的。

       所以我们就证明了当tf(A)=(tf(A_0),tf(A_0)+tf(A_1))时,tf(C=A\times B)=tf(A)*tf(B)是成立的。

        接下来我们就可以用O(n\log n)的时间处理出来tf(A),tf(B)再逐位相乘得到tf(C)

        要求C还不简单吗?

        考虑我们怎么从A变为tf(A)的,是不是先从小的开始,然后将左边的值加到右边。

        我们要变回去,就从最大的区间开始,用用右边的左边的即可。

    and

        那么与运算其实和或运算差不多。tf(A)=(tf(A_0)+tf(A_1),tf(A_1))

        求逆运算的过程类似。关于证明可以自己推导,十分简单。

    xor

         抑或运算有点复杂,但是有公式就很好推导tf(A)=(tf(A_0)+tf(A_1),tf(A_0)-tf(A_1))

         是不是很鬼畜,然而求逆也是非常简单的,自己推导。

     最后,求逆的证明有的博客写得很复杂,其实根本不需要,就相当于将加的减回去即可。

     

#include
#include
#include
#include
using namespace std;

const int maxn=17;
long long a[1<=2;l>>=1)
		for(int i=0;i=2;l>>=1)
		for(int i=0;i=2;l>>=1)
		for(int i=0;i

 

你可能感兴趣的:(学习笔记)