对于最长连续重复子串伤感经历

最新推荐文章于 2023-03-19 19:58:43 发布
最新推荐文章于 2023-03-19 19:58:43 发布
阅读量802 收藏
点赞数
分类专栏： 字符串 字符串—后缀数组
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/My_ACM_Dream/article/details/44090705
版权

最长连续重复子串问题，一般是问重复的最大数，更强的就是要求输出字典序最小的一个答案。

对于这种问题，我们没有什么好办，只能枚举重复子串的长度L，在纸上写些数据会发现，重复子串会覆盖s[0],s[L],s[2*L],...,S[i*L]其中相邻的两个。这样求s[i*L]和s[(i+1)*L]的最长公共前缀为k，通过测试数据会发现这样的规律,重复次数r=k/L+1；这样做只能说肯定能找到对应的连续重复子串的长度和重复次数，但是不一定能找出重复最多的联系子串，原因可通过在纸上模拟L的枚举过程，会发现其实如果先前延伸某个长度在继续执行s[t]和s[t+L]的最长公共子串算法，如果得到的结果比之前的更优，那么之前的k要加1。那么问题来了，对于这样求最长公共子串复杂度最坏是O(n)，这样算下总的复杂度O(n*n*logn)显然太慢。其实可以通过rmq问题对s[t]和s[t+L]最长公共子串进行优化。预处理好，然后就可以通过rmq O(1)的时间复杂度来找到答案。

对于rmq问题，特地去网上找了下，原来就是一种dp而已，相当于利用区间dp来求解问题。状态：dp[i][j]表示区间[i,i+2^j-1]这个区间的最有解。这样状态方程就可以出来了：

dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1])；仔细思考下这个方程。

询问时通过这样找到答案：k=log2(j-i+1), ans=min(dp[i][k],dp[j-(1<<k)][k])；i，j分别是rank值，这个dp方程是以rank的大小顺序进行转移的。可惜悲伤的是自己手残把dp[j-(1<<k)][k]写成了dp[i-(1<<k)][k]；搞了一天，- -！各种找错，瞬间怀疑自己能力了。哎，说多都是泪。

奉上两题：poj 3693(好题) spoj(687)

题解：

poj 3697:

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=100005;
int rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],t4[maxn];
char str[maxn];
int Log2[maxn];
int ans[maxn];

void Swap(int*& x,int*& y){

    int *temp=x;
    x=y;
    y=temp;
}

bool cmp(int t[],int a,int b,int l){

    return t[a]==t[b]&&t[a+l]==t[b+l];
}

void build_SA(char s[],int len,int up){

    int *k1=t1,*k2=t2,*r=t3,*cnt=t4;
    for(int i=0;i<up;i++)cnt[i]=0;
    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;

    for(int d=1,p=0;p<len;d<<=1,up=p){

        p=0;
        for(int i=len-d;i<len;i++)k2[p++]=i;
        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
        for(int i=0;i<len;i++)r[i]=k1[k2[i]];

        for(int i=0;i<up;i++)cnt[i]=0;
        for(int i=0;i<len;i++)cnt[r[i]]++;
        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];

        Swap(k1,k2);
        k1[SA[0]]=0;
        p=1;
        for(int i=1;i<len;i++){
            k1[SA[i]]= cmp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
        }
    }
}

void get_height(char s[],int len){

    //memset(rank,0,sizeof rank);
    for(int i=1;i<=len;i++)rank[SA[i]]=i;
    for(int i=0,p=0;i<len;i++){

        int j=SA[rank[i]-1];
        while(s[i+p]==s[j+p])p++;
        height[rank[i]]=p;
        if(p)p--;
    }
}

void cmin(int& a,int b){

    if(b<a) a=b;
}

int dp[maxn][20];
void rmq(int n){

    int m=Log2[n];
    for(int i=1;i<=n;i++)dp[i][0]=height[i];
    for(int j=1;j<=m;j++){
        for(int i=n;i;i--){
            dp[i][j]=dp[i][j-1];
            if(i+B(j-1)<=n){
                cmin(dp[i][j],dp[i+B(j-1)][j-1]);
            }
        }
    }
}

int query(int l,int r){

    int i=rank[l],j=rank[r];
    if(i>j){ int temp=i; i=j; j=temp; }
    i++;
    int k=Log2[j-i+1];
    return min(dp[i][k],dp[j-B(k)+1][k]);
}

void get_log()
{
    Log2[0]=-1;
    for(int i=1;i<maxn;i++){
        Log2[i]=(i&(i-1))?Log2[i-1]:(Log2[i-1]+1);
    }
}

int main(){

    get_log();
    int cas=1;
    while(scanf("%s",str)!=EOF){
        if(str[0]=='#')break;

        int len=strlen(str);
        build_SA(str,len+1,200);
        get_height(str,len);
        rmq(len);
        int mx=0,tol=0;
        for(int l=1;l<len;l++){

            for(int i=0;i+l<len;i+=l){

                int r=query(i,i+l);///最长公共前缀长度
                int k=r/l+1;///循环次数
                int t=i-(l-r%l);

                if(t>=0&&r%l&&query(t,t+l)>=r) k++;
                if(k>mx){
                    mx=k;
                    tol=0;
                    ans[tol++]=l;
                }else if(k==mx){
                    ans[tol++]=l;
                }
            }
        }
        int L=-1,s=0;
        for(int i=1;i<=len&&L==-1;i++){
            for(int j=0;j<tol;j++){
                int l=ans[j];
                if(query(SA[i],SA[i]+l)>=(mx-1)*l){
                    L=l;
                    s=SA[i];
                    break;
                }
            }
        }
        if(L==-1){
            printf("Case %d: %c\n",cas++,str[SA[1]]);
            continue;
        }
        printf("Case %d: ",cas++);
        for(int i=s;i<s+L*mx;i++)printf("%c",str[i]);
        puts("");
    }
    return 0;
}
/**
baccdbaccdbacbdbacbd
xbcabcab
edbea
abcabcabcab
accdaccaccddcacacaacad
babbabaabaabaabab
abaabaabaaba
ccabababc
daabbccaa
bcbcbaba
pedabacewqpz
ababaaa
baccdbaccdbacbdbacbd
asdfghj
asasasaaa
asasasaaaa
bbbcccaa

Case 1: accdbaccdb
Case 2: bcabca
Case 3: a
Case 4: abcabcabc
Case 5: cacaca
Case 6: abaabaabaaba
Case 7: abaabaabaaba
Case 8: ababab
Case 9: aa
Case 10: baba
Case 11: a
Case 12: aaa
Case 13: accdbaccdb
Case 14: a
Case 15: aaa
Case 16: aaaa
Case 17: bbb

*/

spoj 687: 

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=50005;
int rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],t4[maxn];
char str[maxn];
int Log2[maxn];

void Swap(int*& x,int*& y){

    int *temp=x;
    x=y;
    y=temp;
}

bool cmp(int t[],int a,int b,int l){

    return t[a]==t[b]&&t[a+l]==t[b+l];
}

void build_SA(char s[],int len,int up){

    int *k1=t1,*k2=t2,*r=t3,*cnt=t4;
    for(int i=0;i<up;i++)cnt[i]=0;
    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;

    for(int d=1,p=1;p<len;d<<=1,up=p){

        p=0;
        for(int i=len-d;i<len;i++)k2[p++]=i;
        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
        for(int i=0;i<len;i++)r[i]=k1[k2[i]];

        for(int i=0;i<up;i++)cnt[i]=0;
        for(int i=0;i<len;i++)cnt[r[i]]++;
        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];

        Swap(k1,k2);
        k1[SA[0]]=0;
        p=1;
        for(int i=1;i<len;i++){
            k1[SA[i]]= cmp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
        }
    }
}

void get_height(char s[],int len){

    for(int i=1;i<=len;i++)rank[SA[i]]=i;
    for(int i=0,p=0;i<len;i++){

        int j=SA[rank[i]-1];
        while(s[i+p]==s[j+p])p++;
        height[rank[i]]=p;
        if(p)p--;
    }
}

int dp[maxn][20];
void rmq(int n){

    int m=Log2[n];
    for(int i=1;i<=n;i++)dp[i][0]=height[i];
    for(int j=1;j<=m;j++){
        for(int i=1;i+B(j)-1<=n;i++){
            dp[i][j]=min(dp[i][j-1],dp[i+B(j-1)][j-1]);
        }
    }
}

int query(int l,int r){

    int i=rank[l],j=rank[r];
    if(i>j){ int temp=i; i=j; j=temp; }
    i++;
    int k=Log2[j-i+1];
    return min(dp[i][k],dp[j-B(k)+1][k]);
}

void get_log()
{
    Log2[0]=-1;
    for(int i=1;i<maxn;i++){
        Log2[i]=(i&(i-1))?Log2[i-1]:(Log2[i-1]+1);
    }
}

int main(){

    get_log();
    int T,len;
    scanf("%d",&T);
    while(T--){

        scanf("%d",&len);
        for(int i=0;i<len;i++){
            getchar();
            scanf("%c",&str[i]);
        }
        str[len]='\0';
        build_SA(str,len+1,128);
        get_height(str,len);
        rmq(len);
        int ans=0;
        for(int l=1;l<len;l++){

            for(int i=0;i+l<len;i+=l){

                int r=query(i,i+l);
                int k=r/l+1;
                int t=i-(l-r%l);
                if(t>=0&&r%l&&query(t,t+l)>=r) k++;
                if(k>ans) ans=k;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}
/**
*/



My_ACM_Dream
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
js代码-最长重复子串
07-15
它要求找出一个字符串中最长的连续子序列，这个子序列在原字符串中至少出现两次。这个问题在文本处理、数据压缩、生物信息学等领域都有广泛应用。下面我们将深入探讨如何使用JavaScript解决这个问题。 首先，我们...
Python查找最长不包含重复字符的子字符串算法示例
09-19
对于给定的字符串，如何找到其中最长的不包含重复字符的子字符串是一项典型的算法挑战。这类问题不仅能够帮助我们深入理解字符串操作的基本原理，还能够提升我们在实际开发中的编程技巧。 本篇文章将详细介绍一种...
[算法]最大连续子数组和，最长重复子串，最长无重复字符子串
weixin_30633949的博客
08-10 246
这几道题是我在面试中亲身经历的，在面试滴滴的过程中，我遇到过最大子数组和；在面试阿里的过程中，我遇到过最长重复子串；在面试头条过程中，我遇到过最长无重复字符子串。 1. 最大子数组和 比如，给定一个数组， 1, -2, 3, -4, 5, 6, -7 应该输出， 11。 public static int maxSubArray(int[] arr) { i...
动态规划——最长公共子序列
milu_ELK的博客
09-30 493
比如X总共8位，Y总共5位，那么当j=5，j再增长也没有意义了，比如C[6][6]=C[6][5],不可能比C[5][5]（即C[5][6]，二者同样)要更小,因此C[6][6]=max{C[6][5],C[5][6]},C[i+1][j+1]=max{C[i+1][j],C[i][j+1]}化为C[i][j]=max{C[i][j-1],C[i-1][j]}=Yi,C[i+1][j+1]=max{C[i+1][j],C[i][j+1]}化为C[i][j]=max{C[i][j-1],C[i-1][j]}
算法：字符串全部子序列、子串、全排列
OceanStar的博客
03-29 6716
我们先来看看字符串的子串和子序列有什么区别？ 字符串的子串：必须是连续的一段 字符串的子序列：可以不连续，但是相对次序不能乱（每个字符 要跟不要 所有的路都走一遍）——深度优先遍历 打印出字符串的全部子串 #include <string> #include <vector> using namespace std; void helper(std::string str, std::vector<std::string> & ans){ for(i
最快方法求最长上升子序列（LIS）+最长公共子序列（LCS）模板（C/C++）
m0_71934846的博客
03-19 2956
最快方法求最长上升子序列+最长公共子序列
后缀树与最长回文子串
ice110956的专栏
12-13 5008
后缀树与最长回文子串 首先引入求最长回文子串的问题: 回文串:正反一样的字符串. 如ABCDCBA,反过来是还是ABCDCBA.   回文子串:DABCAEA中有两个回文子串(不包括长度为1的),即ABCA,AEA. 长度ABCA>AEA.   那么字符串DABCAEA的最长回文子串就是ABCA.   解决方法: 1.      遍历所有的子串.  长度N的字符串,共
字符串相关处理kmp，前缀数，后缀树，后缀数组，最长回文串，最长重复字串，最长非重复字串
leileicaocao的专栏
06-04 672
1. 最长回文串 一般用后缀数组或者后缀树可以解决， 用此方法：http://blog.csdn.net/v_july_v/article/details/6897097 预处理后缀树，使得查询LCA的复杂度为O(1)。这步的开销是O(N)，N是单词S的长度 ；对单词的每一位置i(也就是从0到N-1)，获取LCA(S(i), S‘(N-i-1)) 以及LCA(S(i),
《1024程序节-算法修行》无重复字符的最长子串的《四种解决方案》
专注JAVA、算法、数据库等干货分享，同时维护【IT技术资源免费分享】圈子
10-24 539
暴力破解、滑动窗口....，《无重复字符的最长子串的四种解决方案！》
聪明的美食家（最长上升子序列）
qq_51537085的博客
04-17 313
如果有人认为吃东西只需要嘴巴，那就错了。 都知道舌头有这么一个特性，“由简入奢易，由奢如简难”（据好事者考究，此规律也适合许多其他情况）。具体而言，如果是甜食，当你吃的食物不如前面刚吃过的东西甜，就很不爽了。 大宝是一个聪明的美食家，当然深谙此道。一次他来到某小吃一条街，准备从街的一头吃到另一头。为了吃得爽，他大费周章，得到了各种食物的“美味度”。他拒绝不爽的经历，不走回头路而且还要爽歪歪（爽的次数尽量多）。 输入 两行数据。 第一行为一个整数n，表示小吃街上小吃的数量 第二行为n个整数，分别表..
算法题目二：寻找最长重复子序列(2) trie树与后缀树的应用
1ouis的专栏
05-06 1120
本文是转载！原博客地址为：http://www.cnblogs.com/aiyelinglong/archive/2012/04/09/2439777.html 仅作为个人算法学习的记录之用，不用做任何其他用途。   Trie树，即字典树，又称单词查找树或键树，是一种树形结构，是一种哈希树的变种。典型应用是用于统计和排序大量的字符串（但不仅限于字符串），所以经常被搜索引擎系统用于
Leetcode无重复字符的最长子串_leetcode_C++_direct8eo_最长重复子串_
09-29
在本项目中，我们主要探讨的是一个经典的编程问题——寻找字符串中最长的无重复字符的子串。这个问题在LeetCode平台上被广泛讨论，是许多程序员面试时常见的算法题。题目要求设计一个算法，从给定的字符串中找到最长...
Python实现针对给定字符串寻找最长非重复子串的方法
09-20
在示例中，分别对字符串"abcabcd"，"abcabcabd"和"bbbbbbb"进行了处理，结果分别是"abc"，"ab"和"b"，这些都是各自字符串中最长的非重复子串。 这个解决方案的关键在于巧妙地结合了滑窗切片和计数器来快速有效地...
hdu 4513 吉哥系列故事――完美队形II (Manacher)
ACM之梦
02-13 1455
题意： 给出一个队列，求出这个队列满足既是回文并且左断点到终点是成递增的最长子队列。 题解： dp，kmp，Manacher，都可以解决。Manacher最高效。 #include #include #include #include #include #include #include using namespace std; //typedef long long lld; const
yy浅谈后缀数组
ACM之梦
03-07 1063
后缀数组是一种很优美的数组结构，基本上可以完全替代后缀数，后缀数组无论是在内存的占用或者是时间复杂度都和后缀树差不多，但是操作性比后缀树好很多。 后缀数组可以说是字符串处理的神奇啊！想的到的，想不到的都可以处理。 后缀数组的应用也是十分的广泛。后缀数组普遍可以通过O(nlogn)的复杂度解决很多问题。 一下是一些例题和后缀数组的一些应用。 1、求两个串后缀的最长公共前缀，height数组
hihoCoder 1084 扩展KMP(后缀数组)
ACM之梦
04-27 1057
题意： 给出两个串，A和B，问A串中有多少个子串和B是‘匹配’的。 这里的匹配要满足：将一个字符串转化为另一个字符串所需要改变字符的最小数目。对于A中的一个子串，如果它的长度和B的相同且它们之间的海明码距离不超过k，我们认为它们是匹配的. 题解： 用后最数组匹配A的正串的每个后缀和B的正串的最大匹配数，相反同样用后缀数组匹配A的反串的后缀和B的反串的最大匹配。分别用f[i][0]和f[i]
poj 2758 Checking the Text (后缀数组应用)
ACM之梦
03-09 858
题意： 给出一个串，对这个串进行插入，查询，查询是按照原来串的位置进行查询a,b为开头的后缀的最长公共前缀。 题解： 的确是一道好题。参考大牛博客 #include #include #include #include #include #include #include #include #include using namespace std; #define B(x) (1<
uva 719 Glass Beads (后缀自动机)
ACM之梦
03-10 851
题意： 给出一个串，如bcdaf长度n=5，可以将串头尾相连，在这个循环串中求字典序最小且长度为n的串的开头点位置。 题解： 后缀自动机中遍历dfs找到第n个点对应的肯定是长度为n字典序最小的子串，并且结尾点就是遍历的结尾点。 #include #include #include #include #include #include #include #include #inclu
动态规划与最长重复子串
最新发布
11-08
而最长重复子串问题是指在一个字符串中找到最长的重复子串，即至少出现两次的子串中最长的一个。这个问题可以通过动态规划来解决。 具体来说，可以使用一个二维数组dp[i][j]来表示以第i个字符结尾的子串和以第j个...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值