[URAL1519]-Formula 1-插头DP

说在前面

啊…人生中第一道插头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() ;
}