uva 1404 求很大数的素数算法

#include<cstdio>
#include<cstring>
#define MAX(x,y) ((x)>(y)?(x):(y))
using namespace std;
int prime[50050],cnt,vis[50050],ans[1000000];
bool flag[200000000];
void init(){
	for(int i=2;i<=50000;i++){
		if(!vis[i]){
			prime[cnt++]=i;
			for(int j=i+i;j<=50000;j+=i)
			vis[j]=1;
		}
	}
}
int solve(){
	int a,b,k,s;
	scanf("%d%d%d%d",&a,&b,&k,&s);
	if(a>b)	return 0;
	for(int i=0;i<b-a;i++)
	flag[i]=0;	
	for(int i=0;i<cnt&&prime[i]*prime[i]<=b;i++){
		int x=a/prime[i]*prime[i];
		if(x<a)
		x+=prime[i];
		if(x==prime[i])
		x+=prime[i];
		while(x>=a&&x<=b){
			flag[x-a]=1;
			x+=prime[i];
		}
	}
	int len=0;
	for(int i=0;i<=b-a;i++)
	if(!flag[i]){
		ans[len++]=i;
	}
	int res=0;
	for(int i=0;i+k-1<len;i++)
	if(ans[i+k-1]-ans[i]==s)
	res++;
	return res;
}
int main()
{
	init();
	int T;
	scanf("%d",&T);
	while(T--)
	printf("%d\n",solve());
	return 0;
}