URAL1519：Formula 1——题解

http://acm.timus.ru/problem.aspx?space=1&num=1519

https://vjudge.net/problem/URAL-1519

题目大意：给一个网格，有些网格有障碍，问有多少条哈密顿回路。

——————————————————-

插头dp模板题，详见：http://www.cnblogs.com/luyouqi233/p/8256778.html

参考代码和思路于：http://www.yhzq-blog.cc/%E6%8F%92%E5%A4%B4dp-%E4%BB%8E%E4%B8%8D%E4%BC%9A%E5%88%B0%E5%B4%A9%E6%BA%83/

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int INF=;
const int mod=;
const int M=;
struct node{
    int to,nxt;
}edge[M];
int head[M],cnt;
int n,m,e1,e2;
bool mp[][]; //记录是否有障碍
int cur,pre;  //滚动数组
int state[][M];  //记录状态,滚动
ll ans[][M],cntt; //记录答案
int tot[]; //记录状态总数
int bit[];  //提取状态使用
inline void getbit(){
    for(int i=;i<;i++)bit[i]=i<<;
    return;
}
inline void add(int u,int v){
    cnt++;
    edge[cnt].to=v;
    edge[cnt].nxt=head[u];
    head[u]=cnt;
    return;
}
void insert(int now,ll num){
    int u=now%mod;
    for(int i=head[u];i;i=edge[i].nxt){
    int v=edge[i].to;
    if(state[cur][v]==now){
        ans[cur][v]+=num;
        return;
    }
    }
    add(u,++tot[cur]);
    state[cur][tot[cur]]=now;
    ans[cur][tot[cur]]=num;
    return;
}
void plugdp(){
    cur=;
    tot[cur]=;
    ans[cur][]=;
    state[cur][]=;
    for(int i=;i<=n;i++){
    for(int j=;j<=tot[cur];j++){
        state[cur][j]<<=;
    }
    for(int j=;j<=m;j++){
        memset(head,,sizeof(head));cnt=;
        pre=cur,cur^=;
        tot[cur]=;
        for(int k=;k<=tot[pre];k++){
        int now=state[pre][k];
        ll num=ans[pre][k];
        int is_down=(now>>bit[j-])%;
        int is_right=(now>>bit[j])%;
        if(!mp[i][j]){
            if(!is_down&&!is_right)
            insert(now,num);
        }
        else if(!is_down&&!is_right){
            if(mp[i][j+]&&mp[i+][j])
            insert(now+(<<bit[j-])+*(<<bit[j]),num);
        }
        else if(!is_down&&is_right){
            if(mp[i][j+])
            insert(now,num);
            if(mp[i+][j])
            insert(now-is_right*(<<bit[j])+is_right*(<<bit[j-]),num);
        }
        else if(is_down&&!is_right){
            if(mp[i+][j])
            insert(now,num);
            if(mp[i][j+])
            insert(now-is_down*(<<bit[j-])+is_down*(<<bit[j]),num);
        }
        else if(is_down==&&is_right==){
            int count=;
            for(int l=j+;l<=m;l++){
            if((now>>bit[l])%==)count++;
            if((now>>bit[l])%==)count--;
            if(!count){
                insert(now-(<<bit[l])-(<<bit[j-])-(<<bit[j]),num);
                break;
            }
            }
        }
        else if(is_down==&&is_right==){
            int count=;
            for (int l=j-;l>=;l--){
            if((now>>bit[l])%==) count--;
            if((now>>bit[l])%==) count++;
            if(!count){
                insert(now+(<<bit[l])-*(<<bit[j-])-*(<<bit[j]),num);
                break;
            }
            }
        }
        else if(is_down==&&is_right==)
            insert(now-*(<<bit[j-])-(<<bit[j]),num);
        else if(is_down==&&is_right==){
            if(i==e1&&j==e2)
            cntt+=num;
        }
        }
    }
    }
    return;
}
int main(){
    getbit();
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++){
    for(int j=;j<=m;j++){
        char ch=getchar();
        while(ch!='*'&&ch!='.')ch=getchar();
        if(ch=='.')mp[i][j]=,e1=i,e2=j;
    }
    }
    plugdp();
    printf("%lld\n",cntt);
    return ;
}