POJ 1191 棋盘分割 G++ 均方差 dp 没掌握

#include <iostream>
#include <cstdio>
#include <cstring> 
#include <iomanip>
#include <cmath> 
using namespace std;
//英语       看博友分析       抄博友程序       均方差        dp       没掌握
int da[10][10];
int sum[10][10];
int dp[20][10][10][10][10]; 
int cal(int k,int x1,int y1,int x2,int y2)
{
	int ans=0x3f3f3f3f;
	for(int mid=x1+1;mid<x2;mid++)//横切 
	{
		ans=min(ans,dp[k-1][mid][y1][x2][y2]+dp[1][x1][y1][mid][y2]);
		ans=min(ans,dp[k-1][x1][y1][mid][y2]+dp[1][mid][y1][x2][y2]);
	}
	for(int mid=y1+1;mid<y2;mid++)//背 
	{
		ans=min(ans,dp[k-1][x1][mid][x2][y2]+dp[1][x1][y1][x2][mid]);
		ans=min(ans,dp[k-1][x1][y1][x2][mid]+dp[1][x1][mid][x2][y2]);
	}
	return ans;
}
int main()
{
	int n;
	while(cin>>n)
	{
		memset(da,0,sizeof(da));
		memset(sum,0,sizeof(sum));
		memset(dp,0,sizeof(dp));
		for(int i=1;i<=8;i++)//背 
		{
			for(int j=1;j<=8;j++)
			{
				cin>>da[i][j];
			}
		}
		for(int i=1;i<=8;i++)
		{
			for(int j=1;j<=8;j++)
			{
				sum[i][j]=sum[i][j-1]+sum[i-1][j]+da[i][j]-sum[i-1][j-1];
			}
		}
		for(int x1=0;x1<=8;x1++)//背 
		{
			for(int y1=0;y1<=8;y1++)
			{
				for(int x2=x1+1;x2<=8;x2++)
				{
					for(int y2=y1+1;y2<=8;y2++)
					{
						int t=sum[x2][y2]-sum[x2][y1]-sum[x1][y2]+sum[x1][y1];
						dp[1][x1][y1][x2][y2]=t*t;
					}
				}
			}
		}
		for(int k=2;k<15;k++)
		{
			for(int x1=0;x1<=8;x1++)//没掌握 
			{
				for(int y1=0;y1<=8;y1++)
				{
					for(int x2=x1+1;x2<=8;x2++)
					{
						for(int y2=y1+1;y2<=8;y2++)
						{
							dp[k][x1][y1][x2][y2]=cal(k,x1,y1,x2,y2);
						}
					}
				}
			}
		}
		cout<<fixed<<setprecision(3)<<sqrt(1.0*dp[n][0][0][8][8]/n-1.0*sum[8][8]*sum[8][8]/(n*n))<<endl;//背 		
	} 
	return 0;
}