PKU1038 Bugs Integrated, Inc. 状压DP

周伟的论文没能看懂感觉挺羞愧的然而现在我也没能理解。然后滚去看黑书，P138

首先我们先没枚举上一行的状态，dfs这一行的状态，用三进制来表示从上向下有连续的多少个0，然后我们就可以愉快地转移了。

虽然仅仅用时1266ms，但是我并没有明白为什么会这么快，原本以为会T的很惨的。

显然我们先枚举行，然后枚举上一行状态。此时极限数据会达到150*50000+*50000+，然后又5组数据……根本看不到希望

不过实际上会因为这一行很多不能放的和之前的限制把后面的50000+降低了非常多= =，感觉放到UOJ上绝对被hack狂魔折磨的不成人形

周伟的做法dfs同时枚举上一次的状态和这一次的状态，跑的会快很多，然而还没看懂= =

Problem: 1038 User: BPM136
Memory: 1160K Time: 1266MS
Language: G++ Result: Accepted

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#define LL long long
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define down(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
inline LL read()
{
LL d=0,f=1;char s=getchar();
while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
return d*f;
}
#define N 155
#define M 15
#define K 60000
int f[K][2];
bool map
[M];
int b[M],s[M];
int n,m,nm;
int x,y;

void init()
{
memset(map,0,sizeof(map));
memset(f,-1,sizeof(f));
}

void chu()
{
b[0]=1;
fo(i,1,10)b[i]=b[i-1]*3;
}

void dfs(int xu,int t,int k,int cnt,int num)
{
if(k>m)
{
f[cnt][xu&1]=max(f[cnt][xu&1],f[t][(xu-1)&1]+num);
return;
}
int tt=0;
if(map[xu][k])tt=2;
else if(s[k])tt=s[k]-1;else tt=0;
dfs(xu,t,k+1,cnt+tt*b[k-1],num);

if(k+1<=m&&s[k]==0&&s[k+1]==0&&!map[xu][k]&&!map[xu][k+1])
dfs(xu,t,k+2,cnt+2*b[k-1]+2*b[k],num+1);

if(k+2<=m&&s[k]<=1&&s[k+1]<=1&&s[k+2]<=1)
if(!map[xu][k]&&!map[xu][k+1]&&!map[xu][k+2])
dfs(xu,t,k+3,cnt+2*b[k-1]+2*b[k]+2*b[k+1],num+1);
}

int main()
{
chu();
int ttt=read();
while(ttt--)
{
init();
n=read(),m=read(),nm=read();
fo(i,1,nm)
{
x=read(),y=read();
map[x][y]=1;
}
f[b[m]-1][0]=0;
fo(i,1,n)
{
fo(j,0,b[m]-1)f[j][i&1]=-1;
fo(j,0,b[m]-1)
if(f[j][(i-1)&1]!=-1)
{
s[0]=1;
for(int t=j,k=1;k<=m;t/=3,k++)s[k]=t%3;
dfs(i,j,1,0,0);
}
}
int ans=-1;
fo(i,0,b[m]-1)ans=max(ans,f[i][n&1]);
printf("%d\n",ans);
}
return 0;
}