POJ 1088 滑雪 DP

POJ 1088 滑雪 DP

滑雪

Time Limit: 1000MS
Memory Limit: 65536K
Total Submissions: 92862
Accepted: 35121
Description

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更长。事实上，这是最长的一条。
Input

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

输出最长区域的长度。
Sample Input

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
Sample Output

25
Source

SHTSC 2002

下面是我的 32ms 892K 代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
unsigned int f[105][105],a[105][105];//
struct point
{
    unsigned int x;
    unsigned int y;
    unsigned int high;
};
bool cmp(point a,point b)
{
    return a.high<b.high;
}
point p[13001];
int n,m,maxx,b,temp;
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            p[m*(i-1)+j].high=a[i][j];
            p[m*(i-1)+j].x=i;
            p[m*(i-1)+j].y=j;
            f[i][j]=1;
        }
    }
    sort(p+1,p+m*n+1,cmp);
    for (int i=1;i<=n*m;i++)
    {
        b=0;
        if (p[i].x-1>0&&p[i].high>a[p[i].x-1][p[i].y])//down
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x-1][p[i].y]+1);
        }
        if (p[i].x+1<n+1&&p[i].high>a[p[i].x+1][p[i].y])//up
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x+1][p[i].y]+1);
        }
        if (p[i].y-1>0&&p[i].high>a[p[i].x][p[i].y-1])//left
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x][p[i].y-1]+1);
        }
        if (p[i].y+1<m+1&&p[i].high>a[p[i].x][p[i].y+1])//right
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x][p[i].y+1]+1);
        }
        temp=f[p[i].x][p[i].y];
        maxx=max(maxx,temp);
    }
    printf("%d\n",maxx);
    return 0;
}

下面是内存更小，时间更短的：

#include<cstdio>
#include<algorithm>
using namespace std;
int f[102][102],a[102][102];
struct point
{
    int x;
    int y;
};
bool cmp(point c,point d)
{
    return a[c.x][c.y]<a[d.x][d.y];
}
point p[10005];
int n,m,maxx,b,temp;
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            p[m*(i-1)+j].x=i;
            p[m*(i-1)+j].y=j;
            f[i][j]=1;
        }
    }
    sort(p+1,p+m*n+1,cmp);
    for (int i=1;i<=n*m;i++)
    {
        b=0;
        if (p[i].x-1>0&&a[p[i].x][p[i].y]>a[p[i].x-1][p[i].y])//down
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x-1][p[i].y]+1);
        }
        if (p[i].x+1<n+1&&a[p[i].x][p[i].y]>a[p[i].x+1][p[i].y])//up
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x+1][p[i].y]+1);
        }
        if (p[i].y-1>0&&a[p[i].x][p[i].y]>a[p[i].x][p[i].y-1])//left
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x][p[i].y-1]+1);
        }
        if (p[i].y+1<m+1&&a[p[i].x][p[i].y]>a[p[i].x][p[i].y+1])//right
        {
            f[p[i].x][p[i].y]=max(f[p[i].x][p[i].y],f[p[i].x][p[i].y+1]+1);
        }
        temp=f[p[i].x][p[i].y];
        maxx=max(maxx,temp);
    }
    printf("%d\n",maxx);
    return 0;
}