51nod 欧拉函数之和

题目在这里呀！

跟上面一题一样，也是杜教筛的基础题。

推导也跟上面一题基本一样。

for(int i=1;i<=n;i++) ans+=phi[d] (d|i)

φ*1=id 所以呢ans=n*(n+1)/2。

同样的，

for(int i=1;i<=n;i++)
    for(int j=1;j<=n/i;j++) ans+=phi[j];

后者的ans也等于n*(n+1)/2，不再解释了。

同样把i=1的单独拿出来通过移项就可以得到

Σi=1~n φ(i) = n*(n+1)/2 - Σi=2~n Σj=1~⌊n/i⌋ φ(j)

（公式真的好难看qwq

//Suplex
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 5000000+5000
#define ll long long
using namespace std;
const int mod=1000007;
const int P=1000000007;
int cnt,tot,prime[N],head[mod+10],Next[N];
bool vis[N];
ll n,phi[N],t[500000],hashx[500000];

inline void insert(ll k,ll val){
	int key=k % mod;
	for(int i=head[key];i;i=Next[i])
		if(hashx[i]==k) return;
	++cnt;Next[cnt]=head[key];head[key]=cnt;
	hashx[cnt]=k;t[cnt]=val;
}

inline ll find(ll k){
	int key=k % mod;
	for(int i=head[key];i;i=Next[i])
		if(hashx[i]==k) return t[i];
	return -1;
}

void pre(){
	phi[1]=1;
	for(int i=2;i<N;i++){
		if(!vis[i]){
			prime[++tot]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=tot && (ll)i*prime[j]<N;j++){
			vis[i*prime[j]]=1;
			if(i % prime[j]==0){
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}else phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
	}
	for(int i=2;i<N;i++) (phi[i]+=phi[i-1]) %= P;
}

ll solve(ll n){
	if(n<=N) return phi[n];
	ll k=find(n);
	if(k != -1) return k;
	ll next,ans=(n % P*((n+1) % P)/2) % P;
	for(ll i=2;i<=n;i=next+1){
		ll now=n/i;
		next=n/now;
		ans=((ans-(ll)(next-i+1)*solve(n/i) % P) % P+P) % P;
	}
	insert(n,ans);
	return ans;
}

int main()
{
	scanf("%lld",&n);
	pre();
	printf("%lld\n",solve(n));
	return 0;
}