51nod 2026 Gcd and Lcm【杜教筛】

题目描述：

已知：$f(x)={\sum }_{d|x}\mu (d)\ast d$
求$\sum _{i=1}^n\sum _{j=1}^{n}f(\gcd (i,j))\ast f(\text{lcm}(i,j))\mathrm{m}\mathrm{o}\mathrm{d}10^9+7$
$n\le 10^9$

题目分析：

这个题。。认真你就输了。。
根据$\mu (d)\ast d$是积性函数，它的卷积也是积性函数。
所以有：$$f(\gcd (i,j))\ast f(\text{lcm}(i,j))=f(i)\ast f(j)$$
因为$\gcd$和$\text{lcm}$一个是取质因子中较小的，一个是取质因子中较大的，按照质因子拆开后重新组合就可以得到$f(i)\ast f(j)$

于是问题就变成了求$(\sum _{i=1}^{n}f(i){)}^2$
简单变换一下就是$$(\sum _{d=1}^n\mu (d)\ast d\ast \lfloor \frac{n}{d}\rfloor {)}^2$$
杜教筛筛$\mu (d)\ast d$的前缀和即可，$(\mu \cdot id)\ast id=ϵ$

Code：

#include<bits/stdc++.h>
#define maxn 1000005
using namespace std;
const int mod = 1e9+7;
int n,s[maxn],p[maxn/5],F[maxn];
bool v[maxn];
void Prime(const int N){
	s[1]=1;int m=0;
	for(int i=2;i<=N;i++){
		if(!v[i]) p[++m]=i,s[i]=-1;
		for(int j=1,k;j<=m&&(k=i*p[j])<=N;j++){
			v[k]=1; if(i%p[j]==0) break; s[k]=-s[i];
		}
	}
	for(int i=1;i<=N;i++) s[i]=(s[i-1]+s[i]*i)%mod;
}
int solve(int n){
	if(n<=maxn-5) return s[n];
	if(~F[::n/n]) return F[::n/n];
	int ret=1;
	for(int i=2,j;i<=n;i=j+1)
		j=n/(n/i),ret=(ret-1ll*(i+j)*(j-i+1)/2%mod*solve(n/i))%mod;
	return F[::n/n]=(ret+mod)%mod;
}
int main()
{
	Prime(maxn-5),memset(F,-1,sizeof F);
	scanf("%d",&n);
	int ans=0;
	for(int i=1,j;i<=n;i=j+1)
		j=n/(n/i),ans=(ans+1ll*(n/i)*(solve(j)-solve(i-1)))%mod;
	printf("%d\n",1ll*ans*ans%mod);
}