ppt上的题目,这题主要是练习把每个状态的|right|求出来(还是不理解sam的性质,自己没想到咋求,看了爱神的代码),然后利用每个状态s的|right|更新f(s.len)即f(s.len) = max(|s.right|, f(s.len)),ppt上题解还要用f(i)更新f(i-1),虽然这一步实际上是不必要的,但这样简化了思维的复杂度,题解上用f(i)更新f(i-1)是考虑到如果长度为i的子串一定包含一个长度为i-1的子串,但是根据parent(即自动机里的失配)指针的性质,所有长度为i的子串的|right|必定会更新到一个长度为i-1的子串上的,但可能会出现这种情况即一个状态包含多个长度不同的子串,例如 一个状态表示abb和bb这俩个子串,根绝我们的更新办法,看上去bb的|right|将不会被更新,但注意到abb中还包含ab这个子串,所以只要更新ab的|right|即可,实际上ab也有可能会被忽略掉,但递归的想最终总会找到一个len为2的状态的,所以就不需要再用f(i)更新f(i-1)了。
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <queue>
#include <algorithm>
#include <vector>
#include <cstring>
#include <stack>
#include <cctype>
#include <utility>
#include <map>
#include <string>
#include <climits>
#include <set>
#include <string>
#include <sstream>
#include <utility>
#include <ctime>
using std::priority_queue;
using std::vector;
using std::swap;
using std::stack;
using std::sort;
using std::max;
using std::min;
using std::pair;
using std::map;
using std::string;
using std::cin;
using std::cout;
using std::set;
using std::queue;
using std::string;
using std::istringstream;
using std::make_pair;
using std::getline;
using std::greater;
using std::endl;
using std::multimap;
using std::deque;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int, int> PAIR;
typedef multimap<int, int> MMAP;
const int MAXN(250010);
const int SIGMA_SIZE(26);
const int MAXM(110);
const int MAXE(300010);
const int MAXH(18);
const int INFI((INT_MAX-1) >> 1);
const int MOD(2520);
const ULL BASE(31);
const ULL LIM(1000000000000000ull);
struct SAM
{
struct NODE
{
NODE *par, *ch[SIGMA_SIZE];
int val, count;
};
NODE *root, *last;
NODE pool[MAXN << 1];
int size;
void init()
{
last = root = pool;
root->par = 0;
root->val = root->count = 0;
memset(root->ch, 0, sizeof(root->ch));
size = 1;
}
NODE *newnode(int tv)
{
pool[size].val = tv;
pool[size].count = 0;
memset(pool[size].ch, 0, sizeof(pool[size].ch));
return pool+size++;
}
void extend(int id)
{
NODE *p = last, *np = newnode(last->val+1);
while(p && p->ch[id] == 0)
p->ch[id] = np, p = p->par;
if(p == 0)
np->par = root;
else
{
NODE *q = p->ch[id];
if(p->val+1 == q->val)
np->par = q;
else
{
NODE *nq = newnode(p->val+1);
memcpy(nq->ch, q->ch, sizeof(nq->ch));
nq->par = q->par;
q->par = nq;
np->par = nq;
while(p && p->ch[id] == q)
p->ch[id] = nq, p = p->par;
}
}
last = np;
}
};
SAM sam;
int cnt[MAXN];
int buc[MAXN << 1];
int table[MAXN];
char str[MAXN];
int main()
{
while(~scanf("%s", str))
{
sam.init();
int len = strlen(str);
for(int i = 0; i < len; ++i)
sam.extend(str[i]-'a');
memset(cnt, 0, sizeof(cnt[0])*(len+1));
for(int i = 0; i < sam.size; ++i)
++cnt[sam.pool[i].val];
for(int i = 1; i <= len; ++i)
cnt[i] += cnt[i-1];
for(int i = 0; i < sam.size; ++i)
buc[--cnt[sam.pool[i].val]] = i;
SAM::NODE *tp = sam.root;
for(int i = 0; i < len; ++i)
(tp = tp->ch[str[i]-'a'])->count = 1;
memset(table, 0, sizeof(table[0])*(len+1));
for(int i = sam.size; i >= 0; --i)
{
tp = sam.pool+buc[i];
table[tp->val] = max(table[tp->val], tp->count);
if(tp->par)
tp->par->count += tp->count;
}
/* for(int i = len-1; i >= 0; --i)
table[i] = max(table[i], table[i+1]);*/
for(int i = 1; i <= len; ++i)
printf("%d\n", table[i]);
}
return 0;
}