manachar算法小结

1.hdu--4513 吉哥系列故事——完美队形II

http://acm.hdu.edu.cn/showproblem.php?pid=4513

题意：中文题不解释

思路：数字型的manachar算法。将模板中的比较字符改为比较数字就行了。

AC代码：

#include <iostream>
#include<cstdio>
#include<cstring>
#include <algorithm>
#include <cmath>
using namespace std;
int a[200005],p[200005],k[200005];
int manachar(int len)
{
    int maxn,mx,id,i;
    maxn=mx=id=0;
    memset(p,0,sizeof(p));
    for(i=1; i<len; i++)
    {
        if(mx>i)
            p[i]=min(p[id*2-i],mx-i);
        else
            p[i]=1;
        for(; a[i-p[i]]==a[i+p[i]]&&p[i]<=k[i]; p[i]++)
            if(p[i]+i>mx)
            {
                mx=p[i]+i;
                id=i;
            }
        if(p[i]>maxn)
            maxn=p[i];
    }
    printf("%d\n",maxn-1);
    return 0;
}
int main()
{
    int i,t,n;
    while(~scanf("%d",&t))
    {
        while(t--)
        {
            scanf("%d",&n);
            memset(a,0,sizeof(a));
            memset(k,0,sizeof(k));
            a[0]=-1;
            for(i=0; i<n; i++)
            {
                scanf("%d",&a[i*2+2]);
                if(a[i*2+2]>=a[i*2])
                {
                    k[i*2+1]=k[i*2]+1;
                    k[i*2+2]=k[i*2]+2;
                }
                else
                    k[i*2+2]++;
                //printf("%d %d\n",i*2+2,k[i*2+2]);
            }
            manachar(2*n+2);
        }
    }
    return 0;
}

2.hdu--3294 Girls' research

http://acm.hdu.edu.cn/showproblem.php?pid=3294

题意：先输入一个字符ch表示该字符的real是‘a’，循环对应的ch-1的real是‘z’，再输入一个字符串，求这个字符串real值的最大回文串，如果有多个结果输出第一个。输出最大回文串的左端点值和右端点值，如果回文串长度是1的话打印“No solution!”

思路：先把给的字符串改成real串，再求最大回文串，并记录real串的左右端点。

AC代码：

#include <iostream>
#include<cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
char ch,s[200005],str[400005];
int p[400005];
int manachar(int len)
{
    int maxn,mx,id,pi;
    maxn=mx=id=0;
    memset(p,0,sizeof(p));
    for(int i=1;i<len;i++)
    {
        if(mx>i)
            p[i]=min(p[id*2-i],mx-i);
        else
            p[i]=1;
        for(;str[i-p[i]]==str[i+p[i]];p[i]++)
        {
            if(p[i]+i>mx)
            {
                mx=p[i]+i;
                id=i;
            }
        }
        if(maxn<p[i])
        {
            pi=i;
            maxn=p[i];
        }
    }
    maxn--;
    if(maxn<3)
        printf("No solution!\n");
    else
    {
        printf("%d %d\n",(pi-maxn)/2,(pi+maxn-2)/2);
        for(int i=pi-maxn+1;i<=pi+maxn-1;i=i+2)
            printf("%c",str[i]);
        printf("\n");
    }
    return 0;
}
int main()
{
    while(~scanf("%c",&ch))
    {
        getchar();
        gets(s);
        str[0]='$',str[1]='#';
        int len=strlen(s);
        for(int i=0;i<len;i++)
        {
            if(s[i]>=ch)
            str[i*2+2]=s[i]-(ch-'a');
            else
                str[i*2+2]=s[i]+26-(ch-'a');
            str[i*2+3]='#';
        }
        str[len*2+2]='\0';
        manachar(len*2+2);
    }
    return 0;
}

3.hdu--3613 Best Reward

http://acm.hdu.edu.cn/showproblem.php?pid=3613

题意：先输入‘a’~‘z’的价值，然后给出一个字符串，问把这个字符串分成两部分，如果这一部分是回文串则把总价值加上这部分串的价值，求最大的总价值

思路：先把给出的串的总价值进行累加放入数组中，然后对模式串给出的查询串用manachar算法进行计算，如果当前查询到的i<len/2并且p[i]==i时，就max（ans，num[i-1]）,如果i>len/2并且p[i]==len-i时，就max（ans，num[n]-num[n-p[i]+1]）,如果ans<num[n]在统筹算一下

AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=500005;
char s[maxn],str[maxn<<1];
int num[maxn],k[maxn<<1],p[maxn<<1];
int manachar(int len,int n)
{
    int ans,id,mx;
    id=ans=mx=0;
    memset(p,0,sizeof(p));
    memset(k,0,sizeof(k));
    for(int i=1; i<len; i++)
    {
        if(mx>i)
            p[i]=min(p[id*2-i],mx-i);
        else
            p[i]=1;
        for(; str[i-p[i]]==str[i+p[i]]; p[i]++)
            if(p[i]+i>mx)
            {
                mx=p[i]+i;
                id=i;
            }
        if(i<len/2&&p[i]==i)
        {
            ans=max(ans,num[i-1]);
            k[i]=i*2-1;
        }
        else if(i>len/2&&p[i]==len-i)
        {
            ans=max(ans,num[n]-num[n-p[i]+1]);
            k[i]=i*2-len+1;
        }
    }
    for(int i=1; i<len/2; i++)
        if(k[i]==k[len/2+i-1]&&k[i]!=0)
        {
            ans=max(ans,num[n]);
            break;
        }
    printf("%d\n",ans);
    return 0;
}
int main()
{
    int a[26],t;
    while(~scanf("%d",&t))
    {
        while(t--)
        {
            for(int i=0; i<26; i++)
                scanf("%d",&a[i]);
            scanf("%s",s);
            int len=strlen(s);
            memset(str,'#',sizeof(str));
            memset(num,0,sizeof(num));
            str[0]='$';
            for(int i=0; i<len; i++)
            {
                str[i*2+2]=s[i];
                num[i+1]=num[i]+a[s[i]-'a'];
            }
            manachar(len*2+2,len);
        }
    }
    return 0;
}

4.poj--3974 Palindrome

http://poj.org/problem?id=3974

题意：给出一系列字符串以“END”结束，问每个字符串的最大回文串的长度

思路：manachar算法直接算。

AC代码：

#include <iostream>
#include<cstdio>
#include<cstring>
#include <algorithm>
#include <cmath>
using namespace std;
int p[2000005];
char s[1000005],str[2000005];
int manachar(int len)
{
    int maxn,mx,id,i;
    maxn=mx=id=0;
    memset(p,0,sizeof(p));
    for(i=1; i<len; i++)
    {
        if(mx>i)
            p[i]=min(p[id*2-i],mx-i);
        else
            p[i]=1;
        for(; str[i-p[i]]==str[i+p[i]]; p[i]++)
            if(p[i]+i>mx)
            {
                mx=p[i]+i;
                id=i;
            }
        maxn=max(maxn,p[i]);
    }
    printf("%d\n",maxn-1);
    return 0;
}
int main()
{
    int t;
    t=1;
    while(~scanf("%s",s))
    {
        if(strcmp(s,"END")==0)
        break;
        int len=strlen(s);
        memset(str,'#',sizeof(str));
        str[0]='$';
        for(int i=0; i<len; i++)
        {
            str[i*2+2]=s[i];
        }
        str[len*2+2]='\0';
        printf("Case %d: ",t++);
        manachar(len*2+2);
    }
    return 0;
}