POJ2034 线性筛+DFS（下午第二节课给小裴和诗洋同学表演手打一遍过...）

cnyali

于 2016-09-23 16:05:23 发布

阅读量650

点赞数 1

CC 4.0 BY-SA版权

分类专栏：程序代码 POJ 数论相关

本文链接：https://blog.csdn.net/cnyali/article/details/52639280

程序代码同时被 3 个专栏收录

148 篇文章

订阅专栏

POJ

9 篇文章

订阅专栏

数论相关

7 篇文章

订阅专栏

本文介绍了一种生成特定条件下的素数序列算法，通过预先筛选出一定范围内的素数，并利用回溯法来寻找符合条件的素数序列组合，适用于解决特定数学问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
const int maxm=600000;
using namespace std;
bool ispr[600010],p[1010];
int pr[100000],e,s,t,d,list[1010];
void init(){
	for(int i=2;i<=maxm;i++){
		if(!ispr[i])pr[++e]=i;		
		for(int j=1;j<=e;j++){
			if(i*pr[j]>maxm)break;
			ispr[i*pr[j]]=1;			
			if(i%pr[j]==0)break; 
		}
	}
}
bool pd(int x,int y){
	int sum=y;
	if(x<2)return true;
	for(int i=x-1;i>=x-d+1 && i>=1;i--){
		sum+=list[i];
		if(!ispr[sum])return false;
	}
	return true;
}
bool dfs(int x){
	if(x>t-s+1){
		for(int i=1;i<x-1;i++)printf("%d,",list[i]);
		printf("%d\n",list[x-1]);
		return true;
	}		
	for(int i=s;i<=t;i++)
		if(!p[i] && pd(x,i)){
			list[x]=i;p[i]=1;
			if(dfs(x+1))return true;
			p[i]=0;
		}
	return false;
}
int main(){
	init();
	while(1){		
		memset(p,0,sizeof(p));
		scanf("%d%d%d",&s,&t,&d);
		if(s+t+d==0)break;
		if(dfs(1));
		else printf("No anti-prime sequence exists.\n");
	}	
	return 0;
}