noip2007 矩阵取数游戏 （动态规划+高精度）

P1378矩阵取数游戏

Accepted

标签： 动态规划 高精度 NOIP提高组2007

描述

帅帅经常更同学玩一个矩阵取数游戏：对于一个给定的n*m的矩阵，矩阵中的每个元素aij据为非负整数。游戏规则如下：
1. 每次取数时须从每行各取走一个元素，共n个。m次后取完矩阵所有的元素；
2. 每次取走的各个元素只能是该元素所在行的行首或行尾；
3. 每次取数都有一个得分值，为每行取数的得分之和；每行取数的得分 = 被取走的元素值*2^i，其中i表示第i次取数（从1开始编号）；
4. 游戏结束总得分为m次取数得分之和。

帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。

格式

输入格式

包括n+1行；
第一行为两个用空格隔开的整数n和m。

第2~n+1行为n*m矩阵，其中每行有m个用单个空格隔开
l<=n，m<=80，0<=aij<=1000

输出格式

仅包含1行，为一个整数，即输入矩阵取数后的最大的分。

样例1

样例输入1[复制]

2 3
1 2 3
3 4 2

样例输出1[复制]

82

限制

1s

提示

60%的数据满足：1<=n,m<=30,答案不超过10^16
100%的数据满足：1<=n,m<=80,0<=a<=1000

解析：用 f[i][j] 表示左边取 i 个，右边取 j 个的最优值，则：

            f[i][j]=max{f[i-1][j],f[i][j-1]}

        对于每一行，其最大得分即为 ：max{f[i][m-i]}；

代码：

#include<cstdio>
#include<algorithm>
#define maxn 80
#define maxn2 100
#define len 10000
using namespace std;

int n,m;
int a[maxn+20];
int ans[maxn2];
int f[maxn+20][maxn+20][maxn2];
int b[maxn+20][maxn2];
int p[maxn2],q[maxn2];

void multi(int x[],int y[],int z)
{
  int i,j,k,last=0;
  x[0]=y[0];
  for(k=1;k<=x[0];k++)
    {
      x[k]=y[k]*z+last;
      last=x[k]/len,x[k]%=len;
    }
  if(last>0)x[++x[0]]=last;  
}

void add(int x[],int y[])
{
  int i,j,k,last=0;
  k=min(x[0],y[0]);
  for(i=1;i<=k;i++)
    {
      x[i]+=y[i]+last;
	  last=x[i]/len,x[i]%=len;
    }
    
  if(x[0]>k)
    for(i=k+1;i<=x[0] && last>0;i++)
      x[i]+=last,last=x[i]/len,x[i]%=len;
  
  if(y[0]>k)
    for(x[0]=y[0],i=k+1;i<=y[0];i++)
	  x[i]=y[i]+last,last=x[i]/len,x[i]%=len;

  if(last>0)x[++x[0]]=last;	      
}

bool dayu(int x[],int y[])
{
  int i,j,k;
  if(x[0]>y[0])return 1;
  if(x[0]<y[0])return 0;
  for(i=x[0];i>=1;i--)
    {
      if(x[i]>y[i])return 1;
      if(x[i]<y[i])return 0;
    }
  return 0;  
}

void MAX(int x[],int y[],int z[])
{
  int i,j,k;
  if(dayu(y,z))
    for(i=0;i<=y[0];i++)x[i]=y[i];
  else 
    for(i=0;i<=z[0];i++)x[i]=z[i];	   
}

int main()
{  
  int i,j,k;
  scanf("%d%d",&n,&m);
  b[0][0]=1,b[0][1]=1;
  for(i=1;i<=m;i++)multi(b[i],b[i-1],2);
  
  for(k=1;k<=n;k++)
    {
      for(i=1;i<=m;i++)scanf("%d",&a[i]);
      
      f[0][0][0]=0;
      for(i=0;i<=m;i++)
        for(j=0;i+j<=m;j++)if(i>0 || j>0)
          {
          	p[0]=q[0]=0;
   	        
   	        if(i>0)
			  {
			    multi(p,b[i+j],a[i]);
                add(p,f[i-1][j]);
              }
            if(j>0)
			  {  
			    multi(q,b[i+j],a[m-j+1]);
                add(q,f[i][j-1]);
              }
			MAX(f[i][j],p,q);
          }
           
       for(p[0]=i=0;i<=m;i++)MAX(p,p,f[i][m-i]);
	   add(ans,p);  
    }
  
  printf("%d",ans[ans[0]]);  
  for(i=ans[0]-1;i>=1;i--)printf("%04d",ans[i]);
  printf("\n");  
  return 0;
}