HDU 4511 小明系列故事——女友的考驗 （AC自動機+DP）

小明系列故事——女友的考驗

Time Limit: 500/200 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 654    Accepted Submission(s): 136

Problem Description

　　終於放寒假了，小明要和女朋友一起去看電影。這天，女朋友想給小明一個考驗，在小明正准備出發的時候，女朋友告訴他，她在電影院等他，小明過來的路線必須滿足給定的規則：
　　1、假設小明在的位置是1號點，女朋友在的位置是n號點，則他們之間有n-2個點可以走，小明每次走的時候只能走到比當前所在點編號大的位置；
　　2、小明來的時候不能按一定的順序經過某些地方。比如，如果女朋友告訴小明不能經過1 -> 2 -> 3，那么就要求小明來的時候走過的路徑不能包含有1 -> 2 -> 3這部分，但是1 -> 3 或者1 -> 2都是可以的，這樣的限制路徑可能有多條。
　　這讓小明非常頭痛，現在他把問題交給了你。
　　特別說明，如果1 2 3這三個點共線，但是小明是直接從1到3然后再從3繼續，那么此種情況是不認為小明經過了2這個點的。
　　現在，小明即想走最短的路盡快見到女朋友，又不想打破女朋友的規定，你能幫助小明解決這個問題嗎？

 

 

Input

　　輸入包含多組樣例，每組樣例首先包含兩個整數n和m，其中n代表有n個點，小明在1號點，女朋友在n號點，m代表小明的女朋友有m個要求；
　　接下來n行每行輸入2個整數x 和y（x和y均在int范圍），代表這n個點的位置（點的編號從1到n）；
　　再接着是m個要求，每個要求2行，首先一行是一個k，表示這個要求和k個點有關，然后是順序給出的k個點編號，代表小明不能走k1 -> k2 -> k3 ……-> ki這個順序的路徑；
　　n 和 m等於0的時候輸入結束。

　　 [Technical Specification]
　　2 <= n <= 50
　　1 <= m <= 100
　　2 <= k <= 5

 

 

Output

　　對於每個樣例，如果存在滿足要求的最短路徑，請輸出這個最短路徑，結果保留兩位小數；否則，請輸出”Can not be reached!” （引號不用輸出）。

 

 

Sample Input

3 1 1 1 2 1 3 1 2 1 2 2 1 0 0 1 1 2 1 2 5 3 0 0 5 3 1 2 1 22 5 21 3 1 2 3 2 4 5 2 1 5 0 0

 

 

Sample Output

2.00 Can not be reached! 21.65

 

 

一些路徑不能走，通過AC自動機得到狀態轉移。

之后dp[i][j] 表示在i點，狀態在j的距離。

 

網上看到一個floyed的做法，最后發現可以 cha掉，路徑不能只看兩個點的。如1->3  1->2->3  這樣1是不能到3的。

 

 

坐標int相減會溢出int,最好直接用double.

 

/* ***********************************************
Author        :kuangbin
Created Time  :2014/3/2 22:07:17
File Name     :E:\2014ACM\專題學習\AC自動機\HDU4511.cpp
************************************************ */

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
const double INF = 1e20;
int n;

pair<int,int> p[100];
double dis(pair<int,int>a,pair<int,int>b)//坐標相減爆int
{
    return sqrt((double)(1.0 * a.first - b.first) * (1.0 * a.first - b.first) + (double)(1.0 * a.second - b.second)*(1.0 * a.second - b.second));
}
double dp[55][1000];
void CheckMin(double &a,double b)
{
    a = min(a,b);
}
struct Trie
{
    int next[1000][55],fail[1000],end[1000];
    int root,L;
    int newnode()
    {
        for(int i = 1; i <= n;i++)
            next[L][i] = -1;
        end[L++] = 0;
        return L-1;
    }
    void init()
    {
        L = 0;
        root = newnode();
    }
    void insert(int a[],int cnt)
    {
        int now = root;
        for(int i  = 0;i < cnt;i++)
        {
            if(next[now][a[i]] == -1)
                next[now][a[i]] = newnode();
            now = next[now][a[i]];
        }
        end[now] = 1;
    }
    void build()
    {
        queue<int>Q;
        fail[root] = root;
        for(int i = 1;i <= n;i++)
            if(next[root][i] == -1)
                next[root][i] = root;
            else 
            {
                fail[next[root][i]] = root;
                Q.push(next[root][i]);
            }
        while(!Q.empty())
        {
            int now = Q.front();
            Q.pop();
            end[now] |= end[fail[now]];
            for(int i = 1;i <= n;i++)
                if(next[now][i] == -1)
                    next[now][i] = next[fail[now]][i];
                else
                {
                    fail[next[now][i]] = next[fail[now]][i];
                    Q.push(next[now][i]);
                }
        }

    }
    void solve()
    {
        for(int i = 1;i <= n;i++)
            for(int j = 0;j < L;j++)
                dp[i][j] = INF;
        dp[1][next[root][1]] = 0;
        for(int i = 1;i < n;i++)
            for(int j = 0;j < L;j++)
                if(dp[i][j] < INF)
                {
                    for(int k = i+1;k <= n;k++)
                    {
                        int ss = next[j][k];
                        if(end[ss])continue;
                        CheckMin(dp[k][ss],dp[i][j] + dis(p[i],p[k]));
                    }
                }
        double ans = INF;
        for(int i = 0;i < L;i++)
            if(dp[n][i] < INF)
                CheckMin(ans,dp[n][i]);
        if(ans == INF)printf("Can not be reached!\n");
        else printf("%.2f\n",ans);
    }
}ac;

int a[10];
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int m;
    while(scanf("%d%d",&n,&m) == 2)
    {
        if(n == 0 && m == 0)break;
        for(int i = 1;i <= n;i++)
            scanf("%d%d",&p[i].first,&p[i].second);
        ac.init();
        int k;
        while(m--)
        {
            scanf("%d",&k);
            for(int i = 0;i < k;i++)
                scanf("%d",&a[i]);
            ac.insert(a,k);
        }
        ac.build();
        ac.solve();
    }
    return 0;
}