fzu1977之插头DP
Problem 1977 Pandora adventureAccept: 354 Submit: 1177
Time Limit: 1000 mSec Memory Limit : 32768 KB
Problem Description
The pollution of the earth is so serious that people can not survive any more. Fortunately, people have found a new planet that maybe has life, and we call it "Pandora Planet".
Leonardo Da Vinci is the only astronaut on the earth. He will be sent to the Pandora Planet to gather some plant specimens and go back. The plant specimen is important to the people to decide whether the planet is fit to live or not.
Assuming that Da Vinci can only move in an N×M grid. The positions of the plant specimens he wants to collect are all marked by the satellite. His task is to find a path to collect all the plant specimens and return to the spaceship. There are some savage beasts in the planet. Da Vinci can not investigate the grid with the savage beast. These grids are also marked by the satellite. In order to save time Da Vinci could only visit each grid exactly once and also return to the start grid, that is, you can not visit a grid twice except the start grid. You should note that you can choose any grid as the start grid.
Now he wants to know the number of different paths he can collect all the plant specimens. We only care about the path and ignore where the start grid is, so the two paths in Figure 1 are considered as the same.
Input
The first line of the input contains an integer T (T≤100), indicating the number of cases. Each case begins with a line containing two integers N and M (1≤N, M≤12), the size of the planet is N×M. Each of the following N lines contains M characters Gij(1≤i≤N, 1≤j≤M), Gij denotes the status of the grid in row i and column j, where 'X' denotes the grid with savage beast, '*' denotes the safe grid that you can decide to go or not, 'O' denotes the plant specimen you should collect. We guarantee that there are at least three plant specimens in the map.
Output
For each test case, print a line containing the test case number (beginning with 1) and the number of different paths he can collect all the plant specimens. You can make sure that the answer will fit in a 64-bit signed integer.
Sample Input
Sample Output
分析：因为须要求回路个数，所以能想到是插头DP,本题在最简单的求回路个数上添加一些必走的点和可走的点
对于可走的点仅仅要在插头DP的时候假设p=q=0且该点能够不走，则就添加一种不走的状态：
if(!p && !q){
if(mp[i][j] == '*')HashCalState(s,num);
if(mp[i][j+1] == 'X' || mp[i+1][j] == 'X')continue;
s=s+(1<<bit[j-1])+2*(1<<bit[j]);
HashCalState(s,num);
}而问题在于假设求经过全部'O'的回路？
如果求经过全部可经过点的回路时碰到p=1,q=2合并连通块
则发生这样的情况必须是最后一个可走的点。这样全部可走的点才会连通
在这里也是一样，仅仅要记录最后一个必走的点'O就可以
在碰到p=1,q=2合并连通块的时候。这时候仅仅需推断这个点是不是最后一个“O'点而且能够连成环
推断连成环既已决策的点没有插头连着未决策的点，这样全部已决策的点连着必成环
if(p == 1 && q == 2){//成环
if(i == ex && j<ey)continue;
if(i<ex)continue;
s=s-(1<<bit[j-1])-2*(1<<bit[j]);//求解其它已决策的点十是否和未决策的点相连
if(!s)sum+=num;
}其它就是插头DP模版了
採用邻接表寻找状态时head数组能够开小点，head的大小和效率关系非常大
另外本题好坑爹！
。！用long long WA,用__int64 AC,平白无故查错一晚上+一中午
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <queue>
#include <algorithm>
#include <map>
#include <cmath>
#include <iomanip>
#define INF 99999999
typedef __int64 LL;
using namespace std;
const int MAX=30000+10;
const int maxn=50000+10;
const int N=12+10;
int n,m,index,size[2],bit[N];
int ex,ey;
int head[MAX],next[maxn];
LL dp[2][maxn],state[2][maxn],sum;
char mp[N][N];
void HashCalState(LL s,LL num){
int pos=s%MAX;
for(int i=head[pos];i != -1;i=next[i]){
if(state[index][i] == s){
dp[index][i]+=num;
return;
}
}
state[index][size[index]]=s;
dp[index][size[index]]=num;
//头插法
next[size[index]]=head[pos];
head[pos]=size[index]++;
}
void DP(){
sum=0;
index=0;
size[index]=1;
state[index][0]=0;
dp[index][0]=1;
for(int i=1;i<=n;++i){
for(int k=0;k<size[index];++k)state[index][k]<<=2;
for(int j=1;j<=m;++j){
index=index^1;
memset(head,-1,sizeof head);
size[index]=0;
for(int k=0;k<size[index^1];++k){
LL s=state[index^1][k];
LL num=dp[index^1][k];
int p=(s>>bit[j-1])%4;
int q=(s>>bit[j])%4;
if(mp[i][j] == 'X'){//须要绕过
if(!p && !q)HashCalState(s,num);
}else if(!p && !q){
if(mp[i][j] == '*')HashCalState(s,num);
if(mp[i][j+1] == 'X' || mp[i+1][j] == 'X')continue;
s=s+(1<<bit[j-1])+2*(1<<bit[j]);
HashCalState(s,num);
}else if(!p && q){
if(mp[i][j+1] != 'X')HashCalState(s,num);
if(mp[i+1][j] != 'X'){
s=s+q*(1<<bit[j-1])-q*(1<<bit[j]);
HashCalState(s,num);
}
}else if(p && !q){
if(mp[i+1][j] != 'X')HashCalState(s,num);
if(mp[i][j+1] != 'X'){
s=s-p*(1<<bit[j-1])+p*(1<<bit[j]);
HashCalState(s,num);
}
}else if(p == 1 && q == 1){
int b=1;
for(int t=j+1;t<=m;++t){
int v=(s>>bit[t])%4;
if(v == 1)++b;
if(v == 2)--b;
if(!b){
s=s+(1<<bit[t])-2*(1<<bit[t]);
break;
}
}
s=s-(1<<bit[j-1])-(1<<bit[j]);
HashCalState(s,num);
}else if(p == 2 && q == 2){
int b=1;
for(int t=j-2;t>=0;--t){
int v=(s>>bit[t])%4;
if(v == 2)++b;
if(v == 1)--b;
if(!b){
s=s-(1<<bit[t])+2*(1<<bit[t]);
break;
}
}
s=s-2*(1<<bit[j-1])-2*(1<<bit[j]);
HashCalState(s,num);
}else if(p == 1 && q == 2){//成环
if(i == ex && j<ey)continue;
if(i<ex)continue;
s=s-(1<<bit[j-1])-2*(1<<bit[j]);
if(!s)sum+=num;
}else if(p == 2 && q == 1){
s=s-2*(1<<bit[j-1])-(1<<bit[j]);
HashCalState(s,num);
}
}
}
}
}
int main(){
for(int i=0;i<N;++i)bit[i]=i<<1;
int t,num=0;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i)scanf("%s",mp[i]+1);
ex=ey=-1;
for(int i=1;i<=n;++i){
for(int j=1;j<=m;++j)if(mp[i][j] == 'O')ex=i,ey=j;
mp[i][m+1]='X';
}
for(int j=1;j<=m;++j)mp[n+1][j]='X';
DP();
if(ex == -1)++sum;
printf("Case %d: %I64d\n",++num,sum);
}
return 0;
}
/*
100
12 12
************
************
************
************
************
************
************
************
************
************
************
************
*/
