小明系列故事——女友的考验 Time Limit: 500/200 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 863 Accepted Submission(s): 192
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
Source 2013腾讯编程马拉松初赛第二场（3月22日）
Recommend liuyiding | We have carefully selected several similar problems for you: 5209 5208 5207 5206 5205
告诉你某些路径不能存在，这是个很明显的提示，就是个自动机dp，，dp[i][j]表示在自动机节点ｉ，并且在点ｊ时的最短路径，一开始我的循环方式可能不对，后来换了下就可以了，注意坐标还是直接用double 读入吧
/*************************************************************************> File Name: hdu4511.cpp> Author: ALex> Mail: zchao1995@gmail.com> Created Time: 2015年04月16日 星期四 18时24分11秒 ************************************************************************/;const double pi = acos(-1.0);const int inf = 0x3f3f3f3f;const double eps = 1e-15;LL;typedef pair <int, int> PLL;const int MAX_NODE = 510;int n;double dp[MAX_NODE][55];struct node{double x, y;}point[55];double dist(node a, node b){double x = (a.x – b.x);double y = (a.y – b.y);return sqrt(x * x + y * y);}struct AC_Automation{int next[MAX_NODE][55];int fail[MAX_NODE];int end[MAX_NODE];int root, L;int newnode(){for (int i = 0; i < n; ++i){next[L][i] = -1;}end[L++] = 0;return L – 1;}void init (){L = 0;root = newnode();}void Build_Trie (vector <int> buf){int now = root;int len = buf.size();for (int i = 0; i < len; ++i){if (next[now][buf[i]] == -1){next[now][buf[i]] = newnode();}now = next[now][buf[i]];}end[now] = 1;}void Build_AC (){queue <int> qu;fail[root] = root;for (int i = 0; i < n; ++i){if (next[root][i] == -1){next[root][i] = root;}else{fail[next[root][i]] = root;qu.push (next[root][i]);}}while (!qu.empty()){int now = qu.front();qu.pop();if (end[fail[now]]){end[now] = 1;}for (int i = 0; i < n; ++i){if (next[now][i] == -1){next[now][i] = next[fail[now]][i];}else{fail[next[now][i]] = next[fail[now]][i];qu.push(next[now][i]);}}}}void solve(){for (int i = 0; i < L; ++i){for (int j = 0; j < n; ++j){dp[i][j] = (1LL << 60);}}dp[next[root][0]][0] = 0.0;for (int i = 0; i < n; ++i){for (int j = 0; j < L; ++j){if (dp[j][i] < (1LL << 60)){for (int k = i + 1; k < n; ++k){int nxt = next[j][k];if (end[nxt]){continue;}dp[nxt][k] = min(dp[nxt][k], dp[j][i] + dist(point[k], point[i]));}}}}double ans = (1LL << 60);for (int i = 0; i < L; ++i){if (end[i] == 0){ans = min(ans, dp[i][n – 1]);}}if (ans == (1LL << 60)){printf(“Can not be reached!\n”);}else{printf(“%.2f\n”, ans);}}}AC;vector <int> buf;int main(){int m;while (~scanf(“%d%d”, &n, &m)){if (!n && !m){break;}AC.init();int tmp;for (int i = 0; i < n; ++i){scanf(“%lf%lf”, &point[i].x, &point[i].y);}for (int i = 1; i <= m; ++i){int k;buf.clear();scanf(“%d”, &k);while (k–){scanf(“%d”, &tmp);–tmp;buf.push_back(tmp);}AC.Build_Trie(buf);}AC.Build_AC();AC.solve();}return 0;}
只要你扬帆，便会有八面来风。启程了，人的生命才真正开始。
