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 博客园(转载请注明出处)--------------------------------------------------------------- ...
随机推荐
- oracle生成随机数
UPDATE busi_chance_info t SET t.exp_amount = (SELECT floor(dbms_random.value(1, 10000000)) ...
- css 梯形标签页
html 代码 略 css : nav > a{ position: relative; display: inline_block; padding: .3em 1em 0; } nav &g ...
- java集合使用——HashMap
在map中插入.删除和定位元素时,HashMap是最好的选择.如果要按照自然顺序或自定义顺序遍历(获取所有元素),那么treemap更好一些. 第一:构造和添加元素 HashMap map = new ...
- 火狐浏览器,hostadmin hosts文件访问权限不足
开始->附件->以管理员身份运行. cacls %windir%\system32\drivers\etc\hosts /E /G Users:W
- javascript中写不写$(function() {});的区别
原地址 $(document).ready() 里的代码是在页面内容都加载完才执行的,如果把代码直接写到script标签里,当页面加载完这个script标签就会执行里边的代码了,此时如果你标签里执行的 ...
- 关于 jQuery中 function( window, undefined ) 写法的原因
今天在读 jQuery 源码的时候,发现下面的写法: (function(window,undefined){ ...// code goes here })(window); window 作为参数 ...
- mac Word 怎样放大缩小文档结构图文字大小
在文档结构图的侧栏里按住control+option,然后滑动鼠标滚轮/双指上下滚动触摸板.
- NET Core开发-获取所有注入(DI)服务
NET Core开发-获取所有注入(DI)服务 获取ASP.NET Core中所有注入(DI)服务,在ASP.NET Core中加入了Dependency Injection依赖注入. 我们在Cont ...
- linux----ulimit 限制
ulimit -a 显示当前用户的各种限制. ulimit -n 的数值表示每个进程可以打开的文件数目. 一般情况下, ulimit -n 的数值是1024. 当进程打开的文件数目超过此限制时,该进程 ...
- 编译安装zabbix
1. 说明:本例用源码包来安装zabbix,但是zabbix的后台数据库在这里选择mysql,然而mysql的安装方式不在选择源码了,而是选择已经编译好的通用linux包(tar包) zabbix源码 ...