[动态规划] 滑雪

描述

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。

输出

输出最长区域的长度。

样例输入

5 5
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

样例输出

25

解题分析

这是一道很经典的动态规划题目，我们采取了特殊的一个递推方式，而这个递推方式的决定由递推公式决定的，本题有两种方法，一种叫我为人人，一种叫人人为我，这两种方式都可以解决本题。首先，我们可以将dp[i][j]设定为从(i,j)位置一直到它所能达到的最低点的总体长度，那么dp[i][j]=max(dp[i][j],dp[i+dx][j+dy]+1);此处的x+dx,y+dy表示的是(i,j)位置附近的位置，且(i+dx,j+dy)位置雪的高度要小于当前位置雪的高度。之后我们遍历整个dp数组表，寻找一个最大值即可。那么，又有一个问题来了，就是我们怎么去确保我们更新的时候，周围比(i,j)位置雪的高度要小的位置的dp数组是已经更新过的呢？解决方案很简单，我们读入数据后对它们进行一个排序，从小到大，然后我们从最小的那个位置开始用这种方式去递推即可。同理，如果我们设置dp[i][j]表示从其他位置到该位置的最长下降长度，那么我们就要从大到小进行排序，并且递推公式变为了dp[i][j]=max(dp[i][j],dp[i+dx][j+dy]+1);此处的x+dx,y+dy表示的是(i,j)位置附近的位置，且(i+dx,j+dy)位置雪的高度要大于当前位置雪的高度。这就需要确保位置高的雪要先更新了。前面所述两种都是人人为我模式，还有一种我为人人，那么大概就是按第一种思路的话，如果(i,j)位置的雪的高度要小于(i+dx,j+dy)位置的雪的高度，那么，dp[i+dx][j+dy]=max(dp[i+dx][j+dy],dp[i][j]+1);然后，按第二种思路的话就是(i,j)位置的雪的高度要大于(i+dx,j+dy)位置的雪的高度，dp[i+dx][j+dy]=max(dp[i+dx][j+dy],dp[i][j]+1);。

代码实现1

#include <iostream>
#include <cmath>
#include <iomanip>
#include <string>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>
#include <list>
#include <bitset>
using namespace std;

int lot[105][105];

struct node{
	int x,y;
	int num;
};

node a[105*105];

int dp[105][105];

int dx[]={-1,1,0,0};

int dy[]={0,0,-1,1};

int main(){
	int R,C;
	int k=0;
	scanf("%d%d",&R,&C);
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			scanf("%d",&lot[i][j]);
			a[k++]={i,j,lot[i][j]};
			dp[i][j]=1;
		}
	
	sort(a,a+k,[](node a,node b){
		return a.num<b.num;
	});
	
	for(int i=0;i<k;i++){
		for(int c=0;c<4;c++){
			int x1=a[i].x+dx[c];
			int y1=a[i].y+dy[c];
			if(x1>=1 && x1<=R && y1>=1 && y1<=C){
				if(lot[a[i].x][a[i].y]>lot[x1][y1]){
					dp[a[i].x][a[i].y]=max(dp[a[i].x][a[i].y],dp[x1][y1]+1);
				}
			}
		}
	}
	
	int res=1;
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			res=max(res,dp[i][j]);
		}
	cout<<res<<endl;
	return 0;
}

代码实现2

#include <iostream>
#include <cmath>
#include <iomanip>
#include <string>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>
#include <list>
#include <bitset>
using namespace std;

int lot[105][105];

struct node{
	int x,y;
	int num;
};

node a[105*105];

int dp[105][105];

int dx[]={-1,1,0,0};

int dy[]={0,0,-1,1};

int main(){
	int R,C;
	int k=0;
	scanf("%d%d",&R,&C);
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			scanf("%d",&lot[i][j]);
			a[k++]={i,j,lot[i][j]};
			dp[i][j]=1;
		}
	
	sort(a,a+k,[](node a,node b){
		return a.num<b.num;
	});
	
	for(int i=0;i<k;i++){
		for(int c=0;c<4;c++){
			int x1=a[i].x+dx[c];
			int y1=a[i].y+dy[c];
			if(x1>=1 && x1<=R && y1>=1 && y1<=C){
				if(lot[a[i].x][a[i].y]>lot[x1][y1]){
					dp[a[i].x][a[i].y]=max(dp[a[i].x][a[i].y],dp[x1][y1]+1);
				}
			}
		}
	}
	
	int res=1;
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			res=max(res,dp[i][j]);
		}
	cout<<res<<endl;
	return 0;
}

代码实现3

#include <iostream>
#include <cmath>
#include <iomanip>
#include <string>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>
#include <list>
#include <bitset>
using namespace std;

int lot[105][105];

struct node{
	int x,y;
	int num;
};

node a[105*105];

int dp[105][105];

int dx[]={-1,1,0,0};

int dy[]={0,0,-1,1};

int main(){
	int R,C;
	int k=0;
	scanf("%d%d",&R,&C);
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			scanf("%d",&lot[i][j]);
			a[k++]={i,j,lot[i][j]};
			dp[i][j]=1;
		}
	
	sort(a,a+k,[](node a,node b){
		return a.num>b.num;
	});
	
	for(int i=0;i<k;i++){
		for(int c=0;c<4;c++){
			int x1=a[i].x+dx[c];
			int y1=a[i].y+dy[c];
			if(x1>=1 && x1<=R && y1>=1 && y1<=C){
				if(lot[a[i].x][a[i].y]>lot[x1][y1]){
					dp[x1][y1]=max(dp[a[i].x][a[i].y]+1,dp[x1][y1]);
				}
			}
		}
	}
	
	int res=1;
	
	for(int i=1;i<=R;i++)
		for(int j=1;j<=C;j++){
			res=max(res,dp[i][j]);
		}
	cout<<res<<endl;
	return 0;
}