O(n log log n)实现FGT和FLT(Fast GCD/LCM Transformation)

本文是作者看不懂分治FFT之后开始娱乐一下自己写的
看到一道题时候询问了正解后，推出了一个奇怪的变换，发现这个很Transformation，我和正解推出来的奇怪的东西是一样的，但还是想写一下思路。。。
(希望原题作者不要来D我啊）

min卷积和max卷积

考虑两个卷积
(C_{min(i,j)} = sum_{i = 1}^{n}sum_{j = 1}^{n} A_{i}B_{j})
(C_{max(i,j)} = sum_{i = 1}^{n}sum_{j = 1}^{n} A_{i}B_{j})
普通的卷积，即我们对下标进行加法操作后，将两个数的乘积加到该下标，在这里对下标的操作即使取min和取max

和暴力卷积一样，这两个卷积是可以暴力(n^2)出解的，但是有没有更优秀的做法呢

一种变换的出现，往往都是要把卷积变成点积，再用点积反解多项式的系数
例如FFT是插入函数值，利用两个函数相乘的结果，某个x坐标上的值即对应两个函数值的乘积，再利用单位根的特殊性质将系数解出来
例如FWT的异或卷积，我们要达成的成就就是(tf(A)* tf(B) = tf(C))，按最高位是0或1来分组，那么就是(C = (A_{0}*B_{0} + A_{1}*B_{1},A_{0} * B_{1} + A_{1} * B_{0}))，按照点积等于卷积后的答案，构造两个乘法使得可以得到其中的每一项，即为((A_{0} + A_{1})(B_{0} + B_{1}) = A_{0}B_{0} + A_{0}B_{1} + A_{1}B_{0} + A_{1}B_{1})和((A_{0} - A_{1})(B_{0} - B_{1}) = A_{0}B_{0} - A_{1}B_{0} - A_{0}B_{1} + A_{1}B_{1})
两项相加除二即为前半段，相减除二即为后半段，两个变换即为((tf(A_{0}) + tf(A_{1}),tf(A_{0}) - tf(A_{1})))复杂度(O(n log n))由于(n)往往是(2^{k})所以复杂度会被写成(O(2^{k}k))
所以我们尝试把这个卷积进行同样的操作
（然而很多变换的方式都是不知道为啥试出来的）

对于min卷积
我们考虑变换(C(i))为(SC(i)),(SC(i))表示(x>=i)的(C(x))的和，(SA(i))与(SB(i))同理
然后即可发现(SA(i) * SB(i) = SC(i))，最后差分即可还原C数组，复杂度做到了(O(n))

对于max卷积
我们考虑变换(C(i))为(SC(i)),(SC(i))表示(x<=i)的(C(x))的和，(SA(i))与(SB(i))同理
然后即可发现(SA(i) * SB(i) = SC(i))，最后差分即可还原C数组，复杂度做到了(O(n))

数组的标准分解意义下多维后缀/前缀和

我们如果把一个数的标准分解(即质因数分解)作为这个数的坐标，例如
(p_{1}^{k_{1}}p_{2}^{k_2}p_{3}^{k_{3}}...p_{n}^{k_{n}})
每个数即可看做一个无穷维度里，每一个维度的下标为该质因数的指数的一个坐标

数组的多维前缀和即
(f(n) = sum_{d|n}a(d))
数组的多维后缀和即
(f(n) = sum_{n|d}a(d))
朴素的做法可以达到(O(n log n))，可以过很多题目了，即对于一个数枚举它的倍数（用调和级数算复杂度也是一个重要的技巧）

我们考虑如何将这个算法优化到(O(n log log n))

我们以前缀和为例，从最小的质因子2开始，固定其他所有质数的指数，开始考虑以下这个过程
1:a(1)
2:a(2)
3:a(3)
4:a(4)
5:a(5)
6:a(6)
7:a(7)
8:a(8)
9:a(9)
10:a(10)
初始化即每个数组位置存储对应的a(i)
第一次操作2,从后往前扫，我们对每个2的倍数进行累加，具体看下面
1:a(1)
2:a(2) a(1)
3:a(3)
4:a(4) a(2) a(1)
5:a(5)
6:a(6) a(3)
7:a(7)
8:a(8) a(4) a(2) a(1)
9:a(9)
10:a(10) a(5)

第二次操作3
1:a(1)
2:a(2) a(1)
3:a(3) a(1)
4:a(4) a(2) a(1)
5:a(5)
6:a(6) a(3) a(2) a(1)
7:a(7)
8:a(8) a(4) a(2) a(1)
9:a(9) a(3) a(1)
10:a(10) a(5)

第三次操作5
1:a(1)
2:a(2) a(1)
3:a(3) a(1)
4:a(4) a(2) a(1)
5:a(5) a(1)
6:a(6) a(3) a(2) a(1)
7:a(7)
8:a(8) a(4) a(2) a(1)
9:a(9) a(3) a(1)
10:a(10) a(5) a(2) a(1)

第四次操作7
1:a(1)
2:a(2) a(1)
3:a(3) a(1)
4:a(4) a(2) a(1)
5:a(5) a(1)
6:a(6) a(3) a(2) a(1)
7:a(7) a(1)
8:a(8) a(4) a(2) a(1)
9:a(9) a(3) a(1)
10:a(10) a(5) a(2) a(1)

下一个质数是11，已经超出了范围，所以操作完成了，至此，我们已经完成了数组的多维前缀和

下面以相同的例子理解数组的多维后缀和
第一次操作2，从后往前扫
1:a(1) a(2) a(4) a(8)
2:a(2) a(4) a(8)
3:a(3) a(6)
4:a(4) a(8)
5:a(5) a(10)
6:a(6)
7:a(7)
8:a(8)
9:a(9)
10:a(10)

第二次操作3
1:a(1) a(2) a(4) a(8) a(3) a(6) a(9)
2:a(2) a(4) a(8) a(6)
3:a(3) a(6) a(9)
4:a(4) a(8)
5:a(5) a(10)
6:a(6)
7:a(7)
8:a(8)
9:a(9)
10:a(10)

第三次操作5
1:a(1) a(2) a(4) a(8) a(3) a(6) a(9) a(5) a(10)
2:a(2) a(4) a(8) a(6) a(10)
3:a(3) a(6) a(9)
4:a(4) a(8)
5:a(5) a(10)
6:a(6)
7:a(7)
8:a(8)
9:a(9)
10:a(10)

第四次操作7
1:a(1) a(2) a(4) a(8) a(3) a(6) a(9) a(5) a(10) a(7)
2:a(2) a(4) a(8) a(6) a(10)
3:a(3) a(6) a(9)
4:a(4) a(8)
5:a(5) a(10)
6:a(6)
7:a(7)
8:a(8)
9:a(9)
10:a(10)

由于对每个质数枚举倍数（根据埃拉托斯特尼筛法的复杂度）是(O(n log log n))的
所以求数组的标准分解意义下的多维前缀和，后缀和可以做到(O(n log log n))

Fast GCD Transformation 快速最大公约数变换

由于两个数的GCD，即为两个数标准分解下的每个质因数质数取min，即为
(p_{1}^{min(a_{1},b_{1}}p_{2}^{min(a_{2},b_{2})}...p_{k}^{min(a_{k},b_{k})})
类似min卷积，将数组处理成多维后缀和，即为正变换

考虑如何将变换后的结果还原回去
sc数组中点积相乘的结果即为
1:c(1) c(2) c(3) c(4) c(5) c(6) c(7) c(8) c(9) c(10)
2:c(2) c(4) c(6) c(8) c(10)
3:c(3) c(6) c(9)
4:c(4) c(8)
5:c(5) c(10)
6:c(6)
7:c(7)
8:c(8)
9:c(9)
10:c(10)
考虑如何将他们变换回去，我们发现可以采取我们累加后缀和的步骤反序进行

即，从大到小枚举质数，从前往后枚举下标

第一次枚举7
sc数组中点积相乘的结果即为
1:c(1) c(2) c(3) c(4) c(5) c(6) c(8) c(9) c(10)
2:c(2) c(4) c(6) c(8) c(10)
3:c(3) c(6) c(9)
4:c(4) c(8)
5:c(5) c(10)
6:c(6)
7:c(7)
8:c(8)
9:c(9)
10:c(10)

第二次枚举5
1:c(1) c(2) c(3) c(4) c(6) c(8) c(9)
2:c(2) c(4) c(6) c(8)
3:c(3) c(6) c(9)
4:c(4) c(8)
5:c(5) c(10)
6:c(6)
7:c(7)
8:c(8)
9:c(9)
10:c(10)

第三次枚举3
1:c(1) c(2) c(4) c(8)
2:c(2) c(4) c(8)
3:c(3) c(6)
4:c(4) c(8)
5:c(5) c(10)
6:c(6)
7:c(7)
8:c(8)
9:c(9)
10:c(10)

第四次枚举2
1:c(1)
2:c(2)
3:c(3)
4:c(4)
5:c(5)
6:c(6)
7:c(7)
8:c(8)
9:c(9)
10:c(10)

复杂度同样为(O(n log log n))，同时我们可以发现，数组标准分解意义下的多维后缀差分即使前缀和的逆步骤

同时，如果了解莫比乌斯反演，会列出等式
(sc(n) = sum_{n|d}c(d))
(c(n) = sum_{n|d}mu(frac{d}{n})sc(d))
即，该差分本质上是对(sc(i))进行莫比乌斯反演，我们又得到了一种(O(n log log n))求解莫比乌斯反演的办法，而并非(O(n log n))

Fast LCM Transformation 快速最小公倍数变换

我们发现LCM即对两个下标标准分解的每个质因数取max，即
(p_{1}^{max(a_{1},b_{1})}p_{2}^{max(a_{2},b_{2})}...p_{k}^{max(a_{k},b_{k})})
但是由于LCM的性质，我们下标的范围会达到(n^{2})，比暴力好不到哪里，所以讨论起来就没有太多意义

如果我们强制限定长度在N范围内，即变换成数组标准分解的多维前缀和，最后多维前缀差分还原回去，类似FGT

FGT的代码

for(int i = 1 ; i <= tot ; ++i) {
    for(int j = N / prime[i] ; j >= 1 ; --j) {
        a[j] += a[j * prime[i]];
        b[j] += b[j * prime[i]]; 
    }    
}
for(int i = 1 ; i <= n ; ++i) {
    c[i] = a[i] * b[i];
}
for(int i = tot ; i >= 1 ; --i) {
    for(int j = 1 ; j <= N / prime[i] ; ++j) {
        c[j] -= c[j * prime[i]];
    }
}