【NOI2015】品酒大会 SAM

论没有看完题目的危害，以及预处理不做号的危害。两个小时过去了…………………………………………

不难。把字符串逆序建后缀自动机，利用fail边得到后缀树。（理解上的，实际上没有在树上跑。）

由于后缀树上每一个叶子节点所代表的串是原串的后缀。现在逆序之后，就变成了前缀。那

因为一个节点x的sz数组表示这个节点所代表的所有的串都出现了sz[x]次。而且该串所代表的字符串长度为 dep[fa[x]] + 1 ~ dep[x] 所以理论上考虑一个节点对答案的贡献是是要对 dep[fa[x]+1] ~　dep[x]都加一次的。但是很显然有前缀和的思想，那么就不用了。

总而言之。后缀自动机的思想就是用一个长串已有的信息，和它的后缀的信息汇总。进而得出全部的信息。

（特别地，对于任意的 1≤p,q≤n1≤p,q≤n，p≠qp≠q，第 pp 杯酒和第 qq杯酒都是“0相似”的。）不知道是否有人和我一样没看题。而且有些maxl，minl没初始化，WA了好久。心好累。

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 #define rep(i,j,k) for(register int i = j; i <= k; i++)
 5 #define dow(i,j,k) for(register int i = j; i >= k; i--)
 6 #define maxn 302333
 7 #define ll long long
 8 #define inf 0x7fffffff
 9 using namespace std;
10
11 inline int read() {
12     int s = 0, t = 1; char c = getchar();
13     while( !isdigit(c) ) { if( c == ‘-‘ ) t = -1; c = getchar(); }
14     while( isdigit(c) ) s = s * 10 + c - 48, c = getchar();
15     return s * t;
16 }
17
18 int u, p, last = 1, tot = 1, v, nv, dep[maxn<<1], fa[maxn<<1], sz[maxn<<1];
19 int maxl[maxn<<1][2], minl[maxn<<1][2], son[maxn<<1][26];
20 inline void extend(int c,int key) {
21     u = last; last = p = ++tot; dep[p] = dep[u] + 1; sz[p] = 1;
22     maxl[p][0] = minl[p][0] = key; //cout<<"B "<<p<<" "<<key<<endl;
23     maxl[p][1] = -2e9; minl[p][1] = 2e9;
24     while( u && !son[u][c] ) son[u][c] = p, u = fa[u];
25     if( !u ) fa[p] = 1;
26     else {
27         v = son[u][c];
28         if( dep[v] == dep[u] + 1 ) fa[p] = v;
29         else {
30             nv = ++tot; dep[nv] = dep[u] + 1;
31             fa[nv] = fa[v]; fa[v] = fa[p] = nv;
32             maxl[tot][0] = maxl[tot][1] = -2e9;
33             minl[tot][0] = minl[tot][1] = 2e9;
34             memcpy(son[nv],son[v],sizeof(son[v]));
35             while( son[u][c] == v ) son[u][c] = nv, u = fa[u];
36         }
37     }
38 }
39
40 int q[maxn<<1], tong[maxn<<1] = {0}, key[maxn];
41 inline void pre() {
42     rep(i,1,tot) tong[dep[i]]++;
43     rep(i,1,tot) tong[i] += tong[i-1];
44     dow(i,tot,1) q[tong[dep[i]]--] = i;
45     int x;
46     dow(i,tot,1) x = q[i], sz[fa[x]] += sz[x];
47     dep[0] = -1;
48 }
49 ll cnt[maxn], ans[maxn]; char c[maxn];
50 inline void updatemax(int x,int v) {
51 //    cout<<"M "<<x<<" "<<v<<" "<<maxl[x][1]<<" ";
52     maxl[x][1] = max(maxl[x][1],v);
53 //    cout<<maxl[x][1]<<endl;
54     if( maxl[x][1] > maxl[x][0] ) swap(maxl[x][1],maxl[x][0]);
55 }
56 inline void updatemin(int x,int v) {
57 //    cout<<"I "<<x<<" "<<v<<" "<<minl[x][1]<<" ";
58     minl[x][1] = min(minl[x][1],v);
59 //    cout<<minl[x][1]<<endl;
60     if( minl[x][1] < minl[x][0] ) swap(minl[x][1],minl[x][0]);
61 }
62
63 int main() {
64     int n = read(); scanf("%s", c+1);
65     maxl[1][0] = maxl[1][1] = -2e9; minl[1][0] = minl[1][1] = 2e9;
66     rep(i,1,n) key[i] = read();
67     dow(i,n,1) extend(c[i]-‘a‘,key[i]);
68     pre(); int x;
69     memset(ans,128,sizeof(ans));
70     memset(cnt,0,sizeof(cnt));
71     dow(i,tot,1) {
72         x = q[i];
73         if( sz[x] > 1 ) {
74             cnt[dep[x]] += 1ll * (sz[x] - 1) * sz[x] / 2;
75             if( dep[fa[x]] >= 0 ) cnt[dep[fa[x]]] -= 1ll * (sz[x] - 1) * sz[x] / 2;
76             ans[dep[x]] = max(ans[dep[x]],max(1ll*maxl[x][0]*maxl[x][1],1ll*minl[x][0]*minl[x][1]));
77         }
78         updatemax(fa[x],maxl[x][0]);
79         updatemax(fa[x],maxl[x][1]);
80         updatemin(fa[x],minl[x][0]);
81         updatemin(fa[x],minl[x][1]);
82     }
83     dow(i,n-1,0) cnt[i] += cnt[i+1];
84     dow(i,n-1,0) ans[i] = max(ans[i],ans[i+1]);
85     rep(i,0,n-1) if( !cnt[i] ) ans[i] = 0;
86     rep(i,0,n-1) printf("%lld %lld\n", cnt[i],ans[i]);
87  //   x = read();
88     return 0;
89 }