(POJ - 2739)Sum of Consecutive Prime Numbers(尺取或二分)

题目链接：Sum of Consecutive Prime Numbers - POJ 2739 - Virtual Judge

中文题意：

小 X 非常喜欢素数，所以经常研究素数的性质，这天他想知道一个数字是否能用若干个（或许是一个）连续的素数之和表示，并且想知道有多少种方法。例如，53 有两种表示方法 5 + 7 + 11 + 13 + 17 和 53。其中素数v必须是连续的并且一个素数只能使用一次。 现在给你一个整数 n，你要告诉他这个数字有多少种不同的连续素数和表示方式。

Input

多组数据输入，每行一个数字 n，当 n=0 时结束。 n 在 2 与 10000之间。

Output

每个输出占一行，输出对应的答案。

Sample Input

2
3
17
41
20
666
12
53
0

Sample Output

1
1
2
3
0
0
1
2

分析：这道题目有两种做法，下面就两种做法分别给出分析

无论是哪种做法都需要预处理出来素数表，10000的素数表用埃氏筛预处理即可

第一种做法：尺取法

我们用l和r表示当前选中的素数和，若小于n，就使r加一，若大于n就使l减1，若当前选中的素数和等于n，就使得ans+1，这样就可以解决这道题目了，这种方法比较简单，复杂度是o(nlogn)（因为埃氏筛复杂度就是nlogn），也比较容易实现

第二种做法：二分

我们预处理出来素数的前缀和，每次固定右边界r，通过二分左边界来判断是否有满足题意的左边界使得l到r内的素数和为n，这样总体的复杂度也是o（logn），为了优化一点代码，我们可以提前二分出有边界的最大值，因为我们所选择的连续素数中的最大值一定是小于等于n的，所以我们就可以选第一个小于等于n的素数作为右边界，这个过程是用二分实现的，总体复杂度也是o(nlogn)，实现较尺取法复杂

下面是两种方法的代码：

尺取法：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<map>
#include<cmath>
#include<queue>
using namespace std;
const int N=10003;
int prime[N],cnt,s[N];
bool vis[N];
void init()//预处理出来素数
{
	for(int i=2;i<=10000;i++)
	{
		if(vis[i]) continue;
		prime[++cnt]=i;
		for(int j=2;j*i<=10000;j++)
			vis[i*j]=true;
	}
}
int main()
{
	init();
	int n;
	while(1)
	{
		scanf("%d",&n);
		if(!n) break;
		int ans=0,p=0;
		for(int i=1,j=1;i<=cnt&&j<=cnt;)
		{
			while(p<n&&j<=cnt)	p+=prime[j++];
			if(p==n) ans++,p-=prime[i++];
			else while(p>n) p-=prime[i++];
			if(p==n) ans++,p-=prime[i++];
			if(prime[j]>n) break;
		}
		printf("%d\n",ans);
	}
	return 0;
}

二分法：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<map>
#include<cmath>
#include<queue>
using namespace std;
const int N=10003;
int prime[N],cnt,s[N];
bool vis[N];
void init()//预处理出来素数前缀和 
{
	for(int i=2;i<=10000;i++)
	{
		if(vis[i]) continue;
		prime[++cnt]=i;
		for(int j=2;j*i<=10000;j++)
			vis[i*j]=true;
	}
	for(int i=1;i<=cnt;i++)
		s[i]=s[i-1]+prime[i];
}
bool check(int x,int rr)//看x能否被以第rr个素数结尾的连续素数和表示 
{
	int l=0,r=rr-1,mid;
	while(l<r)
	{
		mid=l+r+1>>1;
		if(s[rr]-s[mid]>=x) l=mid;
		else r=mid-1;
	}
	if(s[rr]-s[l]==x) return true;
	return false;
}
int main()
{
	init();
	int n;
	while(1)
	{
		scanf("%d",&n);
		if(!n) break;
		int l=1,r=cnt,mid;
		while(l<r)//二分找出小于等于x的第一个素数 
		{
			mid=l+r+1>>1;
			if(prime[mid]<=n) l=mid;
			else r=mid-1; 
		}
		int ans=0;
		for(int i=1;i<=r;i++)
			ans+=check(n,i);
		printf("%d\n",ans);
	}
	return 0;
}