基础数论练习1——持续更新

（1）f(n)
这题是素数筛的基础上的找规律的题。
规律:
$Gcd(n)=\{\begin{array}{ll}p& n=p^x\\ 1& n=p_1^{x_1}\ast p_2^{x_2}...p_k^{x_k}\end{array}$(其中p及其p下标元素都是素数)

简单描述规律的逻辑证明：
我们要明确一点$C_n^k=\frac{n\ast (n-1)...(n-k+1)}{1\ast 2...k}=n\ast \frac{n-1}{2}...\frac{n-k+1}{k}$
而且$C_n^n=1$,每一项分母都有分子和它结合成1
而在Gcd(n)中，最后一项为$C_n^{n-1}$，意味着n是不可能被完全消除到1,理论上只有最小为n的数能将n变成1
1)当$n=p^x$，时很容易根据上面的分析知道，n始终都会在$C_n^i$中以$p^j的形式存在$,最多被分解到只有p，那么很显然Gcd(n)=p。
2)当$n=p_1^{x_1}\ast p_2^{x_2}...p_k^{x_k}$时，按照上面的逻辑必定有i使得$C_n^i$没有$p_j项$,故任何有关质因数都无法对Gcd做贡献,Gcd的值就为1.

这里使用埃筛，埃筛的利用数因子重复筛数，所以可以。

#include<iostream>
using namespace std;
typedef long long ll;
const int maxn = 1000000+5;
int factor[maxn] = {0},N;
ll ans[maxn];
void e(){
    for(int i = 2;i <= 1000000;i++){
        if(factor[i]==0){
            factor[i] = i;
            for(int j = i*2;j <= 1000000;j+=i){
                if(factor[j]==0)factor[j] = i;
                else factor[j] = 1;
            }
        }
    }
}
void init(){
    ans[2] = 0;
    for(int i = 3;i <= 1000000;i++)ans[i] =ans[i-1] + factor[i];
}
int main(){
    e();
    init();
    while(cin >> N){
        cout << ans[N] <<endl;
    }
    return 0;
}