2019-8-6训练日志欧拉

今天还在看数学知识这部分，已经快看完了，主要都是一些定理，性质，推论等，还是要多看看，记住。  

最大公约数：在所有a和b的公约数中最大的一个，记为gcd（a，b）；最小公倍数：在所有a和b的公倍数中最小的一个，记为Icm（a，b）。 定理：对任意的正整数a，b， gcd（a，b）*Icm（a，b）=a*b。

更相减损术：对任意的正整数a，b。a>=b,有gcd（a，b）=gcd（b，a-b）=gcd（a，a-b）；

欧几里得算法：对任意的正整数a，b，b0，gcd（a，b）=gcd（b，amodb）；

int gcd(int a,int b){
    return b ? gcd(b,a%b) : a;
}

使用欧几里得算法时复杂度为O（log（a+b））。欧几里得算法是我们最常用的方法求最大公约数，但高精度取模不容易实现，需要做高精度运算时，可考虑用更相减损术。

互质即是两个正整数的最大公约数为1.

欧拉函数：1~N中与N互质的数个数为欧拉函数，记为(N).欧拉函数有很多性质，我感觉应该很少用到。我们只要掌握怎么求欧拉函数应该就行了。最快的应该就是利用线性筛法的思想在O（N）内求出1~N每个数的欧拉函数。

int v[510000],prime[510000],phi[510000],m;
void euler(int n)
{
    memset(v,0,sizeof(v));
    m=0;
    for(int i=2;i<=n;i++){
        if(v[i]==0){
            v[i]=i,prime[++m]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=m;j++){
            if(prime[j]>v[i]||prime[j]>n/i) break;
            v[i*prime[j]]=prime[j];
            phi[i*prime[j]]=
               phi[i] * (i%prime[j] ? prime[j]-1 : prime[j]);
        }
    }
phi[1]=1;
}
int main()
{
    int n;
    cin>>n;
    euler(n);
    for(int i=2;i<=n;i++)
    cout<<phi[i]<<' ';
    return 0;
}

还有一个就是通过欧拉函数性质延伸的积性函数了，如果a，b互质，有f（ab）=f（a）*f（b），即为积性函数。