题目链接:点击查看

题目大意:给出一个只含有 ' a ' 和 ' b ' 的字符串,再给出 B 数组的构造方法如下:对于每个位置 i 来说

  1. 如果存在一个位置 j ,使得 j < i 且 s[ j ] == s[ i ] ,则 b[ i ] = i - j
  2. 否则 s[ i ] = 0

现在需要对字符串的每个后缀构造 B 数组,并按照字典序排序

题目分析:参考博客:https://blog.csdn.net/weixin_43965698/article/details/107304525

一个例子:

通过观察,可以将每个后缀的 B 数组分成两段,分别记为段 A 和段 B(上图中表示为段 A 和段 D )

段 A 的特点是:第一个数字和最后一个数字同为 0,且段 A 中只有这两个 0,段 B 的特点是:是字符串 s 的 B 数组的其中一个后缀

因为只有段 A 中含有 0 ,所以在排序时,段 A 的优先级要高于段 B 的,所以段 A 长度更小的,显然应该排在前面

当段 A 的长度相等时,再考虑段 B ,因为此时的段 B 都是字符串 s 的 B 数组的后缀,所以可以利用后缀数组得到 rk 数组,rk[ i ] 当然就是后缀 i 的排名了

当然存在特殊情况,就是如果段 B 为空时,此时显然段 B 为空的字典序更小

代码:

#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=1e5+100;char str[N];int sa[N]; //SA数组,表示将S的n个后缀从小到大排序后把排好序的
//的后缀的开头位置顺次放入SA中
int t1[N],t2[N],c[N];int rk[N],height[N],len;int s[N];void build_sa(int s[],int n,int m)//n为添加0后的总长
{int i,j,p,*x=t1,*y=t2;for(i=0;i<m;i++) c[i]=0;for(i=0;i<n;i++) c[x[i]=s[i]]++;for(i=1;i<m;i++) c[i]+=c[i-1];for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;for(j=1;j<=n;j<<=1) {p=0;for(i=n-j;i<n;i++) y[p++]=i;for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;for(i=0;i<m;i++) c[i]=0;for(i=0;i<n;i++) c[x[y[i]]]++;for(i=1;i<m;i++) c[i]+=c[i-1];for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];swap(x,y);p=1,x[sa[0]]=0;for(i=1;i<n;i++) x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;if(p>=n) break;m=p;}
}void get_height(int s[],int n)//n为添加0后的总长
{int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;}
}void solve(int base=128)
{build_sa(s,len+1,base);get_height(s,len);
}struct Node
{int id,len_a,pos_b;//id:编号,len_a:段A的长度,pos_b:段B的起始位置 Node(int id,int len_a,int pos_b):id(id),len_a(len_a),pos_b(pos_b){}bool operator<(const Node& t)const{if(len_a!=t.len_a)//如果段A的长度不相等,较短的字典序更短 return len_a<t.len_a;if(pos_b!=-1&&t.pos_b!=-1)//如果都有段B,直接根据后缀数组比较大小 return rk[pos_b]<rk[t.pos_b];return pos_b==-1;//否则没有段B的显然更小 }
};int main()
{
#ifndef ONLINE_JUDGE
//  freopen("input.txt","r",stdin);
//  freopen("output.txt","w",stdout);
#endif
//  ios::sync_with_stdio(false);while(scanf("%d",&len)!=EOF){scanf("%s",str);int pre_a=-1,pre_b=-1;for(int i=0;i<len;i++)//根据题目要求初始化s数组 {if(str[i]=='a'){if(pre_a==-1)s[i]=1;elses[i]=i-pre_a+1;pre_a=i;}else{if(pre_b==-1)s[i]=1;elses[i]=i-pre_b+1;pre_b=i;}}s[len]=0;solve(len+1);pre_a=pre_b=-1;vector<Node>ans;for(int i=len-1;i>=0;i--){if(str[i]=='a'){if(pre_b==-1)//如果a之后没再有b的话,说明不存在段B,设置起初始位置为-1,设置段A的长度为len-i+1 ans.push_back(Node(i,len-i+1,-1));else//如果同时存在段A和段B,记录段A的长度和段B的初始位置 ans.push_back(Node(i,pre_b-i+1,pre_b+1));pre_a=i;}else{if(pre_a==-1)ans.push_back(Node(i,len-i+1,-1));elseans.push_back(Node(i,pre_a-i+1,pre_a+1));pre_b=i;}}sort(ans.begin(),ans.end());for(auto v:ans)printf("%d ",v.id+1);puts("");}return 0;
}

牛客多校 - B-Suffix Array(后缀数组)相关推荐

  1. 2019牛客多校第四场 I题 后缀自动机_后缀数组_求两个串de公共子串的种类数

    目录 求若干个串的公共子串个数相关变形题 对一个串建后缀自动机,另一个串在上面跑同时计数 广义后缀自动机 后缀数组 其他:POJ 3415 求两个串长度至少为k的公共子串数量 @(牛客多校第四场 I题 ...

  2. 2020 牛客多校第一场

    2020 牛客多校第一场 A. B-Suffix Array 后缀数组的思想:倍增+桶排序的方式找出一串连续序列后缀的大小.虽说正常使用的时候都是字典序,但是只要修改排序方式,也能够达到一个类似的&q ...

  3. python字符串去重及排序 牛客_2018牛客多校第一场 D.Two Graphs

    题意: n个点,m1条边的图E1,n个点,m2条边的图E2.求图E2有多少子图跟图E1同构. 题解: 用STL的全排列函数next_permutation()枚举映射.对于每一种映射枚举每一条边判断合 ...

  4. 牛客多校三 B Black and white

    牛客多校三 B Black and white 在n*m的棋盘上,每个格子有一个数,初始可以选一定的格子标记为黑色,在任意四个形如(i1, j1)(i1, j2)(i2, j1)(i2, j2)的格子 ...

  5. 2019牛客多校第一场

    2019牛客多校第一场 题号 题目 知识点 A Monotonic Matrix B Symmetric Matrix C Fluorescent 2 D Two Graphs E Removal F ...

  6. LCS(2021牛客多校4)

    LCS(2021牛客多校4) 题意: 让你构造三个字符串s1,s2,s3,长度均为n,要求LCS(s1,s2)=a,LCS(s2,s3)=b,LCS(s1,s3)=c 题解: 先考虑三个串互相LCS为 ...

  7. 24dian(牛客多校第三场)

    24dian(牛客多校第三场) 题意: 给你n张牌,每张牌的大小为1 ~ 13,问这些牌与加减乘除任意组合(可以使用括号),且但所有的有效解在计算过程中都涉及到分数,即非整数,能否组成答案m,如果可以 ...

  8. 2018年牛客多校算法寒假训练营练习比赛(第一场)C. 六子冲

    2018年牛客多校算法寒假训练营练习比赛(第一场)C. 六子冲 题目链接 做法:模拟即可 #include <bits/stdc++.h> #define P pair<int,in ...

  9. 2022牛客多校(十)

    2022牛客多校(十) 一.比赛小结 比赛链接:"蔚来杯"2022牛客暑期多校训练营10 二.题目分析及解法(基础题) F.Shannon Switching Game? 题目链接 ...

  10. 2019牛客多校训练第十场F Popping Balloons

    2019牛客多校训练第十场F Popping Balloons 题意:二维平面内给你若干个点,然后你可以在x轴和y轴分别射三枪(每一枪的间隔是R),问最多能射掉多少气球. 题解:贪心.这个应该只能算作 ...

最新文章

  1. 修改点击cell时显示的颜色
  2. 四十七、面试前,必须搞懂Java中的线程池ThreadPoolExecutor(上篇)
  3. WEKA The workbench for machine learning
  4. 海量数据的topK问题
  5. Mysql学习总结(38)——21条MySql性能优化经验
  6. 高性能javascript读书笔记(三.DOM 编程2)
  7. layUI数据表格(table)
  8. matlab 双音多频 接收端检测到的号码,实验一数字信号处理在双音多频拨号系统中的应用.doc...
  9. python入门之一运算符和表达式第二关:华氏度转摄氏度和进制转换
  10. python aks_使用环回aks和terraform构建基于打字稿的游戏后端
  11. C++游戏编程--基于EasyX图形库--实现走迷宫小游戏
  12. OSError: exception: access violation writing 0x000000003F800000
  13. 搭建GPU服务器(容器篇)
  14. java中级工程师所需的技能_Java 工程师要具备哪些技能或技术?
  15. 给编程一个你热爱它的机会
  16. 学习并掌握结构化写作方法,提高写作能力 ——结构化写作学习笔记(4)
  17. springboot集成mybatis
  18. 网管软件的几种主动监测方式
  19. oracle数据库测评
  20. win10平板读写TF卡异常的解决记录

热门文章

  1. 把linux制作成内存系统,把内存当硬盘,提速你的linux系统
  2. 探究Jvm源码实现-MarkWord
  3. AOP 中必须明白的概念-切面(Aspect)
  4. Redis 到底有多快?
  5. 微信支付 - 提供支付中心商户订单查询
  6. MybatisPlus实现自动填充
  7. 基于注解的方式配置bean
  8. sqoop导入-hdfs
  9. 函数的返回值-改造求和函数
  10. 建造者模式源码解析(jdk-guava+mybatis)