1220 约数之和(杜教筛与莫比乌斯）

题目：求${\sum }_{i=1}^n{\sum }_{j=1}^n{\sigma }_1(ij)$
这题推了一整个下午，留下了不学无术的眼泪orz；

前置公式：
(1)${\sum }_{i=1}^n\sigma (i)={\sum }_{i=1}^n{\sum }_{j=1}^n[j|i]\cdot j={\sum }_{i=1}^{n}i\cdot {\sum }_{j=1}^n[i|j]={\sum }_{i=1}^{n}i\cdot \lfloor \frac{n}{i}\rfloor$

（2）${\sigma }_1(ij)={\sum }_{x|i}{\sum }_{y|j}[gcd(x,y)=1]\frac{xj}{y}={\sum }_{x|i}{\sum }_{y|i}x\cdot y[(\frac{i}{x},y)=1]$

（3）$\sum _{i=1}^n\sum _{j=1}^{\lfloor \frac{n}{i}\rfloor }j=\sum _{j=1}^{n}j\cdot \lfloor \frac{n}{j}\rfloor$

（4）${\sum }_{d|n}\mu （d)=[n==1]$

（5）$\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{j}\rfloor =\lfloor \frac{n}{dj}\rfloor$

至于这些公式的证明
第一个已经顺带证明在那了
第二个公式跟${\sigma }_0(ij)={\sum }_{x|n}{\sum }_{y|n}[(x,y)=1]$ 的出发差不多
以自己的理解解释一下：简单说明，${\sigma }_1(n)={\prod }_i(1+p_i+p_i^2+p_i^3+\cdots +p_i^{a_i})$

${\sigma }_0(n)$ 等于上式连乘后产生项的个数，即每一项选一个质数（指数从0到$a_i$ ）构成一个因数，而对于${\sigma }_0(ij)$ 取每一个$i$的因子$x$与$j$的每一个因子$y在x与y$互质条件下配对（不是指相乘)，假设i可分解的连乘部分包括（$1+p_1+p_1^2+p_1^3$ ），$j$ 可出（$1+p_1+p_1^2$ ） 那对于$ij$便可以分出$1+p_1+p_1^2+p_1^3+p_1^4+p_1^5$ 那对$i$在给其因子配质数$p_1$ 的方式包括(指数)0,1,2,3，同理j可以0,1,2，$i,j$可以0,1,2,3,4,5 ,那么问题来了，5等于什么?，等于2+3,那在保证x,y互质的情况下，即保证了其两种选到的因子里包括$p_1$的指数至少有一个一个为零，就相当于选了都为0，和$i$中的1,2,3和$j$种的1,2加起来选到$p_1$的方式便跟$ij $ 一样了，所以等式成立
$那么对于{\sigma }_1(ij)来说也是同个道理，\frac{xj}{y}相当于你i选了p_1指数加上j的p_1指数，那这里就可以保证ij中p_1指数每一个都能被选到匹配成因子，而互质保证了选择上的不重复$
第三个根据贡献值来推
后面两个百度一下有很多

接下来正式推公式

$ans=\sum _{i=1}^n\sum _{j=1}^n{\sigma }_1(ij)$
$=\sum _{i=1}^n\sum _{j=1}^n\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]\frac{xj}{y}$
$=\sum _{i=1}^n\sum _{j=1}^n\sum _{x|i}\sum _{y|j}\frac{xj}{y}\sum _{d|(x,y)}\mu （d）$
到这一步就是把前置公式往里套，接下来可以换成枚举d的值，对于每一个$(x,y)$d为其的约数，即对于取定一个d，$x和y$都是$d$的倍数，又$x,y$分别为$i,j$的倍数，所以$i,j$只需要枚举d的倍数即可
$=\sum _{d=1}^n\mu （d）\sum _{d|i}^n\sum _{d|j}^n\sum _{x|i,d|x}\sum _{y|i,d|y}\frac{xj}{y}$
$=\sum _{d=1}^n\mu （d）\sum _{d|x}^n\sum _{d|y}^n\sum _{x|i}^n\sum _{y|j}^n\frac{xj}{y}$
到这一步,令$x=a\cdot d,y=b\cdot d,则i=k_{1}ad,j=k_{2}bd$代入有
$=\sum _{d=1}^n\mu （d）\sum _{a=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{b=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k_1=1}^{\lfloor \frac{n}{ad}\rfloor }\sum _{k_2=1}^{\lfloor \frac{n}{bd}\rfloor }{k}_{2}ad$
我们将与$k_{2}ad无关的k_1放下来$
$=\sum _{d=1}^{n}d\mu （d）\sum _{a=1}^{\lfloor \frac{n}{d}\rfloor }a\lfloor \frac{n}{ad}\rfloor \sum _{b=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k_2=1}^{\lfloor \frac{n}{bd}\rfloor }{k}_2$
$=\sum _{d=1}^{n}d\mu （d）\sum _{a=1}^{\lfloor \frac{n}{d}\rfloor }a\cdot \lfloor \frac{\lfloor \frac{n}{d}\rfloor }{a}\rfloor \sum _{b=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{k_2=1}^{\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{b}\rfloor }{k}_2$
代入前置公式（3）得
$=\sum _{d=1}^{n}d\mu （d）\sum _{a=1}^{\lfloor \frac{n}{d}\rfloor }a\cdot \lfloor \frac{\lfloor \frac{n}{d}\rfloor }{a}\rfloor \sum _{b=1}^{\lfloor \frac{n}{d}\rfloor }b\cdot \lfloor \frac{\lfloor \frac{n}{d}\rfloor }{b}\rfloor$
代入前置公式(1)得
$=\sum _{d=1}^{n}d\mu （d）{\sigma }_1^2(\lfloor \frac{n}{d}\rfloor )$
$d\cdot \mu (d)卷积后的前缀和满足MU(n)=1-\sum _{i=2}^{n}i\cdot MU\lfloor \frac{n}{i}\rfloor$
接下来杜教筛分块打表即可
代码：

#include <iostream>
#include<cstring>
#include<algorithm>
#include<string>
#include<cmath>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
const int MX= 1e6;
vector<int>pri;
bool isp[MX];
ll mu[MX];
map<ll, ll>mp;
void init()
{
	isp[0] = isp[1] = false;
	mu[1] = 1;
	for (int i = 2; i < MX; i++) {

		if (isp[i]==0) {
			pri.push_back(i);
			mu[i] = -1;// cout <<i << endl;
		}
		for (int j = 0; j < pri.size() && i * pri[j] < MX; j++) {
			 int cur = i * pri[j];
			isp[cur] =1;
			if (i % pri[j]) {
				mu[cur] = -mu[i];
			}
			else {
				mu[cur] = 0;
				break;
			}
		}
	}
	mu[0] = 0;
	for (int i = 1; i < MX; i++)
	{
		mu[i] = (mu[i - 1] + mu[i] * i)%mod;
	}
}
ll cacl(ll n)
{
	return ((n*(n + 1)) / 2) % mod;
}
ll MU(ll x)
{

	if (x < MX)
	{
		return mu[x];
	}
	if (mp[x])
		return mp[x];
	ll ans = 1, j;
	for (int i = 2; i <= x; i = j + 1)
	{
		j = x / (x / i);
		ans = (ans - 1ll * MU(x / i)*(cacl(j)-cacl( i -1)) %mod+ mod) % mod;

	}
	return mp[x] = ans;
}

ll sigma(ll n)
{
	ll ans = 0;
	ll j;
	for (ll i = 1; i <= n; i++)
	{
		j = n / (n / i);
		ans = (ans + (cacl(j) - cacl(i - 1))*(n/i)) % mod;
		i = j;
	}
		return ans*ans%mod;
}
int main()
{
	init();
	ll n;
	cin >> n;
	ll j;
	ll ans = 0;
	for (int i = 1; i <= n; i++)
	{
		j = n / (n / i);
		ans = (ans + (MU(j) - MU(i - 1))*sigma(n / i)%mod+mod) % mod;
		i = j;
	}
	cout << ans << endl;
	
}在这里插入代码片