P1005 [NOIP2007 提高组] 矩阵取数游戏区间dp

最新推荐文章于 2024-11-14 19:06:32 发布

漫漫人生路^_^

最新推荐文章于 2024-11-14 19:06:32 发布

阅读量102

点赞数

分类专栏： noip 文章标签：动态规划算法

本文为博主原创文章，未经博主允许不得转载

本文链接：https://blog.csdn.net/qq_46198929/article/details/115307032

版权

noip 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

题目链接
用num[i][L]代表行首元素，num[i][R]代表行末元素，f[L][R]代表区间[L,R]取数得到的最大得分。
转移方程：
f[L][R]=max(num[i][L]*p[k]+dp(i,L+1,R),dp(i,L,R-1)+num[i][R]*p[k])

错误代码
#include<bits/stdc++.h> 
using namespace std;
int n,m,i,j,a[100];
#define bll __int128
#define k m-(R-L)
int num[80][80];
bll ans,p[80],f[80][80];
void print(bll x)
{
    if(!x) return;
    if(x) print(x/10);
    putchar(x%10+'0');
}
bll dp(int i,int L,int R){
    if(f[L][R]!=-1) return f[L][R];
    if(R-L>=1) f[L][R]=max(num[i][L]*p[k]+dp(i,L+1,R),dp(i,L,R-1)+num[i][R]*p[k]);
    else f[L][R]=num[i][R]*p[k];
    return f[L][R];
}
int main(){
	p[0]=1;
	cin>>n>>m;
	for(j=1;j<=m;j++){
		p[j]=p[j-1]*2;
	}
	 for(i=1;i<=n;i++)
    {
    	for(j=1;j<=m;j++) 
		cin>>num[i][j];
    }
	for(i=1;i<=n;i++){
		memset(f,-1,sizeof(f));
		ans+=dp(i,1,m);
	}
	if(!ans) 
	printf("0");
    else 
	print(ans);
	return 0;
}

在这里插入图片描述
AC代码

#include<bits/stdc++.h> 
using namespace std;
int n,m,i,j,a[100];
#define bll __int128   //2^80超出了long long的范围,可以用__int128来存储
#define k m-(R-L)//每次都从每一行取出一个，所以共有m次取数，在区间[L+1,R]和区间[L,R-1]均取数R-L次
int num[90][90];//数组要开大
bll ans,p[90],f[90][90];
void print(bll x)
{
    if(!x) return;
    if(x) print(x/10);
    putchar(x%10+'0');
}
//记忆化
bll dp(int i,int L,int R){
    if(f[L][R]!=-1) return f[L][R];
    if(R-L>=1) f[L][R]=max(num[i][L]*p[k]+dp(i,L+1,R),dp(i,L,R-1)+num[i][R]*p[k]);   //取首个或末尾那个
    else f[L][R]=num[i][R]*p[k];
    return f[L][R];
}
int main(){
	p[0]=1;
	cin>>n>>m;
	for(j=1;j<=m;j++){
		p[j]=p[j-1]*2;
	}
	 for(i=1;i<=n;i++)
    {
    	for(j=1;j<=m;j++) 
		cin>>num[i][j];
    }
	for(i=1;i<=n;i++){
		memset(f,-1,sizeof(f));
		ans+=dp(i,1,m);//对每行取得的最大的值求和
	}
	if(!ans) 
	printf("0");
    else 
	print(ans);
	return 0;
}