spoj694/705 Distinct Substrings - 后缀数组

2024-05-12 00:18

本文主要是介绍spoj694/705 Distinct Substrings - 后缀数组,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接:http://acm.hust.edu.cn/vjudge/problem/19282

题目大意:求不同子串的个数


解题思路:后缀数组..
suffix(i)对子串个数所做的贡献为len-sa[i]+1,因为要求要不同的,所以减去与它串重复的子串个数h[i]。即每个后缀对答案的贡献为len-sa[i]+1-h[i];(***)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
#define maxn 1100int sa[maxn],rk[maxn],y[maxn],len;
int Rsort[maxn],wr[maxn],h[maxn];
char s[maxn];
bool cmp(int x,int y,int ln){return (wr[x]==wr[y])&&(wr[x+ln]==wr[y+ln]);}
void get_h()
{int i,k=0;for (i=1;i<=len;i++){int j=sa[rk[i]-1];if (k>0) k--;while (s[i+k]==s[j+k]) k++;h[rk[i]]=k;}
}
void get_sa()
{int i,k,p,ln,m=130;for (i=1;i<=len;i++) rk[i]=s[i];for (i=0;i<=m;i++) Rsort[i]=0;for (i=1;i<=len;i++) Rsort[rk[i]]++;for (i=1;i<=m;i++) Rsort[i]+=Rsort[i-1];for (i=len;i>=1;i--) sa[Rsort[rk[i]]--]=i;p=0;ln=1;while (p<len){for (k=0,i=len-ln+1;i<=len;i++) y[++k]=i;for (i=1;i<=len;i++) if (sa[i]-ln>0) y[++k]=sa[i]-ln;for (i=1;i<=len;i++) wr[i]=rk[y[i]];for (i=0;i<=m;i++) Rsort[i]=0;for (i=1;i<=len;i++) Rsort[wr[i]]++;for (i=1;i<=m;i++) Rsort[i]+=Rsort[i-1];for (i=len;i>=1;i--) sa[Rsort[wr[i]]--]=y[i];memcpy(wr,rk,sizeof(wr));p=1;rk[sa[1]]=1;for (i=2;i<=len;i++){if (!cmp(sa[i],sa[i-1],ln)) p++;rk[sa[i]]=p;}m=p;ln*=2;}sa[0]=s[0]=0;
}
int main()
{int T,i;scanf("%d\n",&T);while (T--){scanf("%s",s+1);len=strlen(s+1);get_sa();get_h();LL ans=0;h[1]=0;for (i=1;i<=len;i++){LL ls=sa[i];ans+=len-ls+1-h[i];}printf("%lld\n",ans);}return 0;
}


这篇关于spoj694/705 Distinct Substrings - 后缀数组的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/981067

相关文章

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

C语言:柔性数组

数组定义 柔性数组 err int arr[0] = {0}; // ERROR 柔性数组 // 常见struct Test{int len;char arr[1024];} // 柔性数组struct Test{int len;char arr[0];}struct Test *t;t = malloc(sizeof(Test) + 11);strcpy(t->arr,

C 语言基础之数组

文章目录 什么是数组数组变量的声明多维数组 什么是数组 数组,顾名思义,就是一组数。 假如班上有 30 个同学,让你编程统计每个人的分数,求最高分、最低分、平均分等。如果不知道数组,你只能这样写代码: int ZhangSan_score = 95;int LiSi_score = 90;......int LiuDong_score = 100;int Zhou

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

C# double[] 和Matlab数组MWArray[]转换

C# double[] 转换成MWArray[], 直接赋值就行             MWNumericArray[] ma = new MWNumericArray[4];             double[] dT = new double[] { 0 };             double[] dT1 = new double[] { 0,2 };

PHP7扩展开发之数组处理

前言 这次,我们将演示如何在PHP扩展中如何对数组进行处理。要实现的PHP代码如下: <?phpfunction array_concat ($arr, $prefix) {foreach($arr as $key => $val) {if (isset($prefix[$key]) && is_string($val) && is_string($prefix[$key])) {$arr[

Go 数组赋值问题

package mainimport "fmt"type Student struct {Name stringAge int}func main() {data := make(map[string]*Student)list := []Student{{Name:"a",Age:1},{Name:"b",Age:2},{Name:"c",Age:3},}// 错误 都指向了最后一个v// a

码蹄集部分题目(2024OJ赛9.4-9.8;线段树+树状数组)

1🐋🐋配对最小值(王者;树状数组) 时间限制:1秒 占用内存:64M 🐟题目思路 MT3065 配对最小值_哔哩哔哩_bilibili 🐟代码 #include<bits/stdc++.h> using namespace std;const int N=1e5+7;int a[N],b[N],c[N],n,q;struct QUERY{int l,r,id;}que

【第0006页 · 数组】寻找重复数

【前言】本文以及之后的一些题解都会陆续整理到目录中,若想了解全部题解整理,请看这里: 第0006页 · 寻找重复数         今天想讨论的一道题在 LeetCode 上评论也是颇为“不错”。有一说一,是道好题,不过我们还是得先理解了它才算真正的好题。这里我们展示一种使用二进制的做法,希望能帮到你哟! 【寻找重复数】给定一个包含 n + 1 个整数的数组 nums ,其数字都