欧拉函数--打表法与直接求值法

#include <iostream>

#include <cstdio>

#include <cstring>

#include <ctime>

#include <cmath>

using namespacestd;

const int maxn =1e6 + 1;

int prime[maxn /10],primesz = 0;

bool isprime[maxn];//0表示是素数

int phi[maxn];//直接求，不用先打素数

//f(a*b) = f(a) * f(b) a,b互素

void allphi(int sz)//打表

{

    for (int i =2; i < sz; i ++) {

        if(!isprime[i]) {phi[i] = i -1;prime[primesz ++] = i;}

        for (int j =0; j < primesz && i *prime[i] < maxn; j ++) {

            prime[i *prime[j]] = true;

            if(i %prime[j] == 0) {

                phi[i *prime[j]] = phi[i] *phi[prime[j]];//**

                break;

            }

            phi[i *prime[j]] = phi[i] * (phi[prime[j]] - 1);//**

        }

    }

}

void getphi(int n)

{

    //没有打表质数，直接求欧拉值

    int ans = n;

    int m =sqrt(n + 0.5);

    for (int i =2; i <= m; i ++) {

        if(n % i ==0) {

            ans = ans / i * (i - 1);

            while (n % i ==0) n /= i;

        }

        if(n ==1) break;

    }

    if(n >1) ans = ans / n * (n - 1);

}

//void getprime(int sz)//线性筛，可以和欧拉打表比较下

//{

//memset(isprime, 0, sizeof(isprime));

//    for (int i = 2; i < sz; i ++) {

//        if(!isprime[i]) prime[primesz ++] = i;

//        for (int j = 0; j < primesz &&  prime[j] * i < maxn; j ++) {//而且不用long long，但是要%

//            isprime[prime[j] * i] = true;

//            if(i % prime[j] == 0) break;

//        }

//    }

//}

int main()

{

    allphi(maxn);

    return0;

}