【jzoj4889】【最长公共回文子序列】【字符串】

题目大意

给出一长一短两个字符串，求最长公共回文子序列

解题思路

枚举短串所有字串判断是否回文，在长串中查找，可以用二分，也可以倒着预处理最近的转移位置。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define min(a,b) ((a<b)?a:b)
#define max(a,b) ((a>b)?a:b)
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const maxn=1e5+1,maxm=20+1,mod=998244353,size=1048576;
int n,m,gra,ans,time,f[30][maxn+10];
char a[maxn+10],b[maxm+10],c[maxm+10];
int main(){
    //freopen("lcps.in","r",stdin);
    //freopen("lcps.out","w",stdout);
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%s",a+1);n=strlen(a+1);
    fo(i,1,n)f[a[i]-'a'][++f[a[i]-'a'][0]]=i;
    scanf("%s",b+1);m=strlen(b+1);
    int maxs=(1<<m)-1,ans=0;
    fo(s,1,maxs){
        int ss=s;c[0]=0;LL tmp=0;
        fo(i,1,m){
            if(ss&1)c[++c[0]]=b[i];
            ss/=2;
        }
        int t=c[0]/2,ok=1;
        fo(i,1,t)
            if(c[i]!=c[c[0]-i+1]){ok=0;break;}
        if(ok){
            int j=1;
            fo(i,1,c[0]){
                int l=1,r=f[c[i]-'a'][0];
                for(;l<r;){
                    int mid=(l+r)/2;
                    if(f[c[i]-'a'][mid]<j)l=mid+1;
                    else r=mid;
                }
                if(f[c[i]-'a'][l]<j)j=n+2;
                else j=f[c[i]-'a'][l]+1;
                if(j>n+1)break;
            }
            if(j<=n+1)ans=max(ans,c[0]);
        }
    }
    printf("%d",ans);
    return 0;
}