欧拉函数(例:HDU1286 找新朋友)

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

通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数(小于等于1)就是1本身)。 (注意:每种质因数只一个。比如12=2*2*3那么φ(12)=12*(1-1/2)*(1-1/3)=4
若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
设n为正整数,以 φ(n)表示不超过n且与n互
素的正整数的个数,称为n的欧拉函数值,这里函数
φ:N→N,n→φ(n)称为欧拉函数。
欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。
特殊性质:当n为奇数时,φ(2n)=φ(n), 证明与上述类似。
若n为质数则φ(n)=n-1。

以上来自百度

欧拉函数伪代码

int euler[maxn];    //用于记录每个数字的欧拉函数值
void make()
{
    euler[1]=0;
    int i;
    for(i=2;i<=maxn;++i)
        euler[i]=i;

    for(i=2;i<=maxn;++i)
        if(euler[i]==i)
            for(int j=i;j<=maxn;j+=i)
                euler[j]=euler[j]/i*(i-1);  //核心循环
}

HDU1286 找新朋友 题解
思路:输入一个数,输出其欧拉函数值即可
代码:

#include<iostream>
using namespace std;
const int maxn=32790;
int euler[maxn+2];
void make()
{
    euler[1]=0;
    int i;
    for(i=2;i<=maxn;++i)
        euler[i]=i;
    for(i=2;i<=maxn;++i)
        if(euler[i]==i)
            for(int j=i;j<=maxn;j+=i)
                euler[j]=euler[j]/i*(i-1);
}
int main()
{
    make();
    int n,t,sum;
    cin>>t;
    while(t--)
    {
        cin>>n;
        cout<<euler[n]<<endl;
    }
    return 0;
}

你可能感兴趣的:(欧拉函数(例:HDU1286 找新朋友))