[51nod] 1088 最长回文子串 #Hash+二分

1088 最长回文子串

基准时间限制：1 秒 空间限制：131072 KB 分值: 0 难度：基础题

 

回文串是指aba、abba、cccbccc、aaaa这种左右对称的字符串。

输入一个字符串Str，输出Str里最长回文子串的长度。

 

Input

输入Str（Str的长度 <= 1000)

Output

输出最长回文子串的长度L。

Input示例

daabaac

Output示例

5

Analysis分析

Hash解决

具体措施：用秦九韶算法做一个可以用前缀和计算字符串哈希码的hash，然后枚举中点，二分回文字串的长度qwq

具体看代码qwq

 

Code代码

#include<stdio.h>
#include<iostream>
#include<cstring>
#define maxn 421357
#define base 21707
#define ULL unsigned long long
using namespace std;

ULL stad[maxn],anti[maxn];

int n,ans = 1,lens; char s[maxn];

ULL Pow(ULL x,int k){
    if(!k) return 1;
    if(k == 1) return x;
    ULL b = 1;
    while(k){
        if(k&1) b = b*x;
        x = x*x;
        k >>= 1;
    }return b;
}

ULL getcode(int L,int R){
    if(!L) return stad[R];
    else return stad[R]-stad[L-1]*Pow(base,R-L+1);
}

ULL getanti(int L,int R){
    if(R == n) return anti[L];
    else return anti[L]-anti[R+1]*Pow(base,R-L+1);
}

int find(int pos){
    int L = 0,R = min(lens-1-pos,pos),mid;
    while(L < R){
        mid = (L+R+1)/2;
//        printf("Checking... [%d , %d]\n",pos-mid,pos+mid);
//        ULL cntcode = (getanti(pos+1,pos+mid)*base+s[pos])*Pow(base,mid)+getcode(pos+1,pos+mid);
        if(getcode(pos-mid,pos-1) == getanti(pos+1,pos+mid)) L = mid;
        else R = mid-1;
    }return L*2+1;
}

int exfind(int pos){
    if(pos == lens-1) return 0;
    int L = 0,R = min(lens-1-pos,pos+1),mid;
    while(L < R){
        mid = (L+R+1)/2;
//        printf("Checking... [%d , %d]\n",pos-mid,pos+mid);
//        ULL cntcode = (getanti(pos+1,pos+mid)*base+s[pos])*Pow(base,mid)+getcode(pos+1,pos+mid);
        if(getcode(pos-mid+1,pos) == getanti(pos+1,pos+mid)) L = mid;
        else R = mid-1;
    }return L*2;
}

int main(){
    scanf("%s",s);
    lens = strlen(s);
    
    ULL cntcode = 0;
    for(int i = 0;i < lens;i++)
        cntcode = cntcode*base+s[i],
        stad[i] = cntcode;
        
    cntcode = 0;
    for(int i = lens-1;i >= 0;i--)
        cntcode = cntcode*base+s[i],
        anti[i] = cntcode;
        
    for(int pos = 0;pos < lens;pos++)
        ans = max(ans,find(pos)),
        ans = max(ans,exfind(pos));
    cout << ans;
    
    return 0;
}