BZOJ 3790 神奇项链 (Manacher 贪心)

本文介绍了一个有趣的算法问题——神奇项链,通过Manacher算法求解如何使用最少次数的机器操作来生成特定样式的回文串项链。文章详细展示了Manacher算法的具体实现过程,并采用贪心策略确定最小覆盖区间的解决方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

3790: 神奇项链

Time Limit: 10 Sec Memory Limit: 64 MB
Description

母亲节就要到了,小 H 准备送给她一个特殊的项链。这个项链可以看作一个用小写字
母组成的字符串,每个小写字母表示一种颜色。为了制作这个项链,小 H 购买了两个机器。第一个机器可以生成所有形式的回文串,第二个机器可以把两个回文串连接起来,而且第二个机器还有一个特殊的性质:假如一个字符串的后缀和一个字符串的前缀是完全相同的,那么可以将这个重复部分重叠。例如:aba和aca连接起来,可以生成串abaaca或 abaca。现在给出目标项链的样式,询问你需要使用第二个机器多少次才能生成这个特殊的项链。
Input

输入数据有多行,每行一个字符串,表示目标项链的样式。
Output

多行,每行一个答案表示最少需要使用第二个机器的次数。
Sample Input

abcdcba

abacada

abcdef
Sample Output

0

2

5
HINT

每个测试数据,输入不超过 5行

每行的字符串长度小于等于 50000

思路:
显然若要选择以一个位置为中心的回文串,那么这个回文串取较长的不会比较短的差。
于是一遍Manacher过后,求出了以每个点为中心的最长回文串,所以问题就变成了选择最少的区间覆盖整个区间。经典贪心。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;

const int N = 200010;
int n, p[N];
char s[N], ss[N];

struct Edge{ 
    int l, r;
}a[N];

bool cmp(Edge a, Edge b){
    return a.l < b.l || (a.l == b.l && a.r > b.r);
}

void manacher(){
    memset(p, 0, sizeof(p));
    int id = 0, mx = 0;
    for(int i=1; i<=n; i++){
        if(mx > i) p[i] = min(p[id*2-i], mx-i+1);
        else p[i] = 1;
        while(s[i+p[i]] == s[i-p[i]]) p[i]++;
        if(i+p[i]-1 > mx){
              mx = i+p[i]-1;
              id = i;
        }
        if(i % 2) a[i].l = i/2 - p[i]/2 + 1, a[i].r = i/2 + p[i]/2;
        else a[i].l = i/2 - p[i]/2 + 1, a[i].r = i/2 + p[i]/2 - 1;
        //原串上的覆盖位置 
     }
}

int main(){
    while( scanf("%s", ss) != EOF ){
        int tot = strlen( ss );
        n = 1; s[0] = '+'; s[1] = '#';
        for(int i=0; i<tot; i++){
            s[++n] = ss[i];
            s[++n] = '#';
        }
        s[n+1] = '-';//在首尾加上两个不相同且原串中没有的字符,防止越界 
        manacher();
        sort(a+1, a+n+1, cmp);
        int rg = a[1].r, ans=1, now=2;
        while(rg < tot){//rg为能覆盖到的最右节点 
            int maxx = rg;
            while(now <= n && a[now].l <= rg + 1){//一个可行的区间 
                maxx = max(maxx, a[now].r);
                now++;
            }//区间是单调排序的,找到最优解 
            ans++; rg = maxx;
        }
        printf("%d\n", ans-1);//拼接数是回文串数减1 
     }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值