Michael喜欢滑雪百这并不奇怪,因为滑雪的确很刺激。可是为了获得速度,滑的区域必须向下倾斜,而且当你滑到坡底,你不得不再次走上坡或者等待升降机来载你。Michael想知道载一个区域中最长底滑坡。区域由一个二维数组给出。数组的每个数字代表点的高度。下面是一个例子 1 2 3 4 516 17 18 19 615 24 25 20 714 23 22 21 813 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。后面是下一组数据;

输出

输出最长区域的长度。

样例输入


1

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,写完后查阅了下资料,发现使用递归算法居多。特此记录下,使用动态规划的算法解决滑雪问题。

      这里的动态规划,可以理解为牺牲存储空间,换取时间资源。在本题中,先初始化所有点的各种信息值,如沿着某点最多可以滑行的长度count。初始化后,就可以进行动态规划。从最低的点A开始,记录它周围比它低的点的个数count。然后找到数值比A大的最近的点B,找到B周围所有比它低的点,然后比较所有点的count值,将最大的count加1作为B的count。

      上面的描述只是一个大概的思想,实施起来,还需要解决一些问题。比如,如何找到点A大的最近的点,找到该点后,又需要对它周围的点进行类似的处理。这种思想,给人使用递归的冲动。实际上,我们只需要把所有的位置点从小到大排序起来,依次记录该点的坐标,高度,该点起最大的滑行长度等信息。

      这样就很自然地构造出了一个数据结构

struct MyPoint
{
    int i,j;
    int height;
    int val;
};



上面提到了按照高度进行排序,再进行其他的处理。在C++的STL中,有一个sort函数可提供排序功能,并且支持自定义的函数比较。我们只需要定义一个比较高度的函数即可。

bool less_height(const MyPoint & m1, const MyPoint & m2)
{
    return m1.height< m2.height;
}




方向的遍历

      当对一个点的四周进行遍历时,虽然可以直接一个一个地引用下标,来读取周围点的信息,但我们有一种更好的实现方式。在这里我们定义一个数组存储自定义的结构体,表示周围点的方向矢量。

      为了实现这一点,我们首先定义点的结构体

struct Pt
{
    int x, y;
}




 

定义完结构体后,我们就按上右下左的顺序,把周围点的方向矢量放到Direction数组中


Pt direction[4]={{-1,0},{0,1},{1,0},{0,-1}};




 

完整代码:


#include <iostream>
#include <iterator>
#include <vector>
#include <algorithm>
using namespace std;

struct MyPoint
{
    int i,j;
    int height;
    int val;
};
struct Pt
{
    int x;
    int y;
};
bool less_height(const MyPoint & m1, const MyPoint & m2)
{
    return m1.height< m2.height;
}
int main()
{
    MyPoint mp;
    vector<MyPoint> v;
    int max;
    int i ,j ,m ,n ,N;
    int x,y;
    //定以矩阵,0表示高度,1表示路径长度
    int data[100][100][2];
    Pt direction[4]={{-1,0},{0,1},{1,0},{0,-1}};    
    cin>>N;
    while(N--)
    {
        v.clear();
        cin>>m>>n;
        //初始化矩阵
        for(i=0; i<m; ++i)
 for(j=0; j<n; ++j)
            {
                data[i][j][0]=0;
                data[i][j][1]=0;
            }
        for(i=0; i<m; ++i)
            for(j=0; j<n; ++j)
            {
                mp.i=i,mp.j=j;
                mp.val=0;
                cin>>mp.height;
                data[i][j][0]=mp.height;
                v.push_back(mp);
            }
        int k=0;
        sort(v.begin(), v.end(), less_height);
        for(i=0; i<v.size(); ++i)
        {
            int t;
            max=0;
            for(t=0; t<4; ++t)
            {
                x=v[i].i + direction[t].x;
                y=v[i].j + direction[t].y;
                //越界检查
                if(x<0 || x>=m || y<0 || y>=n)
                    continue;
                if(data[x][y][0]<v[i].height && data[x][y][1]>max)
                {
                    max = data[x][y][1];
                }
            }
            x=v[i].i; y=v[i].j;
            v[i].val = max+1;
            data[x][y][1]=v[i].val;
 }
        max=0;
        for(i=0; i<m; ++i)
        {
            for(j=0; j<n; ++j)
                if(data[i][j][1]>max)
                    max=data[i][j][1];
        }
        cout<<max<<endl;

    }

    return 0;
}