bzoj 3238 Ahoi2013 后缀数组+单调栈
来源:互联网 发布:工资数据分析 编辑:程序博客网 时间:2024/06/08 15:57
题目链接
思路:先不看lcp部分,若只有len[i] 则
然后就变成了一个跑后缀数组得到height值+RMQ的题目
等等……枚举的时间复杂度为O(n^2)明显超时,那么我们需要在O(n)或者O(nlog)的时间内处理出区间的最小值,注意我们并不需要知道具体的每一个区间是多少,只用知道总值就可以了~~
我们用单调栈来维护这个问题,对于任意一个i答案为向左向右拓展长度的乘积
#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>#include<stack>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fod(i,a,b) for(int i=a;i>=b;i--)using namespace std;const int N=5e5+10;typedef long long ll;int n,p,q,k;char ch[N];int a[N],v[N],h[N];int sa[2][N],rk[2][N]; ll dp[N];void calsa(int sa[N],int rk[N],int SA[N],int RK[N]){ for(int i=1;i<=n;i++)v[rk[sa[i]]]=i; for(int i=n;i;i--) if(sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k; for(int i=n-k+1;i<=n;i++)SA[v[rk[i]]--]=i; for(int i=1;i<=n;i++) RK[SA[i]]=RK[SA[i-1]]+(rk[SA[i-1]]!=rk[SA[i]]||rk[SA[i]+k]!=rk[SA[i-1]+k]);}void getsa(){ memset(v,0,sizeof(v)); p=0,q=1; for(int i=1;i<=n;i++)v[a[i]]++; for(int i=1;i<=30;i++)v[i]+=v[i-1]; for(int i=1;i<=n;i++) sa[p][v[a[i]]--]=i; for(int i=1;i<=n;i++) rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(a[sa[p][i-1]]!=a[sa[p][i]]); for(k=1;k<n;k<<=1,swap(p,q)) calsa(sa[p],rk[p],sa[q],rk[q]);}void geth(){ k=0; for(int i=1;i<=n;i++) if(rk[p][i]==1)h[rk[p][i]]=0; else { int j=sa[p][rk[p][i]-1]; while(a[i+k]==a[j+k])k++; h[rk[p][i]]=k;if(k>0)k--; }}stack<int>s;int main(){ freopen("b.in","r",stdin); freopen("b.out","w",stdout); scanf("%s",ch+1); n=strlen(ch+1); fo(i,1,n) a[i]=ch[i]-'a'+1; getsa(); geth(); ll ans=ll(n+1)*(n-1)*ll(n)>>1; for(int i=n;i;i--) { while(!s.empty()&&h[s.top()]>h[i])s.pop(); int tmp; if(!s.empty()) tmp=s.top(); else tmp=n+1; dp[i]=(ll)h[i]*(tmp-i)+dp[tmp]; ans-=dp[i]*2; s.push(i); } cout<<ans<<endl; return 0;}
阅读全文
0 0
- bzoj 3238 Ahoi2013 后缀数组+单调栈
- BZOJ 3238 [Ahoi2013]差异 后缀数组+单调栈
- bzoj 3238: [Ahoi2013]差异(后缀数组+单调栈)
- 【bzoj3238】[Ahoi2013]差异 后缀数组+单调栈
- 【bzoj3238】[Ahoi2013]差异 后缀数组+单调栈
- 【bzoj3238】【AHOI2013】【差异】【后缀数组+单调栈】
- 后缀数组+单调栈 【Ahoi2013】bzoj3238 差异
- bzoj 3238: [Ahoi2013]差异 后缀数组
- BZOJ 3238: [Ahoi2013]差异 后缀数组
- bzoj 3238 ahoi2013差异 后缀数组
- BZOJ 3238 后缀数组+单调栈
- BZOJ 3238: [Ahoi2013]差异|后缀数组|乘法原理
- bzoj 3238 [Ahoi2013]差异 后缀数组 并查集
- BZOJ 3238 [Ahoi2013]差异 后缀数组+st表
- [BZOJ3238][Ahoi2013]差异(后缀数组+单调栈||后缀自动机+树形dp)
- BZOJ 3238 AHOI 2013 差异 后缀数组+单调栈
- [BZOJ 3238][AHOI 2013]差异(后缀数组+单调栈)
- BZOJ 3238 AHOI2013 差异 后缀自动机
- js中的假值及==、===
- 0-1背包问题
- 2017.08.15【NOIP提高组】模拟赛B组 生日聚餐
- Leetcode-Longest Common Prefix-Python
- 谁说90后不靠谱?95后都做CEO了
- bzoj 3238 Ahoi2013 后缀数组+单调栈
- HDU 6130 Kolakoski 打表
- 2017年8月12日(模拟四(smoj2193,2194,2195)强连通分量、勾股定理探索、动态规划+容斥原理)
- 常用正则表达式大全,手机、电话、邮箱、身份证(最严格的验证)、IP地址、网址、日期等,一般前台js验证
- 【已解决】ubuntu插网线无法联网
- android相对布局中@id和@+id的区别(原理)
- 关于ETCD和Zookeeper的一些简单了解(转)
- Linux下安装Elasticsearch5.5
- Jmeter运行原理