正题

题目链接:https://www.luogu.com.cn/problem/P4070

题目大意

长度为nnn的字符串，对于每个iii求字符串1∼i1\sim i1∼i部分有多少个不同的子串。

解题思路

对于整个串ans=∑i=1nn−i+1−heightians=\sum_{i=1}^nn-i+1-height_ians=∑i=1nn−i+1−heighti，考虑如何维护heightheightheight
往末尾加上一个字符会对heightheightheight造成很大影响，但是如果在字符串前面加上一个heightheightheight的话就不会对其他heightheightheight造成影响了。

所以先翻转字符串，问题变为每次在前面加上一个字符。

首先离线后缀排好序，然后每次加入一个字符时，我们用平衡树找已经插入的前驱和后继然后用STSTST表来求到他们的LCPLCPLCP（也就是更新后的heightheightheight）即可。

codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
using namespace std;
const int N=1e5+10;
int n,m,a[N],b[N],height[N],lg[N];
int c[N],sa[N],rk[N],x[N],y[N];
int f[20][N];long long ans;
set<int> s;
void Qsort(){for(int i=1;i<=m;i++)c[i]=0;for(int i=1;i<=n;i++)c[x[i]]++;for(int i=1;i<=m;i++)c[i]+=c[i-1];for(int i=n;i>=1;i--)sa[c[x[y[i]]]--]=y[i],y[i]=0;return;
}
void SA(){for(int i=1;i<=n;i++)x[i]=a[i],y[i]=i;Qsort();for(int w=1;w<=n;w<<=1){int p=0;for(int i=n-w+1;i<=n;i++)y[++p]=i;for(int i=1;i<=n;i++)if(sa[i]>w)y[++p]=sa[i]-w;Qsort();swap(x,y);p=x[sa[1]]=1;for(int i=2;i<=n;i++)x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+w]==y[sa[i-1]+w])?p:++p;if(p==n)break;m=p;}return;
}
void Get_Height(){int k=0;for(int i=1;i<=n;i++)rk[sa[i]]=i;for(int i=1;i<=n;i++){if(rk[i]==1)continue;if(k)k--;int j=sa[rk[i]-1];while(i+k<=n&&j+k<=n&&a[i+k]==a[j+k])k++;height[rk[i]]=k;}return;
}
void RMQ(){lg[0]=-1;for(int i=1;i<=n;i++)lg[i]=lg[i/2]+1,f[0][i]=height[i];for(int i=1;i<19;i++)for(int j=1;j+(1<<i-1)<=n;j++)f[i][j]=min(f[i-1][j],f[i-1][j+(1<<i-1)]);return;
}
int LCP(int l,int r){l++;int z=lg[r-l+1];return min(f[z][l],f[z][r+1-(1<<z)]);
}
int main()
{scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%d",&a[i]);b[i]=a[i];}sort(b+1,b+1+n);m=unique(b+1,b+1+n)-b-1;for(int i=1;i<=n;i++)a[i]=lower_bound(b+1,b+1+m,a[i])-b;for(int i=1;i<=n/2;i++)swap(a[i],a[n-i+1]);SA();Get_Height();RMQ();for(int i=n;i>=1;i--){s.insert(rk[i]);set<int>::iterator it=s.find(rk[i]);int k=0;if(it!=s.begin()){int p=*(--it);k=LCP(p,rk[i]);++it;}++it;if(it!=s.end()){int p=*it;k=max(k,LCP(rk[i],p));}ans+=n-i+1-k;printf("%lld\n",ans);}
}

P4070-[SDOI2016]生成魔咒【SA,平衡树】相关推荐

P4070 [SDOI2016]生成魔咒
P4070 [SDOI2016]生成魔咒题意: 有n个字符xi,每次在S的末尾加入一个字符,(一开始S为空),每次加入xi后的不相同字串有多少个题解: 做这个题首先要会后缀数组P3809 [模板] ...
洛谷 P4070 [SDOI2016]生成魔咒解题报告
P4070 [SDOI2016]生成魔咒题目描述魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 $1$.$2$ 拼凑起来形成一个魔咒串 $[1,2]$. 一个魔咒 ...
bzoj4516 / P4070 [SDOI2016]生成魔咒
P4070 [SDOI2016]生成魔咒后缀自动机每插入一个字符,对答案的贡献为$len[last]-len[fa[last]]$ 插入字符范围过大,所以使用$map$存储. (去掉第35行就是裸 ...
P4070 [SDOI2016]生成魔咒（SAM len数组的含义）
[SDOI2016]生成魔咒题目描述魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1 , 2 1,2 1,2 拼凑起来形成一个魔咒串 [ 1 , 2 ] [1,2] [1, ...
洛谷 P4070 [SDOI2016]生成魔咒后缀自动机
题目描述魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔咒串 [1,2]. 一个魔咒串 S 的非空字串被称为魔咒串 S 的生成魔咒. 例如 S=[1, ...
P4070 [SDOI2016]生成魔咒(SAM)
传送门用后缀数组做的真的心累后缀数组+ST表+反串思维+set(或平衡树)动态维护 h e i g h t height height数组然后用 S A M SAM SAM就是模板题- 要求本质 ...
luogu P4070 [SDOI2016]生成魔咒
首先有两个结论: 1.后缀自动机具有最简性,即每种不同的子串只会在sam上体现一次,体现形式是sam上一条由root出发的路径. 2.一个字符串不同子串的个数等于所有关键节点的max(x)-min(x ...
BZOJ4516: [Sdoi2016]生成魔咒
BZOJ4516: [Sdoi2016]生成魔咒 Description 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示. 例如可以将魔咒字符 1.2 拼凑起来形成一个魔咒串 [1,2]. 一个魔咒 ...
[SDOI2016]生成魔咒
4516: [Sdoi2016]生成魔咒 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 1460 Solved: 835 [Submit][Sta ...
BZOJ 4516: [Sdoi2016]生成魔咒 [后缀自动机]
4516: [Sdoi2016]生成魔咒题意:询问一个字符串每个前缀有多少不同的子串做了一下SDOI2016R1D2,题好水啊随便AK 强行开map上SAM 每个状态的贡献就是\(Max(s)-M ...

P4070-[SDOI2016]生成魔咒【SA,平衡树】

正题

题目大意

解题思路

codecodecode

P4070-[SDOI2016]生成魔咒【SA,平衡树】相关推荐

最新文章

热门文章