poj 1191 棋盘分割

刘汝佳黑书上的题目。没什么好说的，就是记忆化搜索，懒人必备的动态规划算法。

#include<iostream>
#include<math.h>
#include<cmath>
#include<stdio.h>
using namespace std;
int s[9][9][9][9];
int dp[15][9][9][9][9];
int map[9][9];
int n;
int mini(int x,int y)
{
if(x<y)
return x;
return y;
}
int search(int left,int x1,int y1,int x2,int y2)
{
if(left==1)
{
dp[left][x1][y1][x2][y2]=s[x1][y1][x2][y2];
return dp[left][x1][y1][x2][y2];
}
if(dp[left][x1][y1][x2][y2]>=0)
return dp[left][x1][y1][x2][y2];
int minx=100000000;
for(int a=y1;a<y2;a++)
{
minx=mini(minx,search(left-1,x1,y1,x2,a)+s[x1][a+1][x2][y2]);
minx=mini(minx,search(left-1,x1,a+1,x2,y2)+s[x1][y1][x2][a]);
}
for(int a=x1;a<x2;a++)
{
minx=mini(minx,search(left-1,x1,y1,a,y2)+s[a+1][y1][x2][y2]);
minx=mini(minx,search(left-1,a+1,y1,x2,y2)+s[x1][y1][a][y2]);
}
dp[left][x1][y1][x2][y2]=minx;
return dp[left][x1][y1][x2][y2];
}
int main()
{
int i,j,k,p,q,sum;
double ans;
cin>>n;
sum=0;
for(i=1;i<=8;i++)
for(j=1;j<=8;j++)
{
cin>>map[i][j];
sum=sum+map[i][j];
}
for(i=1;i<=n;i++)
for(j=1;j<=8;j++)
for(k=1;k<=8;k++)
for(p=1;p<=8;p++)
for(q=1;q<=8;q++)
dp[i][j][k][p][q]=-1;
for(i=1;i<=8;i++)
for(j=1;j<=8;j++)
for(k=1;k<=8;k++)
for(p=1;p<=8;p++)
{
if(k>=i&&p>=j)
{
s[i][j][k][p]=0;
for(int x=i;x<=k;x++)
for(int y=j;y<=p;y++)
s[i][j][k][p]=s[i][j][k][p]+map[x][y];
s[i][j][k][p]=s[i][j][k][p]*s[i][j][k][p];
s[k][p][i][j]=s[i][j][k][p];
s[i][p][k][j]=s[i][j][k][p];
s[k][j][i][p]=s[i][j][k][p];
}
}
ans=sqrt(1.0*search(n,1,1,8,8)/n-(sum*1.0/n)*(sum*1.0/n));
printf("%.3f\n",ans);
return 0;
}