欧拉函数
欧拉函数就是指:对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。
φ(n)=n * (1-1/p1) * (1-1/p2) * (1-1/p3) * (1-1/p4)……(1-1/pn),其中pi为n的质因数。
求单个欧拉函数
long long euler(long long n){
long long ans = n;
for(int i=2; i*i<=n; i++){
if(n%i==0){
ans-=ans/i; //等价于通项,把n乘进去
while(n%i==0) //确保下一个i是n的素因数
n/=i; //算术基本定理,每个数都可以分解成几个素数相乘的形式
}
}
if(n>1)
ans-=ans/n; //n应该等于1,如果大于1,表示最后可能还剩下一个素因数没有除
return ans;
}
打表求欧拉函数
void euler() {
e[1]=1;
for(int i=2; i<maxn; i++){
if(!e[i]) {
for(int j=i; j<maxn; j+=i){
if(!e[j])
e[j]=j;
e[j]=e[j]/i*(i-1); //e[j]=e[j]*((i-1)/i);
}
}
}
}
POJ2478 Farey Sequence
可以看出后面的输出的值都是前面所有数用欧拉函数求出来的值得和,所以这样就可以提前跑一遍欧拉函数。预处理一下,后面直接输出就行了。
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
ll e[1000005];
ll num[1000005];
void euler(){
e[1]=1;
for(int i=2; i<1000005; i++){
if(!e[i]) {
for(int j=i; j<1000005; j+=i){
if(!e[j])
e[j]=j;
e[j]=e[j]/i*(i-1); //e[j]=e[j]*((i-1)/i);
}
}
}
}
int main(){
euler();
for (int i=2; i<=1000005; i++){
num[i]=num[i-1]+e[i];
}
int n;
while (scanf("%d", &n) && n)
cout<<num[n]<<endl;
return 0;
}