矩阵取数游戏
题目描述
帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n×m的矩阵,矩阵中的每个元素ai,j均为非负整数。游戏规则如下:
- 每次取数时须从每行各取走一个元素,共n个。经过m次后取完矩阵内所有元素;
- 每次取走的各个元素只能是该元素所在行的行首或行尾;
- 每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值 ×2^i,其中i表示第i次取数(从1开始编号);
- 游戏结束总得分为m次取数得分之和。
帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出取数后的最大得分。
输入输出格式
输入格式:
输入文件包括n+1行:
1行为两个用空格隔开的整数nn和mm。
第2∽n+1行为n×m矩阵,其中每行有m个用单个空格隔开的非负整数。
输出格式:
输出文件仅包含1行,为一个整数,即输入矩阵取数后的最大得分
一开始用了错误的贪心……
还是区间DP妙啊
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,a[85][85];
struct Big
{
static const int BASE=100000000;
long long s[100],len;
Big()
{
memset(s,0,sizeof(s));len=0;
}
Big operator +(const Big &b)
{
Big x;
for(int i=0;i<=len||i<=b.len;i++)
{
if(i<=len)x.s[i]+=s[i];
if(i<=b.len)x.s[i]+=b.s[i];
x.s[i+1]+=x.s[i]/BASE;
x.s[i]%=BASE;
}
x.len=max(len,b.len);
if(x.s[x.len+1])x.len++;
return x;
}
Big operator +(const int &b)
{
Big x=*this;x.s[0]+=b;
for(int i=0;i<=len;i++)
{
x.s[i+1]+=x.s[i]/BASE;
x.s[i]%=BASE;
}
if(x.s[x.len+1])x.len++;
return x;
}
Big operator *(int b)
{
Big x;
for(int i=0;i<=len;i++)
{
x.s[i]+=s[i]*b;
x.s[i+1]+=x.s[i]/BASE;
x.s[i]%=BASE;
}
x.len=len;
if(x.s[x.len+1])x.len++;
return x;
}
Big operator *(const Big &b)
{
Big x;
for(int i=0;i<=len;i++)
for(int j=0;j<=b.len;j++)
{
x.s[i+j]+=s[i]*b.s[j];
x.s[i+j+1]+=x.s[i+j]/BASE;
x.s[i+j]%=BASE;
}
x.len=len+b.len;
if(x.s[x.len+1])x.len++;
return x;
}
bool operator < (const Big &b)
{
if(len<b.len)return 1;
if(len>b.len)return 0;
for(int i=len;i>=0;i--)
{
if(s[i]<b.s[i])return 1;
if(s[i]>b.s[i])return 0;
}
return 0;
}
void print()
{
printf("%lld",s[len]);
for(int i=len-1;i>=0;i--)printf("%08lld",s[i]);
}
}tot,w[81],f[81][81];;
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)scanf("%d",&a[i][j]);
w[0].s[0]=1;
for(int i=1;i<=m;i++)w[i]=w[i-1]*2;
for(int t=1;t<=n;t++)
{
for(int i=0;i<=m;i++)
for(int j=0;j<=m;j++)
memset(f[i][j].s,0,sizeof(f[i][j].s));
for(int j=0;j<m;j++)
for(int i=1;i+j<=m;i++)
{
if(j==0)
{
f[i][i].s[0]=a[t][i]*2;continue;
}
Big r;
f[i][i+j]=f[i+1][i+j]*2+a[t][i]*2;
r=f[i][i+j-1]*2+a[t][i+j]*2;
if(f[i][j+i]<r)f[i][j+i]=r;
}
tot=tot+f[1][m];
}
tot.print();
return 0;
}