迭代启发式搜索 IDA*

本章聚集了一些做了的迭代启发式搜索的题目

为什么只打了迭代启发式搜索?

因为它很好打,有些类似迭代的时候加的最优化剪枝

[因为这个最优化剪枝其实就是你算的估价函数了...]

BZOJ 1085 骑士精神

比较经典的一题,特别是题目中还说了在15步内走完...看上去就是迭代的样子

估价函数->指棋盘上不符合规定的骑士数目...每次移动最多让一名骑士回到自己的位置,由此来剪枝

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4
 5 using namespace std;
 6
 7 int maxd,sx,sy;
 8 char ch[5][5];
 9 int x[8]={1,1,-1,-1,2,2,-2,-2};
10 int y[8]={2,-2,2,-2,1,-1,1,-1};
11 int ans[5][5]={{1,1,1,1,1},{0,1,1,1,1},{0,0,2,1,1},{0,0,0,0,1},{0,0,0,0,0}};
12
13 struct Node{
14     int a[5][5];
15     int illegal(){
16     int cnt=0;
17         for(int i=0;i<5;i++)
18             for(int j=0;j<5;j++)
19                 if(a[i][j]!=ans[i][j]) cnt++;
20         return cnt;
21     }
22 };
23
24 inline bool judge(int a[5][5]){
25     for(int i=0;i<5;i++)
26         for(int j=0;j<5;j++)
27             if(a[i][j]!=ans[i][j]) return false;
28     return true;
29 }
30
31 inline bool inmap(int x,int y){
32     if(x>=5 || x<0 || y>=5 || y<0) return false;
33     return true;
34 }
35
36 bool dfs(int level,Node s,int bx,int by){
37     if(!level){
38         if(judge(s.a)) return true;
39         return false;
40     }
41     int newx,newy;
42     for(int i=0;i<8;i++){
43         newx=bx+x[i],newy=by+y[i];
44         if(inmap(newx,newy)){
45             swap(s.a[bx][by],s.a[newx][newy]);
46             if(s.illegal()<=level)
47                 if(dfs(level-1,s,newx,newy)) return true;
48             swap(s.a[bx][by],s.a[newx][newy]);
49         }
50     }
51     return false;
52 }
53
54 int main(){
55 #ifndef ONLINE_JUDGE
56     freopen("1085.in","r",stdin);
57     freopen("1085.out","w",stdout);
58 #endif
59     int kase;
60     Node s;
61     scanf("%d",&kase);
62     while(kase--){
63         for(int i=0;i<5;i++){
64             scanf("%s",ch[i]);
65             for(int j=0;j<5;j++)
66             if(ch[i][j]==‘*‘) s.a[i][j]=2,sx=i,sy=j;
67             else if(ch[i][j]==‘1‘) s.a[i][j]=1;
68             else s.a[i][j]=0;
69         }
70         if(!s.illegal()) {puts("0");continue;}
71         for(maxd=1;maxd<=15;maxd++)
72             if(dfs(maxd,s,sx,sy)) break;
73         if(maxd<=15)
74             printf("%d\n",maxd);
75         else
76             puts("-1");
77     }
78     return 0;
79 }

codevs1288 埃及分数

强行把从某个状态之后的最小值看做估价函数也是可以的...紫书例题

不过有一些坑点,看codevs题解的第一篇吧(...我写的)。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4
 5 using namespace std;
 6
 7 typedef long long ll;
 8
 9 const int maxn=10010;
10 const int INF=0x7f7f7f7f;
11
12 int a,b,d;
13 ll res[maxn],rec[maxn];
14
15 ll gcd(ll a,ll b){
16     ll t;
17     while(b)
18         t=b,b=a%b,a=t;
19     return a;
20 }
21
22 inline bool better(ll a[]){
23     if(a[1]>=res[1]) return false;
24     else return true;
25     for(int i=1;i<=d;i++)
26         if(a[i]<res[i]) return true;
27         else if(a[i]>res[i]) return false;
28     return false;
29 }
30
31 bool dfs(int level,ll last,ll aa,ll bb){
32     if(level==1){
33         if(bb%aa) return false;
34         rec[1]=bb/aa;
35         if(better(rec))
36             for(int i=1;i<=d;i++) res[i]=rec[i];
37         return true;
38     }
39     bool ok=false;
40     last=max(last+1,(bb/aa)+!(bb%aa==0));
41     for(int i=last;;i++){
42         if(aa*i>=bb*level) break;
43         ll b1=bb*i,a1=aa*i-bb,Gcd;
44         Gcd=gcd(a1,b1);
45         a1/=Gcd,b1/=Gcd;rec[level]=i;
46         if(dfs(level-1,i,a1,b1)) ok=true;
47     }
48     return ok;
49 }
50
51 int main(){
52     scanf("%d%d",&a,&b);
53     memset(res,0x7f,sizeof(res));
54     for(d=1;!dfs(d,0,a,b);d++);
55     for(int i=d;i>=1;i--)
56         printf("%lld ",res[i]);
57
58     return 0;
59 }

codevs2495 水叮当的舞步

这题估价函数居然就是一个颜色?...我其实有想过什么连通块个数的估价然而没有什么用...

其实这题关键还是在于给左上角染色的过程,每次给变色的地方染1,变色的周围染2,然后每次的颜色只要扩展2中相同的就好...比每次dfs快上很多的样子...

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4
 5 using namespace std;
 6
 7 int maxd,n;
 8 int mp[9][9],mark[9][9];
 9 int xx[4]={1,-1,0,0};
10 int yy[4]={0,0,-1,1};
11 bool used[6];
12 bool ans;
13
14 int get(){
15     int cnt=0;
16     memset(used,0,sizeof(used));
17     for(int i=1;i<=n;i++)
18         for(int j=1;j<=n;j++)
19             if(!used[mp[i][j]] && mark[i][j]!=1)
20                 used[mp[i][j]]=1,cnt++;
21     return cnt;
22 }
23
24 void dfs(int a,int b,int x){
25     mark[a][b]=1;
26     for(int i=0;i<4;i++){
27         int nowx=a+xx[i],nowy=b+yy[i];
28         if(nowx<1||nowy<1||nowx>n||nowy>n||mark[nowx][nowy]==1) continue;
29         mark[nowx][nowy]=2;
30         if(mp[nowx][nowy]==x)dfs(nowx,nowy,x);
31     }
32 }
33
34 int fill(int x){
35     int cnt=0;
36     for(int i=1;i<=n;i++)
37         for(int j=1;j<=n;j++)
38             if(mark[i][j]==2&&mp[i][j]==x)
39                 cnt++,dfs(i,j,x);
40     return cnt;
41 }
42
43 bool search(int level){
44     int v=get();
45     if(!v) return true;
46     if(v>level) return false;
47     int temp[9][9];
48     for(int i=0;i<=5;i++){
49         memcpy(temp,mark,sizeof(mark));
50         if(fill(i)) if(search(level-1)) return true;
51         memcpy(mark,temp,sizeof(mark));
52     }
53     return false;
54 }
55
56 int main(){
57
58     while(~scanf("%d",&n)){
59         if(!n) break;
60         memset(mark,0,sizeof(mark));
61         ans=0;
62         for(int i=1;i<=n;i++)
63             for(int j=1;j<=n;j++)
64                 scanf("%d",&mp[i][j]);
65         dfs(1,1,mp[1][1]);
66         for(maxd=0;!search(maxd);maxd++);
67         printf("%d\n",maxd);
68     }
69     return 0;
70 }

UVa11212 编辑书稿

估价函数比较精妙,每次最多让三个元素的后继正确。

为什么是后继呢?因为这个操作的本质就是改变后继...

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cstdlib>
 4 #include<algorithm>
 5
 6 using namespace std;
 7
 8 const int maxn=12;
 9
10 struct Node{
11     int v[maxn];
12 };
13
14 int n,maxd;
15
16 inline int get_h(int a[]){
17     int cnt=0;
18     for(int i=1;i<n;i++)
19         if(a[i+1]!=a[i]+1) cnt++;
20     return cnt;
21 }
22
23 void move(Node &Start,Node &Next,int I,int J,int K){
24     for(int i=0;i<=K-J-1;i++)
25         Next.v[i+I]=Start.v[J+1+i];
26     for(int i=0;i<=J-I;i++)
27         Next.v[K-J+I+i]=Start.v[i+I];
28 }
29
30 bool dfs(int level,Node s){
31     if(get_h(s.v)>3*level) return false;
32     if(get_h(s.v)==0) return true;
33
34     Node next;
35     for(int i=1;i<=n;i++)
36         for(int j=i;j<=n;j++)
37             for(int k=j+1;k<=n;k++){
38                 memcpy(next.v,s.v,sizeof(s.v));
39                 move(s,next,i,j,k);
40                 if(dfs(level-1,next)) return true;
41             }
42     return false;
43 }
44
45 int main(){
46 #ifndef     ONLINE_JUDGE
47     freopen("11212.in","r",stdin);
48     freopen("11212.out","w",stdout);
49 #endif
50
51     int kase=0;
52     Node a0;
53
54     while(~scanf("%d",&n)){
55         if(!n) break;
56         for(int i=1;i<=n;i++)
57             scanf("%d",&a0.v[i]);
58         if(!get_h(a0.v)){
59             printf("Case %d: 0\n",++kase);
60             continue;
61         }
62
63         for(maxd=1;!dfs(maxd,a0);maxd++);
64         printf("Case %d: %d\n",++kase,maxd);
65     }
66
67     return 0;
68 }

时间: 2024-08-24 19:03:57

迭代启发式搜索 IDA*的相关文章

【启发式搜索】A*与IDA*学习笔记

搞了这么久发现自己到现在还不会启发式搜索ヾ(?`Д′?)所以今天正好趁着搜索练习题的风去搞了启发式搜索 A*搜索算法,俗称A星算法.这是一种在图形平面上,有多个节点的路径,求出最低通过成本的算法.常用于游戏中的NPC的移动计算,或在线游戏的BOT的移动计算上. 该算法像Dijkstra算法一样,可以找到一条最短路径:也像BFS一样,进行启发式的搜索. 在此算法中,如果以 g(n)表示从起点到任意顶点n的实际距离,h(n)表示任意顶点n到目标顶点的估算距离,那么 A*算法的公式为:f(n)=g(n

uva 11212 - Editing a Book(迭代加深搜索 IDA*) 迭代加深搜索

迭代加深搜索 自己看的时候第一遍更本就看不懂..是很水,但智商捉急也是没有办法的事情. 好在有几个同学已经是做过了这道题并且对迭代加深搜索的思路有了一定的了解,所以在某些不理解的地方询问了一下他们的见解, 真的是很有帮助,也许自己想要想很久才能想明白,还会很痛苦,稍微问一下别人的想法,点上一个方向,剩下的自己就能想得明白了. 迭代加深. 把answer(需要的步数或其他)在主函数里面从零往上递加,此之谓 "层数",亦可谓之"深度".用书上的话就是: 从小到大枚举深度

UVA 11212 Editing a Book [迭代加深搜索IDA*]

11212 Editing a Book You have n equal-length paragraphs numbered 1 to n. Now you want to arrange them in the order of 1, 2, . . . , n. With the help of a clipboard, you can easily do this: Ctrl-X (cut) and Ctrl-V (paste) several times. You cannot cut

UVA-12558 Egyptian Fractions (HARD version) (IDA* 或 迭代加深搜索)

题目大意:经典的埃及分数问题. 代码如下: # include<iostream> # include<cstdio> # include<cstring> # include<algorithm> using namespace std; # define LL long long int num[5],a,b,k; LL ans[10000],v[10000]; LL gcd(LL a,LL b) { return (b==0)?a:gcd(b,a%b)

【迭代加深搜索】埃及分数问题

谢谢阿苏~http://blog.csdn.net/urecvbnkuhbh_54245df/article/details/5856756 [迭代加深搜索(ID,iterative deepening)]:从小到大枚举上限maxd,每次执行只考虑深度不超过maxd的结点. ------对于可以用回溯法求解但解答树的深度没有明显上限的题目,可以考虑ID算法: ------优点:它主要是在递归搜索函数的开头判断当前搜索的深度是否大于预定义的最大搜索深度,如果大于,就退出这一层的搜索,如果不大于,就

【日常学习】【IDA*】codevs2449 骑士精神题解

题目描述 Description 在一个5×5的棋盘上有12个白色的骑士和12个黑色的骑士, 且有一个空位.在任何时候一个骑士都能按照骑士的走法(它可以走到和它横坐标相差为1,纵坐标相差为2或者横坐标相差为2,纵坐标相差为1的格子)移动到空位上. 给定一个初始的棋盘,怎样才能经过移动变成如下目标棋盘: 为了体现出骑士精神,他们必须以最少的步数完成任务. 输入描述 Input Description 第一行有一个正整数T(T<=10),表示一共有N组数据.接下来有T个5×5的矩阵,0表示白色骑士,

A*和IDA*介绍

\(A*\)算法是一种很神奇的搜索方法,它属于启发式搜索中的一种.A最主要的功能当然就是用来剪枝,提高搜索的效率.A主要的实现方法是通过一个估价函数,每次对下一步进行一个估价,根据估价出的值来决定下一步往哪走.因此,**估价函数的好坏决定了A*算法的优劣**. 估价函数分为三种情况: 如果估计的距离小于到目标状态的实际距离,这种情况下,搜索的点数多,搜索范围大,效率低.但能得到最优解. 如果估计的距离等于最短距离,那么搜索将严格沿着最短路径进行, 此时的搜索效率是最高的. 如果估计的距离大于最短

人工智能08 启发式搜索

启发式搜索 [这一章在某些地方笔者自己也没完全弄清楚,比如在递归最优搜索处没有找到一个很好的例子来理解,比如如何选择启发式函数等等一系列的问题,希望有大神能指明讲解.所以本章重要的只是介绍A*算法流程和简单优化并介绍引出一些改进的A*算法] 使用评估函数 除了搜索过程不是从开始节点统一向外扩展外,本章描述的搜索过程有点像广度优先搜索,不同的是,它会优先顺着有启发性和具有特定信息的节点搜索下去,这些节点可能是到达目标的最好路径.我们称这个过程为最优(best-first)或启发式搜索.下面是其基本

八数码问题强化版:十五数码问题idA*版本

---恢复内容开始--- 上一次介绍过dbfs版本,这次来介绍idA*版本. 首先要理解idA*算法的思想,是将迭代加深与A*的结合,将估价函数h(n)作为迭代的限制值,进行dfs. (A*和迭代加深的介绍等有时间再写出来吧) 对所有点(除0以外的)进行曼哈顿距离计算(目标状态到初始状态),h(n)为当前节点的各点的曼哈顿距离和. 在代码中看: #include<cstdio> #include<algorithm> #include<cstring> #define