题目描述

将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的两部分中的任意一块继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)

原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的平方和最小。

请编程对给出的棋盘及n,求出平方和的最小值。

解析

\(dp[k][i][j][x][y]\)表示当前在切第k刀,当前切\((i,j)\)为左上角,\((x,y)\)为右下角的矩形的最小平方和。

考虑状态转移,一个状态转移到下一个时,它可以继续往两个部分切下去,也就是往两个方向转移。

采用前缀和优化,其中\(calc(i,j,x,y)\)表示以\((i,j)\)为左上角,\((x,y)\)为右下角的矩形的和的平方。

\[dp[k][i][j][x][y]=min(dp[k][i][j][x][y],dp[k-1][i'+1][j][x][y]+calc(i,j,i',y),
\\dp[k-1][i][j][x'][y]+calc(x'+1,j,x,y))
\]

\[dp[k][i][j][x][y]=min(dp[k][i][j][x][y],dp[k-1][i][j'+1][x][y]+clac(i,j,x,j'),
\\dp[k-1][i][j][x][y']+calc(i,y',x,y))
\]

初始化显然是切第0刀时,切任意矩形的最小平方和就是不切,就是该矩形的和的平方。

之前的思路和这个稍稍有些不同,我定义dp数组为当前切剩下以\((i,j)\)为左上角,\((x,y)\)为右下角的矩形的最小平方和,那么就需要四个转移,但是我转移写炸了(摊。

参考代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<cstdlib>
#include<queue>
#include<vector>
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define N 11
#define MOD 2520
#define E 1e-12
using namespace std;
int dp[16][N][N][N][N],mp[N][N],n;
inline int read()
{
int f=1,x=0;char c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
inline int calc(int x1,int y1,int x2,int y2)
{
if(x1>x2) swap(x1,x2);
if(y1>y2) swap(y1,y2);
return (mp[x2][y2]-mp[x1-1][y2]-mp[x2][y1-1]+mp[x1-1][y1-1])*(mp[x2][y2]-mp[x1-1][y2]-mp[x2][y1-1]+mp[x1-1][y1-1]);
}
int main()
{
n=read();
memset(dp,0x3f,sizeof(dp));
memset(dp[0],0,sizeof(dp[0]));
for(int i=1;i<=8;++i)
for(int j=1;j<=8;++j){
dp[0][i][j][i][j]=mp[i][j]=read();
mp[i][j]+=mp[i-1][j]+mp[i][j-1]-mp[i-1][j-1];
}
for(int i=1;i<=8;++i)
for(int j=1;j<=8;++j)
for(int x=i;x<=8;++x)
for(int y=j;y<=8;++y)
dp[0][i][j][x][y]=calc(i,j,x,y);
for(int k=1;k<n;++k){
for(int i=1;i<=8;++i)
for(int j=1;j<=8;++j)
for(int x=i;x<=8;++x)
for(int y=j;y<=8;++y){
for(int dx=i;dx<=x;++dx)
dp[k][i][j][x][y]=min(dp[k][i][j][x][y],min(dp[k-1][i][j][dx][y]+calc(dx+1,j,x,y),dp[k-1][dx+1][j][x][y]+calc(i,j,dx,y)));
for(int dy=j;dy<=y;++dy)
dp[k][i][j][x][y]=min(dp[k][i][j][x][y],min(dp[k-1][i][j][x][dy]+calc(i,dy+1,x,y),dp[k-1][i][dy+1][x][y]+calc(i,j,x,dy)));
}
}
cout<<dp[n-1][1][1][8][8]<<endl;
return 0;
}

P1436 棋盘分割[dp]的更多相关文章

  1. 洛谷 P1436 棋盘分割 解题报告

    P1436 棋盘分割 题目描述 将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的两部分中的任意一块继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共 ...

  2. Luogu P1436 棋盘分割 暴力DP

    我的天,,,,,n=8,k<=15,,,这怕不是暴力DP+高维数组.... 开一个五维数组f[k][i][j][p][q]表示从(i,j)到(p,q)中分成k个矩形最小的平方和. 然后初始化时用 ...

  3. 【Luogu】P1436 棋盘分割 题解

    嗯,点开题目,哇!是一道闪亮亮的蓝题! 不要被吓到了,其实,这道题就是一个简单的DP啦! 我们设 \(f[x1][y1][x2][y2][c]\) 为以 \((x1,y1)\) 为左上角,以 \((x ...

  4. POJ 1191 棋盘分割(DP)

    题目链接 大体思路看,黑书...其他就是注意搞一个in数组,这样记忆化搜索,貌似比较快. #include <cstdio> #include <cstring> #inclu ...

  5. NOI1999 JZYZOJ1289 棋盘分割 dp 方差的数学结论

    http://172.20.6.3/Problem_Show.asp?id=1289 除了下标一坨一坨屎一样挺恶心其他都还挺容易的dp,这道题才发现scanf保留小数位是四舍五入的,惊了. f[k][ ...

  6. 洛谷P1436 棋盘分割

    洛谷题目链接 动态规划: 我们设状态$f[i][j][o][p][k]$表示一个矩形,左上角顶点坐标为$(i,j)$,右下角顶点坐标为$(o,p)$时分割了$k$次,也就是说现在是$k+1$块 我们考 ...

  7. poj 1191 棋盘分割(dp + 记忆化搜索)

    题目:http://poj.org/problem?id=1191 黑书116页的例题 将方差公式化简之后就是 每一块和的平方 相加/n , 减去平均值的平方. 可以看出来 方差只与 每一块的和的平方 ...

  8. [POJ] 1191 [LUOGU] P1436 棋盘分割

    那个均方差,可以通过展开.合并Σ,发现最终只有Xi^2会对答案造成影响,其他都是定值,所以求出最小的和的平方就行. 其实这才是这题最难的部分,以下都是码农部分. f[x1][y1][x2][y2][k ...

  9. HDU 2517 / POJ 1191 棋盘分割 区间DP / 记忆化搜索

    题目链接: 黑书 P116 HDU 2157 棋盘分割 POJ 1191 棋盘分割 分析:  枚举所有可能的切割方法. 但如果用递归的方法要加上记忆搜索, 不能会超时... 代码: #include& ...

随机推荐

  1. SpringBoot 跨域 Access-Control-Allow-Origin

    跨域(CORS)是指不同域名之间相互访问. 跨域,指的是浏览器不能执行其他网站的脚本,它是由浏览器的同源策略所造成的,是浏览器对于JavaScript所定义的安全限制策略. 只要协议,子域名,主域名, ...

  2. windows server 守护进程nssm

    安装nssm 下载程序 https://nssm.cc/ci/nssm-2.24-101-g897c7ad.zip 安装nssm,以管理员运行cmd.exe,执行如下: C:\Users\prolds ...

  3. 18 Ajax、Json以及jackson框架解析json的基本应用

    1. Ajax (1)概念:ASynchronous JavaScript And XML 异步的JavaScript 和 XML 异步和同步:客户端和服务器端相互通信的基础上 * 客户端必须等待服务 ...

  4. 17 JQuery高级----学习笔记

    1. 动画(1) 三种方式显示和隐藏元素 <1> 默认显示和隐藏方式 1. show([speed,[easing],[fn]]) 1. 参数: 1. speed:动画的速度.三个预定义的 ...

  5. 【转帖】处理器史话 | 服务器CPU市场的战役, AMD、Intel和ARM的厮杀

    处理器史话 | 服务器CPU市场的战役, AMD.Intel和ARM的厮杀 https://www.eefocus.com/mcu-dsp/377300   说完了个性鲜明的消费类电子,接下来聊一聊通 ...

  6. c++基础(四)—— 泛型算法

    1.find(first, last, value) 头文件:algorithm 参数:前两个参数是“表示元素范围的迭代器”,第三个是一个值 说明:find 将范围中进行寻找.搜索失败:如果范围中无匹 ...

  7. 19 IO流(十六)——Commons工具包,FileUtils(一)

    Commons包的API:自己查吧懒得传云 Commons包的导入方法 Commons是一个java的IO开源工具,导入方法: 从apache.org下载commons包 解压 copy其中的comm ...

  8. 数组中重复的数字(Python)

    #!/usr/bin/env python # -*- coding: utf-8 -*- # @Time : 2019-08-13 22:35 # @Author : daryl # @File : ...

  9. Layui + thymeleaf org.thymeleaf.exceptions.TemplateProcessingException: Could not parse as expression: " 错误解决

    解决方法: 1.将layui的代码移动到新的js文件当中,用template模板引擎的方式引入: <script th:src="@{/static/js/facility/mover ...

  10. docker查看容器日志

    原文:docker查看容器日志 前言 $ sudo docker logs -f -t --tail 行数 容器名 1 2 1.命令查看 root@c68d4b5dd583c4f4ea30da2989 ...