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 
#include 
#include 
#include 
#include 
#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



你可能感兴趣的:(51nod,OI,杜教筛,分块)