P3804 【模板】后缀自动机 (SAM)

这个玩意 看了一天说白了 就是我们想要建立一个 和trie一样可以识别 所有子串的 东西 但是显然直接 trie n2 是不行的
所以它其实就是 trie的“缩点” 每个点代表的都是自动机里的一个状态
它代表了 所有endpos相同的 一个集合
分为两类边 一类是 类似 trie
一类是 后缀链接树上的 边
满足什么性质呢？

1. 前一个集合的所代表的最长子串+1=后缀链接树上下一个集合的最小串
   长度
2. 子节点的 endpos集合是父节点 的一个划分

#include <bits/stdc++.h>
using namespace  std;
//#define  int long long
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef vector<int> vi;
#define fi first
#define se second
#define pb  push_back
#define inf 1ll<<62
#define endl "\n"
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define de_bug(x) cerr << #x << "=" << x << endl
#define all(a) a.begin(),a.end()
#define IOS   std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define  fer(i,a,b)  for(int i=a;i<=b;i++)
#define  der(i,a,b)  for(int i=a;i>=b;i--)
const int mod = 1e9 + 7;
const int N = 2e6 + 10;
int n, m , k;
string s;
struct node {
	int len, fa;
	int son[26];
} tr[N];
ll f[N], ans;
vi g[N];
int cnt = 1;
int last = 1;
void extend (int c) {
	int p = last, np = last = ++cnt;
	f[cnt] = 1;
	tr[np].len = tr[p].len + 1;
	for(; p && !tr[p].son[c];  p = tr[p].fa) tr[p].son[c] = np;
	if(!p) tr[np].fa = 1;
	else {
		int q = tr[p].son[c];
		if(tr[q].len == tr[p].len + 1)tr[np].fa = q;
		else {
			int nq = ++cnt;
			tr[nq] = tr[q];
			tr[nq].len = tr[p].len + 1;
			tr[q].fa = tr[np].fa = nq;
			for(; p && tr[p].son[c] == q; p = tr[p].fa) tr[p].son[c] = nq;
		}
	}
}
void dfs(int u) {
	for(auto v : g[u]) {
		dfs(v);
		f[u] += f[v];
	}
	if(f[u] > 1)ans = max(ans, f[u] * tr[u].len);
}
void solve() {
	cin >> s;
	n = s.size();
	s = " " + s;
	for(int i = 1; i <= n; i++)extend(s[i] - 'a');
	/*for(int i=2;i<=cnt;i++){
		cout<<tr[i].fa<<" "<<i<<endl;
	}*/
	for(int i = 2; i <= cnt; i++) {
		g[tr[i].fa].push_back(i);
	}
	dfs(1);
	cout << ans << endl;
}
int  main() {
	IOS;
	int _ = 1;
	//cin>>_;
	while( _-- )
		solve();
}