马拉车算法(求回文子串)

先放上一个常用的字符串函数substr()
用法是substr(st,len);
返回的是以st为下标起始的长度为len的一个字符串
比如
假设有一个string类 比如 string str;
初始化为str=“wuhanjiayou”;
然后 string tmp=str.substr(0,5);
那么tmp就为0开始的长度为5的字符串“wuhan"
如果省略了长度len,那么就代表处理到字符串结尾
即tmp=str.substr(0); 那么tmp就为“wuhanjiayou”
--------------------------------------------
--------------------------------------------
关于马拉车算法,我是看这篇博客学的点这里
我大概是模仿那位博主改成自己习惯的书写方式吧

int st,mxlen;//起始点与最大长度
string ans;//最大回文子串
void manacher(string str)
{
    string tmp="$#";
    int len=str.size();
    for(int i=0;i<len;i++)
    {
        tmp+=str[i];
        tmp+='#';
    }
    int mx=0,id=0,tlen=0,mid=0;
    len=tmp.size();
    vector<int>p(len,0);
    for(int i=1;i<=len-1;i++)
    {
        p[i]= mx>i ? min(p[2*id-i],mx-i) : 1;

        while(tmp[i+p[i]]==tmp[i-p[i]]){p[i]++;}
		此处和下面那个小于号 代表找的是第一个最长的回文子串比如
		abab出来是aba 
		如果改成小于等于就是出现的最后一个
		abab出来是bab
        if(mx<i+p[i]){
            mx=i+p[i];
            id=i;
        }
        if(tlen<p[i]){
            tlen=p[i];
            mid=i;
        }
    }
    mxlen=tlen-1;
    st=(mid-tlen)/2;
    ans=str.substr(st,mxlen);
}
void solve()
{
    string str;cin>>str;
    manacher(str);
    cout<<"最大回文长度子串长度为  "<<mxlen<<endl;
    cout<<"起始下标为  "<<st<<endl;
    cout<<"最大回文子串为  "<<ans<<endl;
}

在这里插入图片描述
例题 HDU3068
在这里插入图片描述


//#include<bits/stdc++.h>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<sstream>
#include<algorithm>
#include<queue>
#include<stack>
#include<vector>
#include<cmath>
#include<map>
#include<set>
#include<utility>
#define de(x)       cout<<"\n----debug----"<<x<<"-------\n"
#define rep(i,a,n)  for(int i=a;i<=n;i++)
#define per(i,n,a)  for(int i=n;i>=a;i--)
#define ms(arr,x)   memset((arr),x,sizeof(arr))
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define ll long long
#define MP make_pair
#define pb push_back
#define inf 0x3f3f3f3f
#define INF 0x7fffffff
#define fi  first
#define se  second
//#define ls  (rt<<1)
//#define rs  (rt<<1|1)
//#define md  (l+r>>1)
using namespace std;
typedef pair<int,int> pii;
inline int read()
{
    int ans=0,sign=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            sign=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans=ans*10+(ch-'0');
        ch=getchar();
    }
    return sign*ans;
}
const ll mod=1e9+7;
const double PI=acos(-1.0);
const double eps=1e-5;
const int maxn=1e5+5;
const int maxm=1e5+5;
inline int gcd(int a,int b){return b==0 ? a:gcd(b,a%b);}
inline int lcm(int a,int b){return a*b/gcd(a,b);}
int manacher(string str){
    string tmp="$#";
    int len=str.size();
    rep(i,0,len-1){
        tmp+=str[i];
        tmp+='#';
    }
    len=tmp.size();
    vector<int>p(len,0);
    int mx=0,id=0,mid=0,mxlen=0;
    for(int i=1;i<=len-1;i++){
        p[i]= mx>i ? min(p[id*2-i],mx-i):1;
        while(tmp[i+p[i]]==tmp[i-p[i]]){p[i]++;}
        if(mx<i+p[i]){
            mx=i+p[i];
            id=i;
        }
        if(mxlen<p[i]){
            mxlen=p[i];
            mid=i;
        }
    }
    return mxlen-1;
}
void solve()
{
    string str;
    while(cin>>str){
        cout<<manacher(str)<<endl;
    }
}
int main(){
    int test=1;
    //test=read();
    while(test--)
    {
        solve();
    }
    #ifdef _LOCALE_DEBUG_PAUSE
    system("pause");
    #endif //_LOCALE_DEBUGE_PAUSE
    return 0;
}

猜你喜欢

转载自blog.csdn.net/leoxe/article/details/105538707