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

#include <iostream>

#include <cstdio>

#include <cstring>

#include <ctime>

#include <cmath>

using namespace std;

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);

    return 0;

}