zoukankan      html  css  js  c++  java
  • 双向BFS—>NOIP2002 字串变换

    如果目标也已知的话,用双向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;
    }



  • 相关阅读:
    第六周总结&第四次实验报告
    课程总结
    第十四周课程总结&实验报告
    第十三周课程总结
    第十二周课程总结
    第十一周课程总结
    第十周课程总结
    第九周课程总结&实验报告(七)
    第八周课程总结&实验报告(六)
    第七周课程总结&实验报告(五)
  • 原文地址:https://www.cnblogs.com/tham/p/6827166.html
Copyright © 2011-2022 走看看