说在前面
啊…人生中第一道插头DP题
之前看远航休息栈的「从不会到崩溃」的插头dp学习笔记,me被吓懵了…一个DP的代码居然可以那么那么长???也就一直没敢写。今天A了某道状压dp题之后,丧心病狂的卡了一波常数妄想上榜,结果是失败的。于是闲的没事去看discuss,然后看见了插头dp…
最终还是逼迫自己把这道题敲了hhhhh。成就感++,知识漏洞- -
题目
URAL1519传送门
顺便放一个Vjudge的传送门Vjudge`s URAL - 1519
题目大意
给出一个N行M列的方格矩阵(N,M均不超过12),相邻格子之间有连边,而且矩阵中有一些方格是坏的。
现在询问不经过坏掉的格子的哈密顿回路的条数(保证答案不会超过long long范围)
输入输出格式
输入格式:
第一行两个整数N,M,含义如题
接下来是一个包含N*M个字符的矩阵
如果字符为「*」,则表示当前格子坏掉了;如果为「 . 」,则表示当前格子是好的
输出格式:
输出一行一个整数表示答案
解法
这题的解法..要是真的要详细的写,可能篇幅就太长了,并且me能力也有限…
不过还是安利几个资料吧=w=
陈丹琦(CDQ)的『基于连通性状态压缩的动态规划问题』,PPT做的很精美,看一遍可以大概了解插头DP是个什么东西
远航休息栈的「从不会到崩溃」,讲的很详细,不过代码略长
另外,虽然插头DP写起来真的很烦躁,但是请务必强迫自己写一写,想不通了就在纸上画一画,只要敲一遍就不会觉得很难了!!!
下面是自带大常数的代码
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;
int N , M , bit[15] ;
bool acce[15][15] ;
struct Points{
short i , j ;
}ed ;
long long final_ans , ans[2][300005] ;
int head[10007] , tp[2] , preE[300005] , S[2][300005] , now , pre ;
void Insert( int s , long long num ){
int id = s%10007 ;
for( int i = head[id] ; i ; i = preE[i] )
if( S[now][i] == s ){
ans[now][i] += num ;
return ;
}
preE[ ++tp[now] ] = head[id] ;
head[id] = tp[now] ;
S[now][ tp[now] ] = s ;
ans[now][ tp[now] ] = num ;
}
void solve(){
now = 0 , pre = 1 ;
Insert( 0 , 1 ) ;
//state 计数规则 : 高位靠右, 低位靠左
for( int i = 1 ; i <= N ; i ++ ){
for( int k = 1 ; k <= tp[pre] ; k ++ )
S[now][k] <<= 2 ;//边上的插头没有了...
for( int j = 1 ; j <= M ; j ++ ){
swap( now , pre ) ;
tp[now] = 0 ;
memset( head , 0 , sizeof( head ) ) ;
for( int k = 1 ; k <= tp[pre] ; k ++ ){
int state = S[pre][k] , left = ( state >> bit[j-1] )&3 , up = ( state >> bit[j] )&3 ;
long long tmp = ans[pre][k] ;
// ban
if( !acce[i][j] ){
if( !left && !up ) Insert( state , tmp ) ;
continue ;
}
// accessable
if( !left && !up ){ // 无 , 新建插头
if( acce[i][j+1] && acce[i+1][j] )
Insert( state | ( 1 << bit[j-1] ) | ( 2 << bit[j] ) , tmp ) ;
}
else if ( left && up ) { // 有两个 , 消除插头
state ^= ( left << bit[j-1] ) ^ ( up << bit[j] ) ; //先消除
if( left == 1 && up == 2 ) // ( ) -> # #
final_ans += ( i == ed.i && j == ed.j ) * tmp ;
else if( left == 2 && up == 1 ) // ( ) ( ) -> ( # # )
Insert( state , tmp ) ;
else if( left == 1 && up == 1 ) // ( ( ) ) -> # # ( )
for( int t = j + 1 , count = 1 ; t <= M ; t ++ ){
if( ( state >> bit[t] )&2 ) count -- ;
if( ( state >> bit[t] )&1 ) count ++ ;
if( !count ){
Insert( state ^ ( 3 << bit[t] ) , tmp ) ;
break ;
}
}
else if( left == 2 && up == 2 ) // ( ( ) ) -> ( ) # #
for( int t = j - 2 , count = -1 ; t >= 0 ; t -- ){
if( ( state >> bit[t] )&2 ) count -- ;
if( ( state >> bit[t] )&1 ) count ++ ;
if( !count ){
Insert( state ^ ( 3 << bit[t] ) , tmp ) ;
break ;
}
}
}
else if( left ){ // 延续左插头
if( acce[i+1][j] ) Insert( state , tmp ) ;
if( acce[i][j+1] ) Insert( state ^ ( left << bit[j-1] ) ^ ( left << bit[j] ) , tmp ) ;
}
else if( up ){ // 延续上插头
if( acce[i+1][j] ) Insert( state ^ ( up << bit[j-1] ) ^ ( up << bit[j] ) , tmp ) ;
if( acce[i][j+1] ) Insert( state , tmp ) ;
}
}
}
}
printf( "%lld" , final_ans ) ;
}
int main(){
scanf( "%d%d" , &N , &M ) ;
for( short i = 1 ; i <= N ; i ++ )
for( short j = 1 ; j <= M ; j ++ ){
char ch = getchar() ;
while( ch != '.' && ch != '*' ) ch = getchar() ;
if( ch == '.' ) acce[i][j] = true , ed = (Points){ i , j } ;
}
for( int i = 1 ; i <= 12 ; i ++ )
bit[i] = i * 2 ;
solve() ;
}