写了一半 没了啊啊啊     重新写的
思路:
先不考虑后缀自动机   (我不会啊)

那这道题只能用后缀数组了
先把原串倒一下  后缀->前缀
相当于每回在前面加了一个字母
求不同的子串个数 

首先 正常的求子串个数我们是会的 SPOJ 705
但是这道题比较坑 它让你每回都输出一下
那只好 维护一个前驱 一个后继 求LCP 取max
ans=ans+n-i+1-max(LCP(),LCP())
用set维护前驱和后继就好啦~

//By SiriusRen
#include <set>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=100050;
int n,s[N],sa[N],tsa[N],cntA[N],cntB[N],A[N],B[N],rk[N],ht[N],f[N][20],Log[N];
set<int>SET;long long ans;
void SA(){
    for(int i=1;i<=n;i++)cntA[s[i]]++;
    for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1];
    for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
    rk[sa[1]]=1;
    for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
    for(int l=1;rk[sa[n]]<n;l<<=1){        memset(cntA,0,sizeof(cntA));
        memset(cntB,0,sizeof(cntB));
        for(int i=1;i<=n;i++)
            cntA[A[i]=rk[i]]++,
            cntB[B[i]=(i+l<=n?rk[i+l]:0)]++;
        for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];
        for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
        for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
        rk[sa[1]]=1;
        for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);
    }
    for(int i=1,j=0;i<=n;i++){        j=j?j-1:0;
        while(s[i+j]==s[sa[rk[i]-1]+j])j++;
        f[rk[i]][0]=j;
    }
    for(int j=1;j<=19;j++)
        for(int i=1;i+(1<<(j-1))<=n;i++)
            f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
int LCP(int x,int y){
    int t=Log[y-x+1];
    return min(f[x][t],f[y-(1<<t)+1][t]);
}
int main(){
    scanf("%d",&n),Log[0]=-1;
    for(int i=1;i<=n;i++)Log[i]=Log[i>>1]+1;
    for(int i=1;i<=n;i++)scanf("%d",&s[i]),sa[i]=s[i];
    sort(sa+1,sa+1+n);int u=unique(sa+1,sa+1+n)-sa-1;
    for(int i=1;i<=n;i++)s[i]=lower_bound(sa+1,sa+1+u,s[i])-sa;
    for(int i=1;i<=n/2;i++)swap(s[i],s[n-i+1]);
    SA();
    for(int i=n;i;i--){        int jy=0;
        set<int>::iterator it=SET.upper_bound(rk[i]);
        if(it!=SET.end())jy=LCP(rk[i]+1,*it);
        if(it!=SET.begin())jy=max(jy,LCP((*(--it))+1,rk[i]));
        ans=ans+n-i+1-jy,SET.insert(rk[i]);
        printf("%lld\n",ans);
    }
}

转载于:https://www.cnblogs.com/SiriusRen/p/6532049.html

BZOJ 4516 后缀数组+ST+set相关推荐

  1. 后缀数组 ---- 2018~2019icpc焦作H题[后缀数组+st表+二分+单调栈]

    题目链接 题目大意: 给出nnn个数,定义f[l,r]f[l,r]f[l,r]表示 区间[l,r][l,r][l,r]的最大值,求所有 子区间的最大值的和,要求相同的子区间只能算一次 比如数列 5 6 ...

  2. 【BZOJ4310】跳蚤,后缀数组+ST表求LCP+二分答案

    Time:2016.05.26 Author:xiaoyimi 转载注明出处谢谢 传送门 思路: 首先要求出不同子串的个数 有这样一个性质 一个串中不同子串的总数=∑(len-height[i]-sa ...

  3. [luoguP2463] [SDOI2008]Sandy的卡片(后缀数组 + st表)

    传送门 很容易想到,题目中的相同是指差分数组相同. 那么可以把差分数组连起来,中间加上一个没有出现过的且字典序小的数 双指针移动,用st表维护height数组中的最小值. 当然用单调队列应该也可以且更 ...

  4. NOI 2016 优秀的拆分 (后缀数组+差分)

    题目大意:给你一个字符串,求所有子串的所有优秀拆分总和,优秀的拆分被定义为一个字符串可以被拆分成4个子串,形如$AABB$,其中$AA$相同,$BB$相同,$AB$也可以相同 作为一道国赛题,95分竟 ...

  5. 【BZOJ4566】找相同字符(后缀数组)

    题面 BZOJ 题解 后缀数组的做法,应该不是很难想 首先看到两个不同的串,当然是接在一起求 SA,height SA,height 那么,考虑一下暴力 在两个串各枚举一个后缀,他们的 lcp lcp ...

  6. UOJ #219 BZOJ 4650 luogu P1117 [NOI2016]优秀的拆分 (后缀数组、ST表)

    UOJ #219 BZOJ 4650 luogu P1117 [NOI2016]优秀的拆分 (后缀数组.ST表) 连NOI Day1T1都不会做...看了题解都写不出来还要抄Claris的代码.. 题 ...

  7. [BZOJ 3238] [AHOI 2013] 差异 【后缀数组 + 单调栈】

    题目链接:BZOJ - 3238 题目分析 显然,这道题就是求任意两个后缀之间的LCP的和,这与后缀数组的联系十分明显. 求出后缀数组后,求出字典序相邻两个后缀的LCP,即 Height 数组. 那么 ...

  8. BZOJ3473:字符串(后缀数组,主席树,二分,ST表)

    Description 给定n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串? Input 第一行两个整数n,k. 接下来n行每行一个字符串. Output 一 ...

  9. 【BZOJ】1031: [JSOI2007]字符加密Cipher(后缀数组)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1031 很容易想到这就是将字符串复制到自己末尾然后后缀数组搞出sa然后按区间输出即可. 然后换了下模板 ...

  10. bzoj 1031 [JSOI2007]字符加密Cipher 后缀数组

    题面 题目传送门 解法 后缀数组模板题吧-- 将字符串两倍,然后求一遍sa数组即可 时间复杂度:\(O(n\ log\ n)\) 代码 #include <bits/stdc++.h> # ...

最新文章

  1. 关于Action返回结果类型的事儿(下)
  2. 小鹏发布会上硬刚特斯拉AutoPilot:有视频有真相,特斯拉不懂中国路况
  3. boost::multi_array模块调整 multi_arrays 大小的测试
  4. python按hash分组_Python操作redis系列以 哈希(Hash)命令详解(四)
  5. mysqli模块的安装
  6. hdfs dfs -cat 出现No such file or directory
  7. [透析] 卷积神经网络CNN究竟是怎样一步一步工作的?
  8. 详解MySQL中EXPLAIN解释命令(转)
  9. 【python绘图】Matplotlib绘图及设置(使用python制图)
  10. 60.(leaflet之家)leaflet虚线
  11. 为什么会出现锟斤拷?
  12. adb模拟按键home_adb命令中 模拟按键 KeyCode 部分详解
  13. ps2键盘测试软件,51单片机驱动PS2键盘完整程序
  14. Cordova徽章插件
  15. 30K~65K,春节前最后一批热门技术岗位,快到碗里来
  16. 工作中如何应对那些不讲武德的人?
  17. [Python爬虫] 5-爬虫进阶(多线程爬虫/动态网页抓取/图形验证码识别)
  18. ASM 系列详细教程-02-ASM 整体概览
  19. 计算机exce常用功能,电脑办公实用小技巧
  20. 解决badboy录制时弹框提示“当前页面的脚本发生错误”

热门文章

  1. DML、DDL的概念与区别
  2. PyTorch:Encoder-RNN|LSTM|GRU
  3. flutter弹框怎么封装_Flutter初始化时 弹出弹出框
  4. C++ char,int 之间的相互转换
  5. Git 常见问题总结
  6. linux创建线程pthread_create函数
  7. 用正则表达式生成文章的目录大纲
  8. 服务站: WCF 消息传递基础 -- MSDN Magazine, April 2007
  9. centos上nginx转发tcp请求
  10. P3194 [HNOI2008]水平可见直线