推式子
ans=(n+m−1)∑i=1n+m−1Ai 2−(∑i=1n+m−1Ai)2
a
n
s
=
(
n
+
m
−
1
)
∑
i
=
1
n
+
m
−
1
A
i
2
−
(
∑
i
=
1
n
+
m
−
1
A
i
)
2
令 fi,j,k f i , j , k 表示到格子 (i,j) ( i , j ) 且路径权值和为 k k <script type="math/tex" id="MathJax-Element-67">k</script> 时,权值平方和最小是多少
这样DP就行了
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N=35;
int f[N][N][2*N*N];
int t,n,m,a[N][N];
int main(){
freopen("1.in","r",stdin);
freopen("1.out","w",stdout);
scanf("%d",&t); int Case=0;
while(t--){
memset(f,0x7f>>1,sizeof(f)); int inf=***f;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
f[1][1][a[1][1]]=a[1][1]*a[1][1];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
if(i==1 && j==1) continue;
int cur=a[i][j];
for(int k=cur;k<=30*(i+j);k++)
if(f[i][j-1][k-cur]!=inf || f[i-1][j][k-cur]!=inf)
f[i][j][k]=min(f[i][j-1][k-cur],f[i-1][j][k-cur])+cur*cur;
}
int ans=1<<30;
for(int i=1;i<=(n+m-1)*30;i++)
if(f[n][m][i]!=inf)
ans=min(ans,(n+m-1)*f[n][m][i]-i*i);
printf("Case #%d: %d\n",++Case,ans);
}
return 0;
}