洛谷1005 矩阵取数游戏（dp+高精）

  题目描述 Description

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

NOIP 2007 提高第三题

2 3
1 2 3
3 4 2

82

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <cstdio>
#include <cctype>
#include <algorithm>
#include <memory>
using namespace std;
const int maxn=82;
int ans[maxn],p[maxn],n,m,f[maxn][maxn][32],er[maxn][maxn];
int fpi[maxn],spi[maxn];
int fsum[maxn],ssum[maxn];
inline int get(){
    char c;while(!isdigit(c=getchar()));
    int v=c-48;while(isdigit(c=getchar()))v=v*10+c-48;
    return v;
}
inline void mul(int *a,int *b,int x){
    int c[maxn];
    memset(c,0,sizeof(c));
    c[0]=b[0];
    for(int i=1;i<=c[0];++i){
	    c[i]+=b[i]*x;
	    c[i+1]+=c[i]/10;
	    c[i]%=10;
	}
	while(c[c[0]+1])++c[0],c[c[0]]+=c[c[0]-1]/10,c[c[0]-1]%=10;
	for(int i=0;i<=30;++i)a[i]=c[i];
}
inline void add(int *q,int *a,int *b){
    int c[maxn];
    memset(c,0,sizeof(c));
    c[0]=max(a[0],b[0]);
    for(int i=1;i<=c[0];++i){
	    c[i]+=a[i]+b[i];
	    c[i+1]+=c[i]/10;
	    c[i]%=10;
	}
	if(c[c[0]+1])++c[0];
	for(int i=0;i<=30;++i)q[i]=c[i];//不知道为什么要到30（或>）而不是c[ 0 ]，挂了一个点
}
inline bool mx(int *a,int *b){
    if(a[0]>b[0])return 1;
    else{
	    if(a[0]<b[0])return 0;
	    else
		    for(int i=a[0];i>=1;--i){
			    if(a[i]>b[i])return 1;
			    if(a[i]<b[i])return 0;
			}
	}
	return 1;
}
int main(){
    n=get();m=get();
    memset(ans,0,sizeof(ans));
    memset(er,0,sizeof(er));
    er[0][0]=er[0][1]=1;
    ans[0]=1;
    for(int i=1;i<=m+1;++i)mul(er[i],er[i-1],2);
    for(int i=1;i<=n;++i){
	    for(int j=1;j<=m;++j){
		    p[j]=get();
		    mul(f[j][j],er[m],p[j]);
		}
		for(int k=2;k<=m;++k){
		    for(int l=1;l<=m-k+1;++l){
			    int r=l+k-1;
			    memset(fpi,0,sizeof(fpi));
			    memset(spi,0,sizeof(spi));
			    memset(fsum,0,sizeof(fsum));
			    memset(ssum,0,sizeof(ssum));
			    mul(fpi,er[m-k+1],p[l]);
			    add(fsum,f[l+1][r],fpi);
			    mul(spi,er[m-k+1],p[r]);
			    add(ssum,f[l][r-1],spi);
			    if(mx(fsum,ssum))memcpy(f[l][r],fsum,sizeof(int)*31);
			    else memcpy(f[l][r],ssum,sizeof(int)*31);
			}
		}
		add(ans,ans,f[1][m]);
	}
	for(int i=ans[0];i>=1;--i)printf("%d",ans[i]);printf("\n");
	return 0;
}

思路：刚开始想的贪心是错的，本来认为每次取头尾二者中的最小值，可结果不一定最优。已知最终结果等于每行最大结果的和，想到dp，即f[ i ][ j ]表示从i到j的数最大解=max（f[ i ][ j-1] +a[j]*2^x，f[ i+1 ][ j ]+a[ i ]*2^x）；f[ 1 ][ m ]即为该行最大解。再用高精。