boj 16

dfs,另外用dp[k][x1][y1][x2][y2]进行标记

#include<iostream>
#include<cmath>
using namespace std;

#define N 20
double G[10][10],dp[20][10][10][10][10];//棋盘；
double e;

double f(int x1,int y1,int x2,int y2)
{
	return (G[x2][y2]-G[x1-1][y2]-G[x2][y1-1]+G[x1-1][y1-1])*(G[x2][y2]-G[x1-1][y2]-G[x2][y1-1]+G[x1-1][y1-1]);
}
double Div(int k,int x1,int y1,int x2,int y2)
{
	if(k==0)
		return f(x1,y1,x2,y2);
	if(fabs(dp[k][x1][y1][x2][y2])>1e-6)return dp[k][x1][y1][x2][y2];
	double mins=300000000;
	for(int i=x1;i<x2;i++)
		mins=min(mins,min(Div(k-1,x1,y1,i,y2)+f(i+1,y1,x2,y2),Div(k-1,i+1,y1,x2,y2)+f(x1,y1,i,y2)));
	for(int i=y1;i<y2;i++)
		mins=min(mins,min(Div(k-1,x1,y1,x2,i)+f(x1,i+1,x2,y2),Div(k-1,x1,i+1,x2,y2)+f(x1,y1,x2,i)));
	dp[k][x1][y1][x2][y2]=mins;
	return mins;
}
int main()
{
	int n;
	int sum=0;
	scanf("%d",&n);
	for(int i=1;i<=8;i++)
		G[0][i]=0;
	for(int i=1;i<=8;i++)
	{
		int tmp=0;
		for(int j=1;j<=8;j++)
		{
			int a;
			scanf("%d",&a);
			sum+=a;
			tmp+=a;
			G[i][j]=G[i-1][j]+tmp;
		}
	}
	e=sum/n;
	printf("%.3lf\n",sqrt(Div(n-1,1,1,8,8)/n-pow(e,2)));
}