文档章节

sicily 1002 Anti-prime Sequences

Ciel
 Ciel
发布于 2012/12/17 09:47
字数 737
阅读 64
收藏 0

Description

Given a sequence of consecutive integers n,n+1,n+2,...,m, an anti-prime sequence is a rearrangement of these integers so that each adjacent pair of integers sums to a composite (non-prime) number. For example, if n = 1 and m = 10, one such anti-prime sequence is 1,3,5,4,2,6,9,7,8,10. This is also the lexicographically first such sequence. We can extend the definition by defining a degree danti-prime sequence as one where all consecutive subsequences of length 2,3,...,d sum to a composite number. The sequence above is a degree 2 anti-prime sequence, but not a degree 3, since the subsequence 5, 4, 2 sums to 11. The lexicographically .rst degree 3 anti-prime sequence for these numbers is 1,3,5,4,6,2,10,8,7,9.

Input

Input will consist of multiple input sets. Each set will consist of three integers, n, m, and d on a single line. The values of n, m and d will satisfy 1 <= n < m <= 1000, and 2 <= d <= 10. The line 0 0 0 will indicate end of input and should not be processed.

Output

For each input set, output a single line consisting of a comma-separated list of integers forming a degree danti-prime sequence (do not insert any spaces and do not split the output over multiple lines). In the case where more than one anti-prime sequence exists, print the lexicographically first one (i.e., output the one with the lowest first value; in case of a tie, the lowest second value, etc.). In the case where no anti-prime sequence exists, output No anti-prime sequence exists.

Sample Input

1 10 2
1 10 3
1 10 5
40 60 7
0 0 0

Sample Output

1,3,5,4,2,6,9,7,8,10
1,3,5,4,6,2,10,8,7,9
No anti-prime sequence exists.
40,41,43,42,44,46,45,47,48,50,55,53,52,60,56,49,51,59,58,57,54

分析:

本题的思路很直白,因为数据量并不是非常大,只是1000的范围,长度为10,那么和最大也只有10000以内,所以可以暴力搜索,对每一个节点都进行判断并且迭代搜索,明显是用DFS。注意先用筛法做个素数表辅助判断,直接进行素数判断太麻烦,会超时。这里看到别人的代码,说道DFS卡时,其实就是限制DFS的总搜索次数,模糊判断用来节省时间。当然这样答案并不准确,但是对于数据规模较小的时候,可以取一些极限值测试,本题4000次搜索就足够了。

代码:

// Problem#: 1002
// Submission#: 1792032
// The source code is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License
// URI: http://creativecommons.org/licenses/by-nc-sa/3.0/
// All Copyright reserved by Informatic Lab of Sun Yat-sen University
#include<iostream>
#include<cstring>
using namespace std;

#define MAX 10000
#define N 1010

bool isPrime[MAX];
int prime[MAX];
bool visit[N];
int re[N];
int n,m,d,t,sum;
bool flag;

void primeList(){
    memset(isPrime,true,sizeof(isPrime));
    for(int i = 2;i <= MAX;++i){
        if(isPrime[i])  prime[++prime[0]] = i;
        for(int j = 1,k;j <= MAX && (k = i * prime[j]) <= MAX;++j)
        {
            isPrime[k] = false;
            if(i % prime[j] == 0)   break;
        }
    }
}

bool judge(int len){
    int sum;
    if(len <= 1)    return true;
    for(int i = 0;i < len;++i){
        sum = 0;
        if(i <= len - d){
            for(int j = i;j - i + 1<= d ;++j){
                sum += re[j];
                if(j-i+1 > 1 && isPrime[sum])   return false;
            }
        }else{
            for(int j = i;j < len;++j){
                sum += re[j];
                if(j-i+1 > 1 && isPrime[sum])   return false;
            }
        }
    }
    return true;
}
            
void dfs(int num){
    if(++t > 4000)  return;
    if(flag)    return;
    if(!judge(num)) return;
    if(num== m - n + 1){
        flag = 1;
        return;
    }
    for(int i = n;i <= m && !flag;++i){
        if(visit[i])    continue;
        re[num] = i;
        visit[i] = 1;
        dfs(num+1);
        visit[i] = 0;
    }
}

int main(){
    primeList();
    while(cin>>n>>m>>d && n && m && d){
        flag = false;
        t = 0;
        memset(visit,0,sizeof(visit));
        dfs(0);
        if(flag){
            cout << re[0];
            for( int i=1 ; i<m-n+1 ; ++i ) cout << "," << re[i];
            cout << endl;
        }else cout << "No anti-prime sequence exists." << endl;
    }
    return 0;
}

© 著作权归作者所有

上一篇: sicily 1004 I Conduit!
下一篇: sicily 1031 Campus
Ciel
粉丝 3
博文 24
码字总数 18532
作品 0
程序员
私信 提问
Oracle 11g R2 ADG 搭建

--============Oracle ADG搭建============== --==========准备阶段========= 1.检查primary为archivelog模式。 select log_mode from v$database; 如果为noarchivelog模式,切换到archivelo......

UltraSQL
2018/07/23
0
0
Python中如何使用yield,对于庞大迭代的优化处理

一直困扰于yield的使用,今天看到一篇不错的解释,虽然都是英文,不过不要紧,可以跳开,直接看代码的部分就能懂 Improve Your Python: 'yield' and Generators Explained Posted on Apr 07...

青鸾之旅
2013/08/02
1K
3
ACM Summer Training Warm up

ACM Summer Training Warm up Cover 热身水题 题目 HDU 4500 小Q系列故事——屌丝的逆袭 思路 简单的模拟,一个数组读入数据,一个数组计算维护结果 HDU 2109 Fighting for HDU 思路 简单排序...

SpiffyEight77
2017/08/14
0
0
用Groovy源编程(MOP)动态拦截(AOP)方法(比如记录String的concat和to...

实现AOP的方式有很多种,像Spring的AOP,它只能拦截Spring托管的bean;Groovy AST Transformations、ASM等在编译阶段通过修改字节码也可以做AOP;JAVA HOOK也可以做,但比较麻烦。 Groovy M...

飞飞狐
2012/09/29
155
0
oracle查询问题

两张表 A表 字段1 字段2 字段3 1001 2012-10-12 5 1002 2012-10-10 3 1002 2012-10-11 4 B表 字段1 字段2 字段3 1001 2012-10-11 -2 1002 2012-10-10 -2 1002 2012-10-13 -4 写一个SQL语句、......

mohism
2012/08/13
130
10

没有更多内容

加载失败,请刷新页面

加载更多

插入排序算法

《Java算法总纲目录》 1、定义     元素被分为有序区和无序区两部分。最初有序区只有一个元素。每次从无序区中选择一个元素,插入到有序区的位置,直到无序区变空。 2、代码 public c...

木九天
29分钟前
3
0
ApacheCN 翻译/校对/笔记整理活动进度公告 2019.10.18

注意 请贡献者查看参与方式,然后直接在 ISSUE 中认领。 翻译/校对三个文档就可以申请当负责人,我们会把你拉进合伙人群。翻译/校对五个文档的贡献者,可以申请实习证明。 请私聊片刻(52981...

ApacheCN_飞龙
31分钟前
4
0
Hands-on! 如何给 TiDB 添加新系统表

作者:黄东旭 “TiDB,你已经是一个成熟的数据库了,该学会用自己的 SQL 查自己的状态了。” 对于一个成熟的数据库来说,通过 SQL 来查询系统本身的状态再正常不过,对于 MySQL 来说 INFOMA...

TiDB
37分钟前
4
0
SpringBoot admin+Eureka+钉钉通知

SpringBoot admin+Eureka+钉钉通知 一、效果 登录账号+密码 监控服务 查看实时日志 钉钉通知 二、什么是Spring Boot Admin ? Spring Boot Admin是一个开源社区项目,用于管理和监控SpringB...

小白的成长
40分钟前
22
1
docker-rabbitmq

docker pull rabbitmqmkdir -p /rabbitmqdocker run -d \--name rabbitmq \--hostname rabbitmq \-v /rabbitmq:/var/lib/rabbitmq \-e RABBITMQ_DEFAULT_USER=root \-e RA......

李琼涛
42分钟前
5
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部