HDU4511 小明系列故事——女友的考验（AC自动机 + DP）

题目大概说有平面有n个点，从1点出发走到n点，每一步只能走到序号比当前更大的点且走的序列不能包含给定的m个序列中的任何一个，问1走到n的最短路。

用m个序列建个AC自动机，后缀包含整个序列的结点标记一下，然后用dp[u][S]表示走到u点且走的序列的后缀状态是自动机上第S个结点的最短路，这样在AC自动机上跑着转移就OK了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
using namespace std;
#define MAXN 666
int tn,ch[MAXN][55],fail[MAXN];
bool flag[MAXN];
void insert(int *a,int n){
    int x=0;
    for(int i=0; i<n; ++i){
        int y=a[i];
        if(ch[x][y]==0) ch[x][y]=++tn;
        x=ch[x][y];
    }
    flag[x]=1;
}
int n;
void getfail(){
    memset(fail,0,sizeof(fail));
    queue<int> que;
    for(int i=1; i<=n; ++i){
        if(ch[0][i]) que.push(ch[0][i]);
    }
    while(!que.empty()){
        int x=que.front(); que.pop();
        for(int i=1; i<=n; ++i){
            if(ch[x][i]){
                que.push(ch[x][i]);
                fail[ch[x][i]]=ch[fail[x]][i];
                flag[ch[x][i]]|=flag[ch[fail[x]][i]];
            }else ch[x][i]=ch[fail[x]][i];
        }
    }
}
double x[55],y[55],dist[55][55],d[55][MAXN];
int a[6];
int main(){
    int m,k;
    while(~scanf("%d%d",&n,&m) && (n||m)){
        for(int i=1; i<=n; ++i){
            scanf("%lf%lf",x+i,y+i);
        }
        tn=0;
        memset(ch,0,sizeof(ch));
        memset(flag,0,sizeof(flag));
        while(m--){
            scanf("%d",&k);
            for(int i=0; i<k; ++i) scanf("%d",a+i);
            insert(a,k);
        }
        getfail();
        memset(dist,0,sizeof(dist));
        for(int i=1; i<=n; ++i){
            for(int j=i+1; j<=n; ++j){
                dist[i][j]=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
            }
        }
        for(int i=1; i<=n; ++i){
            for(int j=0; j<=tn; ++j) d[i][j]=-1;
        }
        d[0][0]=0;
        for(int i=0; i<n; ++i){
            for(int j=0; j<=tn; ++j){
                if(d[i][j]<0) continue;
                if(i==0 && !flag[ch[0][1]]){
                    d[1][ch[0][1]]=0;
                    continue;
                }
                for(int k=i+1; k<=n; ++k){
                    if(flag[ch[j][k]]) continue;
                    if(d[k][ch[j][k]]<0 || d[k][ch[j][k]]>d[i][j]+dist[i][k]) d[k][ch[j][k]]=d[i][j]+dist[i][k];
                }
            }
        }
        double res=-1;
        for(int i=0; i<=tn; ++i){
            if(d[n][i]<0) continue;
            if(res<0 || res>d[n][i]) res=d[n][i];
        }
        if(res<0) puts("Can not be reached!");
        else printf("%.2f
",res);
    }
    return 0;
}