(我是不会告诉你我是抄的http://www.cnblogs.com/scau20110726/archive/2013/02/27/2936050.html这个人的)

一开始没有想到要化一下方差的式子 怎么搞都挂……

尴尬

/*
题目固定是8*8,本来想用点的坐标来表示矩形的,但是发现用标号来表示会方便一点
对于最小的小方格,用(i,j)表示,即第i行第j列的小方格,注意不是点的坐标
所以对于一个矩形,我们用它左上角的小方格和右下角的小方格来表示
例如,整个棋盘就是(1,1),(8,8)
另外题目要求,每次分割出一个矩形后,剩下的也必须是矩形
那么其实每次分割只能切一刀,如果是切两刀得到的矩形,那么剩下的就不会是矩形了
只能横着切或者竖着切,而且一切的话要从头切到底(这很容易理解)
另外还有一个东西之前理解错了,就是一刀切下去会得到两个矩形,
选一个为本次切割得到的,以后只能切另一个,选出来的那个以后不能再切了
(如果是两者都能切,感觉复杂很多)
然后动态转移方程觉得还是比较容易想到的,大矩形dp值由小矩形dp值推得来
还要加上次数,dp[n][x1][y1][x2][y2]就是要令当前矩形分出n个小矩形,也就是切割n-1刀
我们要的目标值就是dp[n][1][1][8][8]
一:横着切,当前矩形将会分成上下两份
1.选上面:dp[k][x1][y1][x2][y2]=s[x1][y1][x][y2]+dp[k-1][x+1][y1][x2][y2];
2.选下面:dp[k][x1][y1][x2][y2]=s[x+1][y1][x2][y2]+dp[k-1][x1][y1][x][y2];
x1<=x<x2
二:竖着切,当前矩形将会分成左右两份
1.选左边:dp[k][x1][y1][x2][y2]=s[x1][y1][x2][y]+dp[k-1][x1][y+1][x2][y2];
2.选右边:dp[k][x1][y1][x2][y2]=s[x1][y+1][x2][y2]+dp[k-1][x1][y1][x2][y];
y1<=y<y2
当k=1时也就是不用再继续分割了,dp[1][x1][y1][x2][y2]=s[x1][y1][x2][y2];
显然这个DP用记忆化搜索来做更合适,递推的话感觉很难写
最后强调一个细节问题,dp数组全部初始化为-1,表示还没被计算,
然后 dp[1][x1][y1][x2][y2]=s[x1][y1][x2][y2] 这个也好理解,就是不用切的时候的dp值
最后,比如当前要求的是dp[k][x1][y1][x2][y2],一开始要赋初值为INF,再开始枚举切割方案
因为对于dp[k][x1][y1][x2][y2],要分出k块矩形,但是不一定能分得到,可能根本不够分
所以当前状态如果根本分不出k个小矩形的话,这个状态是一个不可能的状态,为INF
整个枚举过程中它的值也不会更新
*/

这个题解中有两个多余的判断可以去掉。。

//By SiriusRen
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n,a[9][9],f[9][9][9][9][65],s[9][9][9][9],sum;
int dfs(int x1,int y1,int x2,int y2,int k){
if(~f[x1][y1][x2][y2][k])return f[x1][y1][x2][y2][k];
int temp=0x3fffffff;
for(int i=x1+1;i<=x2;i++){
temp=min(temp,dfs(i,y1,x2,y2,k-1)+s[x1][y1][i-1][y2]);
temp=min(temp,dfs(x1,y1,i-1,y2,k-1)+s[i][y1][x2][y2]);
}
for(int i=y1+1;i<=y2;i++){
temp=min(temp,dfs(x1,i,x2,y2,k-1)+s[x1][y1][x2][i-1]);
temp=min(temp,dfs(x1,y1,x2,i-1,k-1)+s[x1][i][x2][y2]);
}
return f[x1][y1][x2][y2][k]=temp;
}
int main(){
memset(f,-1,sizeof(f));
scanf("%d",&n);
for(int i=1;i<=8;i++)
for(int j=1;j<=8;j++)
scanf("%d",&a[i][j]),sum+=a[i][j];
for(int x1=1;x1<=8;x1++)
for(int y1=1;y1<=8;y1++)
for(int x2=x1;x2<=8;x2++)
for(int y2=y1;y2<=8;y2++){
int temp=0;
for(int i=x1;i<=x2;i++)
for(int j=y1;j<=y2;j++)
temp+=a[i][j];
f[x1][y1][x2][y2][1]=s[x1][y1][x2][y2]=temp*temp;
}
double tmp=1.0*sum/n;
printf("%.3lf\n",sqrt(1.0*dfs(1,1,8,8,n)/n-tmp*tmp));
}

POJ 1191 记忆化搜索的更多相关文章

  1. poj 1088(记忆化搜索)

    滑雪 Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 88560   Accepted: 33212 Description ...

  2. Test for Job (poj 3249 记忆化搜索)

    Language: Default Test for Job Time Limit: 5000MS   Memory Limit: 65536K Total Submissions: 9733   A ...

  3. 滑雪(POJ 1088 记忆化搜索)

    滑雪 Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 88094   Accepted: 33034 Description ...

  4. POJ 2704 Pascal's Travels 【DFS记忆化搜索】

    题目传送门:http://poj.org/problem?id=2704 Pascal's Travels Time Limit: 1000MS   Memory Limit: 65536K Tota ...

  5. POJ 1579 Function Run Fun 【记忆化搜索入门】

    题目传送门:http://poj.org/problem?id=1579 Function Run Fun Time Limit: 1000MS   Memory Limit: 10000K Tota ...

  6. 专题1:记忆化搜索/DAG问题/基础动态规划

      A OpenJ_Bailian 1088 滑雪     B OpenJ_Bailian 1579 Function Run Fun     C HDU 1078 FatMouse and Chee ...

  7. POJ 1191 棋盘分割 【DFS记忆化搜索经典】

    题目传送门:http://poj.org/problem?id=1191 棋盘分割 Time Limit: 1000MS   Memory Limit: 10000K Total Submission ...

  8. poj 1191 棋盘切割 (压缩dp+记忆化搜索)

    一,题意: 中文题 二.分析: 主要利用压缩dp与记忆化搜索思想 三,代码: #include <iostream> #include <stdio.h> #include & ...

  9. POJ 1088 滑雪(记忆化搜索)

    滑雪 Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 92384   Accepted: 34948 Description ...

随机推荐

  1. 【Paper Reading】Object Recognition from Scale-Invariant Features

    Paper: Object Recognition from Scale-Invariant Features Sorce: http://www.cs.ubc.ca/~lowe/papers/icc ...

  2. 带入gRPC:gRPC Deadlines

    带入gRPC:gRPC Deadlines 原文地址:带入gRPC:gRPC Deadlines项目地址:https://github.com/EDDYCJY/go... 前言 在前面的章节中,已经介 ...

  3. vue-router query和params传参(接收参数)$router $route的区别

    今天做项目时踩到了vue-router传参的坑(query和params),所以决定总结一下二者的区别. 直接总结干货!!! 1.query方式传参和接收参数 传参: this.$router.pus ...

  4. C语言中数据类型的字节数

    类型 16位 32 位 64位 char 1 1 1 short int 2 2 2 int 2 4 4 unsigned int 2 4 4 float 4 4 4 double 8 8 8 lon ...

  5. HDU 4081 Qin Shi Huang&#39;s National Road System 最小生成树

    点击打开链接题目链接 Qin Shi Huang's National Road System Time Limit: 2000/1000 MS (Java/Others)    Memory Lim ...

  6. CCEditBox/CCEditBoxImplAndroid

    #ifndef __CCEDITBOXIMPLANDROID_H__ #define __CCEDITBOXIMPLANDROID_H__ #include "cocos2d.h" ...

  7. 软件project之软件设计

    英雄是随着历史的时代产生的.软工也不例外.软件project这一门学科主要是为了解决当代软件危机而诞生的, 学习软件project的视频过后,最终让我揭开了它的神奇面纱,让我对软工设计有了一个初步的认 ...

  8. Pleasant sheep and big big wolf

    pid=3046">点击打开链接 题目:在一个N * M 的矩阵草原上,分布着羊和狼.每一个格子仅仅能存在0或1仅仅动物.如今要用栅栏将全部的狼和羊分开.问怎么放,栅栏数放的最少,求出 ...

  9. PermissionError: [Errno 13] in python

    出现该错误,首先尝试以管理员身份运行 cmd.exe 程序,然后关闭所有的与 python 相关的进程. 1. open 打开一个文件夹(目录),而不是文件 这一错误一般发生在使用 open函数对文件 ...

  10. hive安装用mysql作为元数据库,mysql的设置

    mysql的设置 在要作为元数据库的mysql服务器上建立hive数据库: #建立数据库 create database if not exists hive; #设置远程登录的权限 GRANT AL ...