// compile: make data
// run: ./data < data.in
#include <bits/stdc++.h>
using namespace std;
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
#ifdef LOCAL
#include <debug/codeforces.h>
#define debug(x...) _debug_print(#x, x);
#define Debug(x...) _debug_print_format(#x, x);
std::ifstream terminal("/dev/tty");
#define PP cerr<<"\033[1;30mpause...\e[0m",terminal.ignore();
#else
#define debug(x...)
#define Debug(x...)
#define PP
#endif
template<typename...Args> void print_(Args...args){((cout<<args<<" "),...)<<endl;}
#define VI vector<int>
#define VII vector<vector<int>>
#define VIII vector<vector<vector<int>>>
#define rep(i,a,b) for(int i=(a);i<(int)(b);++i)
#define sz(v) ((int)(v).size())
#define print(...) print_(__VA_ARGS__);
#define FIND(a, x) ((find(a.begin(),a.end(),(x))!=a.end())?1:0)
#define cmin(x,...) x=min({(x), __VA_ARGS__})
#define cmax(x,...) x=max({(x), __VA_ARGS__})
#define INTMAX (int)(9223372036854775807)
#define INF (int)(1152921504606846976)
#define NaN (int)(0x8b88e1d0595d51d1)
#define double long double
#define int long long
#define uint unsigned long long
#define MAXN 200010
struct bitree {
int n;
vector<vector<int>> c;
vector<int> sum;
bitree(vector<int> &a) {
n = sz(a); sum.push_back(0);
rep(i, 0, n) sum.push_back(sum[i] + a[i]);
c.resize(2, vector<int>(n+1, 0));
}
int lowbit(int x) {return x & (-x); }
void add(int b, int x, int y) {for (; x <= n; x += lowbit(x)) c[b][x] += y; }
int get(int b, int x) {int ans = 0; for (; x; x -= lowbit(x)) ans += c[b][x]; return ans; }
void modify(int p, int k) {modify(p, p, k); }
void modify(int l, int r, int k) {add(0, l+1, k); add(0, r+2, -k); add(1, l+1, (l+1)*k); add(1, r+2, -(r+2)*k); }
int query(int p) {return query(p, p); }
int query(int l, int r) {int ans = sum[r+1] + (r+2)*get(0, r+1) - get(1, r+1); ans -= sum[l] + (l+1)*get(0, l) - get(1, l); return ans; }
};
struct sam {
struct node {
map<char, node*> next;
node *link;
int minlen, maxlen;
int suf = 0;
int endpos;
int tot = 0, size = 0;
};
vector<node> datalist;
int total = 0;
node *last = nullptr, *root = nullptr;
string s;
node *newnode(int length = 0, node *p = nullptr) {
node *x = &datalist[total++];
x->minlen = 0;
x->maxlen = length;
x->link = p;
return x;
}
sam(string str) {
s = str;
datalist.resize(s.size() * 2);
root = last = newnode();
root->endpos = 0;
for (int i = 0; i < sz(s); ++i) {
char c = s[i];
node *cur = newnode(i+1, root);
cur->endpos = i+1;
cur->tot = 1;
node *parent = last;
while (parent && !parent->next.count(c)) {
parent->next[c] = cur;
parent = parent->link;
}
if (parent) {
node *q = parent->next[c];
if (parent->maxlen + 1 == q->maxlen) cur->link = q;
else {
node *nq = newnode(parent->maxlen + 1, q->link);
nq->endpos = i+1;
nq->next = q->next;
q->link = cur->link = nq;
while (parent && parent->next[c] == q) {
parent->next[c] = nq;
parent = parent->link;
}
}
}
last = cur;
}
// minlen
queue<node*> q; q.push(root);
while (!q.empty()) {
node *x = q.front(); q.pop();
for (auto [c, p]: x->next) {
if (p->minlen) continue;
p->minlen = x->minlen + 1;
q.push(p);
}
}
// is suffix
for (node *x = last; x && x != root; x = x->link) x->suf = 1;
// size of endpos set
vector<vector<node*>> e(total);
rep(i, 0, total) if (datalist[i].link) e[datalist[i].link - &datalist[0]].push_back(&datalist[i]);
function<void(int)> dfs_tot = [&](int x) {
for (auto y: e[x]) {
dfs_tot(y - &datalist[0]);
datalist[x].tot += datalist[y - &datalist[0]].tot;
}
};
dfs_tot(0);
// num of strings
function<int(node*)> dfs_size = [&](node *x) {
if (x->size) return x->size;
x->size = x==root ? 0 : 1;
for (auto [c, p]: x->next) x->size += dfs_size(p);
return x->size;
};
dfs_size(root);
}
string getstr(node *x) {
return s.substr(x->endpos - x->maxlen, x->maxlen);
}
void graphviz_dump(string filename = "graph.dot") {
FILE *f = fopen(filename.c_str(), "w");
fprintf(f, "digraph {\n");
rep(i, 0, total) {
node *x = &datalist[i];
fprintf(f, " %lld[%s shape=\"record\" label=\"", i, x->suf ? "color=\"red\"" : "");
vector<string> desired = {
"min=" + to_string(x->minlen),
"max=" + to_string(x->maxlen),
"endpos=" + to_string(x->endpos),
"sz=" + to_string(x->size)
};
rep(j, 0, sz(desired)) {
fprintf(f, "%s%s", j ? " | " : "", desired[j].c_str());
}
fprintf(f, "\"];\n");
}
rep(i, 0, total) {
for (auto [c, p]: datalist[i].next) {
fprintf(f, " %lld -> %ld ;\n", i, p - &datalist[0], c);
}
}
// rep(i, 0, total) if (datalist[i].link) fprintf(f, " %lld -> %ld [color=\"blue\" style=\"dashed\"];\n", i, datalist[i].link - &datalist[0]);
fprintf(f, "}\n");
}
};
int32_t main() {
ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr);
string s; cin >> s;
sam a(s);
vector<int> cnt(sz(s)+1);
bitree t(cnt);
rep(i, 0, a.total) {
auto *x = &a.datalist[i];
t.modify(x->minlen, x->maxlen, 1);
}
rep(i, 1, sz(s)+1) cout << t.query(i) << " \n"[i==sz(s)];
return 0;
}
No Comments