BSOJ：2909 雷曼兔 DP

2909 -- 【模拟试题】雷曼兔

Description

　　这次，OI山成为了雷曼兔那无尽的冒险传说的新舞台！传说OI山中埋藏着巨大的宝藏，伴随着这个传说的是一个迷题：最瑰丽的舞者将达至精灵世界的彼岸……

　　经过仔细推敲，雷曼兔发现这是一个提示宝藏埋藏位置的谜语，在该谜语中指出了一个特定的路径，只有经过了该路径宝藏才会出现，具体情况如下：

　　OI山的地势图可以看作一个N*N的数字矩阵，由1-N^2的数字组成（每个数字出现且仅出现一次），这些数字表示每个地点的地势高低。雷曼兔的出发点在最高的山顶处，并且每次雷曼兔可以从其当前所在的位置跳跃到任何一个比当前地点高度低的位置，假设雷曼兔该次跳跃从坐标(x1,y1)跳到了坐标(x2,y2)，则这次跳跃的华丽度定义为v=(|x1-x2|+|y1-y2|)^2。而开启宝藏秘密的路径就是从山顶不断跳跃直到山底（高度最低点）的华丽度总和最高的路径，而现在我们想要知道的是这个最高的华丽度总和是多少。

Input

第一行包括一个整数n(n<=50)表示地图的长宽。

接下来n行每行包括n个数表示每个地点的高度。

Output

输出包括一个整数ans，表示从山顶到山底最高华丽度总和。

Sample Input

2

3 2

1 4

Sample Output

9

Hint

【样例解释】

最优路径为4->3->2->1，得分为4+1+4=9。

首先搞清楚，这道题的递推方向是： 按高度递推，因为题目中提到，一个点只有一个高度。

先贴上傻傻的错误DP(有后效性)

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cstdlib>

//#include<cmath>

using namespace std;

inline int read()

{

    int bj=1;

    char ch=getchar();

    while(ch<'0'||ch>'9')

    {

        if(ch=='-')bj=-1;

        ch=getchar();

    }

    int ret=0;

    while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();

    return ret*bj;

}

int n,f[55][55]={0},minn=0x7fffffff,xx,yy,map[55][55]={0};

int squ(int x)

{

    return x*x;

}

int calc(int a,int b,int x,int y)

{

    return squ(abs(a-x)+abs(b-y));

}

void init()

{

    n=read();

    for(int i=1;i<=n;i++)

    {

        for(int j=1;j<=n;j++)

        {

            map[i][j]=read();

            if(minn>map[i][j])

            {

                minn=map[i][j];

                xx=i;yy=j;

            }

        }

    }

}

void DP()

{

    for(int i=1;i<=n;i++)

    {

        for(int j=1;j<=n;j++)

        {

            for(int k=1;k<=n;k++)

            {

                for(int l=1;l<=n;l++)

                {

                    if(i==k&&j==l)continue;

                    if(map[i][j]<map[k][l])f[i][j]=max(f[i][j],f[k][l]+calc(i,j,k,l));

                }

            }

        }

    }

    printf("%d\n",f[xx][yy]);

}

int main()

{

    init();

    DP();

return 0;

}

显然，这个DP的方法不能够让所有能用于更新的点来正确地优化当前点。

改成按高度递推就能A了，方程类似LIS。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
//#include<cmath>
using namespace std;
inline int read()
{
        int bj=1;
        char ch=getchar();
        while(ch<'0'||ch>'9')
        {
                if(ch=='-')bj=-1;
                ch=getchar();
        }
        int ret=0;
        while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
        return ret*bj;
}
struct node
{
        int val,x,y;
}a[1000005];
int f[200005]={0};
int n,maxx=0;
int squ(int x)
{
        return x*x;
}
int calc(int a,int b,int x,int y)
{
        return squ(abs(a-x)+abs(b-y));
}
int main()
{
        n=read();
        for(int i=1;i<=n;i++)
        {
                for(int j=1;j<=n;j++)
                {
                        int p=read();
                        a[p].val=p;
                        a[p].x=i;
                        a[p].y=j;
                }
        }
        for(int i=n*n-1;i;i--)
        {
                for(int j=i+1;j<=n*n;j++)
                {
                        f[i]=max(f[i],f[j]+calc(a[j].x,a[j].y,a[i].x,a[i].y));
                }
        }
        printf("%d\n",f[1]);
return 0;
}