限制条件用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);
}
}