poj 1077 八数码 ida*

IDA*比较好,重复路径也无所谓了,A*存状态判重意义不大,不实用

IDA*就是限制搜索深度的dfs,多加了一个limit限制条件,可以得到最优解

八数码的一些解释http://wenku.baidu.com/view/37a1b7a4284ac850ad024269.html

别人的总结http://www.cnblogs.com/JMDWQ/category/383627.html

这篇写的挺简明http://blog.csdn.net/ysjjovo/article/details/6657008

#include<stdio.h>//2  3  4  1  5  x  7  6  8
#include<string.h>
#include<math.h>
#define size 3
int move[4][2]={{-1,0},{0,-1},{0,1},{1,0}};//上 左 右 下
char op[4]={'u','l','r','d'};
int map[size][size],map2[size*size],limit,path[100];
int flag,length;
//int goal_st[3][3]={{1,2,3},{4,5,6},{7,8,0}};
int goal[9][2]={{2,2},{0,0},{0,1},{0,2},{1,0},{1,1},{1,2},{2,0},{2,1}};//目标位置
int nixu(int a[size*size])//算逆序数,逆序数为偶数状态可达
{
	int i,j,ni;
	ni=0;
	for(i=0;i<size*size;i++)
	{
		if(a[i]==0)             //空格不算,只算八个数块的逆序
			continue;
		for(j=i+1;j<size*size;j++)
		{
			if(a[j]==0)
			continue;
			if(a[i]>a[j])
				ni++;
		}
	}
	if(ni%2==0)
		return 1;
	else
		return 0;
}
int hv(int a[][size])//估价函数,曼哈顿距离,小等于实际总步数
{
	int i,j,cost=0;
	for(i=0;i<size;i++)
	{
		for(j=0;j<size;j++)
		{
			int w=map[i][j];
			cost+=abs(i-goal[w][0])+abs(j-goal[w][1]);
		}
	}
	return cost;
}
void swap(int*a,int*b)
{
	int tmp;
	tmp=*a;
	*a=*b;
	*b=tmp;
}
void dfs(int sx,int sy,int len,int pre_move)//sx,sy是空格的位置
{
	int i,j,nx,ny;
	if(flag)
		return;
	int dv=hv(map);
	if(len==limit)
	{
		if(dv==0)
		{
			flag=1;
			length=len;
			return;
		}
		else
			return;
	}
	else if(len<limit)
	{
		if(dv==0)
		{
			flag=1;
			length=len;
			return;
		}
	}
	for(i=0;i<4;i++)
	{
		if(i+pre_move==3&&len>0)//不和上一次移动方向相反,对第二步以后而言
			continue;    
		nx=sx+move[i][0];
		ny=sy+move[i][1];
		if(0<=nx&&nx<3 && 0<=ny&&ny<3)
		{
			swap(&map[sx][sy],&map[nx][ny]);
			int p=hv(map);
			if(p+len<=limit&&!flag)
			{
				path[len]=i;
				dfs(nx,ny,len+1,i);
				if(flag)
					return;
			}
			swap(&map[sx][sy],&map[nx][ny]);
		}
	}
}
int main()
{
	int i,j,k,l,m,n,sx,sy;
	char c,g;
	flag=0;
	memset(path,-1,sizeof(path));
	i=0;
	while(i<9)
	{
		scanf("%c",&c);
		if(c==32)
			continue;
		if(c=='x')
		{
			map2[i]=0;     
			map[i/3][i%3]=0;
			sx=i/3;sy=i%3;
			i++;
		}
		else
		{
			map2[i]=c-48;
			map[i/3][i%3]=c-48;
			i++;
		}
		
	}
	if(nixu(map2)==1)//该状态可达
	{
		limit=hv(map);
		limit=0;
		while(!flag)//找不到就加大limit,只要状态可达,就总会找到路径,每次搜索的深度都不大于limit,虽然可能错过正确解,但减少了错误路径的搜索长度,总的来说节省了时间
		{
			dfs(sx,sy,0,0);
			if(!flag)
				limit++; //得到的是最小步数
		}
		for(i=0;i<length;i++)
			printf("%c",op[path[i]]);
		printf("\n");
	}
	else
		printf("unsolvable\n");
	return 0;
}
16MS

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这是一道比较经典的计数问题。题目描述如下: 给定一个 $n \times n$ 的网格图,其中一些格子被标记为障碍。一个连通块是指一些被标记为障碍的格子的集合,满足这些格子在网格图中连通。一个格子是连通的当且仅当它与另一个被标记为障碍的格子在网格图中有公共边。 现在,你需要计算在这个网格图中,有多少个不同的连通块,满足这个连通块的大小(即包含的格子数)恰好为 $k$。 这是一道比较经典的计数问题,一般可以通过计算生成函数的方法来解决。具体来说,我们可以定义一个生成函数 $F(x)$,其中 $[x^k]F(x)$ 表示大小为 $k$ 的连通块的个数。那么,我们可以考虑如何计算这个生成函数。 对于一个大小为 $k$ 的连通块,我们可以考虑它的形状。具体来说,我们可以考虑以该连通块的最左边、最上边的格子为起点,从上到下、从左到右遍历该连通块,把每个格子在该连通块中的相对位置记录下来。由于该连通块的大小为 $k$,因此这些相对位置一定是 $(x,y) \in [0,n-1]^2$ 中的 $k$ 个不同点。 现在,我们需要考虑如何计算这些点对应的连通块是否合法。具体来说,我们可以考虑从左到右、从上到下依次处理这些点,对于每个点 $(x,y)$,我们需要考虑它是否能够与左边的点和上边的点连通。具体来说,如果 $(x-1,y)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们就是连通的;同样,如果 $(x,y-1)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们也是连通的。如果 $(x,y)$ 与左边和上边的点都不连通,那么说明这个点不属于该连通块。 考虑到每个点最多只有两个方向需要检查,因此时间复杂度为 $O(n^2 k)$。不过,我们可以使用类似于矩阵乘法的思想,将这个过程优化到 $O(k^3)$ 的时间复杂度。 具体来说,我们可以设 $f_{i,j,k}$ 表示状态 $(i,j)$ 所代表的点在连通块中,且连通块的大小为 $k$ 的方案数。显然,对于一个合法的 $(i,j,k)$,我们可以考虑 $(i-1,j,k-1)$ 和 $(i,j-1,k-1)$ 这两个状态,然后把点 $(i,j)$ 加入到它们所代表的连通块中。因此,我们可以设计一个 $O(k^3)$ 的 DP 状态转移,计算 $f_{i,j,k}$。 具体来说,我们可以考虑枚举连通块所包含的最右边和最下边的格子的坐标 $(x,y)$,然后计算 $f_{x,y,k}$。对于一个合法的 $(x,y,k)$,我们可以考虑将 $(x,y)$ 所代表的点加入到 $(x-1,y,k-1)$ 和 $(x,y-1,k-1)$ 所代表的连通块中。不过,这里需要注意一个细节:如果 $(x-1,y)$ 和 $(x,y)$ 在网格图中没有相邻边,那么它们不能算作连通的。因此,我们需要特判这个情况。 最终,$f_{n,n,k}$ 就是大小为 $k$ 的连通块的个数,时间复杂度为 $O(n^2 k + k^3)$。 参考代码:

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值