Ural 1519. Formula 1 优美的插头DP

今天早上学了插头DP的思想和最基础的应用，中午就开始敲了，岐哥说第一次写不要看别人代码，利用自己的理解一点点得写出来，这样才锻炼代码能力！于是下午慢慢地构思轮廓，一点点地敲出主体代码，其实是很磨蹭的，由于要考虑好多东西，而且昨晚2点睡的有点困，最后终于磨蹭出来了，第一次的代码搓没关系，自己写的才重要。然后果然不出我所料，调试到了晚上才A了（一个郁闷的错误）。。。A的感觉真的是爽呀，虽然搞了差不多一天。当然自己写了自己想的代码后也要把代码优化，不然队友看不懂自己代码就囧了。。。

插头DP，建议大家想学的好好看看陈丹琦的国家集训队论文，这是个优美的DP。

http://www.docin.com/p-46797997.html

#include <stdio.h>
#include <string.h>

#define LL __int64

const int mod = 10007;

//   哈希表
struct HASH{
    int head[mod+10], E, next[80000];
    LL val[80000], cnt[80000];

    void init() {
        memset(head, -1, sizeof(head));
        E = 0;
    }

    int findhash(LL x) {
        return (x%mod + mod)%mod;
    }

    void add(LL x, LL sum) {
        int u = findhash(x);
        for(int i = head[u];i != -1;i = next[i]) if(val[i] == x) {
            cnt[i] += sum;
            return ;
        }
        val[E] = x;
        cnt[E] = sum;
        next[E] = head[u];
        head[u] = E++;
    }

}biao1, biao2;

int c[22], n, m, d[22];
//  编码
void get(LL x) {
    for(int i = m+1;i >= 1; i--) {
        c[i] = x&7;
        x /= 8;
    }
}
//  解码
LL getval() {
    LL ret = 0;
    for(int i = 1;i <= m+1; i++) {
        ret |= d[i];
        ret *= 8;
    }
    ret /= 8;
    return ret ;
}
//   转化成最小表示法
void change() {
    int vis[22];
    memset(vis, 0, sizeof(vis));
    int num = 1;
    for(int i = 1;i <= m+1;i ++) {
        if(!d[i])   continue;
        if(!vis[d[i]]) {
            vis[d[i]] = num;
            d[i] = num++;
        }
        else {
            d[i] = vis[d[i]];
        }
    }
}

void fuzhi() {
    for(int i = 1;i <= m+1;i ++)    d[i] = c[i];
}

char s[22][22];

int main() {
    int i, j, k, l;
    while(scanf("%d%d", &n, &m) != -1) {
        for(i = 1;i <= n; i++)
            scanf("%s", s[i]+1);
        int tot = 0;
        for(i = 1;i <= n; i++)
            for(j = 1;j <= m; j++)
            if(s[i][j] == '.')  tot++;
        if(tot%2==1 || tot < 4) {
            puts("0");
            continue;
        }
        int tox  = -1, toy = -1;
        for(i = 1;i <= n; i++)
        for(j = 1;j <= m; j++) if(s[i][j] == '.') {
            tox = i;
            toy = j;
        }
        biao1.init();
        biao1.add(0, 1);
        LL ans = 0;
        for(i = 1;i <= n; i++){
            for(j = 0;j <= m; j++){
                biao2.init();
                for(l = 0;l < biao1.E; l++)  {
                    get(biao1.val[l]);
                    if(j == m) {
                        for(int ii = 2;ii <= m+1; ii++) d[ii] = c[ii-1];
                        d[1] = 0;
                        change();
                        LL now = getval();
                        biao2.add(now, biao1.cnt[l]);
                        continue;
                    }
                    if(c[j+1] && !c[j+2]) {  //  有左插头无上插头
                        if(s[i][j+1] != '.')  continue;
                        if(j+2 <= m) {
                            fuzhi();
                            d[j+1] = 0;d[j+2] = c[j+1];
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                        if(i < n) {
                            fuzhi();
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                    }
                    else if(!c[j+1] && c[j+2]) {  //  有上插头无左插头
                        if(s[i][j+1] != '.')  continue;
                        if(i < n) {
                            fuzhi();
                            d[j+1] = c[j+2]; d[j+2] = 0;
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                        if(j+2 <= m) {
                            fuzhi();
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                    }
                    else if(!c[j+1] && !c[j+2]) { //   左和上都无插头
                        if(s[i][j+1] != '.') {
                            fuzhi();
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                            continue;
                        }
                        if(j+2 <= m && i < n) {
                            fuzhi();
                            d[j+1] = d[j+2] = 13;
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                    }
                    else { //　　左和上都有插头　，　　要判断左和上插头是否连通
                        if(c[j+2] == c[j+1]) {
                            int tot = 0;
                            for(int ii = 1;ii <= m+1; ii++) if(c[ii])
                                tot++;
                            if(tot == 2 && i == tox && j+1 == toy) ans += biao1.cnt[l];
                        }
                        else {
                            if(s[i][j+1] != '.')    continue;
                            fuzhi();
                            for(int ii = 1;ii <= m+1; ii++) if(ii != j+1 && ii != j+2 && d[ii] == d[j+1]) {
                                d[ii] = d[j+2];
                                break;
                            }
                            d[j+1] = d[j+2] = 0;
                            change();
                            LL now = getval();
                            biao2.add(now, biao1.cnt[l]);
                        }
                    }
                }
                biao1 = biao2;
            }
        }
        printf("%I64d\n", ans);
    }
    return 0;
}