题目描述 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 long
int 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 50
using namespace std;
#define LL long long
struct 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);
}