HDU 3247 Resource Archiver(AC自动机+BFS+状态DP)

转载请注明出处，谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove

题目：给出n个资源，m个病毒，将资源串拼接成一个串，必须包含所有的资源串，可以重叠，但是不能包含病毒

问最小的长度为多少

http://acm.hdu.edu.cn/showproblem.php?pid=3247 

不会做，看题解啊看题解，看了也不会。。。。

将所有的资源串和病毒串都放在Trie树里建立起来，当然作上相应的标记。然后建立fail指针之后

从资源串的结尾出发，BFS，记录能到达其它资源串结尾的步数。得到所有资源串结尾状态的距离邻接阵。

之后是状态压缩DP

dp[i][j]表示资源串的状态为i时，当前处在Trie上的状态j，最短长度。

感觉静态的写得怪怪的，而且不直观。。

#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>#define eps 1e-10#define N 100005#define inf 1<<20#define zero(a) (fabs(a)<eps)#define lson (step<<1)#define rson (step<<1|1)using namespace std;struct Trie{    int next[2],fail;    int virus,source;    void Init(){next[0]=next[1]=0;fail=source=virus=0;}}tree[N];int n,m;int tot,cnt;int path[205][205];int pos[205];int dp[1<<10][205];//dp[i][j]表示状态为i时，以j结尾的最少步数char str[50005];void Insert(char *s,int len,int kind){    int p=0;    for(int i=0;i<len;i++){        if(tree[p].next[s[i]-'0']==0){            tree[++tot].Init();            tree[p].next[s[i]-'0']=tot;        }        p=tree[p].next[s[i]-'0'];    }    if(kind>=0) tree[p].source=1<<kind;    else tree[p].virus=1;}void Bulid_fail(){    int que[N];    int head=0,tail=0;    que[tail++]=0;    while(head<tail){        int p=que[head++];        for(int i=0;i<2;i++){            if(tree[p].next[i]==0) tree[p].next[i]=tree[tree[p].fail].next[i];            else{                int q=tree[p].next[i];                if(p)  tree[q].fail=tree[tree[p].fail].next[i];                tree[q].source|=tree[tree[q].fail].source;                tree[q].virus|=tree[tree[q].fail].virus;                que[tail++]=q;            }        }    }}void Bfs(int u){    int que[N];    int head=0,tail=0;    int dist[N];    memset(dist,-1,sizeof(dist));    que[tail++]=pos[u];    dist[pos[u]]=0;    while(head<tail){        int p=que[head++];        for(int i=0;i<2;i++){            int q=tree[p].next[i];            if(dist[q]<0&&!tree[q].virus){                dist[q]=dist[p]+1;                que[tail++]=q;            }        }    }    for(int i=0;i<cnt;i++)        path[u][i]=dist[pos[i]];}int DP(){    memset(dp,-1,sizeof(dp));    dp[0][0]=0;    for(int i=0;i<(1<<n);i++){        for(int j=0;j<cnt;j++){            if(dp[i][j]<0) continue;            for(int k=0;k<cnt;k++){                if(path[j][k]<0) continue;                if(dp[i|tree[pos[k]].source][k]==-1) dp[i|tree[pos[k]].source][k]=dp[i][j]+path[j][k];                else dp[i|tree[pos[k]].source][k]=min(dp[i|tree[pos[k]].source][k],dp[i][j]+path[j][k]);            }        }    }    int ans=inf;    for(int i=0;i<cnt;i++)        if(dp[(1<<n)-1][i]!=-1) ans=min(ans,dp[(1<<n)-1][i]);    return ans==inf?-1:ans;}int main(){   // freopen("1.in","r",stdin);    while(scanf("%d%d",&n,&m)!=EOF&&n+m){        tree[0].Init();cnt=1;tot=0;        for(int i=0;i<n;i++){            scanf("%s",str);            Insert(str,strlen(str),i);        }        for(int i=0;i<m;i++){            scanf("%s",str);            Insert(str,strlen(str),-1);        }        Bulid_fail();        for(int i=1;i<=tot;i++)            if(tree[i].source)                pos[cnt++]=i;        for(int i=0;i<cnt;i++)            Bfs(i);        printf("%d\n",DP());    }    return 0;}