欧拉函数PHI

欧拉函数

首先什么是欧拉函数，科普一下
在数论，对正整数n，欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名，它又称为Euler’s totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明

1. java版求解方法

public class Main {
// 筛法欧拉函数
static final int MAXN =10000000 ;//要筛选的数
static int phi[] = new int[MAXN + 2];
public static void main(String[] args) {
for (int i = 1; i <= MAXN; i++) {
phi[i] = i;

    }
    for (int i = 2; i <= MAXN; i += 2) {
        phi[i] /= 2;
    }
    for (int i = 3; i <= MAXN; i += 2) {
        if (phi[i] == i) {
            for (int j = i; j <= MAXN; j += i) {
                phi[j] = phi[j] / i * (i - 1);
            }
        }
    }
    System.out.println(phi[MAXN]);


}

}

单独求解

public class Main2 {
// 单独求解法
/**
* 单独求解的本质是公式的应用
*
*/
static final int x = 10000000;

public static void main(String[] args) {
    System.out.println(euler(x));

}

private static int euler(int x) {
    int res = x;
    for (int i = 2; i < (int) Math.sqrt(x * 1.0) + 1; i++) {
        if (x % i == 0) {
            res = res / i * (i - 1);
        }
        while (x % i == 0) {
            x /= i;//保证i一定是素数
        }
    }
    if (x > 1) {
        res = res / x * (x - 1);
    }
    return res;

}

}

---

同时得到素数表和欧拉函数

import java.util.Arrays;

/**
* 同时得到欧拉函数和素数表
*
* @author Administrator
*
*/
public class Main3 {
static int MAXN = 100;
static boolean check[] = new boolean[MAXN + 2];
static int phi[] = new int[MAXN + 10];// 欧拉函数
static int prime[] = new int[MAXN + 10];
static int tot;// 素数个数

public static void main(String[] args) {
    phi_and_prime_table(MAXN);
    System.out.println(phi[MAXN]);
    for (int i = 0; i < prime.length; i++) {
        System.out.print(prime[i]+" ");
    }
}

private static void phi_and_prime_table(int N) {
    Arrays.fill(check, false);
    phi[1] = 1;
    tot = 0;
    for (int i = 2; i <= N; i++) {
        if (!check[i]) {
            prime[tot++] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; j < tot; j++) {
            if (i * prime[j] > N) {
                break;
            }
            check[i * prime[j]] = true;
            if (i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            } else {
                phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            }
        }

    }
}

}

完毕