这道题用深搜,其中有个地方需要注意的是:所求得的序列必须满足任意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;
}