成魔之路−> 算法提高课题解
思路:
-
1.状态表示
集合:区间 [x1,y1,x2,y2] 砍 k 🔪 后的均方差
属性:min -
2. 状态转移
横着切:min(dp(x1,y1,i,y2,k+1)+get(i+1,y1,x2,y2),dp(i+1,y1,x2,y2,k+1)+get(x1,y1,i,y2))
竖着切:min(dp(x1,y1,x2,i,k+1)+get(x1,i+1,x2,y2),dp(x1,i+1,x2,y2,k+1)+get(x1,y1,x2,i))
完整代码
#include<bits/stdc++.h>
using namespace std;
const int N = 9, M = 15;
int n,m=8;
int s[N][N];
double f[N][N][N][N][M];
double X;
//方差
double get(int x1,int y1,int x2,int y2)
{
double sum=s[x2][y2]-s[x2][y1-1]-s[x1-1][y2]+s[x1-1][y1-1]-X;
return sum*sum/n;
}
double dp(int x1,int y1,int x2,int y2,int k)
{
double &v=f[x1][y1][x2][y2][k];
if(v>0) return v; //记忆化搜索,该块已计算过
if(k==n) return get(x1,y1,x2,y2); //最后一块
v=0x3f3f3f3f;
//横着切
for(int i=x1;i<x2;i++)
{
v=min(v,dp(x1,y1,i,y2,k+1)+get(i+1,y1,x2,y2));
v=min(v,dp(i+1,y1,x2,y2,k+1)+get(x1,y1,i,y2));
}
//竖着切
for(int i=y1;i<y2;i++)
{
v=min(v,dp(x1,y1,x2,i,k+1)+get(x1,i+1,x2,y2));
v=min(v,dp(x1,i+1,x2,y2,k+1)+get(x1,y1,x2,i));
}
return v;
}
int main()
{
cin>>n;
for(int i=1;i<=m;i++)
for(int j=1;j<=m;j++)
{
cin>>s[i][j];
s[i][j]+=s[i-1][j]+s[i][j-1]-s[i-1][j-1]; //二维前缀和
}
memset(f,-1,sizeof f);
X=(double)s[m][m]/n; //平均值
printf("%.3lf\n",sqrt(dp(1,1,8,8,1))); //均方差 = sqrt(方差)
return 0;
}
0rz
orz