狄利克雷卷积学习记录

狄利克雷卷积是函数之间的一种运算

定义$o$运算表示狄利克雷卷积

则$h(n)=g(n)of(n)$，可读作"g 卷 f"

其含义为 $h(n)={\sum }_{d\mid n}g(d)f(\frac{n}{d})$

狄利克雷卷积的性质：

狄利克雷卷积满足交换律，结合律，分配率,即 :
$h(n)of(n)=f(n)oh(n)$

$h(n)of(n)og(n)=h(n)o(f(n)og(n))$

$(h(n)+f(n))og(n)=h(n)og(n)+f(n)og(n)$

狄利克雷卷积的另一个性质是，两个积性函数的狄利克雷卷积也是积性函数，即
若$n=p\times q$,且 p 与 q 互质，假如

$f(n)=f(p)\times f(q)$，

$g(n)=g(p)\times g(q)$，

$h(n)=g(n)of(n)$，

那么 $h(n)=h(p)\times h(q)$；

常见的狄利克雷卷积

莫比乌斯函数$\mu (x)$
$$\mu (x)=\{\begin{array}{lr}1& n=1\\ (-1{)}^k& n无平方因数，且n={\prod }_{i=1}^k{p}_i\\ 0& 其他情况\end{array}$$
设函数$I(n)=1,e(n)=[n=1],id(n)=n;$其中[]表示括号内条件成立为1，否则为0
则有四个常见的性质

1：$id=\phi oI$
展开得 $n={\sum }_{d\mid n}\phi (d)$，这是欧拉函数的一个重要性质，因此得证

2：$e=\mu oI$
这个式子ldy大佬讲的时候说，这个不好证，记住就行了

3：任何函数和e做狄利克雷卷积，值都不变
自己展开就得证了

4：$\phi =ido\mu$
证明如下
对于性质二的式子：$e=\mu oI$
两边同时卷上$\phi$，则$\phi oe=\mu oIo\phi$
因为卷积具有交换律，结合律，所以
$\phi oe=(Io\phi )o\mu$
由性质1，性质3得
$\phi =ido\mu$，得证

由以上性质可以得出一些其他性质：

5：比如，因为$e=\mu oI$
所以$[n=1]={\sum }_{d\mid n}\mu (d)$
当 $n=gcd(i,j)$ 的时候
$[gcd(i,j)=1]={\sum }_{d\mid gcd(i,j)}\mu (d)={\sum }_{d\mid i,d\mid j}\mu (d)$

6：同理因为 $id=\phi oI$
所以 $n={\sum }_{d\mid n}\phi (d)$
当$n=gcd(i,j)$时，
$gcd(i,j)={\sum }_{d\mid gcd(i,j)}\phi (d)={\sum }_{d\mid i,d\mid j}\phi (d)$
这也成为了计算 gcd 的另一种方法，在化简式子时也很常用

莫比乌斯反演

若f，F为任意函数，则
若$foI=F$
所以同时卷上$\mu$
$foIo\mu =Fo\mu$
即$foe=Fo\mu$
$f=Fo\mu$

展开就表示为
若$F(n)={\sum }_{d\mid n}f(d)$
则$f(n)={\sum }_{d\mid n}F(d)\mu (\frac{n}{d})$
这就是莫比乌斯反演，若果${\sum }_{d\mid n}F(d)\mu (\frac{n}{d})$很好求，就可以求 f

例题 [NOI2010] 能量采集

洛谷题目传送门
化简一下题意，就是求${\sum }_{i=1}^n{\sum }_{j=1}^m(2\times gcd(i,j)-1)$

那么也就是求出${\sum }_{i=1}^n{\sum }_{j=1}^{m}gcd(i,j)$就可以了

根据上方的性质6，该式子可以化为

${\sum }_{i=1}^n{\sum }_{j=1}^m{\sum }_{d\mid i,d\mid j}\phi (d)$
将最后一个$\sum$提前，得
${\sum }_d^{min(m,n)}\phi (d){\sum }_{d\mid i}^n{\sum }_{d\mid j}^{m}1$
将最后一项化简，得
${\sum }_d^{min(m,n)}\phi (d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$
于是我们预处理处每一个$\phi (d)$，就行了
最后把答案乘2，在减去n*m就行了

#include<bits/stdc++.h>
using namespace std;
const long long N = 1e5+7;
long long phi[N],n,m;
void euler(long long n)
{
	for(long long i=1;i<=n;i++)
	phi[i]=i;
	for(long long i=2;i<=n;i++)
	{
		if(phi[i]==i)
		{
			for(long long j=i;j<=n;j+=i)
			phi[j]=phi[j]/i*(i-1);
		}
	}
}
int main()
{
	cin>>n;
	cin>>m;
	euler(min(n,m));
	long long ans=0;
	for(long long i=1;i<=min(n,m);i++)
	{
		ans=ans+phi[i]*floor(n/i)*floor(m/i);
	}
	cout<<2*ans-(m*n);
	return 0;
}