Light OJ 1402 Assassin`s Creed (状压DP)

参考了Alex » LightOJ1406 – Assassin`s Creed(状压dp)并做了一点简单优化。

通过DFS预处理出所有的可达状态，然后状压，值得注意的是枚举子集的复杂度在n=15是近似1.5e7，所以不会超时。

[code]:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<vector>

#define pb push_back
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 1e5+5;

int n,m,dp[1<<15];
bool link[1<<15],vis[1<<15][15];
vector<int> G[20];

void dfs(int S,int u){
    link[S] = 1;
    int i,v;
    for(i = 0;i < G[u].size();i++){
        v = G[u][i];
        if(vis[S|(1<<v)][v]) continue;
        vis[S|(1<<v)][v] = 1;
        dfs(S|(1<<v),v);
    }
}

int main(){
    int i,j,cas,T;
    scanf("%d",&cas);
    for(T = 1;T <= cas;T++){
        scanf("%d%d",&n,&m);
        for(i = 0;i < n;i++) G[i].clear();
        int u,v;
        while(m--){
            scanf("%d%d",&u,&v);u--,v--;
            G[u].pb(v);
        }

        memset(link,0,sizeof(link));
        memset(dp,63,sizeof(dp));
        link[0] = 1;dp[0] = 0;

        memset(vis,0,sizeof(vis));
        int S,S_,Ed = 1<<n;
        for(i = 0;i < n;i++){
            vis[1<<i][i] = 1;
            dfs(1<<i,i);
        }

        for(S = 1;S < Ed;S++){
            for(S_=S;S_;S_=(S_-1)&S){
                if(link[S_]) dp[S] = min(dp[S],dp[S^S_]+1);
            }
        }
        printf("Case %d: %d\n",T,dp[Ed-1]);
    }

    return 0;
}