PTA 7-43 h0044.连续素数的和

PTA 7-43 h0044.连续素数的和

---

题目描述：

一些正整数能够表示为一个或多个连续素数的和。给出一个正整数，有多少个这样的表示？例如，整数53有两个表示：5+7+11+13+17和53；整数41有三个表示：2+3+5+7+11+13，11+13+17和41；整数3只有一个表示：3；整数20没有这样的表示。注意加法操作数必须是连续的素数，因此，对于整数20，7+13 和 3+5+5+7 都不是有效的表示。
请写一个程序，对于一个给出的正整数，程序给出连续素数的和的表示数。

输入格式：

输入一个正整数序列，每个数一行，在2到10000之间取值。输入结束以0表示。

输出格式：

输出的每一行对应输入的每一行，除了最后的0。输出的每一行对于一个输入的正整数，给出连续素数的和的表示数。输出中没有其他的字符。

输入样例：

2
3
17
41
20
666
12
53
0

输出样例：

1
1
2
3
0
0
1
2

---

思路：

首先用 欧拉筛 （用最小质因子去筛合数）找出所有的素数并存储下来：（欧拉筛法详情见：质数与合数 ）

当 i % primes[j] != 0 时，说明此时遍历到的 primes[j] 不是 i 的质因子，只可能是此时 primes[j] 的最小质因子，所以 primes[j] * i 的最小质因子就是 primes[j]。

当有 i % primes[j] == 0 时，因为我们是从小到大遍历的，说明此时的 prime[j]是满足条件的第一个数，即找到了primes[j] 就是 i 的最小质因子，因此 primes[j] * i 的最小质因子也就是 primes[j]，之后用 st[primes[j + 1] * i] = true 去筛合数时，就不是用最小质因子去更新了，所以此时应该退出循环，避免重复筛选。

接着遍历所有素数，找出匹配的连续素数串并记录答案。

---

代码如下：

#include <iostream>
#include <cstdio>
using namespace std;
const int N = 10010;

int primes[N]; //存素数
int cnt; //记录素数个数
bool st[N]; //判断数组，素数为false否则为true

void get_primes(int n)
{
    for (int i = 2; i <= n; i++){
        if (!st[i]) primes[cnt++] = i;
        for (int j = 0; primes[j] <= n / i; j++)
        {
            //标记，pj一定是pj*i的最小质因子
            st[primes[j] * i] = true;
            
            //从小到大遍历，如果 i%pj=0，则pj一定是i的最小公因子
            if (i % primes[j] == 0) break;
        }
    }
}

int main()
{
    int n;
    
    get_primes(10000);
    
    while (scanf("%d", &n) && n){
        int sum = 0; //存储范围内素数的和
        int ans = 0; //记录答案个数
        int l = 0; //标记指针
        
        //遍历所有素数
        for (int i = 0; i < cnt;){
            if (sum < n){
                sum += primes[i];
                i++;
            }
            else if (sum == n){
                ans++;
                sum = 0; //找到一种方案后归零，再继续查找
                i = ++l;
            }
            else if (sum > n){
                sum = 0;
                i = ++l;
            }
        }
            
        cout << ans << endl;
    }
    
    return 0;
}

---