【POJ 1191】 棋盘分割(DP)
2016-03-22 21:23
621 查看
【POJ 1191】 棋盘分割(DP)
Description
将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)
原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。
均方差
棋盘分割
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
Sample Output
Source
Noi 99
[Submit] [Go Back] [Status] [Discuss]
Home Page
Go Back
To
top
,其中平均值
,xi为第i块矩形棋盘的总分。
请编程对给出的棋盘及n,求出O'的最小值。
Input
第1行为一个整数n(1 < n < 15)。
第2行至第9行每行为8个小于100的非负整数,表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
Output
仅一个数,为O'(四舍五入精确到小数点后三位)。
Sample Input
Sample Output
Source
Noi 99
刷了一圈,刷回到最原始的东西。其实也是学不完的东西。。。DP……花样太多了,烧脑特辑。。。
题目很简单粗暴
8*8的一个正方形,每个点有值(0~100)要求切割n-1次,每次只能从余下部分的边出发切割到另一边,题目中的图画的挺清楚的。然后问最小的方差。。。
方差公式是给了。但如果死盯着这公式。。想再久应该也是想不粗来的……
需要变一下形。就借用别的博客的推导了~~
S^2 = 1/n∑(Xi - X)^2
= 1/n(n*X^2 + ∑Xi^2 - 2*X∑Xi)
= 1/n∑Xi^2 - X^2;
(http://blog.csdn.net/xymscau/article/details/6854410)
这样其实应该就有点小想法了~你会发现X(均值)是固定的(sum/n),只要保证∑Xi^2尽可能小,最终的结果就一定是最小值。并且可以发现S^2不会出现负值,也就是说∑Xi^2再小也比X^2要大,因此就可以放心地dp了。或者说,放心地(记忆化)搜索了,每次切割时,继续搜索的部分继续搜,另一部分需要快速给出一个合值,因此可以预处理一下所有的区间的平方,或者每次用到暴力跑也可以。
需要注意的是切出n个矩形要切n-1刀。
还有,dp数组要记得记录切割刀数。。一开始就是因为没记录,导致混乱,最终WA
另外,dp最先初始化的时候不要初始化为INF,因为可能某些部分最终结果就是INF,也就是无解,如果一开始初始化为INF,并且以这为记忆化返回值的判断,会导致超时。因为到此状态时仍会向下dp。正解应为代码中初始化未-1,当发现是第一次搜索该状态时,立马转变为INF,然后搜出该状态最小值(若不合法,恰巧就是上一步赋值的INF(无穷大))
代码如下:
#include <iostream>
#include <cmath>
#include <vector>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>
#include <list>
#include <algorithm>
#include <map>
#include <set>
#define LL long long
#define Pr pair<int,int>
#define fread() freopen("in.in","r",stdin)
#define fwrite() freopen("out.out","w",stdout)
using namespace std;
const int INF = 0x3f3f3f3f;
const int msz = 10000;
const int mod = 1e9+7;
const double eps = 1e-8;
double arg;
int n;
int mp[8][8];
int dp[16][8][8][8][8];
int sum[8][8][8][8];
void init()
{
for(int x1 = 0; x1 < 8; ++x1)
for(int y1 = 0; y1 < 8; ++y1)
for(int x2 = x1; x2 < 8; ++x2)
for(int y2 = y1; y2 < 8; ++y2)
{
int tmp = 0;
for(int i = x1; i <= x2; ++i)
for(int j = y1; j <= y2; ++j)
tmp += mp[i][j];
sum[x1][y1][x2][y2] = tmp*tmp;
//printf("sum[%d][%d][%d][%d]:%d\n",x1,y1,x2,y2,sum[x1][y1][x2][y2]);
}
}
int dfs(int s,int x1,int y1,int x2,int y2)
{
if(s == n) return sum[x1][y1][x2][y2];
if(x1 == x2 && y1 == y2) return INF;
if(~dp[s][x1][y1][x2][y2]) return dp[s][x1][y1][x2][y2];
dp[s][x1][y1][x2][y2] = INF;
for(int i = x1; i < x2; ++i)
{
int tmp = min(dfs(s+1,x1,y1,i,y2)+sum[i+1][y1][x2][y2],sum[x1][y1][i][y2]+dfs(s+1,i+1,y1,x2,y2));
//printf("(%d,%d),(%d,%d)->(%d,%d),(%d,%d) %d\n",x1,y1,i,y2,i+1,y1,x2,y2,tmp);
dp[s][x1][y1][x2][y2] = min(dp[s][x1][y1][x2][y2],tmp);
}
for(int i = y1; i < y2; ++i)
{
int tmp = min(dfs(s+1,x1,y1,x2,i)+sum[x1][i+1][x2][y2],sum[x1][y1][x2][i]+dfs(s+1,x1,i+1,x2,y2));
//printf("(%d,%d),(%d,%d)->(%d,%d),(%d,%d) %d\n",x1,y1,x2,i,x1,i+1,x2,y2,tmp);
dp[s][x1][y1][x2][y2] = min(dp[s][x1][y1][x2][y2],tmp);
}
return dp[s][x1][y1][x2][y2];
}
int main()
{
//fread();
//fwrite();
while(~scanf("%d",&n))
{
arg = 0;
for(int i = 0; i < 8; ++i)
for(int j = 0; j < 8; ++j)
{
scanf("%d",&mp[i][j]);
arg += mp[i][j];
}
memset(dp,-1,sizeof(dp));
init();
dfs(1,0,0,7,7);
//printf("%d\n",dp[0][0][7][7]);
//printf("%f\n",arg);
printf("%.3f\n",sqrt(dp[1][0][0][7][7]*1.0/n-arg/n*(arg*1.0/n)));
}
return 0;
}
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 13811 | Accepted: 4917 |
将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)
原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。
均方差
棋盘分割
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 13811 | Accepted: 4917 |
将一个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
Source
Noi 99
[Submit] [Go Back] [Status] [Discuss]
Home Page
Go Back
To
top
,其中平均值
,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
Source
Noi 99
刷了一圈,刷回到最原始的东西。其实也是学不完的东西。。。DP……花样太多了,烧脑特辑。。。
题目很简单粗暴
8*8的一个正方形,每个点有值(0~100)要求切割n-1次,每次只能从余下部分的边出发切割到另一边,题目中的图画的挺清楚的。然后问最小的方差。。。
方差公式是给了。但如果死盯着这公式。。想再久应该也是想不粗来的……
需要变一下形。就借用别的博客的推导了~~
S^2 = 1/n∑(Xi - X)^2
= 1/n(n*X^2 + ∑Xi^2 - 2*X∑Xi)
= 1/n∑Xi^2 - X^2;
(http://blog.csdn.net/xymscau/article/details/6854410)
这样其实应该就有点小想法了~你会发现X(均值)是固定的(sum/n),只要保证∑Xi^2尽可能小,最终的结果就一定是最小值。并且可以发现S^2不会出现负值,也就是说∑Xi^2再小也比X^2要大,因此就可以放心地dp了。或者说,放心地(记忆化)搜索了,每次切割时,继续搜索的部分继续搜,另一部分需要快速给出一个合值,因此可以预处理一下所有的区间的平方,或者每次用到暴力跑也可以。
需要注意的是切出n个矩形要切n-1刀。
还有,dp数组要记得记录切割刀数。。一开始就是因为没记录,导致混乱,最终WA
另外,dp最先初始化的时候不要初始化为INF,因为可能某些部分最终结果就是INF,也就是无解,如果一开始初始化为INF,并且以这为记忆化返回值的判断,会导致超时。因为到此状态时仍会向下dp。正解应为代码中初始化未-1,当发现是第一次搜索该状态时,立马转变为INF,然后搜出该状态最小值(若不合法,恰巧就是上一步赋值的INF(无穷大))
代码如下:
#include <iostream>
#include <cmath>
#include <vector>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <queue>
#include <stack>
#include <list>
#include <algorithm>
#include <map>
#include <set>
#define LL long long
#define Pr pair<int,int>
#define fread() freopen("in.in","r",stdin)
#define fwrite() freopen("out.out","w",stdout)
using namespace std;
const int INF = 0x3f3f3f3f;
const int msz = 10000;
const int mod = 1e9+7;
const double eps = 1e-8;
double arg;
int n;
int mp[8][8];
int dp[16][8][8][8][8];
int sum[8][8][8][8];
void init()
{
for(int x1 = 0; x1 < 8; ++x1)
for(int y1 = 0; y1 < 8; ++y1)
for(int x2 = x1; x2 < 8; ++x2)
for(int y2 = y1; y2 < 8; ++y2)
{
int tmp = 0;
for(int i = x1; i <= x2; ++i)
for(int j = y1; j <= y2; ++j)
tmp += mp[i][j];
sum[x1][y1][x2][y2] = tmp*tmp;
//printf("sum[%d][%d][%d][%d]:%d\n",x1,y1,x2,y2,sum[x1][y1][x2][y2]);
}
}
int dfs(int s,int x1,int y1,int x2,int y2)
{
if(s == n) return sum[x1][y1][x2][y2];
if(x1 == x2 && y1 == y2) return INF;
if(~dp[s][x1][y1][x2][y2]) return dp[s][x1][y1][x2][y2];
dp[s][x1][y1][x2][y2] = INF;
for(int i = x1; i < x2; ++i)
{
int tmp = min(dfs(s+1,x1,y1,i,y2)+sum[i+1][y1][x2][y2],sum[x1][y1][i][y2]+dfs(s+1,i+1,y1,x2,y2));
//printf("(%d,%d),(%d,%d)->(%d,%d),(%d,%d) %d\n",x1,y1,i,y2,i+1,y1,x2,y2,tmp);
dp[s][x1][y1][x2][y2] = min(dp[s][x1][y1][x2][y2],tmp);
}
for(int i = y1; i < y2; ++i)
{
int tmp = min(dfs(s+1,x1,y1,x2,i)+sum[x1][i+1][x2][y2],sum[x1][y1][x2][i]+dfs(s+1,x1,i+1,x2,y2));
//printf("(%d,%d),(%d,%d)->(%d,%d),(%d,%d) %d\n",x1,y1,x2,i,x1,i+1,x2,y2,tmp);
dp[s][x1][y1][x2][y2] = min(dp[s][x1][y1][x2][y2],tmp);
}
return dp[s][x1][y1][x2][y2];
}
int main()
{
//fread();
//fwrite();
while(~scanf("%d",&n))
{
arg = 0;
for(int i = 0; i < 8; ++i)
for(int j = 0; j < 8; ++j)
{
scanf("%d",&mp[i][j]);
arg += mp[i][j];
}
memset(dp,-1,sizeof(dp));
init();
dfs(1,0,0,7,7);
//printf("%d\n",dp[0][0][7][7]);
//printf("%f\n",arg);
printf("%.3f\n",sqrt(dp[1][0][0][7][7]*1.0/n-arg/n*(arg*1.0/n)));
}
return 0;
}
相关文章推荐
- 半同步/半异步体系结构模式
- 在gen文件夹下面为什么会有android.support.v7.appcompat这个包
- Spring的JdbcTemplate
- 指针常量
- C++ map 的用法(转载)
- poj2996 Help Me with the Game 结构体多级排序 桶排序
- 南阳题目60-谁获得了最高奖学金
- Matlab中的im2col函数
- [BZOJ2818] gcd - 欧拉函数+筛法
- Spring的AOP自动代理(含AspectJ的AOP开发)
- Class [org.springframework.transaction.config.TxNamespaceHandler] for namespace
- 修改整个app的字体
- 如何改变图片的分辨率
- 258. Add Digits
- 求一个数组中的第二大的数?
- arm_linux_书
- log level
- 深度优先搜索(DFS)
- test4 结对项目
- 且听风吟