如果目标也已知的话,用双向BFS能很大提高速度
单向时,是 b^len的扩展。
双向的话,2*b^(len/2) 快了很多,特别是分支因子b较大时
至于实现上,网上有些做法是用两个队列,交替节点搜索 ×,如下面的伪代码:
while(!empty())
{
扩展正向一个节点
遇到反向已经扩展的return
扩展反向一个节点
遇到正向已经扩展的return
}
但这种做法是有问题的,如下面的图:
求S-T的最短路,交替节点搜索(一次正向节点,一次反向节点)时
Step 1 : S –> 1 , 2
Step 2 : T –> 3 , 4
Step 3 : 1 –> 5
Step 4 : 3 –> 5 返回最短路为4,错误的,事实是3,S-2-4-T
我想,正确做法的是交替逐层搜索,保证了不会先遇到非最优解就跳出,而是检查完该层所有节点,得到最优值。
也即如果该层搜索遇到了对方已经访问过的,那么已经搜索过的层数就是答案了,可以跳出了,以后不会更优的了。
当某一边队列空时就无解了。
优化:提供速度的关键在于使状态扩展得少一些,所以优先选择队列长度较少的去扩展,保持两边队列长度平衡。这比较适合于两边的扩展情况不同时,一边扩展得快,一边扩展得慢。如果两边扩展情况一样时,加了后效果不大,不过加了也没事。
无向图时,两边扩展情况类似。有向图时,注意反向的扩展是反过来的 x->y(NOIP2002G2字串变换)
/* 双向BFS: 1.正向搜索:从初始结点向目标结点方向搜索,按照正向规则(A$->B$)变换。 2.逆向搜索:从目标结点向初始结点方向搜索,按照逆向规则(B$->A$)变换。 当两个方向的搜索生成同一子结点时终止此搜索过程(变换的总步数为此时两个方向BFS的步数总和)。 双向搜索通常有两种方法: 1. 两个方向交替扩展。 2. 选择结点个数较少的那个方向先扩展。 方法2克服了两方向结点的生成速度不平衡的状态,明显提高了效率。本程序使用方法1,两个方向交替BFS,进行正反规则变换。 */ #include <stdio.h> #include <string.h> #include <stdlib.h> #include <iostream> using namespace std; struct Node { char str[41]; int sep; Node() { sep=0; memset(str,sizeof(str),0); } }q1[20000], q2[20000]; int h1, r1, h2, r2; //h1,r1...h2,r2,分别保存起始和目标两个状态的队列头和尾 char s1[6][25], s2[6][25]; //存储变换规则 int n,Min=100; //Min储存最少规则变换次数 void copy2(int start, int use) //逆向搜索,当q2[]搜索到逆向规则匹配的字串B$的时候,进行替换 { int i, j; r2++; //队列中元素增加1个字串 q2[r2].sep = q2[h2].sep + 1; for(i = 0; i < start; i++) //替换匹配规则的字串 { q2[r2].str[i] = q2[h2].str[i]; } for(j = 0; s1[use][j] != ''; j++, i++) { q2[r2].str[i] = s1[use][j]; } for(j = start + strlen(s2[use]); q2[h2].str[j] != ''; j++, i++) { q2[r2].str[i] = q2[h2].str[j]; } //cout<<"q2:"<<q2[r2].str<<endl; for(i = 0; i <= r1; i++) //判断这一次规则替换操作结束后,正向BFS队列中是否会有字串元素与之匹配 { if(strcmp(q1[i].str, q2[r2].str) == 0) { printf("%d ", q1[i].sep + q2[r2].sep); exit(0); // if(q1[i].sep+q2[r2].sep<Min) Min=q1[i].sep+q2[r2].sep; } } } void copy1(int start, int use) //正向搜索,当q1[]搜索到与规则匹配的字串A$的时候,进行替换 { int i, j; r1++; q1[r1].sep = q1[h1].sep + 1; //队列q1[]中增加1个新的元素(一次规则变换后的字串) for(i = 0; i < start; i++) //规则替换操作 { q1[r1].str[i] = q1[h1].str[i]; } for(j = 0; s2[use][j] != ''; j++, i++) { q1[r1].str[i] = s2[use][j]; } for(j = start + strlen(s1[use]); q1[h1].str[j] != ''; j++, i++) { q1[r1].str[i] = q1[h1].str[j]; } //cout<<"q1:"<<q1[r1].str<<endl; for(i = 0; i <= r2; i++) //判断这一次规则替换操作结束后,方向BFS队列中是否会有字串元素与之匹配 { if(strcmp(q2[i].str, q1[r1].str) == 0) { printf("%d ", q2[i].sep + q1[r1].sep); exit(0); // if(q2[i].sep+q1[r1].sep<Min) Min=q2[i].sep+q1[r1].sep; } } } void work(void) { int i, j; while(h1 <= r1 && h2 <= r2) //搜索过程中确保没有一个队列为空,否则搜索不到相交的情况 { if(q1[h1].sep + q2[h2].sep > 10) //正反搜索的步数总和超过了10,说明这样的转换至少要超过10次才能实现,结束 { printf("NO ANSWER! "); exit(0); } for(i = 0; i < strlen(q1[h1].str); i++) { for(j = 0; j < n; j++) //正向搜索,一共n个变换规则 { if(strncmp(s1[j], &q1[h1].str[i], strlen(s1[j])) == 0) { copy1(i, j); } } } h1++; //正向一遍BFS,搜索完所有规则之后,队首元素出队 for(i = 0; i < strlen(q2[h2].str); i++) //加快搜索的速度,同理从目标开始,方向,并根据逆向规则进行BFS { for(j = 0; j < n; j++) { if(strncmp(s2[j], &q2[h2].str[i], strlen(s2[j])) == 0) { copy2(i, j); } } } h2++; } } int main() { freopen("in.in", "r", stdin); // freopen("string.out", "w", stdout); scanf("%s%s", q1[0].str, q2[0].str); while(scanf("%s%s", s1[n], s2[n]) == 2) { n++; } work(); printf("NO ANSWER! "); return 0; }