POJ 1038 - Bugs Integrated, Inc. 三进制状态DP

POJ 1038 - Bugs Integrated, Inc. 三进制状态DP

每个点有三个状态...不放家伙..放横的.放竖的..虽然看上去状态有3^m..最多3^10=59049种...但把自我矛盾的排除后..一行的可能状态最多280种...

由于一个方块最多可能影响到上面两层..所以DP需要三维...dp [ r ] [ x ]  [ y ] 代表第r-1行放状态y...第r行放状态x..

第一次提交爆空间了...把dp的第一维改为滚动的就可以了....

效率不高...3000MS+过的...有些大神不到1000MS过..有些也是三进制状态DP..10000MS~20000MS过...看来状态DP还是写得比较搓 ~~~

Program:

#include#include#include#include#include#include#define oo 1000000007#define ll long long#define pi acos(-1.0)#define MAXN 505using namespace std; int s[155][15],n,m,canuse[280],w[280],tall[280],dp[2][280][280],num;bool f[280][280][2];bool legal(int x){ int s[15],i,j; memset(s,0,sizeof(s)); w[num+1]=0; tall[num+1]=0; for (i=1;i<=m;i++) { if (x%3==1) s[i]++,s[i+1]++,s[i+2]++,tall[num+1]=max(tall[num+1],1); if (x%3==2) s[i]++,s[i+1]++,tall[num+1]=max(tall[num+1],2); if (x%3) w[num+1]++; x/=3; } if (s[m+1]) return false; for (i=1;i<=m;i++) if (s[i]>1) return false; return true; }bool canput(int x,int r){ int i; if (r-tall[x]<1) return false; x=canuse[x]; for (i=1;i<=m;i++) { if (x%3==1 && (s[r][i] || s[r][i+1] || s[r][i+2] || s[r-1][i] || s[r-1][i+1] || s[r-1][i+2])) return false; if (x%3==2 && (s[r][i] || s[r][i+1] || s[r-1][i] || s[r-1][i+1] || s[r-2][i] || s[r-2][i+1])) return false; x/=3; } return true;}bool ok(int x,int y,int tp){ int a[10][15],i,j; memset(a,0,sizeof(a)); x=canuse[x],y=canuse[y]; for (i=1;i<=m;i++) { if (y%3==1) a[0][i]++,a[0][i+1]++,a[0][i+2]++; if (y%3==2) a[0][i]++,a[0][i+1]++; if (x%3==1) a[tp-1][i]++,a[tp-1][i+1]++,a[tp-1][i+2]++; if (x%3==2) a[tp-1][i]++,a[tp-1][i+1]++,a[tp-2][i]++,a[tp-2][i+1]++; y/=3,x/=3; } for (i=0;i<10;i++) for (j=0;j<15;j++) if (a[i][j]>1) return false; return true;}int main(){ int T,K,totol,r,x,y,i,j,ans; scanf("%d",&T); while (T--) { scanf("%d%d%d",&n,&m,&K); memset(s,0,sizeof(s)); while (K--) scanf("%d%d",&y,&x),s[y][x]=1; num=0; totol=1; for (x=1;x<=m;x++) totol*=3; for (x=0;x

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。