HDU 1693 Eat the Trees 插头dp
题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=1693
Eat the Trees
Memory Limit: 32768/32768 K (Java/Others)
样例输出
Case 1: There are 3 ways to eat the trees.
Case 2: There are 2 ways to eat the trees.
题意
求不经过障碍物的哈密顿回路有多少种,可以独立的几个哈密顿回路。
题解
插头dp逐格转移。
对于每个格子:
1、上左都没有插头 --右下有插头
2、只有上方有插头 --上右有插头或上下有插头
3、只有左部有插头 --左下有插头或左右有插头
4、上左部都有插头 --上左有插头
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<ctime>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
using namespace std;
#define X first
#define Y second
#define mkp make_pair
#define lson (o<<1)
#define rson ((o<<1)|1)
#define mid (l+(r-l)/2)
#define sz() size()
#define pb(v) push_back(v)
#define all(o) (o).begin(),(o).end()
#define clr(a,v) memset(a,v,sizeof(a))
#define bug(a) cout<<#a<<" = "<<a<<endl
#define rep(i,a,b) for(int i=a;i<(b);i++)
#define scf scanf
#define prf printf
typedef long long LL;
typedef vector<int> VI;
typedef pair<int,int> PII;
typedef vector<pair<int,int> > VPII;
const int INF=0x3f3f3f3f;
const LL INFL=0x3f3f3f3f3f3f3f3fLL;
const double eps=1e-9;
const double PI = acos(-1.0);
//start----------------------------------------------------------------------
const int maxn=13;
LL dp[2][1<<maxn];
int n,m;
int pre,cur;
int main() {
int tc,kase=0;
scf("%d",&tc);
while(tc--) {
scf("%d%d",&n,&m);
///滚动数组
pre=0,cur=1;
clr(dp[cur],0);
dp[cur][0]=1;
///逐格转移
for(int i=0; i<n; i++) {
for(int j=0; j<m; j++) {
int x;
scf("%d",&x);
swap(pre,cur);
clr(dp[cur],0);
///p1为当前插头的左侧,p2为当前插头的上侧
int p1=(1<<j),p2=(1<<(j+1));
///根据具体情况吧左侧转移到下侧,上侧转移的右侧
for(int k0=0; k0<(1<<(m+1)); k0++) {
int k=k0;
///换行的时候需要位移一位,轮廓线画一画就会有感觉
if(!j) k=k0<<1;
///没有树的情况
if(!x) {
///必须都没有插头才合法
if(!(p1&k)&&!(p2&k)) dp[cur][k]+=dp[pre][k0];
continue;
}
///状压维护的是轮廓线上是否有插头的状态
if(!(p1&k)&&!(p2&k)) {
///右下侧有插头
if(i<n-1&&j<m-1) dp[cur][k^p1^p2]+=dp[pre][k0];
} else if(!(p1&k)&&(p2&k)) {
///上下侧有插头
if(i<n-1) dp[cur][k^p1^p2]+=dp[pre][k0];
///上右侧有插头
if(j<m-1) dp[cur][k]+=dp[pre][k0];
} else if((p1&k)&&!(p2&k)) {
///左右侧有插头
if(j<m-1) dp[cur][k^p1^p2]+=dp[pre][k0];
///左下侧有插头
if(i<n-1) dp[cur][k]+=dp[pre][k0];
} else {
///左上侧有插头
dp[cur][k^p1^p2]+=dp[pre][k0];
}
}
}
}
prf("Case %d: There are %lld ways to eat the trees.\n",++kase,dp[cur][0]);
}
return 0;
}
//end-----------------------------------------------------------------------