51nod 2026 Gcd and Lcm 杜教筛+数论

题意

已知$f(x)=\sum_{d|x}\mu(d)*d$
现在请求出下面式子的值
$\sum\limits_{i=1}^n\sum\limits_{j=1}^nf(gcd(i,j))*f(lcm(i,j))$
答案模1e9+7。
$n<=10^9$

分析

设$x=p_1^{a_1}*...*p_k^{a_k}$，那么有$f(x)=\prod\limits_{i=1}^k(1-p_i)$。这步十分的关键。
然后就可以得到
$f(i)*f(j)=f(i*j)*f(gcd(i,j))$
$f(lcm(i,j))=f(i*j)$
所以$f(gcd(i,j))*f(lcm(i,j))=f(gcd(i,j))*f(i*j)=f(i)*f(j)$
所以说答案就是$(\sum\limits_{i=1}^nf(i))^2$。
然后就可以爱怎么搞怎么搞了。

要求

$$\sum_{i=1}^n\mu(i)*i$$
这个东西显然是个积性函数，考虑用杜教筛来做。将 $f(x)$卷上 $g(x)=x$那么有
$$\sum_{i=1}^n\sum_{d|i}\mu(d)*d*\frac{i}{d}=\sum_{i=1}^n\sum_{d|i}mu(\frac{i}{d})*\frac{i}{d}*d$$
其中左边的部分就等于1，右边的部分有
$$\sum_{i=1}^n\sum_{d|i}mu(\frac{i}{d})*\frac{i}{d}*d$$
$$=\sum_{i=1}^ni*\sum_{d|i}\mu(\frac{i}{d})$$
$$=\sum_{d=1}^nd*\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)*i$$
这就可以用杜教筛来搞了。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;

typedef long long LL;

const int N=10000005;
const int MOD=1000000007;

int n,s[N],tot,prime[N];
bool not_prime[N];
map<int,int> w;

void get_prime(int n)
{
    s[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!not_prime[i]) prime[++tot]=i,s[i]=-1;
        for (int j=1;j<=tot&&i*prime[j]<=n;j++)
        {
            not_prime[i*prime[j]]=1;
            if (i%prime[j]==0) break;
            s[i*prime[j]]=-s[i];
        }
    }
    for (int i=1;i<=n;i++) s[i]=(s[i-1]+s[i]*i)%MOD;
}

int get_sum(int l,int r)
{
    return ((LL)r*(r+1)/2-(LL)l*(l-1)/2)%MOD;
}

int get_s(int n)
{
    if (n<=10000000) return s[n];
    if (w[n]) return w[n];
    int ans=1;
    for (int i=2,last;i<=n;i=last+1)
    {
        last=n/(n/i);
        (ans-=(LL)get_sum(i,last)*get_s(n/i)%MOD)%=MOD;;
    }
    return w[n]=ans;
}

int main()
{
    scanf("%d",&n);
    get_prime(10000000);
    int ans=0;
    for (int i=1,last;i<=n;i=last+1)
    {
        last=n/(n/i);
        (ans+=(LL)(n/i)*(get_s(last)-get_s(i-1))%MOD)%=MOD;
    }
    ans+=ans<0?MOD:0;
    printf("%d",(LL)ans*ans%MOD);
    return 0;
}