闲聊欧拉函数入门

之前看过了欧拉数，想想还是记录一下关键部分，主要方便自己以后忘了可以查看一下，也可以帮助一些朋友学习（虽然不一定有用= =）

n的欧拉(phi)数 的定义 ：小于n且与n互素的整数个数。

两个方法可以求解欧拉值：

一 。初始化phi为n，直接从2遍历sqrt(n)，当遇见素因子x时，phi = phi/x*(x-1)，然后把这个因子除尽，这个式子挺好理解 ：当x能整除n时，把n分成 n/x 份 ， 其中x的倍数一定不是素因子，而x的倍数只占一份，那么  phi/x  就是每份的数的个数，再乘以(x-1)就除去了x的倍数了，最后得到的就是答案。

要注意有些题目可以先预处理出素数表，然后直接遍历素数表，在很多次计算欧拉值得情况下这样可以大大优化效率。

然后如果想求出1~n的所有数的欧拉数也有比较快的做法，详见代码。

求1~n的所有数的phi:复杂度O(nloglogn)

inline void phi_table (int n){
    memset(phi,0,sizeof(phi));
    phi[1]=1;
    for(int i=2;i<=n;i++) if(!phi[i]){
        for(int j=i;j<=n;j+=i){
            if(!phi[j])phi[j]=j;
            phi[j] = phi[j]/i*(i-1);
        }
    }
}

求n的phi：复杂度O(sqrt(n))

inline int get_phi(int n){
    int m = sqrt(n+0.5);
    int phi = n;
    for(int i=2;i<=m&&n>1 ;i++)
    {
        if(n%i==0){
            phi = phi / i * (i-1);
            do{n/=i;}while(n%i==0);
        }
    }if(n>1)phi = phi / n * (n-1);
    return phi ;
}

二  已知n的唯一分解式 然后可以转换成 就可以O(k)次求解了，这个转换有点复杂，建议把第二个式子手动分解一下大概就清楚了~maybe，不然会用就行了，但是这个转换给出了一种解决容斥问题的思想，还是有点用的。

另附POJ2773代码(欧拉+循环节)

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

int prime[1000000];
bool p[1000000+5];


inline int work(int x){
    int num = 1;prime[0]=1;
    memset(p,0,sizeof(p));
    for(int i=2;i<=x;i++)
    {
        if(!p[i]){
            if(x%i==0){
                for(int j=i;j<=x;j+=i)
                    p[j]=true;
            }else
                prime[num++]=i;
        }
    }return num ;
}


int main()
{
    int n,k;
    while(scanf("%d%d",&n,&k)==2)
    {
        int sum = work(n);
        if(k % sum != 0){
            printf("%I64d\n",(long long)k/sum * n + prime[k%sum-1]);
        }else {
            printf("%I64d\n",(long long)(k/sum-1)*n + prime[sum-1]);
        }
    }return 0;
}