POJ 3518 : Prime Gap - 素数筛

题意：

长度为n素数间隔：两个相邻素数p和p+n之间，n-1个连续合数组成的序列

给出一个正整数k，编写一个程序，计算包含k的素数间隔的长度。若没有包含k的素数间隔，则长度为0。

分析：

令ans[k]为 包含k的素数间隔的长度。

显然，若k是素数，则ans[k]=0；若k为合数且k在p1与p2中间，则k所在合数区间内每个合数的ans值都为同一个数。

ans[p1+1]=ans[p1+2]=……=ans[p2-1]=(p2-1)-(p1+1)+2

步骤：

1）用筛选法计算[2…1299709]的素数筛u[]

2）顺序枚举[2…max n-1]中每个数i：若i为素数，则ans[i]=0；否则计算i右邻的素数j（u[i]=…=u[j-1]=false,u[j]=true），置ans[i]=…=ans[j-1]=j-i+2，并设i=j，以提高枚举效率。

输入k，ans[k]就是应该输出的答案了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn=1299710;//整数值的上限
bool u[maxn];//素数筛
int ans[maxn];//即答案的数组了

void prepare(){
	int i,j,k;
	for(i=2;i<maxn;i++)
		u[i]=1;//使用筛法计算[2…1299709]的素数筛u[]
	for(i=2;i<maxn;i++)
		if(u[i])//若i是素数，将其倍数从筛中筛出
			for(j=2;j*i<maxn;j++)
				u[i*j]=0;
	for(i=2;i<maxn;i++){//若i是合数，计算合数区间[i，j]
		if(!u[i]){
			j=i;
			while(j<maxn&&!u[j])
				j++;
			j--;
			for(k=i;k<=j;k++)
				ans[k]=j-i+2;
			i=j;
		}
		else ans[i]=0;//素数ans为0
	}
}

int main(){
	int i,j,k;
	prepare();
	while(~scanf("%d",&k)&&k){
		printf("%d\n",ans[k]);
	}
	return 0;
}