题目描述 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(" "); } 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); }