hdu 4511 小明系列故事——女友的考验(AC自动机+dp)

题目链接：hdu 4511 小明系列故事——女友的考验

题意：

中文不解释。

题解：

考虑将那些限制的条件插进AC自动机里面，然后在AC自动机上dp一下最短路就行了。

考虑dp[i][j]表示当前在第i个点上，在AC自动机的第j个点上的最短路。

然后dp转移看代码。

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
typedef pair<double,double>P;
const int AC_N=1000,tyn=50;
int n,m,s[50];
P a[60];
double dp[55][505],inf=1e18;

struct AC_automation{
    int tr[AC_N][tyn],cnt[AC_N],Q[AC_N],fail[AC_N],tot;
    void nw(){cnt[++tot]=0,fail[tot]=0;memset(tr[tot],0,sizeof(tr[tot]));}
    void init(){tot=-1,fail[0]=-1,nw();}
    void insert(int *s,int x=0){
        for(int i=0,w;s[i];x=tr[x][w],i++)
            if(!tr[x][w=s[i]])nw(),tr[x][w]=tot;
        cnt[x]++;
    }
    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];
                cnt[tr[x][i]]+=cnt[tr[fail[x]][i]];
            }
            else tr[x][i]=tr[fail[x]][i];
    }
    double dis(P &a,P &b)
    {
        return sqrt((a.first-b.first)*(a.first-b.first)
              +(a.second-b.second)*(a.second-b.second));
    }
    void solve()
    {
        F(i,0,n)F(j,0,tot)dp[i][j]=inf;
        dp[1][tr[0][1]]=0;
        F(i,1,n-1)F(j,0,tot)if(dp[i][j]!=inf)
        {
            F(k,i+1,n)
            {
                int nxt=tr[j][k];
                if(cnt[nxt])continue;
                dp[k][nxt]=min(dp[k][nxt],dp[i][j]+dis(a[i],a[k]));
            }
        }
        double ans=inf;
        F(j,0,tot)ans=min(ans,dp[n][j]);
        if(ans==inf)puts("Can not be reached!");
        else printf("%.2f
",ans);
    }
}AC;

int main()
{
    while(scanf("%d%d",&n,&m),n+m)
    {
        F(i,1,n)scanf("%lf%lf",&a[i].first,&a[i].second);
        AC.init();
        F(i,1,m)
        {
            memset(s,0,sizeof(s));
            int num,x;
            scanf("%d",&num);
            F(j,1,num)scanf("%d",s+j-1);
            AC.insert(s);
        }
        AC.build(),AC.solve();
    }
    return 0;
}