贴一下欧拉函数的模版,已备日后使用
int oula(int n) { //O(n)
int rea = n;
for (int i = 2; i <= n; i ++) {
if(n % i == 0) {
rea = rea - rea/i;
do{
n /= i;
}while(n % i == 0);
}
}
return rea;
}
int oula_youhua(int n){//O(√n)
int rea = n;
for (int i = 2; i * i <= n; i ++) {
if(n % i == 0) {
rea = rea - rea / i;
do{
n /= i;
}while(n % i == 0);
}
}
if(n > 1)
rea = rea - rea / n;
return rea;
}
bool Isprime[50000];
int p[maxn];
void prim() {//素数打表
memset(Isprime, 0, sizeof(Isprime));
Isprime[0] = Isprime[1] = 1;
int k = 0;
for (int i = 2; i < 50000; i ++) {
if(!Isprime[i])
p[k++] = i;
for (int j = 0; j < k && i * p[j] < 50000; j ++) {
Isprime[i *p[j]] = 1;
if(!(i % p[j]))
break;
}
}
}
//筛法打表
int phi(int n){
int rea = n;
for (int i = 0; p[i] * p[i] <= n; i ++) {
if(n % p[i] == 0) {
rea = rea - rea / n;
do{
n /= p[i];
}while(n % p[i] == 0);
}
}
if(n > 1)
rea = rea - rea / n;
return rea;
}
void ditui() {//递推打表
int i,j;
for(i = 1; i < maxn; i++)
p[i] = i;
for(i = 2;i < maxn; i++){
if(p[i] == i){
for(j = i; j < maxn; j += i){
p[j] = (p[j] / i) * (i-1);
}
}
}
}