ZOJ3466-The Hive II

题意

有一个六边形格子，共 \(n\) 行，每行有 8 个位置，有一些格子不能走。求用一些环覆盖所有可走格子的方案数。\(n\le 10\) 。

分析

插头dp，只不过是六边形上的，分奇数列和偶数列讨论，转移不太一样，但大同小异。编号需要插空，注意以下就行了。复杂度为 \(O(nm2^{2m})\) 。

代码

cc_hash_table 多次用实力证明，他跑得比其他任何内置的 hash_map 都快！

#include<bits/stdc++.h>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
using namespace std;
using namespace __gnu_pbds;
typedef long long giant;
typedef cc_hash_table<int,giant> Map;
typedef Map::iterator itt;
const int maxn=11;
const int m=8,maxm=m+1;
int n,q;
bool no[maxn][maxm];
Map f,g;
inline int get(int x,int p) {return (x>>p)&1;}
inline int mod(int x,int p,int d) {return (x&(~(1<<p)))+(d<<p);}
inline int mor(int x,int p,int d) {
    int bef=x&((1<<p)-1);
    return ((x-bef)<<d)+bef;
}
inline int mol(int x,int p,int d) { // cover the pth
    int bef=x&((1<<p)-1);
    return (((x-bef)>>(p+2))<<p)+bef;
}
int work() {
    memset(no,0,sizeof no);
    while (q--) {
        static char s[5];
        scanf("%s",s);
        no[s[0]-'A'+1][s[1]-'A'+1]=true;
    }
    f.clear(),g.clear();
    f[0]=1;
    for (int i=1;i<=n;++i) {
        f.swap(g),f.clear();
        for (itt it=g.begin();it!=g.end();++it) {
            const int &d=it->first;
            const giant &w=it->second;
            if (get(d,1)==0 && get(d,17)==0) f[d]=w;
        }
        for (int j=1;j<=m;++j) {
            f.swap(g),f.clear();
            for (itt it=g.begin();it!=g.end();++it) {
                const int &d=it->first;
                const giant &w=it->second;
                if (j&1) {
                    int a=j*2-1,t=mor(d,a+2,2),x=get(d,a),y=get(d,a+1),sum=x+y;
                    if (no[i][j]) {
                        if (!sum) f[t]+=w;
                        continue;
                    }
                    if (!sum) {
                        for (int k=0;k<4;++k) for (int r=k+1;r<4;++r) f[mod(mod(t,a+k,1),a+r,1)]+=w;
                    } else {
                        int v=mod(mod(t,a,0),a+1,0);
                        if (sum==1) for (int k=0;k<4;++k) f[mod(v,a+k,1)]+=w; else
                        if (sum==2) f[v]+=w;
                    }
                } else {
                    int a=j<<1,t=mol(d,a+2,2),u=get(d,a),x=get(d,a+1),y=get(d,a+2),z=get(d,a+3),sum=u+x+y+z;
                    if (no[i][j]) {
                        if (!sum) f[t]+=w;
                        continue;
                    }
                    if (!sum) f[mod(mod(t,a,1),a+1,1)]+=w; else {
                        int v=mod(mod(t,a,0),a+1,0);
                        if (sum==1) for (int k=0;k<2;++k) f[mod(v,a+k,1)]+=w; else
                        if (sum==2) f[v]+=w;
                    }
                }
            }
        }
    }
    printf("%lld\n",f[0]);
}
int main() {
#ifndef ONLINE_JUDGE
    freopen("test.in","r",stdin);
#endif
    while (~scanf("%d%d",&n,&q)) work();
    return 0;
}

转载于:https://www.cnblogs.com/owenyu/p/7511317.html