欧拉函数学习笔记整理 POJ 2407+POJ 1284+POJ 2478+POJ 3090

定义

欧拉函数φ(x)代表1-x中与x互质的数的个数。

基本相关性质

1.设$x={a1}^{p1}*{a2}^{p2}*...*{an}^{pn}$
则φ(x)=$x*(1-\frac{1}{a1})*(1-\frac{1}{a2})*...*(1-\frac{1}{an})$;
2.当gcd(m,n)=1时，φ(mn)=φ(m)*φ(n);
3.当x为质数时，φ(x)=x-1,当x为奇数时，φ(2x)=x-1;
4.若x=${p}^{k}$,则$φ(x)={p}^{k}-{p}^{k-1}=(p-1){p}^{k-1};$
5.当p为质数时，
　(1)若满足p|n且${p}^{2}∤n$，$φ(n)=φ(\frac{n}{p})*(p-1)$;

　(2)若满足p|n且${p}^{2}|n$，$φ(n)=φ(\frac{n}{p})*p$;

6.当n>1时，1-n中与n 互质的数的整数和为$\frac{n*φ(n)}{2}$；

7.$\sum_{d|n}φ(d)=n$

证明：(1)若n==1，φ(n)=1=n;

　　　(2)若n为质数时，φ(n)=n-1，则φ(n)+φ(1)=(n-1)+1=n;

　　　(3)若$n={p}^{k}$，且p为质数,k>1，

　　　　$\sum_{d|n}φ(d)=φ(1)+φ(p)+φ({p}^{2})+...+φ({p}^{k})$

　　　　$=1+(p-1)+p(p-1)+{p}^{2}(p-1)+...+{p}^{k}(p-1)$

　　　　$=1+(p-1)\frac{1-{p}^{k}}{1-p}={p}^{k}=n$;
　　　(4)n含有多个质因子，即$n={a1}^{p1}{a2}^{p2}...{an}^{pn}$,

　　　　$\sum_{d|n}φ(d)=\sum_{i1=0}^{p1}\sum_{i2=0}^{p2}...\sum_{in=0}^{pn}φ({a1}^{i1}{a2}^{i2}...{an}^{in})$

8.欧拉定理：若gcd(a,n)=1,有${a}^{φ(n)}≡1(mod~~n)$

　费马定理：若gcd(a,p)=1,则${a}^{p-1}≡1(mod~~p)~~~$(将φ(p)=p-1带入欧拉定理)。
9.快速求出欧拉函数的值（P为N的质因子）
若(N%P==0&&(N/P)%P==0) 则有：E(N)=E(N/P)*P;
若(N%P==0&&(N/P)%P!=0) 则有：E(N)=E(N/P)*(P-1)

相关题目

POJ 2407 Relatives(基础题)

POJ 2407题目链接

很直白的使用欧拉函数

#include <iostream>
#include<cstdio>
using namespace std;
int n;
void cal()
{
    int ans=n;
    int t=n;
    for(int i=2; t>1; i++)
    {
        if(t%i==0)
        {
            while(t%i==0)
            {
                t/=i;
            }
            ans=ans-ans/i;
        }
    }
    cout<<ans<<endl;
}
int main()
{
    while(cin>>n)
    {
        if(n==0)break;
        cal();
    }
    return 0;
}

POJ 1284 Primitive Roots(+原根)

POJ 1284题目链接

主要涉及一个有关原根的定理，
定理：如果p有原根，则它具有φ(φ(n))个不同的原根，又因为p为质数，所以φ(p)=p-1，所以有φ(p-1)个原根。

#include <iostream>
#include<cstdio>
using namespace std;
#define ll __int64
int n;
void cal()
{
    int ans=n;
    int t=n;
    for(int i=2; t>1; i++)
    {
        if(t%i==0)
        {
            while(t%i==0)
            {
                t/=i;
            }
            ans=ans-ans/i;
        }
    }
    cout<<ans<<endl;
}
int main()
{
    while(cin>>n)
    {
        if(n==0)break;
        n--;
        cal();
    }
    return 0;
}

POJ 2478 Farey Sequence(打表)

POJ 2478题目链接

数据范围太大，先打表统计欧拉数。

#include <iostream>
#include<cstdio>
using namespace std;
#include<cstring>
#define maxn 1000006
#define ll __int64
int x;
ll phi[maxn];
void init()
{
    for(int i=2;i<=maxn;i++)
    {
        if (!phi[i])
        {
            for (int j=i;j<=maxn;j+=i)
            {
                if (!phi[j])
                    phi[j]=j ;
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
}
int main()
{
    memset(phi,0,sizeof(phi));
    init();
    while(cin>>x)
    {
        if(x==0)break;
        ll ans=0;
        for(int i=1; i<=x; i++)
            ans+=phi[i];
        cout<<ans<<endl;
    }
    return 0;
}

POJ 3090 Visible Lattice Points

POJ 3090题目链接

该题主要求从原点看可以看到几个点，覆盖到的忽略掉。
主要是统计斜率，求斜率（0，1）开区间为基准进行求解，最后结果*2+3；
与2478相类似，每次增加的都是该值的素因子个数，如5*5，增加了1/5,2/5,3/5,4/5；6*6会增加1/6,5/6。

#include <iostream>
#include<cstdio>
using namespace std;
#include<cstring>
#define maxn 1000006
#define ll __int64
int x;
ll ans[maxn];
ll phi[maxn];
void init()
{
    for(int i=2; i<=maxn; i++)
    {
        if (!phi[i])
        {
            for (int j=i; j<=maxn; j+=i)
            {
                if (!phi[j])
                    phi[j]=j ;
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
    for(int i=2; i<=maxn; i++)
        ans[i]=ans[i-1]+2*phi[i];
}
int main()
{
    memset(phi,0,sizeof(phi));
    memset(ans,0,sizeof(ans));
    init();
    int T;
    cin>>T;
    for(int num=1; num<=T; num++)
    {
        cin>>x;
        cout<<num<<" "<<x<<" "<<ans[x]+3<<endl;  //+3分别为+斜率=1，斜率=0，斜率=∞
    }
    return 0;
}