杭电ACM算法入门习题(素数与埃氏筛法)

素数判定(HDU 2012)

Problem Description

对于表达式n^2+n+41，当n在（x,y）范围内取整数值时（包括x,y）(-39<=x<y<=50)，判定该表达式的值是否都为素数。

Input

输入数据有多组，每组占一行，由两个整数x，y组成，当x=0,y=0时，表示输入结束，该行不做处理。

Output

对于每个给定范围内的取值，如果表达式的值都为素数，则输出"OK",否则请输出“Sorry”,每组输出占一行。

Sample Input

0 1
0 0

Sample Output

OK

解题思路：试除法模板题
AC代码：

#include <iostream>

using namespace std;

int cal_num(int n) { return n * n + n + 41; }

bool is_prime(int u)
{
    if (u <= 1) return false;//1既不是素数也不是合数
    for (int i = 2;i <= u / i;i++)//u / i防止int溢出
    {
        if (u % i == 0) return false;
    }
    return true;
}

int main()
{
    int x,y;
    while (cin >> x >> y)
    {
        if (x == 0 && y == 0) break;
        
        int flag = 1;//结果是素数
        for (int i = x;i <= y;i++)//枚举区间[x,y]中所有的数
        {
            if (!is_prime(cal_num(i)))
            {
                flag = 0;
                break;//区间内找到一个数字使得结果不为素数,结束循环
            }
        }
        
        if (flag == 0) cout << "Sorry" << endl;
        else cout << "OK" << endl;
    }
    return 0;
}

分拆素数和(HDU 2098)

Problem Description

把一个偶数拆成两个不同素数的和，有几种拆法呢？

Input

输入包含一些正的偶数，其值不会超过10000，个数不会超过500，若遇0，则结束。

Output

对应每个偶数，输出其拆成不同素数的个数，每个结果占一行。

Sample Input

30
26
0

Sample Output

3
2

解题思路：埃氏筛模板题。我们可以直接枚举质数p，n为输入的偶数，令x = n - p。

• 如果x也为质数，则x + p = n - p + p = n，符合题意，此时方案数加1。
• 如果x不是质数，则继续枚举。

需要注意两个细节：

• 题中要求将一个偶数拆分成两个不同的素数之和，故需要排除4 = 2 + 2这种情况，即需要使用语句p != x，防止将一个偶数拆分为两个一样的素数之和。
• 对于偶数10，可以将其拆分为10 = 3 + 7以及10 = 7 + 3，计算机会认为这是两种情况，但其实是一种结果，故需要将总结果除以2或者只枚举区间[2,n / 2]中的素数p即可。如12 = 5 +7或者12 = 7 + 5，其实在枚举质数p的时候，只需要枚举到12 / 2 = 6即可。

AC代码：

#include <iostream>

using namespace std;
const int N = 1e4 + 5;
int prime[N];//prime数组用来存储所有的质数
int vis[N];//vis[i] = 0表示i是质数,vis[i] = 1表示i是合数
int pcnt;//素数的个数

void init()//埃氏筛法模板
{
    vis[0] = vis[1] = 1;//0、1不是素数
    for (int i = 2;i < N;i++)
    {
        if (vis[i] == 0)//i是质数
        {
            prime[pcnt++] = i;//将i存储在数组中
            for (int j = i + i;j < N;j += i)//质数的倍数是合数,筛掉
            {
                vis[j] = 1;
            }
        }
    }
}

int main()
{
    int n;
    init();
    while (cin >> n && n)
    {
        //if (!n) break;
        
        int ans = 0;//表示满足题意的方案数
        for (int i = 0;i < pcnt && prime[i] <= n / 2;i++)//枚举n/2之前所有的质数
        {
            int x = n - prime[i];
            if (vis[x] == 0 && x != prime[i])//拆分为不同的质数
            {
                ans ++;
            }
        }
        
        cout << ans << endl;
    }
    return 0;
}