P1005 矩阵取数游戏

题目描述

帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的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;
}