# | Time | Username | Problem | Language | Result | Execution time | Memory |
---|---|---|---|---|---|---|---|
22821 | ↓우리보다잘하는팀 (#40) | New Ocurrences (KRIII5_NO) | C++98 | 1000 ms | 5552 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include<cstdio>
#include<algorithm>
using namespace std;
int n, C[101000], ord1[101000], ord2[101000], SA[101000], Rank[101000], LCP[101000];
char p[101000];
struct AA{
int a, b;
}w[101000];
void Suffix_Array(){
int i, M = max(n, 200), LL = 1;
for(i=0;i<n;i++){
w[i].a = p[i],w[i].b = p[i+1];
}
while(1){
for(i=0;i<=M;i++)C[i]=0;
for(i=0;i<n;i++)C[w[i].b]++;
for(i=1;i<=M;i++)C[i]+=C[i-1];
for(i=n-1;i>=0;i--)ord1[--C[w[i].b]] = i;
for(i=0;i<=M;i++)C[i]=0;
for(i=0;i<n;i++)C[w[i].a]++;
for(i=1;i<=M;i++)C[i]+=C[i-1];
for(i=n-1;i>=0;i--)ord2[--C[w[ord1[i]].a]] = ord1[i];
int cnt = 0;
for(i=0;i<n;i++){
if(i==0||w[ord2[i]].a!=w[ord2[i-1]].a||w[ord2[i]].b!=w[ord2[i-1]].b)cnt++;
Rank[ord2[i]] = cnt;
}
if(cnt == n)break;
LL<<=1;
for(i=0;i<n;i++){
w[i].a = Rank[i];
if(i+LL >= n)w[i].b = 0;
else w[i].b = Rank[i+LL];
}
}
for(i=0;i<n;i++)SA[Rank[i]] = i;
int t = 0, x, y;
for(i=0;i<n;i++){
x = Rank[i];
if(t)t--;
if(x==n)continue;
y = SA[x+1];
while(p[i+t]==p[y+t])t++;
LCP[x] = t;
}
}
int top, R[101000];
struct TP{
int l, b;
}st[101000];
int main(){
int i, j;
scanf("%s",p);
for(i=0;p[i];i++);
n = i;
for(i=0;i<n/2;i++){
swap(p[i],p[n-i-1]);
}
Suffix_Array();
for(i=1;i<=n;i++){
int b = -1;
while(top && st[top].l >= LCP[i-1]){
b = st[top].b;
top--;
}
if(b!=-1){
st[++top] = {LCP[i-1],b};
}
int pv = i-1;
int r = 0, c = 0;
for(j=top;st[j].l;j--){
while(pv>=1 && st[j].b <= pv){
if(SA[pv] > SA[i])c++;
pv--;
}
r += (st[j].l-st[j-1].l) * c;
}
R[n-SA[i]] += r*2;
st[++top] = {n - SA[i], i};
}
top = 0;
for(i=n;i>=1;i--){
int b = -1;
while(top && st[top].l >= LCP[i]){
b = st[top].b;
top--;
}
if(b!=-1){
st[++top] = {LCP[i],b};
}
int pv = i+1;
int r = 0, c = 0;
for(j=top;st[j].l;j--){
while(pv<=n && st[j].b >= pv){
if(SA[pv] > SA[i])c++;
pv++;
}
r += (st[j].l-st[j-1].l) * c;
}
R[n-SA[i]] += r*2;
st[++top] = {n - SA[i], i};
}
long long sum = 0;
for(i=1;i<=n;i++){
sum = sum + i + R[i];
printf("%lld\n",sum);
}
}
Compilation message (stderr)
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |