POJ 1191

这哪里是动规 明明是暴搜。

万恶的黑皮书。。。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=16;
const int M=8;
const int inf=1<<30;
int ans[maxn][maxn][maxn][maxn][maxn];
bool flag[maxn][maxn][maxn][maxn][maxn];
int s[maxn][maxn][maxn][maxn];
int ar[maxn][maxn];

int doubp(int a)
{
	return a*a;
}

int dp(int k,int x1,int y1,int x2,int y2)
{
	if(flag[k][x1][y1][x2][y2])
	{
		return ans[k][x1][y1][x2][y2];
	}
	int i,sum=inf;
	for(i=x1;i<x2;i++)
	{
		sum=min(sum,min(dp(k-1,x1,y1,i,y2)+doubp(s[i+1][y1][x2][y2]), dp(k-1,i+1,y1,x2,y2)+doubp(s[x1][y1][i][y2])));
	}
	for(i=y1;i<y2;i++)
	{
		sum=min(sum,min(dp(k-1,x1,y1,x2,i)+doubp(s[x1][i+1][x2][y2]), dp(k-1,x1,i+1,x2,y2)+doubp(s[x1][y1][x2][i])));
	}
	ans[k][x1][y1][x2][y2]=sum;
	flag[k][x1][y1][x2][y2]=true;
	return sum;
}

int main()
{
	int n;
	while(~scanf("%d",&n))
	{
		memset(s,0,sizeof(s));
		memset(ans,-1,sizeof(ans));
		memset(flag,0,sizeof(flag));
		int i,j,k,l,p,sum=0;
		for(i=1;i<=M;i++)
		{
			for(j=1;j<=M;j++)
			{
				scanf("%d",&ar[i][j]);
				s[i][j][i][j]=ar[i][j];
				sum+=ar[i][j];
			}
		}
		for(i=1;i<=M;i++)
		{
			for(j=1;j<=M;j++)
			{
				for(k=i;k<=M;k++)
				{
					for(l=j;l<=M;l++)
					{
						s[i][j][k][l]=s[i][j][k-1][l]+s[i][j][k][l-1]-s[i][j][k-1][l-1]
						               +s[k][l][k][l];
						ans[1][i][j][k][l]=doubp(s[i][j][k][l]);
						flag[1][i][j][k][l]=true;
					}
				}
			}
		}
		dp(n,1,1,M,M);
		double dn=n;
		printf("%.3lf\n",sqrt(ans[n][1][1][M][M]/dn-(sum/dn)*(sum/dn)));
	}
	return 0;
}