标签：blog http io ar os sp for on div
http://www.spoj.com/problems/NSUBSTR/
题意：给一个字符串S，令F(x)表示S的所有长度为x的子串中，出现次数的最大值。求F(1)..F(Length(S))
这题做法：
首先建立字符串的后缀自动机。
因为自动机中的每个状态都代表一类子串前缀，且任意状态的最长的|max|所对应的子串是唯一的。
所以我们算出每个子串（即找到的状态是end态），他们的right值为++（即保证长度为当前子串的出现次数为1），然后自底向上在parent树中更新right值（parent树中right值就是出现次数，所以包含子树的要加上子树的right）。
那么可以得到所有长度的right值。right值就是子串的出现次数。
此时可以更新出长度为每个状态长度为|max|的出现次数。
因为我们求出的只是max，可是因为：长子串中一定有短子串那么多个（parent树中就是祖先），那么可以用长串来更新短串。
更新right值时要从长度长的更新（因为长度从长到短，对应了parent树上的叶节点和内节点的关系，状态的right值的数目单调不减），所以我们先基排一次，然后从大到小更新。
这样就求出了所有子串的出现次数。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
struct sam {
static const int N=1000005;
int c[N][26], l[N], f[N], root, last, cnt;
sam() { cnt=0; root=last=++cnt; }
void add(int x) {
int now=last, a=++cnt; last=a;
l[a]=l[now]+1;
for(; now && !c[now][x]; now=f[now]) c[now][x]=a;
if(!now) { f[a]=root; return; }
int q=c[now][x];
if(l[q]==l[now]+1) { f[a]=q; return; }
int b=++cnt;
memcpy(c[b], c[q], sizeof c[q]);
l[b]=l[now]+1;
f[b]=f[q];
f[q]=f[a]=b;
for(; now && c[now][x]==q; now=f[now]) c[now][x]=b;
}
void build(char *s) {
int len=strlen(s);
rep(i, len) add(s[i]-‘a‘);
}
}a;
const int N=250005;
char s[N];
int f[N], r[1000005], t[N], b[1000005];
void getans(int len) {
int *l=a.l, *p=a.f, cnt=a.cnt;
for(int now=a.root, i=0; i<len; ++i) now=a.c[now][s[i]-‘a‘], ++r[now];
for1(i, 1, cnt) ++t[l[i]];
for1(i, 1, len) t[i]+=t[i-1];
for1(i, 1, cnt) b[t[l[i]]--]=i;
for3(i, cnt, 1) r[p[b[i]]]+=r[b[i]];
for1(i, 1, cnt) f[l[i]]=max(f[l[i]], r[i]);
for3(i, len, 1) f[i]=max(f[i+1], f[i]);
for1(i, 1, len) printf("%d\n", f[i]);
}
int main() {
scanf("%s", s);
a.build(s);
getans(strlen(s));
return 0;
}
8222. SubstringsProblem code: NSUBSTR |
You are given a string S which consists of 250000 lowercase latin letters at most. We define F(x) as the maximal number of times that some string with length x appears in S. For example for string ‘ababa‘ F(3) will be 2 because there is a string ‘aba‘ that occurs twice. Your task is to output F(i) for every i so that 1<=i<=|S|.
String S consists of at most 250000 lowercase latin letters.
Output |S| lines. On the i-th line output F(i).
Input:
ababa
Output:
3
2
2
1
1
标签：blog http io ar os sp for on div
原文地址：http://www.cnblogs.com/iwtwiioi/p/4141913.html
