sicily 1002 Anti-prime Sequence

这道题用深搜,其中有个地方需要注意的是:所求得的序列必须满足任意i(i = 1,2,3,...,k)个连续的数的和为合数。一开始自己写,过了,但是时间却是0.65s,效率不高,看了师兄师姐的代码才知道,有个卡时的方法:当深度搜索达到了一定的次数,可以认为继续沿着这条路径搜下去不会有答案了。于是可以用个计数dfs的次数来提前返回,用了卡时之后,效率大大提高,0s

/**
 * Problem: sicily 1002 Anti-not_prime Sequences
 * Description: 
 *    给定一个起始为m,终止为n的连续自然数的序列,重组这个序列,使得这个序列
 * 任意连续的i(i = 1,2,3,...,k)个数的和均为合数
 * Author: Wenzhu Luo @ Loen
 * Date: 12.07.2012
 */

 #include <iostream>
 #include <memory.h>
 #include <math.h>
 #include <cstdio>
 using namespace std;

 
 #define MAXN 10000
 bool not_prime[MAXN];	//用来标记是否为素数
 bool visited[1005]; //用来标记序列的某个数是否已经使用过
 int t;
 /**
  * Description: 筛法求素数
  * @para: None
  * @return: void
  */
 void Generate_prime()
 {
 	memset(not_prime, 0, sizeof(not_prime));//初始化所有数为素数
 	for (int i = 2; i < sqrt(MAXN); i++)
 	{
 		//如果一个数为素数,则这个数的倍数则不是素数
 		if (!not_prime[i])
 		{
 			for (int j = i * i; j < MAXN; j += i)
 				not_prime[j] = true;
 		}
 	}
 }

 /**
  * Description:
  *		给定一个起始为start,终止为end的连续自然数的序列,重组这个序列,使得
  * 这个序列任意连续的i(i = 1,2,3,...,k)个数的和均为合数,采用深度搜索的方法求解
  * @para: start-序列的起始,end-序列的终止,求得的解任意连续k个数的和为素数
  *        res数组-用来储存解,n当前解的长度
  * @return: 如果有解则返回true,如果无解返回false
  */
 bool dfs(int start, int end, int k,int res[], int n)
 {
 	t++;//用于卡时
 	if (t > 4000)//当次数大于4000时,认为搜不出答案了,这个数字根据题目规模而定,可以自己测试下
 		return false;
 	 if (n == end - start + 1)
 	 	return true;
 	 else
 	 {
 	 	for (int i = start; i <= end; i++)
 	 	{
 	 		if (!visited[i])
 	 		{
 	 			res[n] = i;
 	 			if (n == 0)
 	 			{
 	 				visited[i] = true;
 	 				if (dfs(start, end, k, res, n+1))
 	 					return true;
 	 				visited[i] = false;
 	 			}
 	 			else 
 	 			{
 	 				int sum = res[n];
 	 				bool flag = true;
 	 				for (int j = n - 1; j >= 0 && j >= n - k + 1; j--)
 	 				{
 	 					sum += res[j];
 	 					if (!not_prime[sum])
 	 					{
 	 						flag = false;
 	 					}
 	 				}
 	 				if (flag)
 	 				{
 	 					visited[i] = true;
 	 					if (dfs(start, end, k, res, n+1))
 	 						return true;
 	 					visited[i] = false;
 	 				}
 	 			}
 	 		}
 	 	}
 	 	return false;
 	 }
 	 
 }


 int main()
 {

 	int start,end,k;
 	Generate_prime();
    while (scanf ("%d%d%d",&start, &end, &k))
 	{
 		if (start == 0)
 			break;
 		memset(visited, false, sizeof(visited));
 		int res[1005];//用来存储解
 		t = 0;
 		//如果存在合法的序列,则输出结果
 		if (dfs(start, end, k, res, 0))
 		{
 			printf ("%d",res[0]);
 			
 			for (int i = 1; i < end - start + 1; i++)
 				printf (",%d",res[i]);
 			printf ("\n");
 		}
 		//否则输出不存在合法序列的提示信息
 		else
 			printf ("No anti-prime sequence exists.\n");
 	}
 	return 0;
 	
 	
 }


 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值