线性欧拉筛+前缀和

解题思路:
可以看我写的题解   https://blog.csdn.net/weixin_61067952/article/details/131122831?spm=1001.2014.3001.5501

注意事项:

参考代码:

#include<bits/stdc++.h> 
using namespace std;
typedef long long ll;
const int N=1e7+5;
int presum[N];
vector<int>prime;
bool vis[N];
int phi[N];
void get_phi(int n)
{
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            phi[i]=i-1;
            prime.emplace_back(i);
        }
        for(int j=0;j<prime.size()&&i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];break;
            }
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}
int main( )
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    int n;cin>>n;
    int mod=1e9+7;
    get_phi(n);
    for(int i=1;i<=n;i++)
    {
        presum[i]=presum[i-1]+phi[i];presum[i]%=mod;
    }
    ll ans=0;
    for(ll i=1;i<=n;i++)
    {
        ans=(ans+((i*i%mod)*(presum[n/i]*2-1))%mod)%mod;
    }
    cout<<ans;
    return 0;
}