hdu_3247_Resource Archiver(AC自动机+bfs+TSP)

题目链接：hdu_3247_Resource Archiver

题意：

有n个资源串，m个病毒串，现在要将所有的资源串整合到一个串内，并且这个串不能包括病毒串，问最短的串长为多少

题解：

将资源串和病毒串都插入到AC自动机中，分别做好标记，然后用bfs求出0节点和所有资源串互相的最短距离，最后就是一个TSP的状态压缩DP。

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;i++)
using namespace std;

const int AC_N=1e5+7,tyn=2;
struct AC_automation{
    int tr[AC_N][tyn],so[AC_N],v[AC_N],Q[AC_N],fail[AC_N],tot;
    inline int getid(char x){return x-'0';}
    void nw(){so[++tot]=0,fail[tot]=0,v[tot]=0;memset(tr[tot],0,sizeof(tr[tot]));}
    void init(){tot=-1,fail[0]=-1,nw();}
    void insert(char *s,int idx,int fg,int x=0){
        for(int len=strlen(s),i=0,w;i<len;x=tr[x][w],i++)
            if(!tr[x][w=getid(s[i])])nw(),tr[x][w]=tot;
        if(fg)so[x]=1<<idx;else v[x]=1;
    }
    void build(int head=1,int tail=0){
        for(int i=0;i<tyn;i++)if(tr[0][i])Q[++tail]=tr[0][i];
        while(head<=tail)for(int x=Q[head++],i=0;i<tyn;i++)
            if(tr[x][i])fail[tr[x][i]]=tr[fail[x]][i],Q[++tail]=tr[x][i],
            so[tr[x][i]]|=so[tr[fail[x]][i]],v[tr[x][i]]+=v[tr[fail[x]][i]];
            else tr[x][i]=tr[fail[x]][i];
    }
    
}AC;

char s[AC_N];
int path[205][205],pos[205],dp[1<<10][205],cnt,Q[AC_N],dis[AC_N],n,m;

inline void up(int &a,int b){if(a>b)a=b;}

void bfs(int u,int head=1,int tail=0)
{
    memset(dis,-1,sizeof(dis));
    Q[++tail]=pos[u],dis[pos[u]]=0;
    while(head<=tail)
    {
        int now=Q[head++];
        F(i,0,1)
        {
            int q=AC.tr[now][i];
            if(dis[q]<0&&!AC.v[now])dis[q]=dis[now]+1,Q[++tail]=q;
        }
    }
    F(i,0,cnt)path[u][i]=dis[pos[i]];
}

void fuck()//TSP的dp
{
    memset(dp,-1,sizeof(dp)),dp[0][0]=0;
    for(int i=0;i<(1<<n);i++)F(j,0,cnt)if(dp[i][j]>=0)F(k,0,cnt)if(path[j][k]>=0)
    if(dp[i|AC.so[pos[k]]][k]==-1)dp[i|AC.so[pos[k]]][k]=dp[i][j]+path[j][k];
    else up(dp[i|AC.so[pos[k]]][k],dp[i][j]+path[j][k]);
    int ans=INT_MAX;
    F(i,0,cnt)if(dp[(1<<n)-1][i]!=-1)up(ans,dp[(1<<n)-1][i]);
    printf("%d
",ans==INT_MAX?-1:ans);
}

int main()
{
    while(~scanf("%d%d",&n,&m)&&n+m)
    {
        AC.init(),cnt=0;
        F(i,0,n-1)scanf("%s",s),AC.insert(s,i,1);
        F(i,0,m-1)scanf("%s",s),AC.insert(s,i,0);
        AC.build();
        F(i,1,AC.tot)if(AC.so[i])pos[++cnt]=i;
        F(i,0,cnt)bfs(i);//求出0到所以资源串的距离和资源串到资源串之间的距离
        fuck();
    }
    return 0;
}