1. 程式人生 > >Bzoj3926: [Zjoi2015]諸神眷顧的幻想鄉

Bzoj3926: [Zjoi2015]諸神眷顧的幻想鄉

get using zjoi ref amp typedef AR degree truct

題面

傳送門

Sol

求多個串的不同子串的個數

廣義後綴自動機
也就是可以表示所有的串的所有的後綴的自動機

那麽每次建一個串後,另外一個串接在初始節點下面建就好了

葉子節點最多\(20\)
那麽對於每個葉子結點遍歷一遍樹建立\(sam\)
註意回溯時要把\(last\)指回來

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
 
IL int
Input(){ RG int x = 0, z = 1; RG char c = getchar(); for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1; for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48); return x * z; } const
int maxn(1e5 + 5); const int maxm(4e6 + 5); int trans[10][maxm], fa[maxm], tot = 1, last, len[maxm]; int n, first[maxn], cnt, col[maxn], degree[maxn]; ll ans; struct Edge{ int to, next; } edge[maxn << 1]; IL void Add(RG int u, RG int v){ edge[cnt] = (Edge){v, first[u]}, first[u] = cnt++; } IL void
Extend(RG int c, RG int p){ RG int np = ++tot; len[last = np] = len[p] + 1; while(p && !trans[c][p]) trans[c][p] = np, p = fa[p]; if(!p) fa[np] = 1; else{ RG int q = trans[c][p]; if(len[q] == len[p] + 1) fa[np] = q; else{ RG int nq = ++tot; fa[nq] = fa[q], len[nq] = len[p] + 1; for(RG int i = 0; i < 10; ++i) trans[i][nq] = trans[i][q]; fa[q] = fa[np] = nq; while(p && trans[c][p] == q) trans[c][p] = nq, p = fa[p]; } } } IL void Dfs(RG int u, RG int ff, RG int p){ Extend(col[u], p), p = last; for(RG int e = first[u]; e != -1; e = edge[e].next) if(edge[e].to != ff) Dfs(edge[e].to, u, p); } int main(RG int argc, RG char *argv[]){ n = Input(), Input(); for(RG int i = 1; i <= n; ++i) first[i] = -1, col[i] = Input(); for(RG int i = 1; i < n; ++i){ RG int u = Input(), v = Input(); Add(u, v), Add(v, u), ++degree[u], ++degree[v]; } for(RG int i = 1; i <= n; ++i) if(degree[i] == 1) last = 1, Dfs(i, 0, 1); for(RG int i = 1; i <= tot; ++i) ans += len[i] - len[fa[i]]; printf("%lld\n", ans); return 0; }

Bzoj3926: [Zjoi2015]諸神眷顧的幻想鄉