ACM新手DAY 20 欧拉函数&线性筛

题解

A - Prime Gap

题目：读入一个n，是素数就输出0，不是就输出它左右两边最近的两个素数形成的区间的长度。

• 素数线性筛
• 思路：埃式筛法打表（素数表），用结合比较快的素数判定函数。
• 注意点：上限是1299709，可以设成const int maxn = 1299709+100;
• 判定素数函数部分

bool isprime(int n)
{
    if(n < 2)return 0;
    if(n ==2 || n == 3)return 1;
    if(n % 6 != 1 && n % 6 != 5)return 0;
    for(int i = 5; i * i <= n; i += 6)
    {
        if(n % i == 0 || n % (i + 2) == 0)
        return 0;
    }
    return 1;
}

• 打表部分（埃式）

    //构造素数表
    memset(bl, 0, sizeof bl);
    for(int i=2; i<maxn; ++i)
        for(int k=i*2; k<maxn; k+=i)
            bl[k]=1;
    int num=0;
    for(int i=2; i<maxn; ++i)//素数队列
        if(bl[i]==0)
            prime[num++]=i;

• 逻辑部分

int n;
    while(cin>>n)
    {
        if(n==0)
            break;
        if(isprime(n))
        {
            cout<<0<<endl;
            continue;
        }
        for(int i=0; i<num; ++i)
        {
            if(prime[i]<n&&prime[i+1]>n)
                cout<<(prime[i+1]-prime[i])<<endl;
        }
    }

B - Relatives

题目：求单个欧拉函数

• 求单个欧拉函数
• 样题

#include <iostream>
using namespace std;
int eular(int n)           //求单个数的欧拉函数
{
    int ans = n;
    for(int i = 2; i*i <= n; i++)
    {
        if(n % i == 0)
        {
            ans = ans/i*(i-1);       //ans=ans*(1-1/i)=ans*(i-1)/i
            //        =ans/i*(i-1) //先除防止溢出
            while(n % i == 0)      //消除i因子
                n /= i;
        }
    }
    if(n > 1)ans = ans/n*(n-1);    //最后可能还剩下一个质因数没有除
    return ans;
}
int main()
{
    int n;
    while(cin >> n && n)
    {
        cout << eular(n) << endl;
    }
    return 0;
}

C - Sum of Consecutive Prime Numbers

题目：读入一个正整数n，求有多少种方式可以把它表示成几个连续的素数的和。

• 素数线性筛（一个正整数转换多个连续的素数的和）
• 打表（埃式筛法）+暴力枚举
• 注意点：几个连续的素数的和
• 还用A题的打表模板，下面是逻辑部分

int n;
    while(cin>>n&&n)
    {
        int cnt(0);
        for(int i=0; prime[i]<=n; i++)
        {
            int sum(0);
            for(int j=i; ; j++)
            {
                if(sum == n)
                { cnt++; break; }
                else if(sum > n) break;
                else sum += prime[j];
            }
        }
        cout << cnt << endl;
    }

D - Farey Sequence

题目：读入正整数n，求1到n里面有几对互质数。（如果 $gcd(a,b)=1$，则称a，b为互质数）

• 欧拉函数
• 思路：

1. 分析：找到所有的互质数对，从2开始每一个数找它前面的存在的互质数，恰好符合欧拉函数的特点。
2. 做法：构造欧拉表（打表）

• 构造欧拉表（埃式；模板）。直接函数读入maxn的值就可以打表了。

void euler(int n)
{
    for (int i=1; i<=n; i++) phi[i]=i;
    for (int i=2; i<=n; i++)
    {
        if (phi[i]==i)//这代表i是质数
        {
            for (int j=i; j<=n; j+=i)
            {
                phi[j]=phi[j]/i*(i-1);   //把i的倍数更新掉
            }
        }
    }
}

• 另一种方法

    memset(a,0,sizeof a);
    //构造欧拉表
    for(int i=2; i<=mx; i++)
    {
        if(a[i]==0)
        {
            for(int j=i; j<=mx; j+=i)
            {
                if(a[j]==0)
                    a[j]=j;
                a[j]=a[j]/i*(i-1);
            }
        }
    }

E - Visible Lattice Points

题目；方格网中取$(n+1)\ast (n+1)$个坐标点，左下角为（0，0），求从（0，0）点可以看到几个点。（从原点看去，在同一条线上的点只能看见最前面的点）

• 多个数的欧拉函数
• 思路：

1. 分析：可以知道，可见的点的总数，等于按照对角线分开的其中一个三角形的可见点的二倍加上斜率为1的那条线上的一个点。对一个三角形分析可知，n增加的时候，点数增加的数对应着新的斜率。如图。
   (图片上有点错误，n=4时，不是两个1/3，而是一个1/3一个2/3)
2. 做法：欧拉表加输出调整；按照欧拉函数求那些斜率大于0小于1的线对应的点数的二倍，最后加3（一个斜率为0，一个斜率为1，一个没有斜率）

• 对应代码

//n*n的方格，从（0,0）点可以看到几个点
#include <iostream>
using namespace std;
const int mx = 1000;//开大了也可能出现run的问题问题
int phi[mx+10];
void euler(int n)
{
    for (int i=1; i<=n; i++) phi[i]=i;
    for (int i=2; i<=n; i++)
    {
        if (phi[i]==i)//这代表i是质数
        {
            for (int j=i; j<=n; j+=i)
            {
                phi[j]=phi[j]/i*(i-1);   //把i的倍数更新掉
            }
        }
    }
}
int main()
{
    int n, t;
    euler(mx);
    int cnt=1;
    cin >> t;
    for(; cnt<=t; cnt++)
    {
        cin >> n;
        int sum(0);
        for(int i=2; i<=n; i++)
            sum+=phi[i];
        //就是改了个输出，分析可知n的增加，对应着增加了1/n,2/n,3/n````(n-1)/n,而且
        //其中的假分数代表被挡到，看不见了，所以还是欧拉数列，每个n是前面那个数对应的欧拉值的和
        //再加上斜率为0和斜率为1的三个点
        cout << cnt << " " << n << " " <<sum*2+1+2<< endl;
    }
    return 0;
}

F - Almost Prime Numbers

题目：定义almost prime就是除了1和本身仅可以被一个素数整除的非素数，求一个给定区间里的almost prime数量。
如4，8就是符合条件的数

• 素数线性筛
• 思路：打表+暴力枚举
• 对应代码
• 注意点：这里的long long很重要

//自己的按照a题弄的，错在了没用ll，一直没发现
#include<iostream>
#include<string.h>
using namespace std;
#define ll long long
const int maxn = 1000005;
ll prime[maxn],vis[maxn];
int main()
{
    //构造素数表
    memset(vis, 0, sizeof vis);
    for(ll i=2; i<=maxn; i++)
        for(ll k=i*2; k<=maxn; k+=i)
            vis[k]=1;
    ll len=0;
    for(ll i=2; i<=maxn;i++)//素数队列
        if(vis[i]==0)
            prime[len++]=i;
    ll t;
    cin>>t;
    while (t--)
    {
        ll cnt = 0;
        ll l,r;
        cin>>l>>r;
        for (ll i = 0; i < len; i++)
            for ( ll j = prime[i]*prime[i]; j <= r; j*=prime[i] )
                if ( j >= l )
                    cnt++;
        cout<<cnt<<endl;
    }
    return 0;
}