欧拉函数的应用(1)——水题

例题一：Relatives(pku 2407)

给出一个正整数n，小于n且和n互素的正整数有多少个

输入：输入数据有多组，每组数据占一行，输入n(n<=1 000 000 000),当n==0时结束

输出：对应每组数据，输出一行相应的结果。

Sample Input

7
12
0

Sample Output

6
4

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int phi(int n)
{
    int ans=n;
    for(int i=2;i*i<=n;i++)
        if(n%i==0)
    {
        ans=ans-ans/i;
        while(n%i==0)
            n/=i;
    }
    if(n>1)
        ans=ans-ans/n;
    return ans;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)
            break;
        cout<<phi(n)<<endl;
    }
    return 0;
}

例题二：Primitive Roots(原根)(pku 1284)

一个整数x(0<x<p)是奇数p的原根，当且仅当集合{ (x ⁱ mod p) | 1 <= i <= p-1 }与集合{ 1, ..., p-1 }是相同的。例如，3的连续次幂对7取模的结果是3，2，6，4，5，1，所以3是7的原根。

给出一个奇数p(3 <= p < 65536),编程求出p的原根的个数。

输入：输入数据有多组，每行一个素数p

输出：输出原根的个数

Sample Input

23
31
79

Sample Output

10
8
24

分析：本题应用到一个结论：p是素数，则p有phi(p-1)个原根。(虽然不知道为什么，姑且认为它是对的吧)

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int phi(int n)
{
    int ans=n;
    for(int i=2;i*i<=n;i++)
        if(n%i==0)
    {
        ans=ans-ans/i;
        while(n%i==0)
            n/=i;
    }
    if(n>1)
    {
        ans=ans-ans/n;
    }
    return ans;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        cout<<phi(n-1)<<endl;
    }
    return 0;
}

例题三：Farey Sequence

法雷级数Fn(n>=2)由一系列不能约分的分数a/b(0< a < b <= n 且gcd(a,b) = 1)按递增的顺序排列组成

F2 = {1/2}
F3 = {1/3, 1/2, 2/3}
F4 = {1/4, 1/3, 1/2, 2/3, 3/4}
F5 = {1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5}

计算法雷级数Fn中有多少个分数

输入：输入数据有多组，每组测试数据占一行，输入n(2<=n<=1 000 000);当n=0时，输入结束。

输出：对应每组输入数据，输出占一行，输出法雷级数Fn中分数的个数。

Sample Input

2
3
4
5
0

Sample Output

1
3
5
9

分析：一个法雷序列Fn中数的个数就是分别与2，3，4，5.。。。。。n-1，n互素的数的个数和，继而就是求从2到n连续的欧拉函数值得和(以这些数为分母，求出与这些数互素的个数，进而相加即可)。根据n的取值范围，可以直接用递推的方法来求解。再打表处理前n项和即可。

#include<cstdio>
#include<iostream>
#include<cstring>
#define maxn 1000000
using namespace std;
long long phi[maxn+10];
void pretreatment()///预处理，进行打表
{
    int i,j;
    for(i=1;i<=maxn;i++)
        phi[i]=i;
    for(i=2;i<=maxn;i+=2)
        phi[i]/=2;
    for(i=3;i<=maxn;i+=2)
        if(phi[i]==i)
    {
        for(j=i;j<=maxn;j+=i)
            phi[j]=phi[j]/i*(i-1);
    }
    for(i=3;i<maxn+10;i++)
        phi[i]+=phi[i-1];
}
int main()
{
    int n;
    pretreatment();
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)
            break;
        cout<<phi[n]<<endl;
    }
    return 0;
}

例题四:

Visible Lattice Points(可见格点)(pku 3090)

题目描述：找出有多少点，使其与原点的连线不经过其他点；

输入：第一行输入一个整数t(1<=t<=1000),代表有t组测试数据。每组测试数据占一行，每行有一个整数N(1<=N<=1000),代表了区域的大小。

输出：对于每组测试数据占一行，输出包括：测试数据的编号  空格  区域的大小N  空格  该区域可见点数

Sample Input

4
2
4
5
231

Sample Output

1 2 5
2 4 13
3 5 21
4 231 32549

分析：采用递推地方式从1开始到n进行递推。 就是求从2到n连续的欧拉函数值得和(以这些数为分母，求出与这些数互素的个数，进而相加即可)，最后乘以2加1就是最终结果

#include<cstdio>
#include<iostream>
#include<cstring>
#define maxn 1010
using namespace std;
long long phi[maxn];
void pretreatment()
{
    int i,j;
    for(i=1;i<=maxn;i++)
    {
        if(i&1)
            phi[i]=i;
        else
            phi[i]=i/2;
    }
     for(i=3;i<=maxn;i+=2)
        if(phi[i]==i)
    {
        for(j=i;j<=maxn;j+=i)
            phi[j]=phi[j]/i*(i-1);
    }
    for(i=2;i<maxn;i++)
        phi[i]+=phi[i-1];
}
int main()
{
    int t,n;
    cin>>t;
    pretreatment();
    for(int i=1;i<=t;i++)
    {
        scanf("%d",&n);
        cout<<i<<" "<<n<<" "<<phi[n]*2+1<<endl;
    }
    return 0;
}