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 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);
}