HDU 5492 Find a path 2015合肥

N* M(1<=N，M<=30)的矩阵，矩阵的每一格有一个非负权值(<=30)，从(1，1)出发，每次只能向右或向下移动，到达(n，m)时，经过的格子的权值形成一序列，求此序列方差最小值*(n-m+1)。

图片有点大

这样，我们要求的就是平方和，与和的平方，利用dp来编写。dp[i][j][k]表示在i，j这个点，元素的和为k时，元素的平方和的最小值。
dp[i][j][k] = min(dp[i][j][k], dp[i-1][j][k-mp[i][j]] + mp[i][j]*mp[i][j]); //左边的点
dp[i][j][k] = min(dp[i][j][k], dp[i][j-1][k-mp[i][j]] + mp[i][j]*mp[i][j]); //上边的点

#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#define INF 0x3f3f3f3f
#define N 1800	//总和不会超过1800
typedef long long LL;
using namespace std;
const LL maxn=40;

LL dp[maxn][maxn][1850];
LL mp[maxn][maxn];

int main()
{
    LL T;
    scanf("%lld",&T);
    for(LL I=1; I<=T; I++){
        LL n,m;
        scanf("%lld%lld",&n,&m);

        for(LL i=1; i<=n; i++)
            for(LL j=1; j<=m; j++)
                scanf("%lld",&mp[i][j]);

        for(LL i=1; i<=n; i++)
            for(LL j=1; j<=m; j++)
                for(LL k=0; k<=N; k++)
                    dp[i][j][k] = INF;// 初始化

        dp[1][1][ mp[1][1] ] = mp[1][1] * mp[1][1];

        for(LL i=1; i<=n; i++){
            for(LL j=1; j<=m; j++){
                for(LL k=mp[i][j]; k<=N; k++){
                    if(i-1>=1) dp[i][j][k] = min(dp[i][j][k], dp[i-1][j][k-mp[i][j]] + mp[i][j]*mp[i][j]);
                    if(j-1>=1) dp[i][j][k] = min(dp[i][j][k], dp[i][j-1][k-mp[i][j]] + mp[i][j]*mp[i][j]);

                    //printf("%d %d %d - %d\n",i,j,k,dp[i][j][k]);
                }
            }
        }
        LL ans=INF;
        for(LL k=0; k<=N; k++){
            if(dp[n][m][k]!=INF){
                ans =min( (n+m-1)*dp[n][m][k] - k*k , ans); // (n+m-1)*平方和 - 和的平方
            }
        }
        printf("Case #%lld: %lld\n",I,ans);
    }
}