棋盘分割 (POJ - 1191,二维区间 DP)

2024-03-30 13:48

本文主要是介绍棋盘分割 (POJ - 1191,二维区间 DP),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一.题目链接:

POJ-1191

二.题目大意:

中文题...

三.分析:

卧艹,这尼玛也太暴力了!

大佬讲解1   大佬讲解2

一开始没想明白状态表示,是我太菜.

dp[i][x1][y1][x2][y2] 表示 将区域 (x1, y1, x2, y2) 切 i 次后的子矩形平方和最小值.

搞清状态表示后,状态转移就很好写啦.

详见代码(没错,我就是懒得讲)

四.代码实现:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;const int M = 10;
const int inf = 0x3f3f3f3f;int sum[M][M];
int dp[M + 5][M][M][M][M];int cost(int x1, int y1, int x2, int y2)
{int s = sum[x2][y2] - sum[x2][y1 - 1] - sum[x1 - 1][y2] + sum[x1 - 1][y1 - 1];return s * s;
}int main()
{int n;scanf("%d", &n);for(int i = 1; i <= 8; ++i){for(int j = 1; j <= 8; ++j)scanf("%d", &sum[i][j]);}for(int i = 1; i <= 8; ++i){for(int j = 1; j <= 8; ++j)sum[i][j] += sum[i][j - 1] + sum[i - 1][j] - sum[i - 1][j - 1];}memset(dp, inf, sizeof(dp));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){dp[0][x1][y1][x2][y2] = cost(x1, y1, x2, y2);}}}}for(int i = 1; i <= n - 1; ++i){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){for(int x = x1; x < x2; ++x)dp[i][x1][y1][x2][y2]= min(dp[i][x1][y1][x2][y2],min(dp[i - 1][x1][y1][x][y2] + dp[0][x + 1][y1][x2][y2],dp[i - 1][x + 1][y1][x2][y2] + dp[0][x1][y1][x][y2]));for(int y = y1; y < y2; ++y)dp[i][x1][y1][x2][y2]= min(dp[i][x1][y1][x2][y2],min(dp[i - 1][x1][y1][x2][y] + dp[0][x1][y + 1][x2][y2],dp[i - 1][x1][y + 1][x2][y2] + dp[0][x1][y1][x2][y]));}}}}}printf("%.3f\n", sqrt(1.0 * dp[n - 1][1][1][8][8] / n - 1.0 * sum[8][8] * sum[8][8] / n / n));return 0;
}

 

这篇关于棋盘分割 (POJ - 1191,二维区间 DP)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/861651

相关文章

poj2576(二维背包)

题意:n个人分成两组,两组人数只差小于1 , 并且体重只差最小 对于人数要求恰好装满,对于体重要求尽量多,一开始没做出来,看了下解题,按照自己的感觉写,然后a了 状态转移方程:dp[i][j] = max(dp[i][j],dp[i-1][j-c[k]]+c[k]);其中i表示人数,j表示背包容量,k表示输入的体重的 代码如下: #include<iostream>#include<

hdu2159(二维背包)

这是我的第一道二维背包题,没想到自己一下子就A了,但是代码写的比较乱,下面的代码是我有重新修改的 状态转移:dp[i][j] = max(dp[i][j], dp[i-1][j-c[z]]+v[z]); 其中dp[i][j]表示,打了i个怪物,消耗j的耐力值,所得到的最大经验值 代码如下: #include<iostream>#include<algorithm>#include<

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int