NYOJ 10 skiing (深搜和动归)

skiing

时间限制：3000 ms  |  内存限制：65535 KB

难度：5

描写叙述

Michael喜欢滑雪百这并不奇怪。 由于滑雪的确非常刺激。但是为了获得速度。滑的区域必须向下倾斜。并且当你滑到坡底，你不得不再次走上坡或者等待升降机来载你。Michael想知道载一个区域中最长底滑坡。区域由一个二维数组给出。数组的每一个数字代表点的高度。以下是一个样例 

1 2 3 4 5



16 17 18 19 6



15 24 25 20 7



14 23 22 21 8



13 12 11 10 9



一个人能够从某个点滑向上下左右相邻四个点之中的一个，当且仅当高度减小。

在上面的样例中，一条可滑行的滑坡为24-17-16-1。当然25-24-23-...-3-2-1更长。其实，这是最长的一条。

输入

第一行表示有几组測试数据，输入的第二行表示区域的行数R和列数C(1 <= R,C <= 100)。以下是R行，每行有C个整数，代表高度h。0<=h<=10000。

后面是下一组数据；

输出

输出最长区域的长度。

算法分析：

一看本题就感觉这不是才看过的单调递减子序列问题么。但这个难点在于数据并非线性的，每一步都有四个方向可走。于是乎。，深搜解此题必须是手到擒来啊。。注意中间的剪枝问题。利用记忆话数组去除反复的计算类似于递归型DP，以下看一个来自网上的深搜代码：

#include<stdio.h>
#include<string.h>
int a[101][101],visit[101][101];
int dx[4]={0,0,1,-1};
int dy[4]={1,-1,0,0};
int s,t,max,r,c;
int bfs(int x,int y)
{
	if(visit[x][y]>1)//**剪枝。不剪枝应该是TLE的，測试数据有点弱，只是就算没TLE。加了剪枝的时间优化了不少*//
	{
		return visit[x][y];//**假设曾经搜索过这点，就直接返回搜索的这点。不用再进行搜索**//
	}
	for(int k=0;k<4;k++)
	{
		s=x+dx[k];
		t=y+dy[k];
		if(s>=0 && s<r && t>=0 && t<c && a[x][y]<a[s][t])//**搜索的范围。

DP看作单调递增数列**//
		{
			max=bfs(s,t);//**动态规划**//
			if(visit[x][y]<max+1)
			{
				visit[x][y]=max+1;
			}
		}
	}
	return visit[x][y];//**返回这点能取得的最大值**//
}
int main()
{
	int ncases,i,j,ans;
	scanf("%d",&ncases);
	while(ncases--)
	{
		ans=-1;
		scanf("%d %d",&r,&c);
		for(i=0;i<=r-1;i++)
		{
			for(j=0;j<=c-1;j++)
			{
				scanf("%d",&a[i][j]);
				visit[i][j]=1;
			}
		}
		for(i=0;i<=r-1;i++)
		{
			for(j=0;j<=c-1;j++)
			{
				bfs(i,j);
				if(visit[i][j]>ans)
				{
					ans=visit[i][j];
				}
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}        

注意，这里为了保证以随意起点開始搜索，必需要以m*n个起点进行向下深搜，这也自然使得时间复杂度提升了。于是乎我选择的是人人为我型动归，以其它四个方向的值来不断更新当前点的值，然后不断进行递推，但问题来了。为了使得每次递归都以当前最低位置開始，那么每一次递推的起点怎么确定呢？？？这里我的确被坑了一下，想了各种不同的方案，比方，增设一维数组来保存map,j进行排序。但怎么返回坐标呢。之后我有通过双循环不断产生去除上一个最低位置后的最低位置，但这不是比深搜还有麻烦，最后无奈之下，结构体走起。。。问题就直接攻克了

#include<iostream>
#include<algorithm>
using namespace std;
int map[101][101];
int d[101][101];
struct node{
	int x,y;
	int value;
}tmap[10001];
int K;
int dir[2][4]={{1,0,-1,0},{0,1,0,-1}};
int max(int a,int b)
{return a>b?a:b;}
int comp(node a,node b)
{
	return a.value<b.value;
}
int dp(int m,int n)
{
	int i,j,k,x,y,ans=0,p;
	for(p=0;p<K;p++)
	{
		i=tmap[p].x;
		j=tmap[p].y;
		for(k=0;k<4;k++)
		{
			x=i+dir[0][k];
			y=j+dir[1][k];
			if(x>=1&&x<=m&&y>=1&&y<=n&&map[x][y]<map[i][j])
				d[i][j]=max(d[i][j],d[x][y]+1);
			ans=ans>d[i][j]?

ans:d[i][j];
		}
	}
	return ans;
}
int main()
{
	int N,m,n,i,j;
	cin>>N;
	while(N--)
	{
		K=0;
		cin>>m>>n;
		for(i=1;i<=m;i++)
			for(j=1;j<=n;j++)
			{
				cin>>map[i][j];
				d[i][j]=1;
				tmap[K].x=i;
				tmap[K].y=j;
				tmap[K++].value=map[i][j];
			}
			sort(tmap,tmap+K,comp);
			cout<<dp(m,n)<<endl;
	}
	return 0;
}