CCPC网络赛 HDU-6706 huntian oy（莫比乌斯反演+杜教筛+sum(i*phi(i))模板）

链接：http://acm.hdu.edu.cn/showproblem.php?pid=6706

题意：求

官方题解：

 

杜教筛强推博客：https://www.cnblogs.com/peng-ym/p/9446555.html

 

#include <bits/stdc++.h> 
#include <tr1/unordered_map>
#define ll long long
#define ull unsigned long long 
using namespace std;
const int N = 3e6+10;
const ll mod = 1e9+7;
ll phi[N],prime[N],cnt,ni2,ni6;
ll pre1[N];
tr1::unordered_map <ll,ll> mp1; 
template<typename T>inline void read(T &x)
{
    x=0;
    static int p;p=1;
    static char c;c=getchar();
    while(!isdigit(c)){if(c=='-')p=-1;c=getchar();}
    while(isdigit(c)) {x=(x<<1)+(x<<3)+(c-48);c=getchar();}
    x*=p;
}
void Init()
{
	cnt=0;
	phi[1]=1;
	for(int i=2;i<=N-10;i++)
	{
		if(!phi[i])
		{
			prime[cnt++]=i;
			phi[i]=i-1;
		}
		for(int j=0;j<cnt&&(ll)i*prime[j]<=N-10;j++)
		{
			if(i%prime[j])
			{
				phi[i*prime[j]]=phi[i]*(prime[j]-1+mod)%mod;
			}
			else
			{
				phi[i*prime[j]]=phi[i]*prime[j]%mod;
				break;
			}
		}
	}
	for(int i=1;i<=N-10;i++)
	{
		pre1[i]=(pre1[i-1]+i*phi[i]%mod)%mod;	
	}	
	return ;
}
//sum(i*phi(i)) i=1,2,...,n
ll S(ll n)
{
	if(n+10<=N) return pre1[n]; 
	if(mp1[n]) return mp1[n];
	ll ans=n*(n+1)%mod*(2*n%mod+1)%mod*ni6%mod; 
	ll r;
	for(ll l=2;l>=0&&l<=n;l=r+1)
	{
		r=n/(n/l);
		ans-=(r-l+1)*(l+r)%mod*ni2%mod*S(n/l)%mod; 
		ans=(ans%mod+mod)%mod;
	} 
	return mp1[n]=ans; 
}
ll poww(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;	
		b>>=1;
	}
	return ans;
}
int main(void)
{
	Init();
	int t,n,a,b;
	ni2=poww(2,mod-2);
	ni6=poww(6,mod-2);
	scanf("%d",&t);
	while(t--)
	{
		read(n),read(a),read(b); 
		printf("%lld\n",(S(n)-1+mod)%mod*ni2%mod);	
	}
	return 0;	
}