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

最新推荐文章于 2024-08-11 23:00:35 发布

daq0411

最新推荐文章于 2024-08-11 23:00:35 发布

阅读量432

点赞数

分类专栏：题解文章标签： __int128 矩阵状态转移方程高精度计算贪心算法

本文链接：https://blog.csdn.net/yingtingqi/article/details/113913670

版权

题解专栏收录该内容

33 篇文章 0 订阅

订阅专栏

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

1、__int128的使用（不能直接cout 或者 printf，需要写函数）
2、状态转移方程：

if(i==j) dp[i][i] = map[k][i]*sum(2,m);
else{
	dp[i][j] = max(dp[i+1][j]+map[k][i]*sum(2,m-j+i), dp[i][j-1]+map[k][j]*sum(2,m-j+i));
}

主要是高精的__int128的问题，我们也可以用贪心的方法进行空间的压缩，当前行最大，加和之后也就是最大。
完整代码如下：

#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
const int maxn = 210;
int n,m;
__int128 dp[maxn][maxn];
int map[maxn][maxn];
__int128 ans = 0;
__int128 sum(int a,int n){
	if(n==1) return a;
	__int128 Ans = sum(a,n/2);
	if(n&1) return Ans * Ans * a;
	else return Ans * Ans;
}
void print(__int128 x){
    if(!x) return;
    if(x) print(x/10);
    putchar(x%10+'0');
}
int main(){
	cin >> n >> m;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			cin >> map[i][j];
		}
	}
	for(int k=1;k<=n;k++){
		memset(dp,0,sizeof(dp));
		for(int i=m;i>=1;i--){
			for(int j=i;j<=m;j++){
				if(i==j) dp[i][i] = map[k][i]*sum(2,m);
				else{
					dp[i][j] = max(dp[i+1][j]+map[k][i]*sum(2,m-j+i), dp[i][j-1]+map[k][j]*sum(2,m-j+i));
				}
			}
		}
		ans += dp[1][m];
	}
	if(!ans) printf("0");
    else print(ans);
	return 0;	
}