1574 矩阵取数游戏(NOIP2007 LOJ10152 LUOGU1005 提高+/省选-) 从外围到内区间动归ULL60分 高精度替代者int128 90分 补上int128输出0漏洞100分

已于 2022-01-25 09:18:36 修改
阅读量2.2k 收藏 1
点赞数
分类专栏： 信息学奥赛一本通 提高篇 文章标签： 矩阵 动态规划
于 2022-01-21 21:32:16 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mrcrack/article/details/122331203
版权

总目录

在线测评地址(ybt)

在线测评地址(LOJ)

在线测评地址(LUOGU)

看了数据范围，需高精度，不过，比赛中，是不编的，试试int128，先算算该题的数据最大值。

推导过程中，忘记了等比数列求和公式，自个重新推导了一遍。

n=80,m=80,aij=1000

80*1000*(2^1+2^2+2^3+......+2^80)

=8*10^4*(2-2*2^80)/(1-2)

=2^81*8*10^4

=193428131138340667952988160000
=2*10^29

int128最大值
2^127-1=1.7014118346046923173168730371588e+38

完全够用，
那么，该题 高精度靠边站，
用int128

区间动规，在一维层面处理惯了，二维层面极不适应。还是要从样例入手，找出规律。

看了一眼思路，发现完全想不到，那么套用思路，整理样例吧：

  • 求n行最大得分和，每一行取数又不会影响到其他行，那么只要确保每一行得分最大，管好自家孩子就行了。（这个在动规中叫最优子结构）
  • 每次取数是在边缘取，那么每次取数完剩下来的元素一定是在一个完整的一个区间中，又是求最优解，区间DP应运而生。
  • 我们用f{i,j}表示区间变为[i,j]时，获得的最大分数。
  • 当区间变为[i,j]时，上一次取数的时候区间一定是[i-1,j]或[i,j+1]从这两个状态转移即可。

样例2

输入:
1 4
4 5 0 5
输出:
122

位置1 2 3 4
数值4 5 0 5

取了0次
dp[1][4]=0

取了1次
dp[1][3]=5*2^1=10
dp[2][4]=4*2^1=8

取了2次
dp[1][2]=dp[1][3]+0*2^2=10
dp[2][3]=max(dp[1][3]+4*2^2,dp[2][4]+5*2^2)=max(26,28)=28
dp[3][4]=dp[2][4]+5*2^2=28

取了3次
dp[1][1]=dp[1][2]+5*2^3=50
dp[2][2]=max(dp[1][2]+4*2^3,dp[2][3]+0*2^3)=max(42,28)=42
dp[3][3]=max(dp[2][3]+5*2^3,dp[3][4]+5*2^3)=max(68,68)=68
dp[4][4]=dp[3][4]+0*2^3=28

取了4次，因无法找到合适的状态，故沿用取了3次的状态
dp[1][1]=50+4*2^4=114
dp[2][2]=42+5*2^4=122
dp[3][3]=68+0*2^4=68
dp[4][4]=28+5*2^4=108

取最大值，故输出122

先编写unsigned long long范围内的代码

有了上述模拟的基础，随便一遍，3个样例尽然都通过了，先放着，等int128编号后，再提交

以下是unsigned long long范围内的代码

ybt

未通过

测试点结果内存时间
测试点1答案正确620KB2MS
测试点2答案正确612KB2MS
测试点3答案正确620KB2MS
测试点4答案正确644KB2MS
测试点5答案正确644KB2MS
测试点6答案错误664KB2MS
测试点7答案错误664KB2MS
测试点8答案错误704KB3MS
测试点9答案错误716KB3MS
测试点10答案正确616KB2MS

LOJ

LUOGU

1.从外围到内区间动归ULL60分 

#include <bits/stdc++.h>
#define ULL unsigned long long
#define maxn 85
using namespace std;
ULL dp[maxn][maxn],a[maxn][maxn],two[maxn],ans=0,mx;
int main(){
	int n,m,i,j,lt,rt,row,cnt,k;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
			scanf("%llu",&a[i][j]);
	two[0]=1;
	for(i=1;i<=m;i++)two[i]=two[i-1]*2;
	for(row=1;row<=n;row++){//一行一行处理
		for(i=1;i<=m;i++)
			for(j=i;j<=m;j++)
				dp[i][j]=0;//初始化
		for(cnt=0;cnt<m;cnt++){//取数个数为cnt 
			for(k=0;k<=cnt;k++){
				lt=k+1;//左边取k个
				rt=m-(cnt-k);//右边取cnt-k个
				dp[lt][rt]=max(dp[lt-1][rt]+a[row][lt-1]*two[cnt],dp[lt][rt+1]+a[row][rt+1]*two[cnt]); 
			}
		}
		for(i=1;i<=m;i++)//取数个数为m时
			dp[i][i]+=a[row][i]*two[m];
		mx=0;
		for(i=1;i<=m;i++)
			mx=max(mx,dp[i][i]);
		ans+=mx; 
	}
	printf("%llu\n",ans);
	return 0;
}

有了上面代码的基础，int128的代码很快编好，3个样例也都通过。

ybt

未通过

测试点结果内存时间
测试点1答案错误596KB2MS
测试点2答案正确612KB2MS
测试点3答案正确648KB2MS
测试点4答案正确676KB2MS
测试点5答案正确684KB2MS
测试点6答案正确696KB2MS
测试点7答案正确736KB2MS
测试点8答案正确796KB3MS
测试点9答案正确804KB3MS
测试点10答案正确620KB1MS

LOJ

LUOGU

 

2.高精度替代者int128 90分 

#include <bits/stdc++.h>
#define maxn 85
using namespace std;
__int128 dp[maxn][maxn],a[maxn][maxn],two[maxn],ans=0,mx;
inline __int128 read(){//int128读取 
	__int128 x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//读取非数字字符 
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x*f;
}
inline void print(__int128 x){
	if(x<0){
		putchar('-');
		x=-x;
	}
	if(x){
		print(x/10);//递归写法
		putchar(x%10+'0'); 
	}
} 
int main(){
	int n,m,i,j,lt,rt,row,cnt,k;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
			a[i][j]=read();
	two[0]=1;
	for(i=1;i<=m;i++)two[i]=two[i-1]*2;
	for(row=1;row<=n;row++){//一行一行处理
		for(i=1;i<=m;i++)
			for(j=i;j<=m;j++)
				dp[i][j]=0;//初始化
		for(cnt=0;cnt<m;cnt++){//取数个数为cnt 
			for(k=0;k<=cnt;k++){
				lt=k+1;//左边取k个
				rt=m-(cnt-k);//右边取cnt-k个
				dp[lt][rt]=max(dp[lt-1][rt]+a[row][lt-1]*two[cnt],dp[lt][rt+1]+a[row][rt+1]*two[cnt]); 
			}
		}
		for(i=1;i<=m;i++)//取数个数为m时
			dp[i][i]+=a[row][i]*two[m];
		mx=0;
		for(i=1;i<=m;i++)
			mx=max(mx,dp[i][i]);
		ans+=mx; 
	}
	print(ans);
	return 0;
}

才90分，还有什么没考虑到？

仔细对了60分代码，发现第一个点是正确的，那么，还是__int128的读写有问题，需改进。

上述90分代码争对如下输入数据，没有输出

1 1
0

明显，上述争对__int128的输出有问题，0是无法输出的，要补上这个漏洞。此处需特判。

ybt

通过

测试点结果内存时间
测试点1答案正确612KB2MS
测试点2答案正确612KB1MS
测试点3答案正确648KB1MS
测试点4答案正确688KB2MS
测试点5答案正确676KB2MS
测试点6答案正确708KB2MS
测试点7答案正确728KB2MS
测试点8答案正确804KB3MS
测试点9答案正确836KB3MS
测试点10答案正确616KB1MS

LOJ

LUOGU

3.补上int128输出0漏洞100分 

#include <bits/stdc++.h>
#define maxn 85
using namespace std;
__int128 dp[maxn][maxn],a[maxn][maxn],two[maxn],ans=0,mx;
inline __int128 read(){//int128读取 
	__int128 x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//读取非数字字符 
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x*f;
}
inline void print(__int128 x){
	if(x<0){
		putchar('-');
		x=-x;
	}
	if(x){
		print(x/10);//递归写法
		putchar(x%10+'0'); 
	}
} 
int main(){
	int n,m,i,j,lt,rt,row,cnt,k;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
			a[i][j]=read();
	two[0]=1;
	for(i=1;i<=m;i++)two[i]=two[i-1]*2;
	for(row=1;row<=n;row++){//一行一行处理
		for(i=1;i<=m;i++)
			for(j=i;j<=m;j++)
				dp[i][j]=0;//初始化
		for(cnt=0;cnt<m;cnt++){//取数个数为cnt 
			for(k=0;k<=cnt;k++){
				lt=k+1;//左边取k个
				rt=m-(cnt-k);//右边取cnt-k个
				dp[lt][rt]=max(dp[lt-1][rt]+a[row][lt-1]*two[cnt],dp[lt][rt+1]+a[row][rt+1]*two[cnt]); 
			}
		}
		for(i=1;i<=m;i++)//取数个数为m时
			dp[i][i]+=a[row][i]*two[m];
		mx=0;
		for(i=1;i<=m;i++)
			mx=max(mx,dp[i][i]);
		ans+=mx; 
	}
	if(ans==0)putchar('0');//__int128的print函数的漏洞 
	else print(ans);
	return 0;
}

该题习得什么？

学习了从外围到内的区间动归。

进一步强化了unsigned long long的读取%llu

int128的读写更加的熟练。

补上了int128输出是0的漏洞。

mrcrack
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
1574：矩阵取数游戏
weixin_52308504的博客
02-25 1694
1574：矩阵取数游戏 时间限制: 1000 ms 内存限制: 524288 KB 提交数: 618 通过数: 248 【题目描述】 经常和同学玩一个矩阵取数游戏：对于给定的 n×m 的矩阵，矩阵中每个元素 aij 均为非负整数。游戏规则如下： 1、每次取数时必须从每行各取走一个元素，共 n 个，m 次取完所有元素。 2、每次取走的各个元素只能是该元素所在行行首或行尾。 3、每次取数都有一个的分值，为每行取数得分之和，每行取数得分==被取走元素值×2i ，其中 i 表示第 i 次取数，从 1 开始计数。 4
一本通1574矩阵取数游戏
weixin_30904593的博客
02-07 261
1574：矩阵取数游戏 时间限制: 1000 ms 内存限制: 524288 KB 题目描述 原题来自：NOIP 2007 帅帅经常和同学玩一个矩阵取数游戏：对于给定的n×m的矩阵，矩阵中每个元素aij​均为非负整数。游戏规则如下： 每次取数时必须从每行各取走一个元素，共n个，m次取完所有元素。 每次取走的各个元素只能是该元素所在行行首...
矩阵取数游戏
sundujing的博客
07-17 577
矩阵取数游戏
NOIP 2007 矩阵取数游戏 源代码
03-23
NOIP 2007 的题目矩阵取数游戏的pascal源代码，测试数据全部通过。
NOIP2004-2017提高组初赛试题集
06-10
《全国青少年信息学奥林匹克联赛（NOIP）2004-2017提高组初赛试题集》 全国青少年信息学奥林匹克联赛（National Olympiad in Informatics in Provinces，简称NOIP）是由中国计算机学会主办的一项面向全国中学生的...
NOIP2018提高组C++（试题+答案_初赛题目）
10-14
**全国青少年信息学奥林匹克联赛（NOIP）是针对中学生的一项编程竞赛，旨在培养和选拔在计算机编程方面具有潜力的学生。2018年的NOIP提高组比赛使用了C++语言，这是一门广泛应用于算法竞赛和软件开发的高级编程语言...
NOIP2019 提高组初赛 CSP-S C++试题及答案.rar
10-28
《全国青少年信息学联赛NOIP CSP-J/S历年C++试题解析》 全国青少年信息学联赛（NOIP）是一项旨在激发青少年计算机科学兴趣、提升编程能力的竞赛活动，分为初级组（CSP-J）和提高组（CSP-S）。CSP-J/S考试包括两轮，...
1、省选＋NOI-第一部分 动态规划DP_2020.08.29.pdf
11-25
1、省选＋NOI-第一部分 动态规划DP_2020.08.29.pdf
P1005 [NOIP2007 提高组] 矩阵取数游戏（区间dp+__int128）
sancpp的博客
02-25 363
传送门 思路不解释了：大佬的题解 看到这题__int128也能过，想拿__int128试试。 遇到的小坑： f[i][j] = max(f[i - 1][j] + (1 << (m - j + i - 1)) * d[i - 1], f[i][j]); f[i][j] = max(f[i][j + 1] + (1 << (m - j + i - 1)) * d[j + 1], f[i][j]); 上面的写法是错的（不知道怎么解释 ） f[i][j] = max(f[i - 1][j
<区间DP> 【noip 2007】 矩阵取数游戏（60%）
Loi_LXTT
09-28 287
*题目链接http://codevs.cn/problem/1166/ 因为剩下的40%是高精度而且是恶心的高精度，我就先打60%吧。。。 一开始想的转移方程错了，调了一会发现不满足无后效性原则，一开始的想法是这样的： dp[i][j]是第i次取数，取的数是第j个的最大的答案，枚举第i-1次取的数，但是这样做不能满足j之前没有被取过，因此不可以 正确的转移方程： 对于每一行： dp[i]
三分模板+int128模板
我们仍未知道那天所看见的花的名字。
11-16 266
while(l<r) { ll midl=l+(r-l)/3; ll midr=r-(r-l)/3; if(check(midl)<=check(midr)) r=midr-1; else l=midl+1; } #include <bits/stdc++.h> using namespace std; inline __int128 read() { __int128
【2019CSP-J普及组】T1 数字游戏
蒟蒻こんにゃく
04-19 601
P5660 数字游戏 题目传送门 没什么好说的，按照惯例，T1还是这么简单！ 方法一： #include<cstdio> #include<iostream> #include<cmath> #include<cstring> #include<string> #include<algorithm> #define fre(x...
矩阵取数游戏（动态规划）
wowhhhgd的Blog
06-10 2643
小明经常玩一个矩阵取数游戏：对于一个给定的n x n的矩阵，矩阵中的每个元素aij表示该格子的价值，均为非负整数。游戏规则如下： 小明从左上角走到右下角，只能向下向右走，经过某个格子，就能获得格子相应价值的奖励， 请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大奖励值。 输入格式: 输入包括n+1行： 第1行为一个整数n。(2 <= n <= 500) 第2~n+1行为nxn矩阵，其中每行有n个用单个空格隔开的非负整数。（0 <= aij <= 10000) 输出格式: 仅包含
关于CSP2021
Wu_while的博客
10-21 1271
今年（2021年），CCF发布了 NOI Linux 2.0 以及 关于NOI系列活动中编程语言使用限制的补充说明 ，将NOI系列比赛的标准进行了改动。
关于__int128
热门推荐
HumveeA6的博客
04-03 1万+
这种类型只有部分GCC编译器支持。经测试，正常运算与普通int无异，但是输入输出时无论是cin，cout还是printf都会报错，所以必须自己写输入输出函数。 定义时与别的数据类型并没有什么区别 eg: __int128 a,b,c; cin&amp;gt;&amp;gt;a&amp;gt;&amp;gt;b&amp;gt;&amp;gt;c; 方法一：重载输出 std::ostream&amp;amp; operator&amp;l
10、省选＋NOI-第十部分 概率统计与多项式_2020.08.29.pdf
11-25
10、省选＋NOI-第十部分 概率统计与多项式_2020.08.29.pdf

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值