DP49 矩阵取数游戏（dp+高精度）

我要去看得最远的地方

已于 2022-09-22 18:42:06 修改

阅读量615

点赞数

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

于 2022-09-22 18:39:25 首次发布

本文链接：https://blog.csdn.net/qq_42434171/article/details/126997390

版权

动态规划专栏收录该内容

70 篇文章 0 订阅

订阅专栏

1.首先了解什么是高精度：(107条消息) c++ __int128_t使用_荼白777的博客-CSDN博客_c++ int128https://blog.csdn.net/weixin_45724872/article/details/120656141

2.这个结题思路和DP53取数游戏一样，就多加了个循环和高精度；

dp[L][R]:表示从左端取L个，从右端取R个

状态转移： dp[L][R]=max(dp[L][R],max(dp[L-1][R]+a[L]*power[L+R],dp[L][R-1]+a[n-R+1]*power[L+R]));

对每一行进行相同的处理，最后求和。

描述

帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的 n*m 的矩阵，矩阵中的每个元素 a_{ij}aij均为非负整数。游戏规则如下：
1.每次取数时须从每行各取走一个元素，共 n 个。m 次后取完矩阵所有元素；
2.每次取走的各个元素只能是该元素所在行的行首或行尾；
3.每次取数都有一个得分值，为每行取数的得分之和，每行取数的得分 = 被取走的元素值 * 2i，其中i表示第 i 次取数（从1开始编号）；
4.游戏结束总得分为 m 次取数得分之和。

帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。

输入描述：

第一行输入两个正整数 n 和 m ，表示矩阵的长宽。

后续 n 行每行输入 m 个正整数，表示矩阵的元素

输出描述：

输出最大得分

示例1

输入：

2 3
1 2 3
3 4 2

复制输出：

复制说明：

第1次：第1行取行首元素，第2行取行尾元素，本次得分为1 * 21 + 2 * 21 = 6
第2次：两行均取行首元素，本次得分为2 * 22 + 3 * 22 = 20

第3次：得分为3 * 23 + 4 * 23 = 56。

总得分为6 + 20 + 56 = 82

示例2

输入：

1 4
4 5 0 5

复制输出：

复制

示例3

输入：

2 10
96 56 54 46 86 12 23 88 80 43
16 95 18 29 30 53 88 83 64 67

复制输出：

#include <iostream>
using namespace std;
#include<string.h>
using namespace std;
using ll=__int128_t;
ll dp[103][103];
int a[105][105];
ll mod=1e9+7;
ll power2[102];
void power(int m)
{
    power2[0]=1;
    for(int i=1;i<=m;i++)
    {
        power2[i]=(power2[i-1]*2);
    }
    return;
}
int main() {
 int n,m;
//   scanf("%d%d",&n,&m);
  cin>>n>>m;
  for(int i=1;i<=n;i++)
  {
      for(int j=1;j<=m;j++)
      {
        cin>>a[i][j];
      }
  }
  power(m);
  ll ans=0;
  for(int row=1;row<=n;row++)
  {
      ll maxx=0;
      memset(dp,0,sizeof(dp));
      for(int L=0;L<=m;L++)
      {
          for(int R=0;L+R<=m;R++)
          {
              if(L==0 &&R>0)
              {
                  dp[L][R]=max(dp[L][R],(dp[L][R-1]+power2[L+R]*a[row][m-R+1]));
              }
              else if(L>0&&R==0)
              {
                    dp[L][R]=max(dp[L][R],(dp[L-1][R]+power2[L+R]*a[row][L]));
              }
              else if(L>0&&R>0)
              {
                  dp[L][R]=max(dp[L][R],max((dp[L][R-1]+power2[L+R]*a[row][m-R+1]),(dp[L-1][R]+power2[L+R]*a[row][L])));
              }
              dp[L][R]=dp[L][R];
          }
      }
      for(int l=0;l<=m;l++)
      {
        maxx=max(maxx,dp[l][m-l]);
      }
      ans+=maxx;
      ans=ans;
  }
    cout<<int(ans%mod)<<endl;
    return 0;
}
// 64 位输出请用 printf("%lld")