codevs1166 矩阵取数游戏

题目描述 Description

【问题描述】
帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的n*m 的矩阵，矩阵中的每个元素aij均
为非负整数。游戏规则如下：
1. 每次取数时须从每行各取走一个元素，共n个。m次后取完矩阵所有元素；
2. 每次取走的各个元素只能是该元素所在行的行首或行尾；
3. 每次取数都有一个得分值，为每行取数的得分之和，每行取数的得分= 被取走的元素值*2i，
其中i 表示第i 次取数（从1 开始编号）；
4. 游戏结束总得分为m次取数得分之和。
帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。

输入描述 Input Description

第1行为两个用空格隔开的整数n和m。
第2~n+1 行为n*m矩阵，其中每行有m个用单个空格隔开的非负整数。

输出描述 Output Description

输出 仅包含1 行，为一个整数，即输入矩阵取数后的最大得分。

样例输入 Sample Input

2 3
1 2 3
3 4 2

样例输出 Sample Output

82

数据范围及提示 Data Size & Hint

样例解释

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

 

【限制】
60%的数据满足：1<=n, m<=30, 答案不超过1016
100%的数据满足：1<=n, m<=80, 0<=aij<=1000

唉没有高精就是sb题

加了个高精从30行变成110+行QAQ当高精模板用吧

这是不加高精的，60分

#include<cstdio>#include<iostream>#include<cstring>using namespace std;#define LL long longint n,m;LL f[100][100];LL a[100];LL ans;int main(){scanf("%d%d",&n,&m);while (n--){memset(f,0,sizeof(f));memset(a,0,sizeof(a));for (int i=1;i<=m;i++)  scanf("%lld",&a[i]);for (int i=1;i<=m;i++)  f[i][i]=a[i]*(1 << m);for (int i=m;i>=1;i--)  for (int j=i+1;j<=m;j++)  {  int s=m-(j-i);  LL ss=1<<s;    f[i][j]=max(f[i+1][j]+a[i]*ss,f[i][j-1]+a[j]*ss);  }ans+=f[1][m];}printf("%lld",ans);}

然后高精就变成这么长

#include<cstdio>#include<iostream>#include<cstring>#define mx 50using namespace std;#define LL long longstruct gaojing{ int len; int a[mx+10];}f[100][100],a[100],ans,mul[100];int n,m;inline void set0(gaojing &s)  {      s.len=1;    for (int i=1;i<mx+5;i++)s.a[i]=0;  }inline void inputn(gaojing &a){char ch=getchar();while (ch<'0'||ch>'9')ch=getchar();while (ch>='0'&&ch<='9'){a.a[a.len++]=ch-'0';ch=getchar();}a.len--;int change[1000];for (int i=1;i<=a.len;i++)  change[i]=a.a[i];for (int i=1;i<=a.len;i++)  a.a[i]=change[a.len-i+1];}inline void put(gaojing a)  {      for (int i=a.len;i>=1;i--)printf("%d",a.a[i]);      printf("\n");  }inline int cmp(const gaojing &a,const gaojing &b)//a<b:1  a>b:-1  a==b:0{if (a.len!=b.len){if (a.len<b.len)return 1;else return -1;}for (int i=a.len;i>=1;i--)  if(a.a[i]<b.a[i])return 1;  else if (a.a[i]>b.a[i])return -1;return 0;}inline gaojing max(const gaojing &a,const gaojing &b){int opr=cmp(a,b);if (opr==1)return b;else return a;}inline gaojing operator + (const gaojing &a,const gaojing &b){gaojing c;set0(c);    int maxlen=max(a.len,b.len);        for (int i=1;i<=maxlen;i++)        {            c.a[i]=c.a[i]+a.a[i]+b.a[i];            if (c.a[i]>=10)            {                c.a[i+1]+=c.a[i]/10;              c.a[i]%=10;          }      }        c.len=maxlen+4;        while (!c.a[c.len]&&c.len>1) c.len--;      return c;}inline gaojing operator * (const gaojing &a,const gaojing &b){gaojing c;set0(c);      for(int i=1;i<=a.len;i++)        for (int j=1;j<=b.len;j++)          c.a[i+j-1]+=a.a[i]*b.a[j];      int mxlen=a.len+b.len+10;    for (int i=1;i<=mxlen;i++)        {          c.a[i+1]+=c.a[i]/10;          c.a[i]%=10;        }      while (c.a[mxlen]==0)mxlen--;      c.len=mxlen;      return c;}int main(){scanf("%d%d",&n,&m);mul[0].len=1;mul[0].a[1]=1;for (int i=1;i<=m;i++)  mul[i]=mul[i-1]+mul[i-1];while (n--){for (int i=1;i<=m;i++)for (int j=1;j<=m;j++)set0(f[i][j]);for (int i=1;i<=m;i++)set0(a[i]);for (int i=1;i<=m;i++)inputn(a[i]);for (int i=1;i<=m;i++)  f[i][i]=a[i]*mul[m];for (int i=m;i>=1;i--)  for (int j=i+1;j<=m;j++)  {  int s=m-(j-i);  gaojing si=a[i]*mul[s];  gaojing sj=a[j]*mul[s];    f[i][j]=max(f[i+1][j]+si,f[i][j-1]+sj);  }ans=ans+f[1][m];}put(ans);}