Manacher算法、扩展Kmp

已于 2024-08-15 16:16:37 修改
阅读量174 收藏 4
点赞数 3
文章标签： 算法
于 2024-08-14 15:26:51 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/2301_80275822/article/details/141181780
版权

求一段字符串中的最长回文字符串

暴力方法：枚举每个点不停的向左右两边扩展，但这样偶数的回文串不能求出来

解决：在每两个字符之间及字符串的左右两头都插上任意字符(扩展串)，然后枚举每个点向左右两边扩展，得到的扩展半径-1就是回文串长度

   P3805 【模板】manacher - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;

public class Main {

	public static int MAXN = 11000001;

	public static char[] ss = new char[MAXN << 1];

	public static int[] p = new int[MAXN << 1];

	public static int n;

	public static void main(String[] args) throws IOException {
		BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
		PrintWriter out = new PrintWriter(new OutputStreamWriter(System.out));
		out.println(manacher(in.readLine()));
		out.flush();
		out.close();
		in.close();
	}

	public static int manacher(String str) {
		manacherss(str.toCharArray());
		int max = 0;
		for (int i = 0, c = 0, r = 0, len; i < n; i++) {
			len = r > i ? Math.min(p[2 * c - i], r - i) : 1;
			while (i + len < n && i - len >= 0 && ss[i + len] == ss[i - len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			max = Math.max(max, len);
			p[i] = len;
		}
		return max - 1;
	}

	public static void manacherss(char[] a) {
		n = a.length * 2 + 1;
		for (int i = 0, j = 0; i < n; i++) {
			ss[i] = (i & 1) == 0 ? '#' : a[j++];
		}
	}

}

扩展Kmp

 z数组：一个字符串从每个位置开始的字串和原字符串开头匹配上的最长公共前缀长度

 e数组：字符串a的每一个位置和字符串b的开头匹配上的最长公共前缀长度 

P5410 【模板】扩展 KMP/exKMP（Z 函数） - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;

public class Main {

	public static int MAXN = 20000001;

	public static int[] z = new int[MAXN];

	public static int[] e = new int[MAXN];

	public static void main(String[] args) throws IOException {
		BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
		PrintWriter out = new PrintWriter(new OutputStreamWriter(System.out));
		char[] a = in.readLine().toCharArray();
		char[] b = in.readLine().toCharArray();
		zArray(b, b.length);
		eArray(a, b, a.length, b.length);
		out.println(eor(z, b.length));
		out.println(eor(e, a.length));
		out.flush();
		out.close();
		in.close();
	}

	// 非常像Manacher算法
	public static void zArray(char[] s, int n) {
		z[0] = n;
		for (int i = 1, c = 1, r = 1, len; i < n; i++) {
			len = r > i ? Math.min(r - i, z[i - c]) : 0;
			while (i + len < n && s[i + len] == s[len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			z[i] = len;
		}
	}

	// 非常像Manacher算法
	public static void eArray(char[] a, char[] b, int n, int m) {
		for (int i = 0, c = 0, r = 0, len; i < n; i++) {
			len = r > i ? Math.min(r - i, z[i - c]) : 0;
			while (i + len < n && len < m && a[i + len] == b[len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			e[i] = len;
		}
	}

	public static long eor(int[] arr, int n) {
		long ans = 0;
		for (int i = 0; i < n; i++) {
			ans ^= (long) (i + 1) * (arr[i] + 1);
		}
		return ans;
	}

}

   

 3031. 将单词恢复初始状态所需的最短时间 II - 力扣（LeetCode）

从第一个字符开始枚举间隔k长度的字符是否和原字符串的前缀完全相同：如果相同，那么跳过几个间隔就花费几秒；如果到字符串最后还没匹配上，那么字符串的总长度除以向上取整就是花费的时间。所以用z数组就可以解决

public class Solution {

	public static int minimumTimeToInitialState(String word, int k) {
		char[] s = word.toCharArray();
		int n = s.length;
		zArray(s, n);
		for (int i = k; i < n; i += k) {
			if (z[i] == n - i) {
				return i / k;
			}
		}
		return (n + k - 1) / k;
	}

	// leetcode增加了数据量
	// 所以把这个值改成10^6规模
	public static int MAXN = 1000001;

	public static int[] z = new int[MAXN];

	public static void zArray(char[] s, int n) {
		z[0] = n;
		for (int i = 1, c = 1, r = 1, len; i < n; i++) {
			len = r > i ? Math.min(r - i, z[i - c]) : 0;
			while (i + len < n && s[i + len] == s[len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			z[i] = len;
		}
	}

}

 

 5. 最长回文子串 - 力扣（LeetCode）

class Solution {
public:
    static const int maxn=1001;
    int p[maxn<<1];
    int end=0,maxnum=0;
    void get_manacher(string s){
        for(int i=0,r=0,c=0,len;i<s.size();i++){
            len=r>i?min(p[2*c-i],r-i):1;
            while(i+len<s.size()&&i-len>=0&&s[len+i]==s[i-len])
            len++;
            if(i+len>r){
                r=i+len;
                c=i;
            }
            if(len>maxnum){
                maxnum=len;
                end=(i+len-1);
            }
            p[i]=len;
        }
    }
    string longestPalindrome(string s) {
         string cur="#";
         for(int i=0;i<s.size();i++){
            cur+=s[i];
            cur+='#';
         }      
         get_manacher(cur);
         int epos=end/2;
         int bpos=epos-(maxnum-1);
         string ans;
         for(int i=bpos;i<epos;i++)
         ans+=s[i];
         return ans;

    }
};

manacher的过程中记录最长的回文半径的结束位置，然后除2是原始字符串的终止位置，往前推原始字符串的长度截取字串即可

647. 回文子串 - 力扣（LeetCode）

class Solution {
public:
    static const int maxn=1001;
    int p[maxn<<1];
    void get_mancher(string s){
        for(int i=0,r=0,c=0,len;i<s.size();i++){
            len=r>i?min(p[2*c-i],r-i):1;
            while(i+len<s.size()&&i-len>=0&&s[i+len]==s[i-len])
            len++;
            if(i+len>r){
                r=i+len;
                c=i;
            }
            p[i]=len;
        }
    }
    int countSubstrings(string s) {
        string cur="#";
        for(int i=0;i<s.size();i++){
            cur+=s[i];
            cur+="#";
        }
        get_mancher(cur);
        int ans=0;
        for(int i=0;i<cur.size();i++)
        ans+=p[i]/2;
        return ans;

    }
};

回文半径长度/2就是原始字符串中一个回文串中含有的回文串个数

  2472. 不重叠回文子字符串的最大数目 - 力扣（LeetCode）

class Solution {
public:
   static const int maxn=2001;
   int p[maxn<<1];
   string cur;
   int manacher(int l,int k){
      for(int i=l,r=l,c=l,len;i<cur.size();i++){
         len=r>i?min(p[2*c-i],r-i):1;
         while(i+len<cur.size()&&i-len>=l&&cur[i-len]==cur[i+len]){
            if(++len>k)
            return i+k+(cur[i+k]=='#'?0:1);
         }
         if(i+len>r){
            r=len+i;
            c=i;
         }
         p[i]=len;
      }
      return -1;
   }
   int f(int k){
      int ans=0;
      int nex=0;
      while((nex=manacher(nex,k))!=-1)
      ans++;
      return ans;
     
   }
    int maxPalindromes(string s, int k) {
         cur="#";
        for(int i=0;i<s.size();i++){
            cur+=s[i];
            cur+="#";
        }
        return f(k);
    }
};

为了更多的切出大于等于k长度的回文串，只要遇到长度满足的回文串，就加1，且返回距离结束位置最近的‘#’字符的位置，从此位置开始重新进行mancher算法，最终就会得到最多的回文串

 P1659 [国家集训队] 拉拉队排练 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

// 长度前k名的奇数长度回文子串长度乘积
// 给定一个字符串s和数值k，只关心所有奇数长度的回文子串
// 返回其中长度前k名的回文子串的长度乘积是多少
// 如果奇数长度的回文子串个数不够k个，返回-1
// 测试链接 : https://www.luogu.com.cn/problem/P1659

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;

public class Main {

	public static int MOD = 19930726;

	public static int MAXN = 10000001;

	public static String[] mk;

	public static int m, n;

	public static long k;

	public static char[] ss = new char[MAXN << 1];

	public static int[] p = new int[MAXN << 1];

	public static int[] cnt = new int[MAXN];

	public static void main(String[] args) throws IOException {
		BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
		PrintWriter out = new PrintWriter(new OutputStreamWriter(System.out));
		mk = in.readLine().split(" ");
		m = Integer.valueOf(mk[0]);
		k = Long.valueOf(mk[1]);
		out.println(compute(in.readLine()));
		out.flush();
		out.close();
		in.close();
	}

	public static int compute(String s) {
		manacher(s);
		for (int i = 1; i < n; i += 2) {
			cnt[p[i] - 1]++;
		}
		long ans = 1;
		long sum = 0;
		for (int len = (m & 1) == 1 ? m : (m - 1); len >= 0 && k >= 0; len -= 2) {
			sum += cnt[len];
			ans = (ans * power(len, Math.min(k, sum))) % MOD;
			k -= sum;
		}
		return k < 0 ? (int) ans : -1;
	}

	public static void manacher(String str) {
		manacherss(str.toCharArray());
		for (int i = 0, c = 0, r = 0, len; i < n; i++) {
			len = r > i ? Math.min(p[2 * c - i], r - i) : 1;
			while (i + len < n && i - len >= 0 && ss[i + len] == ss[i - len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			p[i] = len;
		}
	}

	public static void manacherss(char[] a) {
		n = a.length * 2 + 1;
		for (int i = 0, j = 0; i < n; i++) {
			ss[i] = (i & 1) == 0 ? '#' : a[j++];
		}
	}

	public static long power(long x, long n) {
		long ans = 1;
		while (n > 0) {
			if ((n & 1) == 1) {
				ans = (ans * x) % MOD;
			}
			x = (x * x) % MOD;
			n >>= 1;
		}
		return ans;
	}

}

manacher算出每个中心的最长回文串，统计以‘#’字符为中心的各回文长度的词频(因为不以‘#’为中心的回文串一定不是奇数串)，奇数串最大长度不超过原子符串长度，从最大的奇数串开始向下依次累乘(乘法快速幂)。注意：长度更大的回文串中一定包含长度较小的各奇数回文串，所以乘的时候词频要向下传递

P4555 [国家集训队] 最长双回文串 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

// 最长双回文串长度
// 输入字符串s，求s的最长双回文子串t的长度
// 双回文子串就是可以分成两个回文串的字符串
// 比如"aabb"，可以分成"aa"、"bb"

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;

public class Main {

	public static int MAXN = 100002;

	public static char[] ss = new char[MAXN << 1];

	public static int[] p = new int[MAXN << 1];

	public static int[] left = new int[MAXN << 1];

	public static int[] right = new int[MAXN << 1];

	public static int n;

	public static void main(String[] args) throws IOException {
		BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
		PrintWriter out = new PrintWriter(new OutputStreamWriter(System.out));
		out.println(compute(in.readLine()));
		out.flush();
		out.close();
		in.close();
	}

	public static int compute(String s) {
		manacher(s);
		for (int i = 0, j = 0; i < n; i++) {
			while (i + p[i] > j) {//最先影响到的一定是最长的，就统计
				left[j] = j - i;
				j += 2;
			}
		}
		for (int i = n - 1, j = n - 1; i >= 0; i--) {
			while (i - p[i] < j) {
				right[j] = i - j;
				j -= 2;
			}
		}
		int ans = 0;
		for (int i = 2; i <= n - 3; i += 2) {
			ans = Math.max(ans, left[i] + right[i]);
		}
		return ans;
	}

	public static void manacher(String str) {
		manacherss(str.toCharArray());
		for (int i = 0, c = 0, r = 0, len; i < n; i++) {
			len = r > i ? Math.min(p[2 * c - i], r - i) : 1;
			while (i + len < n && i - len >= 0 && ss[i + len] == ss[i - len]) {
				len++;
			}
			if (i + len > r) {
				r = i + len;
				c = i;
			}
			p[i] = len;
		}
	}

	public static void manacherss(char[] a) {
		n = a.length * 2 + 1;
		for (int i = 0, j = 0; i < n; i++) {
			ss[i] = (i & 1) == 0 ? '#' : a[j++];
		}
	}

}

预处理得到left数组记录每个‘#’字符左边的最大紧贴回文串长度和right数组记录紧贴‘#’字符右边的最大回文串长度，最后统计最大值

little~钰
关注
  • 3
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Manacher算法和扩展kmp
b1169909203的博客
02-21 243
【代码】Manacher算法和扩展kmp。
浅谈Manacher算法与扩展KMP之间的联系
Keep Moving的专栏
07-03 1123
首先，在谈到Manacher算法之前，我们先来看一个小问题:给定一个字符串S，求该字符串的最长回文子串的长度.对于该问题的求解，网上解法颇多，时间复杂度也不尽相同，这里列述几种常见的解法. 解法一       通过枚举S的子串，然后判断该字串是否为回文，由于S的子串个数大约为，加上每次判断需要的时间，所以总的时间复杂度为，空间复杂度为. bool check(string &S, int l
字符串匹配三大O(n)算法：KMP，Manacher，扩展KMP
偶耶的博客
12-20 362
讲三个玄学算法，我肝要爆了
KMP&扩展KMP&Manacher算法基础与习题(第一更）
weixin_43870697的博客
07-27 534
KMP&扩展KMP&Manacher算法基础与习题(第二更） KMP&扩展KMP&Manacher算法基础与习题(第三更） 目录 KMP算法讲解 NEXT数组详解 模板 例题 A:HDU-1711 Number Sequence B:HDU-1686 Oulipo C:HDU-2087 剪花布条 D:HDU-3746 Cyclic Nack...
KMP 、扩展KMP、Manacher算法 总结
alince20008的博客
11-26 193
一. KMP 1 找字符串x是否存在于y串中，或者存在了几次 HDU1711 Number Sequence HDU1686 Oulipo HDU2087 剪花布条 2.求多个字符串的最长公共子串 POJ3080 Blue Jeans HDU1238 Substrings 3.next数组的应用 1) 求循环节 HDU3746 C...
Manacher算法及其扩展
hch977的博客
01-09 2877
回文字符串： 正着看反着看是一样的 abccba abcba 存在一个轴对称 最长回文字符串问题 在一个字符串中找到最长回文字符串，而manacher算法就是去找这个最长回文字符串。 有啥用？ DNA序列 回文基因序列有一些生理学意义 暴力解最长回文问题（O(N2)O(N^2)O(N2)） ab|ab 按照虚轴对称 ababa 按照实轴中间的a对称 遍历每个实轴与虚轴，分别向两边扩展，存最长的 public static String getmaxstr(String s) {
KMP算法和Manacher算法
BadBadBad的博客
03-03 1042
KMP算法用来解决字符串匹配问题: 找到长串中短串出现的位置.
Manacher算法详解
weixin_53566169的博客
04-26 123
Manacher算法 马拉车算法解决最长回文字串的问题 时间复杂度：O(n) 传统思想：以一个字符为中心，向两边扩大，就能返回最长回文字串，但是长度为偶数的字符串是不成立的，需要分情况讨论，而且算法时间复杂度 o(n^2) ,也可以使用动态规划，时间复杂度仍为o(n^2) 解决办法：在传统方法的基础上将每个字符两边添加字符 # ，然后再向两边扩展，这样奇数和偶数的回文串都找到了！ 问题：添加的字符是不是要求字符串中没出现的字符？ 不是，随便添加，因为无论什么时候都是实际上有的的和有的比较，新添加的和新
KMP&扩展KMP&Manacher算法基础与习题(第二更）
weixin_43870697的博客
08-02 437
KMP&扩展KMP&Manacher算法基础与习题(第一更） KMP&扩展KMP&Manacher算法基础与习题(第三更） 目录 扩展KMP算法讲解 例题 A:HDU-2594 Simpsons’ Hidden Talents B:HDU-3336 Count the string C:HDU-4300 Clairewd’s message D...
算法讲解103【扩展】 Manacher算法、扩展KMP.pptx
最新发布
08-27
代码：https://download.csdn.net/download/m0_73085893/89681787
字符串处理- 单模式匹配- MP 算法与 KMP 算法.rar
09-16
2. 并行匹配：Manacher算法使用两个指针P和R，P指向当前匹配的中心，R记录当前能匹配的最大半径。在匹配过程中，可以同时尝试向两边扩展。 3. 利用回文结构：如果新位置的左侧有已匹配的回文子串，那么可以利用这个...
算法进阶课——7.1 启发式合并、Manacher算法
m0_62021646的博客
01-31 1471
7.1 启发式合并、Manacher算法 启发式合并 启发式合并类似于并查集的按秩合并. 例如，一开始有 nnn 个数，每个数在仅包含自己的单独的集合当中，每次将某个集合内的所有元素，合并到另外一个集合当中，最终要将所有数都合并到一个集合当中. 如果暴力地去实现的话，将包含 aaa 个元素的集合，合并到另外一个集合当中，则其时间复杂度为 O(a×op)O(a \times op)O(a×op)，其中 opopop 为每个元素合并到另一个集合的时间...
代码随想录算法训练营19-回溯1
weixin_53416771的博客
08-22 420
画出来的树是这样记录所有组合，一个变量current装当前的处理结果，一个res装所有的处理的结果下面的代码还可以剪枝优化。
算法日记day 45（单调栈之每日温度|接雨水）
m0_74980681的博客
08-22 1343
如果大于，则取栈顶元素为中间值并弹出栈顶元素，将遍历的元素与新的栈顶元素值比较取较小者，减去中间值得到高度h，然后栈顶元素下标值减去当前遍历的元素下标值再减一得到宽度w，最后高×宽得到面积，并记录，重复上述操作，直到遍历结束。第三种，如果遍历的元素值大于栈顶元素，用该下标减去栈顶元素下标，并将栈顶元素弹出，继续用该元素与新的栈顶元素比较，如果小于等于，则压入栈，如果大于，重复第三步操作。第一种，如果之后遍历的元素值小于栈顶压入的元素，将该元素的下标压入栈中。如果小于栈顶元素值，则压入栈内。
算法025：前缀和
m0_62319039的博客
08-22 636
算法025：前缀和
1087：级数求和
如果没有特别说明，c++都是基于GCC-4.9.2运行的
08-22 655
【代码】：1087：级数求和
算法笔试-编程练习-01-B-23
志远的博客
08-23 1197
小红组内有 n 个人，大家合作完成了一个初版方案，初始时大家的愤怒值都是 0。 但是领导对方案并不满意，共需要修改 m 次方案，每次修改会先让第 l 到 r 个人的愤怒值加 1，然后再修改方案；
代码随想录算法训练营第二十二天| 回溯算法理论基础 77. 组合 216.组合总和III 17.电话号码的字母组合
Mwicogito的博客
08-23 1361
回溯法章节开始，逐渐开始正式复习算法设计课的知识。文章图片来源：代码随想录 (https://programmercarl.com/)
字符串匹配算法：KMP、Manacher与扩展KMP解析
"本文主要介绍了三种字符串匹配算法：KMP算法、Manacher算法以及扩展KMP算法，并结合实例分析了暴力方法的不足之处，提出了更高效的解决方案。" ### KMP算法 KMP算法是由Knuth、Morris和Pratt提出的，主要用于解决...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值