a tutorial on computer science

  C++博客 :: 首页 :: 新随笔 :: 联系 :: 聚合  :: 管理 ::
  21 随笔 :: 0 文章 :: 17 评论 :: 0 Trackbacks
   据说不作此题人生不完整。好吧。很久以前就做过了,写过BFS,A*,和双搜。A*用了200+ms,汗,BFS都比他快。正好这几天在看搜索估价函数之类的东西,就把这道经典题拿出来,再做一遍,突然发现,估价函数+迭代加深搜索就是IDA*算法,好吧。以前傻傻看黑书的时候,理解不了A* ,觉得巨麻烦(现在也觉得挺麻烦),现在写起来IDA*,觉得还挺简洁,并且比较通用,而且这玩意又好写又比较通用,就详细研究了一下。看了别人的一个IDA*的算法,觉得写的很简洁很工整,就参详了一下,然后改造成了自己的,A掉了1077题。楼教主写的那个百度之星的版本的Allyes.com,还没有详细看,觉得有点复杂。有机会要好好研究下。
   我感觉IDA*的主要的价值在于,在深度搜索的时候,可以把问题的条件放宽,只要找到一个接近并且小于实际解的估价函数,就可以非常快的得到解了。而寻求估价函数的过程,非常的有意思,可以转换成各种经典问题。下面是代码,我觉得还算工整吧。
#include <cstdio>
#include <cstring>


char board[3][3];

int dest_goal[9][2] = {{0,0},{0,1},{0,2},{1,0},{1,1},{1,2},{2,0},{2,1},{2,2}}; 

int val[4][2] = {{0,1},{1,0},{-1,0},{0,-1}};  //r,l,d,u

char op[4] ={'r','d','u','l'};

char solution[100];

int find;

int testdata,low_bound;

int abs(int i)
{
  if(i>0) return i; return -i;
}

void swap(char& a,char& b)
{
  char c = a; a= b;b = c;
}

int getH()
{
   int i,j;
   int nRet = 0;
   for(i=0;i<3;i++)
    for(j=0;j<3;j++)
    {
      if(board[i][j] != 'x')
        nRet += abs(i - dest_goal[board[i][j] - '1'][0]) + abs(j - dest_goal[board[i][j] - '1'][1] );
    }
   return nRet;
}

int bound(int x,int y)
{
  if(x<0||y<0||x>=3||y>=3)  return 0;
  return 1;
}

int min(int a,int b)
{
    if(a>b) return b; return a; 
}

int dfs(int x,int y,int step,int maxstep)
{
   if(getH() == 0 || find == 1)
   {
     find = 1;
     return step;
   }

   if(step + getH() > maxstep)
     return step + getH();
  
   int i;
   int now_bound = 100000;
   for(i=0;i<4;i++)
   {
     int nx = x + val[i][0];
     int ny = y + val[i][1];
     if(bound(nx,ny))
     {
       swap(board[x][y],board[nx][ny]);
       solution[step] = op[i]; 
       int next_bound = dfs(nx,ny,step+1,maxstep);
       now_bound = min(now_bound,next_bound);
       if(find == 1)
         return now_bound;
       swap(board[x][y],board[nx][ny]);
     }
   }
   return now_bound;
}
int main()
{
   freopen("in_1077.txt","r",stdin);
   freopen("out.txt","w",stdout);
   int i,j,sx,sy;
   char c;
   for(i=0;i<3;i++)
    for(j=0;j<3;j++)
    {
       while(scanf("%c",&board[i][j]) && board[i][j] == ' ')
         NULL;
       if(board[i][j] == 'x')
       {
        sx = i;
        sy = j;
       }
    }
   
   find = 0;
   low_bound = getH();
   while(low_bound < 100 && !find) //如果启发函数比较好(很靠近结果并且小于结果),基本上很少次迭代就可以出结果。时间浪费在迭代次数上 
   {
     low_bound = dfs(sx,sy,0,low_bound);
   }
   
   if(find)
     {
       solution[low_bound] = '\0';
       printf("%s\n",solution);
     }
   else if(find == 0)
     printf("unsolvable\n");
   return 0;
}
posted on 2012-04-07 22:57 bigrabbit 阅读(3195) 评论(1)  编辑 收藏 引用

评论

# re: IDA*算法-POJ1077 2012-04-10 09:16 tb
学习算法了   回复  更多评论
  

# re: IDA*算法-POJ1077[未登录] 2013-09-07 12:35 Peter
搔年,copy也是要动脑子的,比如说这一句

while(low_bound < 100 && !find) //如果启发函


为什么你的最大上界要设置成 low_bound < 100 呢??  回复  更多评论
  


只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   博问   Chat2DB   管理