题目描述
帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数。游戏规则如下:
1.每次取数时须从每行各取走一个元素,共n个。m次后取完矩阵所有元素;
2.每次取走的各个元素只能是该元素所在行的行首或行尾;
3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值*2^i,其中i表示第i次取数(从1开始编号);
4.游戏结束总得分为m次取数得分之和。
帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出取数后的最大得分。
输入输出格式
输入格式:
输入文件game.in包括n+1行:
第1行为两个用空格隔开的整数n和m。
第2~n+1行为n*m矩阵,其中每行有m个用单个空格隔开的非负整数。
数据范围:
60%的数据满足:1<=n, m<=30,答案不超过10^16
100%的数据满足:1<=n, m<=80,0<=aij<=1000
输出格式:
输出文件game.out仅包含1行,为一个整数,即输入矩阵取数后的最大得分。
输入输出样例
输入样例#1:
2 3
1 2 3
3 4 2
输出样例#1:
82
说明
NOIP 2007 提高第三题
C++:
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int n,m,s[82][82],ans[100],f[82][82][100],l[82][82];
void ccj(int jj,int i,int j){
int y[100],yy[100],l1=0,g=s[jj][j],l2=l[j+1][j+i-1],l3=0;
memset(y,0,sizeof(y));memset(yy,0,sizeof(yy));
for (int a1=0;a1<l[j+1][j+i-1];a1++)
y[a1]=f[j+1][j+i-1][a1];
while(g)
{yy[l1++]=g%10;g/=10;}
int yyy=m-i+1,a[8];
for (int a1=0;a1<8;a1++){if (yyy>10){yyy-=10;a[a1]=1<<10;}else {a[a1]=1<<yyy;yyy=0;}}
for (int a1=0;a1<8;a1++)
{
int k=0;if (a[a1]<2)break;
for (int a2=0;a2<l1||k;a2++,l1=max(a2,l1))
{
yy[a2]=yy[a2]*a[a1]+k;
k=yy[a2]/10;
yy[a2]%=10;
}
}
int k=0;
for (int a1=0;a1<max(l1,l2)||k;a1++,l3=a1)
{y[a1]=y[a1]+yy[a1]+k;k=y[a1]/10;y[a1]%=10;}
int x[100],xx[100],ll1=0,gg=s[jj][j+i-1],ll2=l[j][j+i-2],ll3=0;
memset(x,0,sizeof(x));memset(xx,0,sizeof(xx));
for (int a1=0;a1<l[j][j+i-2];a1++)
x[a1]=f[j][j+i-2][a1];
while(gg)
{xx[ll1++]=gg%10;gg/=10;}
for (int a1=0;a1<8;a1++)
{
int k=0;if (a[a1]<2)break;
for (int a2=0;a2<ll1||k;a2++,ll1=max(ll1,a2))
{
xx[a2]=xx[a2]*a[a1]+k;
k=xx[a2]/10;
xx[a2]%=10;
}
}
k=0;
for (int a1=0;a1<max(ll1,ll2)||k;a1++,ll3=a1)
{x[a1]=x[a1]+xx[a1]+k;k=x[a1]/10;x[a1]%=10;}
bool u=1;
if (l3!=ll3){if (l3>ll3)u=1;else u=0;}
else
{
for (int a1=l3-1;a1>=0;a1--)
{if (y[a1]>x[a1]){u=1;break;}if (y[a1]<x[a1]){u=0;break;}}
}
if (u)
{
l[j][j+i-1]=l3;
for (int a1=0;a1<l3;a1++)
f[j][j+i-1][a1]=y[a1];
}
else
{
l[j][j+i-1]=ll3;
for (int a1=0;a1<ll3;a1++)
f[j][j+i-1][a1]=x[a1];
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
scanf("%d",&s[i][j]);
for (int jj=1;jj<=n;jj++)
{
memset(l,0,sizeof(l));memset(f,0,sizeof(f));
for (int i=1;i<=m;i++)
for (int j=1;j<=m&&j+i-1<=m;j++)
ccj(jj,i,j);
for (int i=0;i<100;i++)
{ans[i]+=f[1][m][i];ans[i+1]+=ans[i]/10;ans[i]%=10;}
}
bool u=0;
for (int i=99;i>=0;i--)
{
if (ans[i])u=1;
if (u)printf("%d",ans[i]);
}
if (!u)printf("0");
return 0;
}