BZOJ2310

2310: ParkII

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 86  Solved: 34
[Submit][Status][Discuss]

Description

Hnoi2007-Day1有一道题目 Park：给你一个 m * n 的矩阵，每个矩阵内有个
权值V(i,j) (可能为负数)，要求找一条回路，使得每个点最多经过一次，并且经过
的点权值之和最大，想必大家印象深刻吧． 
无聊的小 C 同学把这个问题稍微改了一下：要求找一条路径，使得每个点
最多经过一次，并且点权值之和最大，如果你跟小 C 一样无聊，就麻烦做一下
这个题目吧．

Input

第一行 m, n，接下来 m行每行 n 个数即． 
V( i,j)

Output

一个整数表示路径的最大权值之和．

Sample Input

2 3
1 -2 1
1 1 1

Sample Output

5
【数据范围】
30%的数据，n≤6．
100%的数据，m<=100，n ≤ ≤8．
注意：路径上有可能只有一个点．

恶心的杀头DP...
可以看cdq的<基于连通性状态压缩的动态规划问题>
PS:我的代码不知道比mcfx短到哪里去了...

#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#define maxp 2000000#define yjq 32767#define bit(x) ((x)<<1)using namespace std;int a[110][10];int nxt[maxp],head[yjq+10],tot[2],n,m,nn,mm;int status[2][maxp],ptr,nptr;int dp[2][maxp],ans=-1<<30;void add(int x,int s){    int p=x&yjq;    for(int t=head[p];t;t=nxt[t])        if(status[ptr][t]==x){            dp[ptr][t]=max(dp[ptr][t],s);            return ;        }    tot[ptr]++;    status[ptr][tot[ptr]]=x;    dp[ptr][tot[ptr]]=s;    nxt[tot[ptr]]=head[p],head[p]=tot[ptr];}int conn(int x,int j,int k){    if(k==1){        int b=1;        for(int i=j+1;i<=m;++i){            int p=(x>>bit(i))&3;            if(p==1)b++;            else if(p==2)b--;            if(b==0)return i;        }    } else if(k==2){        int b=1;        for(int i=j-2;i>=0;--i){            int p=(x>>bit(i))&3;            if(p==2)b++;            else if(p==1)b--;            if(b==0)return i;        }    }    return -1;}bool check0(int x){//for(int i=0;i<=m;++i)if(((x>>bit(i))&3)==1)return false;//return true;//return ()}int replace(int x,int j,int to){int p=(x>>bit(j))&3;return x-(p<<bit(j))+(to<<bit(j));}int main(){//freopen("in.txt","r",stdin);scanf("%d%d",&n,&m);for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)scanf("%d",&a[i][j]),ans=max(ans,a[i][j]);ptr=1,nptr=0,tot[ptr]=1;for(int i=1;i<=n;++i){for(int j=1;j<=tot[ptr];++j)status[ptr][j]<<=2;for(int j=1;j<=m;++j){swap(ptr,nptr);tot[ptr]=0;//printf("[%d]",tot[nptr]);memset(head,0,sizeof(head));for(int k=1,c1;k<=tot[nptr];++k){int st=status[nptr][k];if(st>=(1<<bit(m+1)))continue;int x=dp[nptr][k];int p=(st>>bit(j-1))&3,q=(st>>bit(j))&3;int xt=st-(p<<bit(j-1))-(q<<bit(j));if(p==0&&q==0){add(st,x);add(st+(3<<bit(j-1)),x+a[i][j]);add(st+(3<<bit(j)),x+a[i][j]);add(st+(1<<bit(j-1))+(2<<bit(j)),x+a[i][j]);} else if(p==0&&q>0){add(st,x+a[i][j]);add(st+(q<<bit(j-1))-(q<<bit(j)),x+a[i][j]);if(q!=3)add(replace(st-(q<<bit(j)),conn(st,j,q),3),x+a[i][j]);else if(!xt)ans=max(ans,x+a[i][j]);} else if(p>0&&q==0){add(st,x+a[i][j]);add(st-(p<<bit(j-1))+(p<<bit(j)),x+a[i][j]);if(p!=3)add(replace(st-(p<<bit(j-1)),conn(st,j,p),3),x+a[i][j]);else if(!xt)ans=max(ans,x+a[i][j]);} else if(p==q&&p!=3){ add(replace(st-(p<<bit(j-1))-(q<<bit(j)),conn(st,j,p),p),x+a[i][j]); } else if(p==2&&q==1){ add(st-(p<<bit(j-1))-(q<<bit(j)),x+a[i][j]); }else if(p==3&&q==3){ if(!xt)ans=max(ans,x+a[i][j]); } else if(p==3&&q!=3){ add(replace(st-(p<<bit(j-1))-(q<<bit(j)),conn(st,j,q),3),x+a[i][j]); }else if(p!=3&&q==3){ add(replace(st-(p<<bit(j-1))-(q<<bit(j)),conn(st,j,p),3),x+a[i][j]); }}}}    printf("%d",ans); }