UOJ.35.[模板]后缀排序(后缀数组 倍增)

题目链接

论找到一个好的教程的正确性。。
后缀数组

下标从1编号:

//299ms 2560kb
#include <cstdio>
#include <cstring>
#include <algorithm>
const int N=1e5+5;int n,sa[N],rk[N],sa2[N],tm[N],ht[N];
char s[N];void Get_SA()
{int *x=rk,*y=sa2,m=200;for(int i=0; i<=m; ++i) tm[i]=0;for(int i=1; i<=n; ++i) ++tm[x[i]=s[i]-'a'+1];for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];for(int i=n; i; --i) sa[tm[x[i]]--]=i;for(int p=0,k=1; k<n; k<<=1,m=p,p=0){for(int i=n-k+1; i<=n; ++i) y[++p]=i;for(int i=1; i<=n; ++i) if(sa[i]>k) y[++p]=sa[i]-k;for(int i=0; i<=m; ++i) tm[i]=0;for(int i=1; i<=n; ++i) ++tm[x[i]];for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];for(int i=n; i; --i) sa[tm[x[y[i]]]--]=y[i];std::swap(x,y), p=1, x[sa[1]]=1;for(int i=2; i<=n; ++i)
//          if(y[sa[i-1]]==y[sa[i]]&&((sa[i-1]+k<=n&&sa[i]+k<=n&&y[sa[i-1]+k]==y[sa[i]+k])||(sa[i-1]+k>n&&sa[i]+k>n))) x[sa[i]]=p;
//          if(y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]) x[sa[i]]=p;
//          else x[sa[i]]=++p;x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p:++p;if(p>=n) break;}
}
void Calc_ht()
{for(int i=1; i<=n; ++i) rk[sa[i]]=i;ht[1]=0;for(int k=0,p,i=1; i<=n; ++i)//ht[rk[i]]>=ht[rk[i-1]]-1{if(rk[i]==1) continue;//排名为0的 ht[rk[i]]=ht[0]=0 if(k) --k;p=sa[rk[i]-1];while(i+k<=n&&p+k<=n&&s[i+k]==s[p+k]) ++k;ht[rk[i]]=k;}
}int main()
{scanf("%s",s+1), n=strlen(s+1);Get_SA();for(int i=1; i<=n; ++i) printf("%d ",sa[i]);putchar('\n'), Calc_ht();for(int i=2; i<=n; ++i) printf("%d ",ht[i]);return 0;
}
//321ms 2552kb
#include <cstdio>
#include <cstring>
#include <algorithm>
const int N=1e5+5;int n,sa[N],rk[N],sa2[N],tm[N],ht[N];
char s[N];void Get_SA()
{int *x=rk,*y=sa2,m=28;for(int i=0; i<m; ++i) tm[i]=0;for(int i=0; i<n; ++i) ++tm[x[i]=s[i]-'a'];for(int i=1; i<m; ++i) tm[i]+=tm[i-1];for(int i=n-1; ~i; --i) sa[--tm[x[i]]]=i;for(int p=0,k=1; k<n; k<<=1,/*m=p,*/m=p+1,p=0){for(int i=n-k; i<n; ++i) y[p++]=i;for(int i=0; i<n; ++i) if(sa[i]>=k) y[p++]=sa[i]-k;for(int i=0; i<m; ++i) tm[i]=0;for(int i=0; i<n; ++i) ++tm[x[i]];for(int i=1; i<m; ++i) tm[i]+=tm[i-1];for(int i=n-1; ~i; --i) sa[--tm[x[y[i]]]]=y[i];//      std::swap(x,y), p=1, x[sa[0]]=0;
//      for(int i=1; i<n; ++i)
//          if(y[sa[i-1]]==y[sa[i]]&&((sa[i-1]+k<n&&sa[i]+k<n&&y[sa[i-1]+k]==y[sa[i]+k])||(sa[i-1]+k>=n&&sa[i]+k>=n))) x[sa[i]]=p-1;
//          else x[sa[i]]=p++;
//      if(p>=n) break;std::swap(x,y), p=0, x[sa[0]]=0;for(int i=1; i<n; ++i)if(y[sa[i-1]]==y[sa[i]]&&((sa[i-1]+k<n&&sa[i]+k<n&&y[sa[i-1]+k]==y[sa[i]+k])||(sa[i-1]+k>=n&&sa[i]+k>=n))) x[sa[i]]=p;else x[sa[i]]=++p;if(p+1>=n) break;}
}
void Calc_ht()
{for(int i=0; i<n; ++i) rk[sa[i]]=i;ht[0]=0;for(int k=0,p,i=0; i<n; ++i)//ht[rk[i]]>=ht[rk[i-1]]-1{if(!rk[i]) continue;//排名为0的 ht[rk[i]]=ht[0]=0 if(k) --k;p=sa[rk[i]-1];while(i+k<n&&p+k<n&&s[i+k]==s[p+k]) ++k;ht[rk[i]]=k;}
}int main()
{scanf("%s",s), n=strlen(s);Get_SA();for(int i=0; i<n; ++i) printf("%d ",sa[i]+1);putchar('\n'), Calc_ht();for(int i=1; i<n; ++i) printf("%d ",ht[i]);return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/8567563.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/427978.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

linux 使cpu使用率升高_关于linux系统CPU篇---gt;CPU使用率升高

1.CPU使用率为单位时间内CPU使用情况的统计&#xff0c;以百分比的方式展示。LINUX作为一个多任务操作系统&#xff0c;将每个CPU的时间划分为很短的时间片&#xff0c;再通过调度器轮流分配给各个任务使用&#xff0c;因此造成多任务同时运行的错觉2.如何查看CPU使用率&#x…

软件测试之逻辑覆盖测试理论总结(白话文)

1.语句覆盖 语句覆盖就是所有可执行的语句都可以得到一次执行。注意可执行那三个字就可以&#xff0c;因为可能有的判定条件比较狗&#xff0c;它下面的分支的语句无论如何都不会执行。 例如这随手就能写出一个: if(x>100&&x<100) {x100; }当然实际上一般不会有…

datatables 的导出button自定义

1.dom 方式 $(#myTable).DataTable( {dom: Bfrtip,buttons: [copy, excel, pdf] } ); 2.手动插入 var table $(#example).DataTable( {buttons: [copy, excel, pdf] } );table.buttons().container().appendTo( $(.col-sm-6:eq(0), table.table().container() ) ); $(#myTable…

excel显著性检验_#如何用excel做anova分析#用excel做显著性分析

如何用Excel做方差分析&#xff1f;&#xff1f;方差分析的步骤&#xff1a;(1)分别计算行与列平方和。行平方与组差相似&#xff0c;是每值与总体均值的离差平方和&#xff0c;列平方和是每列的均值与总体均值的离差平方和。(2)总的平方和的计算与单因素方差分析一样&#xff…

软件测试之控制流图以及环形复杂度独立路径求解问题

首先需要明确的是&#xff0c;控制流图并不等于流程图&#xff0c;可以理解为控制流图的出现是为了后续的环形复杂度的计算和写出独立路径和配以相应的测试用例。 所以控制流图是核心&#xff0c;画图的时候务必谨慎再谨慎&#xff0c;要不然可能你后面的全部崩盘。 控制流图考…

MS Code 使用 TFVC 插件时遇到的问题

使用 TFVC 插件&#xff0c; 在 singin 时提示 &#xff1a; “It appears you have configured a non-English version of the TF executable. Please ensure an English version is properly configured.” 这是 Code 的问题&#xff0c; 官方解释&#xff1a; https://gith…

server的自增主键返回函数 sql_SQL自增主键函数

自动生成编码的主键函数比如CRM00001CRM00002CRM00003create table IntKey(KeyChar char(10))gocreate function GetKey()returns char(10)asbegindeclare KeyValue intdeclare KeyReturn varchar(20)set KeyValue cast(isnull((select max(KeyChar) from IntKey),0) as int) …

软件测试之黑盒测试-等价类划分法

首先&#xff0c;明确等价类分为有效等价类和无效等价类两种&#xff0c;一般无效等价类都是根据有效等价类写出来的。为了严谨&#xff0c;在白话文之前还是附个定义吧&#xff0c;以免误人子弟。 说白了&#xff0c;有效等价类就是有效输入数据的集合&#xff0c;无效等价类…

PM2.5环境检测系统的设计与分析

PM2.5环境检测系统的设计与分析 摘要&#xff1a; 大气颗粒物污染对人类健康和生态环境造成了很大的影响&#xff0c;这让人们逐渐重视起对细颗粒物PM2.5检测技术的研究。本文阐述了PM2.5浓度检测的五种方法&#xff0c;在对上述各方法分析总结的基础上针对日常生活中PM2.5污染…

关于C和C++中for循环对于中间逗号的判定

for&#xff08;&#xff1b;A,B,;&#xff09; 虽然学了好久C语言了&#xff0c;但是自己还真是一直没有注意过for循环中间有逗号表达式的时候判定&#xff0c;之前一直都是以为A和B同时满足了才会进行执行。学编译原理&#xff0c;深入理解一些语言的结构的时候才发现了这一点…

c语言中的取模运算符_C语言除法算法和取模运算的实现(多种算法,多种思路)...

对计算机来说&#xff0c;除法与求模是整数算术运算中最复杂的运算。相对其他运算(如加法与减法)来说&#xff0c;这两种算法的执行速度非常慢。例如&#xff0c;ARM 硬件上不支持除法指令&#xff0c;编译器调用 C 库函数来实现除法运算。直接利用 C 库函数中的标准整数除法程…

【CF603E】Pastoral Oddities cdq分治+并查集

【CF603E】Pastoral Oddities 题意&#xff1a;有n个点&#xff0c;依次加入m条边权为$l_i$的无向边&#xff0c;每次加入后询问&#xff1a;当前图是否存在一个生成子图&#xff0c;满足所有点的度数都是奇数。如果有&#xff0c;输出这个生成子图中边权最大的边的权值最小可能…

关于C/C++中函数参数传递的规则

f&#xff08;A,B,C&#xff09;&#xff1b; 在编译程序中&#xff0c;对函数参数传递的处理是由右向左进行的。 因此&#xff0c;先传递的是最右边的参数的值。这个结论在日常使用的时候其实也很难发现有什么具体的区别&#xff0c;但是在编译原理中或者一些偏向于理论的考试…

sql如何遍历几百万的表_SQL Server遍历表中记录的2种方法(使用表变量和游标)

SQL Server遍历表一般都要用到游标&#xff0c;SQL Server中可以很容易的用游标实现循环&#xff0c;实现SQL Server遍历表中记录。本文将介绍利用使用表变量和游标实现数据库中表的遍历。表变量来实现表的遍历以下代码中&#xff0c;代码块之间的差异已经用灰色的背景标记。DE…

SQL的主键和外键

SQL的主键和外键的作用 SQL的主键和外键的作用&#xff1a; 外键取值规则&#xff1a;空值或参照的主键值。 (1)插入非空值时&#xff0c;如果主键表中没有这个值&#xff0c;则不能插入。 (2)更新时&#xff0c;不能改为主键表中没有的值。 (3)删除主键表记录时&#xff0c;你…

win 11bitlocker恢复密匙一般情况的解决方式(这里指的是你现在使用的微软账户一直没有变更过)

本来没有打算写这一篇解决帖子的&#xff0c;但是最近好多微博的友友都在私我询问解决方法&#xff0c;孩子虽然很热心&#xff0c;但是真的回复不过来了&#xff0c;打字太累了&#xff0c;这里给大家简单指个路&#xff0c;希望能够对大家有所帮助。 当时出现这个问题的时候…

java类验证和装载顺序_java中类的加载顺序介绍(ClassLoader)

转自&#xff1a;http://blog.csdn.net/eff666/article/details/522034061、ClassNotFoundExcetpion我们在开发中&#xff0c;经常可以遇见java.lang.ClassNotFoundExcetpion这个异常&#xff0c;今天我就来总结一下这个问题。对于这个异常&#xff0c;它实质涉及到了java技术体…

[转]ASP.NET MVC4中@model使用多个类型实例的方法

本文转自&#xff1a;http://blog.csdn.net/hulihui/article/details/48199897 有时需要在ASP.NET MVC4的视图的model中使用多个类型的实例&#xff0c;.NET Framework 4.0版本引入的System.Tuple类可以轻松满足这个需求。 假设Person和Product是两个类型&#xff0c;如下是控制…

软件测试之黑盒测试-边界值分析法(理论白话学习/期中期末备考)

这里简单进行一下知识体系的梳理。 文章目录前言一、边界值分析是什么&#xff1f;二、因果图1.分析规格说明中给出的原因&#xff0c;结果和是否可能存在中间态。2.标约束&#xff0c;画图三、画判定表并画出对应的测试用例表。总结前言 主要从三个方面谈一下简单应用 边界值…

window下搭建Python3.7+selenium3.1.1+pycharm环境

1.安装Python3.7 1.1 下载 Python并安装 Python3.5 &#xff08;勾选上 Add Python3.7 to PATH&#xff09; 点击 Install Now&#xff0c;安装完成后将python路径加入环境变量Path 如D:\python35\Python35;D:\python35\Python35\Scripts; 1.2 安装完成之后&#xff0c;winr…