设p为质数,phi(n)为n的欧拉函数值。
- phi( p)=p-1,显然从1到p-1都与p互质
- 若i%p==0,则phi(i*p)=phi(i)*p
证明:首先知道两个条件,①用反证法很容易得出gcd(n,i)=1,则gcd(n+i,i)=1;②n,i不互质,则n+i,i不互质,gcd(n,i)=k(k>1),n=ak,i=bk,n+i=(a+b)k。
[1,i]中与i不互质的数有i-phi(i)个, [i+1,i+i]中与i不互质的数也有i-phi(i)个,依次递推,得出[1,i*p]中与i不互质的数有p*(i-phi(i))个。
因为i%p==0,则若有k整除p,则必有k整除i,因此i*p的因数和i的因数是相同的。
故[1,i*p]中与i*p不互质的数为p*(i-phi(i))个,同时也为i*p-phi(i*p),据此得出结论phi(i*p)=phi(i)*p。
- 若i%p!=0,则phi(i*p)=phi(i)*phi( p),即phi(i*p)=phi(i)*(p-1),积性函数的性质。
模板
int prime[maxn], phi[maxn], tot;
bool vis[maxn];
void init() {
phi[1] = 1;
for (int i = 2; i < maxn; i++) {
if (!vis[i]) {
prime[++tot] = i;
phi[i] = i - 1;
}
for (int j = 1; j <= tot && i * prime[j] < maxn; j++) {
vis[i * prime[j]] = true;
if (i % prime[j] == 0) {
phi[i * prime[j]] = phi[i] * prime[j];
break;
}
phi[i * prime[j]] = phi[i] * (prime[j] - 1);
}
}
}
应用
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn = 1e6 + 5;
int prime[maxn], phi[maxn], tot;
bool vis[maxn];
ll ans[maxn];
void solve() {
phi[1] = 1;
for (int i = 2; i < maxn; i++) {
if (!vis[i]) {
prime[++tot] = i;
phi[i] = i - 1;
}
for (int j = 1; j <= tot && i * prime[j] < maxn; j++) {
vis[i * prime[j]] = true;
if (i % prime[j] == 0) {
phi[i * prime[j]] = phi[i] * prime[j];
break;
}
phi[i * prime[j]] = phi[i] * (prime[j] - 1);
}
}
for (int i = 1; i < maxn; i++)
for (int j = 1; i * j < maxn; j++)
ans[i * j] += 1LL * j * phi[j];
for (int i = 1; i < maxn; i++)ans[i] = (ans[i] + 1) * i / 2;
}
int main(void) {
solve();
int t, n;
scanf("%d", &t);
while (t--){
scanf("%d", &n);
printf("%lld\n", ans[n]);
}
return 0;
}