51nod 1239 欧拉函数之和
链接:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1239
φ∗id1=id0
Sφ(n)=n−∑i=2niSφ(⌊ni⌋)
#include <algorithm>
#include <string.h>
#include <stdio.h>
#include <cmath>
#define MAXN 10000000
#define SQRT 111111
using namespace std;
typedef long long LL;
const LL P=1e9+7;
int phi[MAXN];
int tmp[SQRT];
void clat(int d,LL n)
{
int m=sqrt(n)+1,ans=0;
for(int L=1;L<m;L++)
{
ans+=(LL)phi[L]*((n/L-n/(L+1))%P)%P;
if(ans>=P)ans-=P;
}
for(int i=(int)(n/m);i>1;i--)
{
LL u=(n/i);
if(u<MAXN)
ans+=phi[u];
else
ans+=tmp[i*d];
if(ans>=P)ans-=P;
}
if(n&1)
n=(((n+1)>>1)%P)*(n%P)%P;
else
n=((n>>1)%P)*((n+1)%P)%P;
tmp[d]=(n-ans+P)%P;
}
int solve(LL n)
{
if(n<MAXN)return phi[n];
for(LL i=(n/MAXN);i;i--) clat((int)i,n/i);
return tmp[1];
}
int main ()
{
phi[1]=1;
for(int i=2;i<MAXN;i++)
{
if(phi[i])continue;
phi[i]=i-1;
LL size=(LL)i*i;
for(int j=i<<1;(LL)j<size&&j<MAXN;j+=i) phi[j]=(phi[j])?phi[j]*phi[i]:phi[i];
if(size>=MAXN)continue;
for(int j=i*i,k=i;j<MAXN;j+=i,k++)
if(k%i)
phi[j]=(phi[k])?phi[k]*phi[i]:phi[i];
else
phi[j]=(phi[k])?phi[k]*i:i;
}
for(int i=1;i<MAXN;i++)
{
phi[i]+=phi[i-1];
if(phi[i]>=P)phi[i]-=P;
}
LL n;
scanf("%lld",&n);
printf("%d\n",solve(n));
return 0;
}