USACO 1.3... 虫洞 解题报告(搜索+强大剪枝+模拟)
这题可真是又让我找到了八数码的感觉。。。哈哈。
首先,第一次见题,没有思路,第二次看题,感觉是搜索,就这样写下来了。
这题我几乎是一个点一个点改对的(至于为什么是这样,后面给你看一个神奇的东西),让我发现了许多搜索上的问题。思路非常简
单:搜索出每一种可能的配对方式,然后从每一个点出发
一遍,模拟走的过程,如果到死循环里面就答案加一,那个模拟走的过程话说我敲的还是很爽的,这题难在对配对序列的搜索和配对状态重
复的剪枝上面。
记得配对序列的搜索我敲了5、6个版本,发现AC之后一个都没有用上,但是在这个过程中,你脑袋里面模拟那个搜索的过程的确是很爽
啊,确定了搜索序列的方法之后,就分别对每一个序列执行模拟走的操作,但是,这是一个严重的问题出现了,大量的重叠状态,导致所有
N = 12的数据严重超时,让人更加不理解的是,USACO上原题还有小于12的数呢,但是COGS上的数据好像全部都是12,坑了。但是,剪
枝的方法十分的妙,什么,单调。我给出一个不重复的序列作为说明:
当N等于6时,所有可能的序列如下:
12 34 56 | 12 35 46 | 12 36 45 |
13 24 56 | 13 25 46 | 13 26 45 |
14 23 56 | 14 25 36 | 14 26 35 |
15 23 46 | 15 24 36 | 15 26 34 |
16 24 35 | 16 23 45 | 16 25 34 |
我们可能清楚的发现,如果你的搜索顺序对的话,那么有一个神奇而又能让你通过测试的性质,就是这个序列无论什么样子,他的奇数位上
的数字都是单调递增的,正是因为这个东西,让我剪掉了数量非常大的无用状态。
至于怎么体现这句话,关键地方说完了,直接贴代码,自己看吧。。。有许多大神0.0几秒就过了,不能理解,我的0.6几秒,这效率,天
差地
别啊,,,自愧不如。。。。。
Code:
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <set>
#include <map>
using namespace std; int N;
int Ans = 0;
bool vi[15];
bool a[15][15];
int que[15];
struct data{
int x,y;
}G[15]; int T[30];
int tot = 0; void dfs(int,int,int);
bool cmp(data,data);
bool judge();
bool go(int,int);
void outandin(); int main(){
freopen("wormhole.in","r",stdin);
freopen("wormhole.out","w",stdout);
scanf("%d",&N);
for(int i = 1;i <= N;++ i){
scanf("%d%d",&G[i].x,&G[i].y);
T[++ tot] = G[i].x;
T[++ tot] = G[i].y;
}
outandin();
sort(G+1,G+N+1,cmp); for(int i = 2;i <= N;++ i){
dfs(1,i,2);
}
printf("%d\n",Ans);
return 0;
} void dfs(int now,int next,int depth){
vi[now] = vi[next] = true;
que[depth] = now;que[depth+1] = next;
if(depth == N && N%2==0){
if(judge())
++ Ans;
return;
}
else if(depth == N-1 && N%2 == 1){
if(judge())
++ Ans;
return;
} for(int i = 1;i <= N;++ i){
if(!vi[i] && i > now){//这括号里面很神奇的一句。。。i > now 剪枝的神器。。
for(int j = 1;j <= N;++ j){
if(!vi[j] && j != i && j > i){
dfs(i,j,depth+2);
vi[j] = false;
}
}
vi[i] = false;
}
}
vi[now] = vi[next] = false;
} bool cmp(data a,data b){
if(a.x == b.x)
return a.y < b.y;
return a.x < b.x;
} bool judge(){
for(int i = 2;i <= (N&1?N:N+1);++ i)
if(go(que[i],i))
return true;
return false;
} bool go(int now,int pos){
bool v[15];
memset(v,false,sizeof v); while(1){
if(v[now])
return true;
v[now] = true;
pos = ((pos&1) ? (pos-1):(pos+1));
bool flag = false;
for(int i = 1;i <= N;++ i){
int n = que[pos];
if(i != n && G[i].x > G[n].x && G[i].y == G[n].y){
for(int j = 2;j <= (N&1?N:N+1);++ j){
if(que[j] == i){
now = que[j];
pos = j;
flag = true;
break;
}
}
}
if(flag)
break;
}
if(!flag)
return false;
}
} map <int,int> pos;
map <int,bool> fla;
void outandin(){
int no = 0;
for(int i = 1;i <= tot;++ i){
if(!fla[T[i]]){
fla[T[i]] = true;
pos[T[i]] = ++ no;
}
}
for(int i = 1;i <= N;++ i){
G[i].x = pos[G[i].x];
G[i].y = pos[G[i].y];
}
} 如果你不剪枝,但是时间又特别特别长的话(也就是没有超时这种东西),你的输出就要这样写:
if(N > 11)
printf("%d\n",Ans/120);
else if(N > 9)
printf("%d\n",Ans/24);
else if(N > 7)
printf("%d\n",Ans/6);
else if(N > 5)
printf("%d\n",Ans/2);
else
printf("%d\n",Ans);
正是这个让我一个点一个点的过,也是这个让我发现了大量的重复状态,由那个/120你就知道重复状态有多么可怕了。。。。。
总之,模拟走的过程练码力,搜索的过程练搜索。。一道非常不错的题目。。。
USACO 1.3... 虫洞 解题报告(搜索+强大剪枝+模拟)的更多相关文章
- USACO Section1.3 Wormholes 解题报告
wormhole解题报告 —— icedream61 博客园(转载请注明出处)------------------------------------------------------------- ...
- USACO Section1.2 Transformations 解题报告
transform解题报告 —— icedream61 博客园(转载请注明出处)------------------------------------------------------------ ...
- USACO 1.4 ariprog 解题报告
这是继虫洞之后又让我为难的一个 剪枝题目,无论如何,做的再快,也只能过6个点,最后三个点也TLE.后来参考了一下标答,大概思路是这样的. 朴素算法就不多说了,枚举a,b然后判断就行,网上说这样优化到位 ...
- USACO 6.3 章节 你对搜索和剪枝一无所知QAQ
emmm........很久很久以前 把6.2过了 所以emmmmmm 直接跳过 ,从6.1到6.3吧 Fence Rails 题目大意 N<=50个数A1,A2... 1023个数,每个数数值 ...
- USACO Section1.4 Mother's Milk 解题报告
milk3解题报告 —— icedream61 博客园(转载请注明出处)---------------------------------------------------------------- ...
- USACO Section2.2 Preface Numbering 解题报告 【icedream61】
preface解题报告----------------------------------------------------------------------------------------- ...
- USACO Section2.1 Hamming Codes 解题报告 【icedream61】
hamming解题报告----------------------------------------------------------------------------------------- ...
- USACO Section2.1 Healthy Holsteins 解题报告 【icedream61】
holstein解题报告 --------------------------------------------------------------------------------------- ...
- USACO Section2.1 The Castle 解题报告
castle解题报告 —— icedream61 博客园(转载请注明出处)--------------------------------------------------------------- ...
随机推荐
- xcode 7 添加空模板
文件下载链接: http://pan.baidu.com/s/1pKbyf4R 密码: ppi7 下载完成后将模板放入以下路径应用程序->Xcode->右键显示包内容->/Conte ...
- 使用Spring JDBCTemplate简化JDBC的操作
使用Spring JDBCTemplate简化JDBC的操作 接触过JAVA WEB开发的朋友肯定都知道Hibernate框架,虽然不否定它的强大之处,但个人对它一直无感,总感觉不够灵活,太过臃肿了. ...
- java集合使用——HashMap
在map中插入.删除和定位元素时,HashMap是最好的选择.如果要按照自然顺序或自定义顺序遍历(获取所有元素),那么treemap更好一些. 第一:构造和添加元素 HashMap map = new ...
- K - Ignatius and the Princess IV
Description "OK, you are not too bad, em... But you can never pass the next test." ...
- SeekBar和RatingBar
今天在看一个音乐播放器的源代码时候用到了SeekBar,就翻出来mars老师的视频复习了一下,然后综合使用了一下. 首先先看下运行效果: 下来我们看看布局文件的设计: main.xml: < ...
- TeamTalk自定义IM协议的理解
一.TeamTalk自定义IM协议 TeamTalk自定义IM协议是一种基于protocol buffer的消息传递协议,protocol buffer可以自定义消息格式.protocol buffe ...
- C语言基础10
栈区间:在函数内部声明的变量都存放在栈区间,比如int char 数组 结构体 指针,只管申请,系统会自动帮我们回收,收回的时间是作用域结束之后,遵循的原则是"先进后出". int ...
- js刷新页面
reload 方法,该方法强迫浏览器刷新当前页面. 语法:location.reload([bForceGet]) 参数: bForceGet, 可选参数, 默认为 false,从客户端缓存里取当前页 ...
- typecho路由机制详解
本文介绍的是typecho的路由机制,引自 不烦恼路由机制是typecho的核心,有很多功能都是基于路由功能设计的,理解并熟悉TE的路由机制将非常有助于插件的开发. 完整的路由表如下: array ( ...
- BZOJ 2440 完全平方数(莫比乌斯反演,容斥原理)
http://www.lydsy.com/JudgeOnline/problem.php?id=2440 题意:求第K个没有平方因子的数 思路:首先,可以二分数字,然后问题就转变成x以内有多少无平方因 ...