棋盘分割

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 501    Accepted Submission(s): 248

Problem Description
将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行) 原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。  均方差 ,其中平均值 xi为第i块矩形棋盘的总分。  请编程对给出的棋盘及n,求出O'的最小值。 
Input
第1行为一个整数n(1 < n < 15)。  第2行至第9行每行为8个小于100的非负整数,表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
Output
仅一个数,为O'(四舍五入精确到小数点后三位)。
Sample Input
3 1 1 1 1 1 1 1 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 3
Sample Output
1.633

题解:

re了半天,知道真相的我眼泪流了下来,TMD,步数是15,我开的是10;

记忆化搜索,dp[x1][y1][x2][y2]代表到这个矩形状态的最小s[i]平方和;

下面是大神的详细解释:

n刀切割棋盘

下面是8*8的棋盘,每个数字代表棋盘对应点的权值,问切割n刀后,每一块的和  的均方差最小是多少

均方差的公式需要先化简:

由上式得,均方差最小 显然是要 Xi^2 最小

d[k][x1][y1][x2][y2]代表棋盘从(x1,y1)->(x2,y2)已经切了k刀 获得的最小的平方和

用sum[i][j] 代表 从(1,1)点 到 (i,j)点的权值和

这样答案就是 dp[n][1][1][8][8]/n  -(sum[8][8]/n)^2

用S[ (x1,y1) ,( x2,y2) ] 代表 这两点间的权值和

这里用递归dp

状态转移方程:

d[k][x1][y1][x2][y2]=

Min{  横向切: d[k-1] +剩下未切部分的平方和,纵向切:d[k-1]+剩下未切部分的平方和  }

横向切的最优解就是 Min{ d[ k-1,(x1,y1) , ( i ,y2) ] + S[ (i +1,y1) , ( x2 ,y2) ] ,  d[ k-1,(i+1,y1) , ( x2 ,y2) ] + S[ (x1,y1) , ( i ,y2) ] }  ( x1<= i <x2)

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int INF=0x3fffffff;
#define mem(x,y) memset(x,y,sizeof(x))
#define SI(x) scanf("%d",&x)
#define PI(x) printf("%d",x)
#define SD(x) scanf("%lf",&x)
#define P_ printf(" ")
typedef long long LL;
const int MAXN=10;
double dp[20][MAXN][MAXN][MAXN][MAXN];//卧槽
int s[MAXN][MAXN];
double ss[MAXN][MAXN][MAXN][MAXN];
double SUM(int x1,int y1,int x2,int y2){if(ss[x1][y1][x2][y2]>=0)return ss[x1][y1][x2][y2];int temp=s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1];return ss[x1][y1][x2][y2]=temp*temp;
}
double dfs(int k,int x1,int y1,int x2,int y2){if(k==1){dp[k][x1][y1][x2][y2]=SUM(x1,y1,x2,y2);return dp[k][x1][y1][x2][y2];}if(dp[k][x1][y1][x2][y2]>=0)return dp[k][x1][y1][x2][y2];double  temp = 99999999999;for(int i=x1;i<x2;i++){temp=min(temp,SUM(x1,y1,i,y2)+dfs(k-1,i+1,y1,x2,y2));temp=min(temp,SUM(i+1,y1,x2,y2)+dfs(k-1,x1,y1,i,y2));}for(int i=y1;i<y2;i++){temp=min(temp,SUM(x1,i+1,x2,y2)+dfs(k-1,x1,y1,x2,i));temp=min(temp,SUM(x1,y1,x2,i)+dfs(k-1,x1,i+1,x2,y2));}return dp[k][x1][y1][x2][y2]=temp;
}
int main(){int N;while(~SI(N)){int temp;mem(ss,-1);mem(s,0);mem(dp,-1);for(int i=1;i<=8;i++){for(int j=1;j<=8;j++){scanf("%d",&temp);s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+temp;}}double ave=s[8][8]*1.0/N;double ans=dfs(N,1,1,8,8);//printf("%lf %d %lf\n",ans,s[8][8],ave);printf("%.3lf\n",sqrt(ans*1.0/N-ave*ave));}return 0;
}

棋盘分割(记忆化搜索)相关推荐

  1. pku 1191 棋盘分割 DP / 记忆化搜索

    http://poj.org/problem?id=1191 题意:中文省略. 思路:黑说p116有讲解, 主要的状态转移方程为 横着切: dp[k][x1][y1][x2][y2]  = min(d ...

  2. NYOJ 87-棋盘分割(记忆化搜索)

    棋盘分割 时间限制:1000 ms  |  内存限制:65535 KB 难度:6 描述 将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割 ...

  3. [BZOJ 1048] [HAOI2007] 分割矩阵 【记忆化搜索】

    题目链接:BZOJ - 1048 题目分析 感觉这种分割矩阵之类的题目很多都是这样子的. 方差中用到的平均数是可以直接算出来的,然后记忆化搜索 Solve(x, xx, y, yy, k) 表示横坐标 ...

  4. 棋盘里的数学(qduoj 记忆化搜索)

    Description lhcoder有一个n行m列的棋盘,有一颗棋子从左上角(1,1)开始移动,每次只能往右或者往下移动一格,到右下角(n,m)一共有多少移动方案? Input 有多组测试数据,每组 ...

  5. 【BZOJ1048】分割矩阵(记忆化搜索,动态规划)

    [BZOJ1048]分割矩阵(记忆化搜索,动态规划) 题面 BZOJ 洛谷 题解 一个很简单的\(dp\),写成记忆化搜索的形式的挺不错的. #include<iostream> #inc ...

  6. dfs记忆化搜索(带限制的选择问题) 讲解:LeetCode打家劫舍||| / 蓝桥 地宫取宝/蓝桥 k进制数//剪格子//方格分割

    记忆化递归的必要性: 普通的递归可能会重复求解某一值,类似斐波那契数列.同样的子问题可能会被求解多次,这样就会很慢很慢很慢 解决方法:我们把历史求解(子问题)记录下来,如果下次需要求解子问题,那么直接 ...

  7. bzoj 1048: [HAOI2007]分割矩阵(记忆化搜索)

    1048: [HAOI2007]分割矩阵 Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 1133  Solved: 822 [Submit][Sta ...

  8. 【AcWing】数位统计DP、树形DP、状态压缩DP、记忆化搜索

    [AcWing]数位统计DP.树形DP.状态压缩DP.记忆化搜索 一.数位统计DP 二.状态压缩DP 三.树形DP 四.记忆化搜索 一.数位统计DP 计数问题 给定两个整数 a 和 b,求 a 和 b ...

  9. hdu2067 简单dp或者记忆化搜索

    题意: 小兔的棋盘 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Sub ...

最新文章

  1. 使用搜索引擎的10个搜索技巧
  2. 成员资格、授权和安全性(一)
  3. Android 操作SQLite基本用法
  4. cnetos7 vncserver安装与配置
  5. 毕业即失业?你到底在焦虑什么
  6. 通过IHttpHandler实现图片验证码
  7. 优雅的使用Laravel之phpstorm配置
  8. doxygen问题集锦
  9. css 大于号 标签_CSS设计基础选择器篇
  10. python爬取豆瓣电影信息数据
  11. Http请求返回最外层的模型
  12. 用sql语句也样创建oracle,Oracle中用sql语句创建数据库
  13. python视频教程免费慕课网-python视频教程慕课 | 最好的python视频教程谁有
  14. 解决Linux中Ubuntu16.04声音一直静音的方法
  15. 西门子200smart自创库与说明
  16. ajax用户登录注册
  17. 如何解决apt-get中Unmet dependencies问题
  18. PLC抑制干扰电路的设计
  19. 《Real-Time Rendering》第四版学习笔记——Chapter 9 Physically Based Shading(一)
  20. Python3 juypterlab pyecharts 饼图制作

热门文章

  1. python3.8新特性 逻辑表达式_Python3.8新特性
  2. maven显示1.5版本过期的解决办法
  3. LUOGU P2764 最小路径覆盖问题 (最小路径点覆盖)
  4. Scrapy+Scrapy-redis+Scrapyd+Gerapy 分布式爬虫框架整合
  5. MySQL数据库学习笔记(十二)----开源工具DbUtils的使用(数据库的增删改查)
  6. POI 导出excel带小数点的数字格式显示不对解决方法
  7. springboot使用logback日志,部署到tomcat不生效问题解决
  8. OSS上传图片无法在线预览的解决方案
  9. 日常工作问题解决:配置NTP服务器以及一些常见错误解决
  10. eNSP模拟器路由器无法正常启动一直显示“#”——问题解决方法