bzoj·入门OJ·Sgu167 I-country

最新推荐文章于 2020-12-17 16:22:25 发布
最新推荐文章于 2020-12-17 16:22:25 发布
阅读量247 收藏
点赞数 2
分类专栏： 动态规划 文章标签： dp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_43326267/article/details/91872884
版权

初见安~本题解参考：lyd的《进阶》。

传送门为：入门OJ P4859，如原题那样输出方案……我的代码爆内存了QAQ暂时也想不到如何优化……但是在Acwing上，内存128MB是可以过的，所以也贴出来了【顺便求教如何把内存控制在64MBQwQ】

Description

在 N*M 的矩阵中，每个格子有一个权值，要求寻找一个包含 K 个格子的凸连通块（连通块中间没有空缺，并且轮
廓是凸的，如右图所示），使这个连通块中的格子的权值和最大。求出这个最大的权值和.N,M≤15,K≤225。

Input

第一行给出N，M，K，接下来给出这个矩阵

Output

如题

Sample Input

5 5 12
7 1 0 1 5 
6 8 0 4 10 
6 6 9 4 6 
9 2 9 1 0 
8 0 0 3 7

Sample Output

80

Sol

是个dp题。思路很好想，但是也很难想到——因为连通块的轮廓为凸，所以一定满足：对于连通块里的每一行而言，会有一个左边界和右边界越来越大的过程和逐步缩小的过程。所以易得：

若现在的状态是左右扩大，那么上一行的状态便同样为左右扩大或者没选；

若现在的状态是左扩大右缩小，则上一行应为左扩大右缩小或者左扩大右扩大【表示转折】

若现在的状态是左缩小右扩大，则上一行应为左扩大右扩大或者左缩小右扩大【表示转折】

若现在的状态是左右缩小，那么上一行的状态这里的四种均可。【除了不选】

过度过来的状态决定了枚举的上一行状态的区间左右边界。所以我们的dp里需要存的变量有：dp[i][j][l][r][x][y]，表示第i行，至此共选了j个，这一行选择的是第l~r个，x和y表示状态【纵坐标0增1减，所以左右扩大为[1, 0]】。接下来的递推就很简单了：

例如：.左右扩大

dp[i][j][l][r][1][0] = \sum _{p = l}^{r} {a[p]} + max_{1 <=p <= l <= r <= q <= m} \begin{Bmatrix} dp[i - 1][j - (r - l + 1)][p][q][1][0] \end{Bmatrix}

*若j == r - l + 1，则dp[i][j][l][r][1][0] = dp[i - 1][0][0][0][1][0]，即从这一行才开始选。

同理可得其他的三个表达式了。【懒得打公式QwQ】同时也为了方便处理，我们的a数组都存为前缀和，这样只要作差就可得到l~r的值。

最后的目标状态为：dp[i][K][l][r][x][y]。因为我们求解过程中的任意状态都是属于凸形，并且如果在某行中间就得到了K个方格，是不会地推到第n行的，所以还要枚举上来找。

先上代码了——【有点恶心，其实核心部分都是可以边写边复制粘贴的QwQ】

#include<bits/stdc++.h>
using namespace std;
int read()
{
	int x = 0, f = 1, ch = getchar();
	while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
	while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
	return x * f;
}

int n, m, K;
int w[20][20], dp[20][230][20][20][2][2];
int main()
{
	n = read(), m = read(), K = read();
	for(int i = 1; i <= n; i++) for(int j = 1; j <= m; j++) w[i][j] = w[i][j - 1] + read();
	
	for(register int i = 1; i <= n; i++)
		for(register int j = 0; j <= min(i * m, K); j++)
			for(register int l = 1; l <= m; l++)
				for(register int r = l; r <= m; r++)//0增1减 
				{
					register int len = r - l + 1;
					if(len > j) continue;//这一句不合法判定很重要！！！
					//左右扩大，1，0 
					if(len != j) for(register int p = l; p <= r; p++)//这里特判
						for(register int q = p; q <= r; q++)
							dp[i][j][l][r][1][0] = max(dp[i][j][l][r][1][0], dp[i - 1][j - len][p][q][1][0]);
					else dp[i][j][l][r][1][0] = dp[i - 1][0][0][0][1][0];
					//左扩大右收缩， 1,1
					for(register int p = l; p <= r; p++)
						for(register int q = r; q <= m; q++)
							dp[i][j][l][r][1][1] = max(dp[i][j][l][r][1][1], max(dp[i - 1][j - len][p][q][1][1], dp[i - 1][j - len][p][q][1][0]));
					//左收缩右扩大，0,0
					for(register int p = 1; p <= l; p++)
						for(register int q = l; q <= r; q++)
							dp[i][j][l][r][0][0] = max(dp[i][j][l][r][0][0], max(dp[i - 1][j - len][p][q][0][0], dp[i - 1][j - len][p][q][1][0]));
					//左右收缩，0，1 
					for(register int p = 1; p <= l; p++)
						for(register int q = r; q <= m; q++)
							dp[i][j][l][r][0][1] = max(max(dp[i][j][l][r][0][1], max(dp[i - 1][j - len][p][q][0][0], dp[i - 1][j - len][p][q][0][1])), max(dp[i - 1][j - len][p][q][1][0], dp[i - 1][j - len][p][q][1][1]));
					//上面的一大坨其实可以写个for循环0~1，看起来会简洁很多，后面关于方案输出的代码就是用for简化的。
					dp[i][j][l][r][0][0] += w[i][r] - w[i][l - 1];//最后再来加上
					dp[i][j][l][r][0][1] += w[i][r] - w[i][l - 1];
					dp[i][j][l][r][1][0] += w[i][r] - w[i][l - 1];
					dp[i][j][l][r][1][1] += w[i][r] - w[i][l - 1]; 
				}
				
	
	int ans = 0;
	for(int i = 1; i <= n; i++)
		for(int l = 1; l <= m; l++)
			for(int r = l; r <= m; r++)
				ans = max(max(ans, max(dp[i][K][l][r][0][0], dp[i][K][l][r][0][1])), max(dp[i][K][l][r][1][0], dp[i][K][l][r][1][1]));
        //枚举答案【其实也可以简化一下枚举状态的过程QAQ

	printf("%d\n", ans);
	return 0;
}
/*
5 5 12
7 1 0 1 5 
6 8 0 4 10 
6 6 9 4 6 
9 2 9 1 0 
8 0 0 3 7 
*/

再来就是输出方案的问题了——一般这种问题其实存一下前导状态然后一路递归回去输出就行了。只是说存的时候会有点麻烦，因为量有点多……

希望有大佬能解答一下如何在64MB以内完成存路径，输出路径的操作QwQ下面的代码是128MB可以过的：

#include<bits/stdc++.h>
using namespace std;
int read()
{
	int x = 0, f = 1, ch = getchar();
	while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
	while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
	return x * f;
}

int n, m, K;
long long w[20][20], dp[20][230][20][20][2][2];
int fa[20][230][20][20][2][2][5];//0 = j, 1 = l, 2 = r, 3 = op1, 4 = op2
void dfs(int i, int j, int l, int r, int op1, int op2)
{
	if(!i || !j) return;
	dfs(i - 1, fa[i][j][l][r][op1][op2][0], fa[i][j][l][r][op1][op2][1], fa[i][j][l][r][op1][op2][2], fa[i][j][l][r][op1][op2][3], fa[i][j][l][r][op1][op2][4]);
	for(int x = l; x <= r; x++) printf("%d %d\n", i, x);
}

int main()
{
	n = read(), m = read(), K = read();
	for(int i = 1; i <= n; i++) for(int j = 1; j <= m; j++) w[i][j] = w[i][j - 1] + read();
	for(register int i = 1; i <= n; i++)
		for(register int j = 0; j <= min(i * m, K); j++)
			for(register int l = 1; l <= m; l++)
				for(register int r = l; r <= m; r++)//0增1减 
				{
					register int len = r - l + 1;
					if(len > j) continue;
					//左右扩大，1，0 
					if(len != j) {//这里一定要打大括号……
					for(register int p = l; p <= r; p++)
						for(register int q = p; q <= r; q++)
							if(dp[i - 1][j - len][p][q][1][0] > dp[i][j][l][r][1][0])
								dp[i][j][l][r][1][0] = dp[i - 1][j - len][p][q][1][0], fa[i][j][l][r][1][0][0] = j - len, fa[i][j][l][r][1][0][1] = p, fa[i][j][l][r][1][0][2] = q, fa[i][j][l][r][1][0][3] = 1, fa[i][j][l][r][1][0][4] = 0;
					}//好长的说……后面因为涉及到是哪一个状态前导过来的所以只能枚举状态了。
					else dp[i][j][l][r][1][0] = dp[i - 1][0][0][0][1][0], fa[i][j][l][r][1][0][0] = 0, fa[i][j][l][r][1][0][1] = 0, fa[i][j][l][r][1][0][2] = 0, fa[i][j][l][r][1][0][4] = 0, fa[i][j][l][r][1][0][3] = 1;
					//左扩大右收缩， 1,1
					for(register int p = l; p <= r; p++)
						for(register int q = r; q <= m; q++)
							for(register int y = 0; y <= 1; y++)
								if(dp[i - 1][j - len][p][q][1][y] > dp[i][j][l][r][1][1])
									dp[i][j][l][r][1][1] = dp[i - 1][j - len][p][q][1][y], fa[i][j][l][r][1][1][0] = j - len, fa[i][j][l][r][1][1][1] = p, fa[i][j][l][r][1][1][2] = q, fa[i][j][l][r][1][1][3] = 1, fa[i][j][l][r][1][1][4] = y;
					//左收缩右扩大，0,0
					for(register int p = 1; p <= l; p++)
						for(register int q = l; q <= r; q++)
							for(register int x = 0; x <= 1; x++)
								if(dp[i - 1][j - len][p][q][x][0] > dp[i][j][l][r][0][0])
									dp[i][j][l][r][0][0] = dp[i - 1][j - len][p][q][x][0], fa[i][j][l][r][0][0][0] = j - len, fa[i][j][l][r][0][0][1] = p, fa[i][j][l][r][0][0][2] = q, fa[i][j][l][r][0][0][3] = x, fa[i][j][l][r][0][0][4] = 0;
					//左右收缩，0，1 
					
					for(register int p = 1; p <= l; p++)
						for(register int q = r; q <= m; q++)
							for(register int x = 0; x <= 1; x++)
								for(register int y = 0; y <= 1; y++)
									if(dp[i - 1][j - len][p][q][x][y] > dp[i][j][l][r][0][1])
										dp[i][j][l][r][0][1] = dp[i - 1][j - len][p][q][x][y], fa[i][j][l][r][0][1][0] = j - len, fa[i][j][l][r][0][1][1] = p, fa[i][j][l][r][0][1][2] = q, fa[i][j][l][r][0][1][3] = x, fa[i][j][l][r][0][1][4] = y;
					
					dp[i][j][l][r][0][0] += w[i][r] - w[i][l - 1];
					dp[i][j][l][r][0][1] += w[i][r] - w[i][l - 1];
					dp[i][j][l][r][1][0] += w[i][r] - w[i][l - 1];
					dp[i][j][l][r][1][1] += w[i][r] - w[i][l - 1]; 
					
				}
				
	register long long ans = 0;
	register int is, ls, rs, op1, op2;
	for(register int i = 1; i <= n; i++)
		for(register int l = 1; l <= m; l++)
			for(register int r = l; r <= m; r++)
				for(register int x = 0; x <= 1; x++)//这里也必须手动枚举状态
					for(register int y = 0; y <= 1; y++)
						if(dp[i][K][l][r][x][y] > ans)
							ans = dp[i][K][l][r][x][y], is = i, ls = l, rs = r, op1 = x, op2 = y;

	printf("Oil : %lld\n", ans);
	dfs(is, K, ls, rs, op1, op2);
	return 0;
}

迎评：）
——End——

樱狸❀
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
OJ的入门使用
06-16
ACM入门，OJ系统的使用及简介代码说明、 常见问题
【(伪)Top tree】【BZOJ3153 Sone1 题解】-附件资源
03-05
【(伪)Top tree】【BZOJ3153 Sone1 题解】-附件资源
【BZOJ 入门oj】1103: 那些N位数(dfs)
reverie_mjp的博客
08-01 971
1103: 那些N位数 Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 581  Solved: 232 [Submit][Status][Web Board] Description 一个n位数，只由1，2，3，4...p这几个数字组成。 请编写程序输出这些n位数，先小后大，每行一个。   Input 两个整数n，p(n＜9，p＜＝9，且...
【SGU167】I-country
Dusk_Of_Alchemist
04-27 488
Description在N*M的矩阵中，每个格子有一个权值，要求寻找一个包含K个格子的凸连通块（连通块中间没有空缺，并且轮廓是凸的，如右图所示），使这个连通块中的格子的权值和最大。求出这个最大的权值和，并给出连通块的具体方案。InputOn the first line of input there are 3 integer numbers N,M,K (1&lt;=N,M&lt;=15, 0&...
bzoj-入门OJ·简单的网络游戏
樱狸✨愿携清风归来处
06-12 299
初见安~这里是传送门【我知道这是一个很冷门的题QAQ】入门OJ P4033 Description 在某款极具技术含量的网络游戏中，佳佳靠着他的聪明智慧垄断了游戏中的油田系统。油田里有许多油井，这些油 井排成一个M*N的矩形。每个油井都有一个固定的采油量。每两个相邻的油井之间有一条公路，这些公路是油井与 油井之间唯一的运油方式。佳佳的领地在油田的右方和下方，他需要把采到的油通过这些公路运输到他...
bzoj入门4839 蒲公英
nbl97
10-09 276
题目： n个数m个询问，每个询问是一个[l,r]，输出[l,r]之前的众数。有相同的输出小的。n &lt;=40000,m&lt;=50000,ai&lt;=1e9 分析： 分块。 预处理出f[i][j]表示第i个块到第j个块的众数。对于一个询问[l,r]，中间整块的部分可以直接得到答案ans，然后对于两边的2*sqrt(n)个元素，查询出每个元素的出现次数，看是否能更新ans。至于，查找一个[l...
bzoj 1000-1499
11-26
bzoj部分数据.
[资料共享] BZOJ所有题 OJ 在线评测 题面 PDF版 非常良心-附件资源
03-05
[资料共享] BZOJ所有题 OJ 在线评测 题面 PDF版 非常良心-附件资源
bzoj·入门OJ·统计损失
樱狸✨愿携清风归来处
08-20 219
bzoj-入门OJ P2004 统计损失 树形dp
sgu:OJ-SGU-answer-backup
05-11
ID 标题 交流电 A + B 18881 互质 7697 第3分部 6906 总和 6185 日历 4336 画线 4159 987654321问题 4014 肉饼 3998 几乎素数 3845 a ^ b-b ^ a 3673 骨牌 3621 花小店 3417 数数 3290 电话目录 3282 数字根 2874 盒子 2836 276 安德鲁的烦恼 2658 130 圆圈 2555 154 阶乘 2529 114 广播电台 2476 222 小白嘴鸦 2354 116 超优质指数 2253 133 边界 2159 108 自我数字II 2149 144 会议 2121 175 编码方式 2032年 231 本金 2003年 134 质心 1931年 180 反演 1906年 106 等式 1884
分块入门 bzoj2141排队
Wall-E
08-09 192
To infinity and beyond. ——《Wall-E》   首先可以用离散化+树状数组来计算逆序对个数 时间复杂度 （n log n） 对于有关逆序对的题目，首先可以想到逆序对的特质： 然后对于每个询问(x,y) (x&lt;y) 两侧的数是没有影响的，区间(x,y)的数a[i]讨论如下： a[i]&lt;a[x] --ans a[i]&gt;a[x] ++ans a[i...
数据结构合并算法入门（BZOJ4719 天天爱跑步题解）
hzk_cpp的博客
09-07 411
一.合并算法的引入. 我们先来看一道例题. 题目：BZOJ4719. 题目大意：给定一棵nnn个点的树和树上每个节点的wiw_iwi​，现在给出mmm对sis_isi​和tit_iti​，表示从sis_isi​到tit_iti​会有一个人沿树上的路径走过，并且这个人每秒移动到下一个点.现在每个人都在时刻000走出，询问每一个点iii，请你输出第wiw_iwi​时刻有多少人会在点iii. 1≤n≤3...
bzoj3944 Sum(杜教筛入门/欧拉函数+莫比乌斯函数)
Code92007的博客
08-26 200
题目 T(T<=10)组样例，每次给出一个非负整数n(0<=n<=2的31次方-1) 求和 思路来源 https://blog.csdn.net/skywalkert/article/details/50500009唐老师博客例题 题解 n为2e9，预处理前的值，大概前2e6 以下图片来自唐老师博客： 卡常题目，得把一切能开int的地方都开int才能过 ...
统计损失
Jokercold的博客
09-07 808
统计损失 题目描述 SJY有一天被LLT紧急召去计算一些可能的损失。LLT元首管理的SHB国的交通形成了一棵树，现在将会出现一颗陨石砸在SHB国中，并且陨石砸毁的必定是SHB国构成的交通树上的一条路径。SHB国的损失可表示为被砸毁的路径上的所有城市价值之积。现在还暂时无法确定陨石的掉落路线，所以LLT元首希望SJY能够告诉他SHB国在受到每一种砸毁方式后会受到的损失之和模10086之后的值。注
BZOJ1654 奶牛舞会+tarjan算法（入门题目）+判断图中有几个环
地衣的博客——qdu
10-12 234
题目大意： 判断图中有几个环 #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int num=10010; /******************/ struct node{ int u,v,next; }e[num*5]; int hea...
OJ入门介绍
st_hunter的博客
07-07 805
OJ——online judge要求：正确处理输入格式；正确处理输出格式；输入输出均是标准输入/输出，不用考虑读写文件。1.正确处理输入格式 常见输入格式： 预先不输入数据的组数； 预先输入数据的组数； 只有一组输入数据应对： 读到文件结尾； 读数据组数再循环； 直接读入数据。预先不输入数据的组数 C++while(cin>>a>>b){ cout<<a+b<<endl;
算法学习笔记--OJ入门
santa12138的博客
12-17 354
从零开始 书籍：算法笔记 入门 OJ：online judgement，在线评测系统 单点测试与多点测试 单点测试 pat为单点测试 “一次性的写法”，执行一遍程序return 0即可 多点测试 codeup为多点测试 EOF 判断输入是否结束 while(scanf("%d",&a) != EOF){ } while break 当输入数据满足某个条件时停止 为上一种情况延伸 简洁写法 while(scanf("%d%d"),&a,&b),a||b){ } whi
BZOJ JOI 2013~2014 春季training合宿 系列题解1
最新发布
08-04
BZOJ JOI 2013~2014 春季training合宿 系列题解1

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值