题意:从 1 走到 n,要求所走路径不能出现给定的路径,求最短路
思路:因为要求不能出现给定路径,那么我可以求助ac自动机完成判断。
我们可以在build的时候标记哪些路径不能出现,显然下面这种表示后缀不能出现,那么他也不能出现
if(node[node[u].fail].cnt && u) node[u].cnt = 1; //都不能取
然后再把图建完整。因为如果一个路径不在Trie中有两种情况,一种是他可能是某个不能走的串的前缀,那么我就重新指向这个不能走的串,比如Trie中只有AT,那么我走ATA相当于走AT再走回A;另一种是这个路径根本没出现过,那么相当于从0开始走。
if(!node[u].next[i]){ if(u == 0) node[u].next[i] = 0; else node[u].next[i] = node[node[u].fail].next[i]; }
dp[i][j]表示走到i点在AC自动机上走到j点的最小路径。
代码:
#include<cmath> #include<set> #include<map> #include<queue> #include<cstdio> #include<vector> #include<cstring> #include <iostream> #include<algorithm> using namespace std; typedef long long ll; typedef unsigned long long ull; const int maxn = 500 + 5; const int M = 50 + 5; const ull seed = 131; const double INF = 1e20; const int MOD = 1000000007; double x[M], y[M]; int s[M]; double dp[M][maxn]; int m, n; double dis(int i, int j){ return sqrt((x[i] - x[j]) * (x[i] - x[j]) + (y[i] - y[j]) * (y[i] - y[j])); } struct Aho{ struct state{ int next[M]; int fail, cnt; }node[maxn]; int size; queue<int> q; void init(){ size = 0; newtrie(); while(!q.empty()) q.pop(); } int newtrie(){ memset(node[size].next, 0, sizeof(node[size].next)); node[size].cnt = node[size].fail = 0; return size++; } void insert(int s[], int len){ int now = 0; for(int i = 0; i < len; i++){ int c = s[i]; if(node[now].next[c] == 0){ node[now].next[c] = newtrie(); } now = node[now].next[c]; } node[now].cnt = 1; } void build(){ node[0].fail = -1; q.push(0); while(!q.empty()){ int u = q.front(); q.pop(); if(node[node[u].fail].cnt && u) node[u].cnt = 1; //都不能取 for(int i = 0; i <= n; i++){ if(!node[u].next[i]){ if(u == 0) node[u].next[i] = 0; else node[u].next[i] = node[node[u].fail].next[i]; } else{ if(u == 0) node[node[u].next[i]].fail = 0; else{ int v = node[u].fail; while(v != -1){ if(node[v].next[i]){ node[node[u].next[i]].fail = node[v].next[i]; break; } v = node[v].fail; } if(v == -1) node[node[u].next[i]].fail = 0; } q.push(node[u].next[i]); } } } } void query(){ for(int i = 0; i <= n; i++){ for(int j = 0; j < size; j++){ dp[i][j] = INF; } } dp[1][node[0].next[1]] = 0; for(int i = 1; i <= n; i++){ for(int j = 0; j < size; j++){ if(dp[i][j] == INF) continue; for(int k = i + 1; k <= n; k++){ if(node[node[j].next[k]].cnt) continue; dp[k][node[j].next[k]] = min(dp[k][node[j].next[k]], dp[i][j] + dis(i, k)); } } } double ans = INF; for(int i = 0; i < size; i++) if(dp[n][i] < INF) ans = min(ans, dp[n][i]); if(ans == INF) printf("Can not be reached!\n"); else printf("%.2f\n", ans); } }ac; int main(){ while(~scanf("%d%d", &n, &m) && n + m){ for(int i = 1; i <= n; i++){ scanf("%lf%lf", &x[i], &y[i]); } ac.init(); while(m--){ int k; scanf("%d", &k); for(int i = 0; i < k; i++){ scanf("%d", &s[i]); } ac.insert(s, k); } ac.build(); ac.query(); } return 0; }
原文地址:https://www.cnblogs.com/KirinSB/p/11180002.html
时间: 2024-10-04 03:52:23