poj 1191dp 클래식

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
using namespace std;
int dp[20][10][10][10][10];
int s[10][10][10][10];
int mat[10][10];
int n;

int cal(int k,int x1,int y1,int x2,int y2)
{
	if(k==1 || dp[k][x1][y1][x2][y2]!=-1)
		return dp[k][x1][y1][x2][y2];

	int &t=dp[k][x1][y1][x2][y2];
	for(int a=x1;a<x2;++a)
	{
		int tt1=cal(k-1,a+1,y1,x2,y2);
		if(t==-1 || t>tt1+s[x1][y1][a][y2])
			t=tt1+s[x1][y1][a][y2];
		int tt2=cal(k-1,x1,y1,a,y2);
		if(t==-1 || t>s[a+1][y1][x2][y2]+tt2)
			t=s[a+1][y1][x2][y2]+tt2;
	}
	for(int b=y1;b<y2;++b)
	{
		int tt3=cal(k-1,x1,b+1,x2,y2);
		if(t==-1 || t>tt3+s[x1][y1][x2][b])
			t=tt3+s[x1][y1][x2][b];
		int tt4=cal(k-1,x1,y1,x2,b);
		if(t==-1 || t>s[x1][b+1][x2][y2]+tt4)
			t=s[x1][b+1][x2][y2]+tt4;
	}
	return dp[k][x1][y1][x2][y2];
}
int val[10][10];
int main ()
{
	memset(mat,0,sizeof(mat));
	memset(s,0,sizeof(s));
	while(scanf("%d",&n)!=EOF)
	{
		memset(dp,-1,sizeof(dp));
		int sum=0;
		for(int i=1;i<=8;++i)
			for(int j=1;j<=8;++j)
				scanf("%d",&mat[i][j]),sum+=mat[i][j],val[i][j]=mat[i][j];
		for(int i=1;i<=8;++i)
		{
			for(int j=1;j<=8;++j)
			{
				int tt=mat[i][j-1];
				for(int k=1;k<=i;++k)
					tt+=val[k][j];
				mat[i][j]=tt;
			}
		}
		for(int x2=1;x2<=8;++x2)
			for(int y2=1;y2<=8;++y2)
				for(int x1=1;x1<=x2;++x1)
					for(int y1=1;y1<=y2;++y1)
					{
						s[x1][y1][x2][y2]=mat[x2][y2]-mat[x2][y1-1]-mat[x1-1][y2]+mat[x1-1][y1-1];
						s[x1][y1][x2][y2]*=s[x1][y1][x2][y2];
						dp[1][x1][y1][x2][y2]=s[x1][y1][x2][y2];
					}
		
		printf("%.3lf
",sqrt(cal(n,1,1,8,8)*1.0/n-pow(sum*1.0/n,2)));
	}
	return 0;
}