蓝桥杯 筛数（欧拉函数+线性筛）

问题描述：

炫炫学了筛法之后，很想用筛法求欧拉函数。他决定求 $1$ 到 $N$ 的所有数的欧拉函

数值。

输入格式

输入的第一行包含 $1$ 个整数 $n$ 。

输出格式

输出若干行，每行包含一个整数，第 $i$ 行表示 $i$ 的欧拉函数值

样例输入

2

样例输出

1
1

数据规模与约定

$n\le 500000$

前置知识

欧拉函数的定义

$1\sim N$ 中与 $N$ 互质的数的个数被称为欧拉函数，记为 $\varphi (N)$ 。

若在算数基本定理中，$N=p_1^{a_1}{p}_2^{a_2}\dots p_m^{a_m}$

则： $\varphi (N)=N\times \frac{p_1-1}{p_1}\times \frac{p_2-1}{p_2}\times ...\times \frac{p_m-1}{p_m}$

欧拉函数的证明

解题思路

利用线性筛法求出每个数的欧拉函数，质数的欧拉函数值就是 $P-1$

对于一个数 $i$ ，和一个质数 $p$ ，有以下两种情况：

情况1：

$i\mathrm{mod}p==0$ ，此时 $\phi (i\ast p)=\phi (i)\times p$，证明如下：

如果 $p$ 是 $i$ 的一个因子，$N$ 的某一项可以变成 $p_k^{a_k+1}$，即

$\phi (p_k^{a_k+1})=p_k^{a_k+1}-p_k^{a_k}=p_k^{a_k+1}\times (1-\frac{1}{p_k})=p_k\times p_k^{a_k}\times (1-\frac{1}{p_k})$

那么：

$\phi (i\ast p)=p_k\times N\times (1-\frac{1}{p_1})\times ...\times (1-\frac{1}{p_m})=p\times \phi (i)$

得证。

情况二：

$i\mathrm{mod}p\ne 0$，此时 $\phi (i\ast p)=\phi (i)\times (p-1)$，证明如下：

其实就是 $N$ 多了一个质因数 $N=p_1^{a_1}\times ...\times p_k^{a_k}\times p_{k+1}$

所以 $\phi (i\ast p)=\phi (i)\ast \phi (p)$ ，而 $\phi (p)=p-1$

故：$\phi (i\ast p)=\phi (i)\times (p-1)$，结论得证

代码展示

import java.io.*;
import java.util.*;

public class Main {
	static final int N = 500010;
	static int[] primes = new int[N];
	static int[] phi = new int[N];
	static boolean[] st = new boolean[N];
	static int cnt;
	
	static void get_eulers(int n) {
	    for (int i = 2; i <= n; i ++ ) {
	        if (!st[i]) {
	            primes[cnt++] = i;
	            phi[i] = i-1;
	        }
	        for (int j = 0; primes[j] <= n/i; j ++ ) {
	            st[i*primes[j]] = true;
	            if (i % primes[j] == 0) {
	                phi[i*primes[j]] = primes[j] * phi[i];
	                break;
	            }
	            phi[i*primes[j]] = phi[i] * (primes[j]-1);
	        }
	    }
	}

	public static void main(String[] args) throws Exception {
		BufferedReader bf = new BufferedReader(new InputStreamReader(System.in));
		PrintWriter pw = new PrintWriter(System.out);
		int n = Integer.parseInt(bf.readLine().split(" ")[0]);
		get_eulers(n);
		phi[1] = 1;
		for (int i = 1; i <= n; i ++ ) pw.println(phi[i]);
		pw.close();
	}
}