题目描述
帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的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 提高第三题
题解
从DP角度来说,是一道简单题;
对于每一层,层与层之间的各种状态互不影响。所以,问题简单了,求每一层的解 其sigma为答案
f[i][j] 表示 在区间[i][j]取数的最优解 f[i][j]=max(f[i+1][j]+a[i](2^当前次幂),f[i][j-1]+a[j](2^当前次幂))
可以记忆化搜索
然而本题最终数值过大 需要高精 不用高精的能过60
60分代码
#include <iostream>
#include <cstring>
using namespace std;
unsigned long long n,m,f[1111][1111],map[1111][1111],ans;
unsigned long long dfs(int x,int l,int r,unsigned long long s)
{
if(f[l][r]) return f[l][r];
if(l==r)
{
f[l][r]=map[x][l]*s;
return f[l][r];
}
f[l][r]=max(dfs(x,l,r-1,s*2)+map[x][r]*s,dfs(x,l+1,r,s*2)+map[x][l]*s);
return f[l][r];
}
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
cin>>map[i][j];
}
}
for(int i=1;i<=n;i++)
{
memset(f,0,sizeof(f));
ans+=dfs(i,1,m,2);
}
cout<<ans<<endl;
}
高精的话 则有f[i][j][k] 第三维相当于高精度的每一位
并且我进行了压位(压四位);
并在开始预处理出了2的n次幂
下面是既恶心又恶心而且恶心的高精。。。。。
#include <iostream>
#include <cstring>
using namespace std;
#define pp 10000
#define ll long long
int n,m,f[111][111][111],map[1111][1111],ans[111],z[111][1111];
int gd1[111],gd2[111],gd3[111],gd4[111],gd5[111];
int bo[111][111];
int init()//初始化到 2的m次幂
{
z[1][1]=2;z[1][0]=1;
for(int ii=2;ii<=m;ii++)
{
int tot=z[ii-1][0]+1;
for(int i=1;i<=z[ii-1][0];i++)
{
z[ii][i]=z[ii-1][i]*2;
}
for(int i=1;i<=tot;i++)
{
z[ii][i+1]+=z[ii][i]/pp;
z[ii][i]%=pp;
}
while(!z[ii][tot]&&tot>1) tot--;
z[ii][0]=tot;
}
}
int print(int * a)//输出数组a中的数 可用于debug
{
for(int i=a[0];i>=1;i--)
{
if(i==a[0])
{
//printf("%d",a[i]);
cout<<a[i];
continue;
}
int now=a[i];
int xx[6];
xx[4]=now%10;
now/=10;
xx[3]=now%10;
now/=10;
xx[2]=now%10;
now/=10;
xx[1]=now%10;
now/=10;
cout<<xx[1]<<xx[2]<<xx[3]<<xx[4];
//printf("%d%d%d%d%d",xx[1],xx[2],xx[3],xx[4]);
}
}
int gjj(int *a,int *b,int *c)//把b数组与c数组的数相加 赋给a数组
{
for(int i=1;i<=max(b[0],c[0]);i++)
{
a[i]=b[i]+c[i];
}
a[0]=max(b[0],c[0])+1;
for(int i=1;i<=a[0];i++)
{
a[i+1]+=a[i]/pp;
a[i]%=pp;
}
while(!a[a[0]]&&a[0]>1) a[0]--;/****要a[0]大于1 有一组数值ans==0 不大于1 会无输出。。。。*********/
}
int gjc(int *a,int x,int b)//把x与2的b次方的乘积赋值给a数组
{
//memset()
a[0]=z[b][0]+4;
//for(int i=1;i<=a[0];i++)
// a[i]=0;
//cout<<a[2]<<" a 2 ";
for(int i=1;i<=z[b][0];i++)
{
a[i]=z[b][i]*x;
//cout<<a[i]<<' '<<i<<' '<<z[b][i]<<' '<<x<<' '<<i<<' '<<z[b][0]<<' '<<a[2]<<endl;
}
for(int i=1;i<=a[0];i++)
{
// cout<<a[i]<<' '<<i<<' '<<pp<<' '<<a[i+1]<<endl;
a[i+1]+=a[i]/pp;
a[i]%=pp;
}
while(!a[a[0]]&&a[0]>1) a[0]--;
//cout<<"fdfdfdfdfd"<<endl;
//for(int i=a[0];i>=1;i--)
// cout<<a[i];
//cout<<endl<<x<<' '<<b<<' '<<"fdfdfdf"<<' '<<z[b][1]<<' '<<z[b][1]*x<<endl;
}
int max(int *a,int *b,int *c)//把b数组与c数组中数值大的赋值给a数组
{
int pd=0;
if(b[0]>c[0]) pd=1;
if(b[0]<c[0]) pd=2;
if(b[0]==c[0])
{
int h=b[0];
while(h>=1)
{
if(b[h]>c[h])
{
pd=1;break;
}
if(b[h]<c[h])
{
pd=2;
break;
}
h--;
}
}
if(pd==1)//DEBUG!!!!!!!! pd=1!!!!!!!@!$@!$@!
for(int i=0;i<=b[0];i++)
a[i]=b[i];
else
for(int i=0;i<=c[0];i++)
a[i]=c[i];
while(!a[a[0]]&&a[0]>1) a[0]--;
/*print(b);
cout<<"~~~"<<' '<<b[0]<< ' '<<c[0]<<' '<<pd<<' '<<b[1]<<' '<<c[1]<<endl;
print(a);
cout<<"~~~~~~~~~!!"<< ' ';
print(c);*/
}
int dfs(int x,int l,int r,int s)//记忆化搜索
{
if(bo[l][r]) return 0;//记忆化
bo[l][r]=1;
if(l==r)
{
gjc(f[l][r],map[x][l],s);/****相当于下面一行的状态转移方程*******/
// f[l][r]=map[x][l]*s;
// print(f[l][r]);
//cout<<f[l][r][1]<<' '<<l<<' '<<r<<' '<<map[x][l]<<' '<<s<<' '<<z[s][1]<<endl;
// return f[l][r];
return 0;
}
memset(gd1,0,sizeof(gd1));/*******************/
memset(gd2,0,sizeof(gd2));/*初始化很重要**/
memset(gd3,0,sizeof(gd3));/**!!!!!*******/
memset(gd4,0,sizeof(gd4));/*******************/
dfs(x,l,r-1,s+1);
dfs(x,l+1,r,s+1);
memset(gd1,0,sizeof(gd1));/*******************/
memset(gd2,0,sizeof(gd2));/*初始化很重要->AC**/
memset(gd3,0,sizeof(gd3));/**!!!!!*******/
memset(gd4,0,sizeof(gd4));//不能只在上面初始化 因为中间有dfs 和递归过程 期间 gd1 2 3 4 数组的值早已改变!!!!
gjc(gd1,map[x][r],s); ///***********************/
gjc(gd2,map[x][l],s); // ~~~~**********************/
gjj(gd3,f[l][r-1],gd1); //~~~~**********************/
gjj(gd4,f[l+1][r],gd2); //~~~ *********************/
max(f[l][r],gd3,gd4); //~~~此五行 相当于下面 无高精的状态转移方程
//f[l][r]=max(dfs(x,l,r-1,s*2)+map[x][r]*s,dfs(x,l+1,r,s*2)+map[x][l]*s);
//return f[l][r];
}
int main()
{
cin>>n>>m;
init();//预处理 2的n次方
//gjc(f[1][1],68,7);
/*map[1][1]=99;
dfs(1,1,1,80);*/
/*for(int i=40;i>=1;i--)
{
int now=f[1][1][i];
int xx[6];
xx[4]=now%10;
now/=10;
xx[3]=now%10;
now/=10;
xx[2]=now%10;
now/=10;
xx[1]=now%10;
now/=10;
cout<<xx [1]<<xx [2]<<xx [3]<<xx [4];
}*/
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
cin>>map[i][j];
}
}
for(int i=1;i<=n;i++)
{
memset(f,0,sizeof(f));//
memset(bo,0,sizeof(bo));//对每一层for i的循环 需要初始化 用来记忆的bo数组 以及每层的f数组(滚动)
dfs(i,1,m,1);
memset(gd5,0,sizeof(gd5));
for(int ii=0;ii<=ans[0];ii++)
gd5[ii]=ans[ii];
gjj(ans,gd5,f[1][m]);
/*print(f[1][1]);
print(f[2][2]);
print(f[3][3]);*/
// print(f[1][3]);
// cout<<endl;
}
print(ans);
}
无debug的注释的
#include <iostream>
#include <cstring>
using namespace std;
#define pp 10000
#define ll long long
int n,m,f[111][111][111],map[1111][1111],ans[111],z[111][1111];
int gd1[111],gd2[111],gd3[111],gd4[111],gd5[111];
int bo[111][111];
int init()//初始化到 2的m次幂
{
z[1][1]=2;z[1][0]=1;
for(int ii=2;ii<=m;ii++)
{
int tot=z[ii-1][0]+1;
for(int i=1;i<=z[ii-1][0];i++)
{
z[ii][i]=z[ii-1][i]*2;
}
for(int i=1;i<=tot;i++)
{
z[ii][i+1]+=z[ii][i]/pp;
z[ii][i]%=pp;
}
while(!z[ii][tot]&&tot>1) tot--;
z[ii][0]=tot;
}
}
int print(int * a)//输出数组a中的数 可用于debug
{
for(int i=a[0];i>=1;i--)
{
if(i==a[0])
{
cout<<a[i];
continue;
}
int now=a[i];
int xx[6];
xx[4]=now%10;
now/=10;
xx[3]=now%10;
now/=10;
xx[2]=now%10;
now/=10;
xx[1]=now%10;
now/=10;
cout<<xx[1]<<xx[2]<<xx[3]<<xx[4];
}
}
int gjj(int *a,int *b,int *c)//把b数组与c数组的数相加 赋给a数组
{
for(int i=1;i<=max(b[0],c[0]);i++)
{
a[i]=b[i]+c[i];
}
a[0]=max(b[0],c[0])+1;
for(int i=1;i<=a[0];i++)
{
a[i+1]+=a[i]/pp;
a[i]%=pp;
}
while(!a[a[0]]&&a[0]>1) a[0]--;/****要a[0]大于1 有一组数值ans==0 不大于1 会无输出。。。。*********/
}
int gjc(int *a,int x,int b)//把x与2的b次方的乘积赋值给a数组
{
a[0]=z[b][0]+4;
for(int i=1;i<=z[b][0];i++)
{
a[i]=z[b][i]*x;
}
for(int i=1;i<=a[0];i++)
{
a[i+1]+=a[i]/pp;
a[i]%=pp;
}
while(!a[a[0]]&&a[0]>1) a[0]--;
}
int max(int *a,int *b,int *c)//把b数组与c数组中数值大的赋值给a数组
{
int pd=0;
if(b[0]>c[0]) pd=1;
if(b[0]<c[0]) pd=2;
if(b[0]==c[0])
{
int h=b[0];
while(h>=1)
{
if(b[h]>c[h])
{
pd=1;break;
}
if(b[h]<c[h])
{
pd=2;
break;
}
h--;
}
}
if(pd==1)//DEBUG!!!!!!!! pd=1!!!!!!!@!$@!$@!
for(int i=0;i<=b[0];i++)
a[i]=b[i];
else
for(int i=0;i<=c[0];i++)
a[i]=c[i];
while(!a[a[0]]&&a[0]>1) a[0]--;
}
int dfs(int x,int l,int r,int s)//记忆化搜索
{
if(bo[l][r]) return 0;//记忆化
bo[l][r]=1;
if(l==r)
{
gjc(f[l][r],map[x][l],s);/****相当于状态转移方程*******/
return 0;
}
memset(gd1,0,sizeof(gd1));/*******************/
memset(gd2,0,sizeof(gd2));/*初始化很重要**/
memset(gd3,0,sizeof(gd3));/**!!!!!*******/
memset(gd4,0,sizeof(gd4));/*******************/
dfs(x,l,r-1,s+1);
dfs(x,l+1,r,s+1);
memset(gd1,0,sizeof(gd1));/*******************/
memset(gd2,0,sizeof(gd2));/*初始化很重要->AC**/
memset(gd3,0,sizeof(gd3));/**!!!!!*******/
memset(gd4,0,sizeof(gd4));//不能只在上面初始化 因为中间有dfs 和递归过程 期间 gd1 2 3 4 数组的值早已改变!!!!
gjc(gd1,map[x][r],s); ///***********************/
gjc(gd2,map[x][l],s); // ~~~~**********************/
gjj(gd3,f[l][r-1],gd1); //~~~~**********************/
gjj(gd4,f[l+1][r],gd2); //~~~ *********************/
max(f[l][r],gd3,gd4); //~~~此五行 相当于无高精的状态转移方程
}
int main()
{
cin>>n>>m;
init();//预处理 2的n次方
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
cin>>map[i][j];
}
}
for(int i=1;i<=n;i++)
{
memset(f,0,sizeof(f));
memset(bo,0,sizeof(bo));//对每一层for i的循环 需要初始化 用来记忆的bo数组 以及每层的f数组(滚动)
dfs(i,1,m,1);
memset(gd5,0,sizeof(gd5));
for(int ii=0;ii<=ans[0];ii++)
gd5[ii]=ans[ii];
gjj(ans,gd5,f[1][m]);
}
print(ans);
}