HDU 4511 小明系列故事——女友的考验 【AC自动机+DP】

限制条件用AC自动机找出合法状态

但是这题的DP并不是矩阵加速的那种，转移需要涉及到实际点

所以定义 dp[i][j] ：当前点为 i，在自动机上为 j 的最小路程，枚举点转移

顺带一提，枚举转移是n*(n+自动机节点数)，因为除了自动机的根节点其它点唯一对应一个编号。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 505
#define maxc 50
using namespace std;
const double INF = 1e20;
int n,m,s[5];
double x[maxc],y[maxc];
struct ac_automaton
{
	int ch[maxn][maxc],fail[maxn],id[maxn],tot;
	bool cnt[maxn];
	double f[maxc][maxn];
	inline void init(){
		memset(ch,0,sizeof ch);memset(cnt,0,sizeof cnt);
		tot=0;
	}
	inline void insert(int *s,int len)
	{
		int r=0,v;
		for(int i=0;i<len;r=ch[r][v],i++)
			if(!ch[r][v=s[i]]) id[ch[r][v]=++tot]=s[i];
		cnt[r]=1;
	}
	inline double dis(int i,int j){return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));}
	int Q[maxn],head,tail;
	inline void build()
	{
		Q[head=0]=0,tail=1;
		while(head<tail)
		{
			int r=Q[head++],c;
			for(int i=0;i<n;i++)
			{
				if(c=ch[r][i]) fail[c]=r?ch[fail[r]][i]:0,Q[tail++]=c,cnt[c]|=cnt[fail[c]];
				else ch[r][i]=ch[fail[r]][i];
			}
		}
	}
	inline void solve(int n)
	{
		for(int i=0;i<n;i++) for(int j=0;j<=tot;j++) f[i][j]=INF;
		f[0][ch[0][0]]=0;
		for(int i=0;i<n-1;i++)
			for(int j=0;j<=tot;j++) if(f[i][j]!=INF)
				for(int k=i+1;k<n;k++)
					if(!cnt[ch[j][k]]) f[k][ch[j][k]]=min(f[k][ch[j][k]],f[i][j]+dis(i,k));
		double ans=INF;
		for(int i=0;i<=tot;i++) ans=min(ans,f[n-1][i]);
		if(ans<INF) printf("%.2lf\n",ans);
		else puts("Can not be reached!");
	}
}A;
int main()
{
	while(scanf("%d%d",&n,&m),n)
	{
		A.init();
		for(int i=0;i<n;i++) scanf("%lf%lf",&x[i],&y[i]);
		for(int i=0,k;i<m;i++)
		{
			scanf("%d",&k);
			for(int j=0;j<k;j++) scanf("%d",&s[j]),s[j]--;
			A.insert(s,k);
		}
		A.build();
		A.solve(n);
	}
}