当前位置: 代码迷 >> 综合 >> [caioj]1493: 基于连通性状态压缩的 动态规划问题:Plan
  详细解决方案

[caioj]1493: 基于连通性状态压缩的 动态规划问题:Plan

热度:28   发布时间:2023-10-29 12:55:51.0

题目描述
【题目描述】
给你一个 n * m 的矩阵,矩阵中的每个格子都有一个分数,求从左上角走到右下角的最大分数,每个格子只能经过一次,也可以不经过。

 如图,一个 3*3 的矩阵,最大分数为61 。

【输入】
多组数据,每组数据的第一行有两个整数 n 和 m(1<= n <=8,1<= m <=9) ,下面给出一个 n * m 的矩阵,矩阵中每个分数的值在 -2000 到 2000之间。
【输出】
对于第 i 组数据,输出 Case i: ans ,ans 为最大分数。
【样例输入】
2 2
1 2
3 1
3 3
0 -20 100
1 -20 -20
1 1 1
【样例输出】
Case 1: 5
Case 2: 61

其实和裸题也差不多。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define LL long long
const int mod=100037;
int n,m,map[15][15],now=0,list[2][mod],num[2];
// 状态列表 状态总数 
LL state[2][mod],ans,val[15][15];
int HASH[mod];//每种状态的数目
LL get(int s,int p)//获取状态s从右向左数的第p位 
{
   return (s>>((p-1)*2))&3;}
void change(int &s,int p,int v)//把状态s的第p位改成v 
{s^=get(s,p)<<((p-1)*2);s^=(v)<<((p-1)*2);
}
void update(LL &x,LL y){x=max(x,y);}
void add(int now,int st,LL sum)
{int ss=st%mod;while(HASH[ss]!=-1&&list[now][HASH[ss]]!=st){ss++;ss%=mod;if(ss==0)ss=1;}if(HASH[ss]==-1){HASH[ss]=++num[now];list[now][num[now]]=st;state[now][num[now]]=sum;}else update(state[now][HASH[ss]],sum);
}
void work()
{state[0][1]=val[1][1];num[0]=1;list[0][1]=1;//一开始的状态必须设成这样,确保起点为(1,1) for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){now^=1;num[now]=0;memset(HASH,-1,sizeof(HASH));for(int k=1;k<=num[now^1];k++){int st=list[now^1][k];LL sum=state[now^1][k];int p=get(st,j);int q=get(st,j+1);int d=get(st,j+2);//为了避免重复计算某个格子价值,引入d这个括号 if(!p&&!q){add(now,st,sum);//因为不一定经过所有格子,所以有这种状态 if(map[i][j+1]&&map[i+1][j]){change(st,j,1);change(st,j+1,2);if(!d)add(now,st,sum+val[i][j]+val[i][j+1]+val[i+1][j]);else add(now,st,sum+val[i][j]+val[i+1][j]);//上面两个一定要加上val[i][j]啊!因为它的价值没有被累加过! }}else if(!p&&q){if(i==n&&j==m) {update(ans,sum);continue;}if(map[i][j+1])//从上面到达(n,m){if(!d)add(now,st,sum+val[i][j+1]);else add(now,st,sum);}if(map[i+1][j]){change(st,j,q);change(st,j+1,0);add(now,st,sum+val[i+1][j]);}}else if(p&&!q){if(i==n&&j==m)//从左面到达(n,m) {update(ans,sum);continue;}if(map[i+1][j])add(now,st,sum+val[i+1][j]);if(map[i][j+1]){change(st,j,0);change(st,j+1,p);if(!d)add(now,st,sum+val[i][j+1]);else add(now,st,sum);}}else if(p==2&&q==1){change(st,j,0);change(st,j+1,0);add(now,st,sum);}else if(p==1&&q==1){int top=1;for(int pos=j+2;pos<=m+1;pos++){int temp=get(st,pos);if(temp==1)top++;if(temp==2)top--;if(top==0){change(st,j,0);change(st,j+1,0);change(st,pos,1);add(now,st,sum);break;}}}else if(p==2&&q==2){int top=1;for(int pos=j-1;pos;pos--){int temp=get(st,pos);if(temp==2)top++;if(temp==1)top--;if(top==0){change(st,j,0);change(st,j+1,0);change(st,pos,2);add(now,st,sum);break;}}}}}for(int j=1;j<=num[now];j++)list[now][j]<<=2;}
}
int main()
{int Case=0;while(scanf("%d%d",&n,&m)!=EOF){ans=-(1LL<<31);now=0;memset(map,0,sizeof(map));for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){scanf("%lld",&val[i][j]);map[i][j]=1;}work();printf("Case %d: %lld\n",++Case,ans);}
}
  相关解决方案