最长公共回文子序列 【NOIP2016提高A组集训第14场11.12】

Description

（直接复制了，偷个懒O(∩_∩)O~）

YJC最近在学习字符串的有关知识。今天，他遇到了这么一个概念：最长公共回文子序列。一个序列S，如果S是回文的且分别是两个或多个已知序列的子序列，且是所有符合此条件序列中最长的，则 S 称为已知序列的最长公共回文子序列。YJC很聪明，他很快就学会了如何求最长公共回文子序列。他现在想把问题规模扩大一些，于是他提出了这么一个问题：给一个长度为n(1≤n≤100000)的字符串a和一个长度为m(1≤m≤20)的字符串b，求a和b的最长公共回文子序列的长度。YJC发现他不会做了，于是他来问你这个问题的答案。

Input

第一行包含一个字符串，表示字符串a。
第二行包含一个字符串，表示字符串b。

Output

一行，包含一个整数，表示a和b的最长公共回文子序列的长度。

Sample Input

aaaba
aabba

Sample Output

3

Data Constraint

对于50%的数据，满足1≤n≤100。
对于100%的数据，满足1≤n≤100000，1≤m≤20，字符串只包含小写字母。

Hint

最长公共回文子序列是aaa，长度为3。最长公共子序列是aaba，但它不是回文的。

---

思路

当时真没什么想法，看着m这么小，就是想着暴力，枚举b串的回文串，然后在a串中寻找，问题是如何搜索才能在a串中迅速找出？旁边有人说，直接输出b串中回文串的最大长度理论上有30分
然后我就想到kmp（入魔），自然就想到了next数组，懵逼。觉得好像可以用类似思想记录字符的位置的情况。
用一个模拟链表，记录a串里，每一个字母的位置情况，也就是每一个字母的位置向下一个相同字母的位置连一条边，并记录这个字母在a串中的位置。
然后暴力b串中回文串，在a里找。
注意：由于模拟链表记录时是记录最后一条连出的边，所以在a里找时，必须把当前暴力的回文串倒过来搜。
然后我一个煞笔错误少了10分QwQ。
这个10分的答案是m的最大值20，还好我没有听那个逗逼说直接输出。
时间复杂度$2^m*m$

---

代码

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)

using namespace std;

const int maxn=1e5+5,maxm=150;
int fre[maxm],next[maxn],go[maxn],a[maxn],ans,n,m,num,now[25],way[maxm];
char s1[maxn],s2[26];
bool bz[25];

void add(int x,int y)
{
    go[++num]=y;
    next[num]=fre[x];
    fre[x]=num;
}
void check(int p)
{
    bool bo=true;
    fo(i,1,p/2)
    if (now[i]!=now[p-i+1]) {
        bo=false;
        break;
    }
    if (!bo) return;
    bo=true;
    memset(way,255,sizeof(way));
    now[p+1]=0; way[0]=0;
    for(int i=p;i>=1;i--){
        int val=now[i];
        int u;
        if (way[val]==-1) u=fre[val];
        else u=next[way[val]];
        while (u&&go[u]>=go[way[now[i+1]]]) u=next[u];
        if (u) way[val]=u;
        else {
            bo=false;
            break;
        }
    }
    if (bo) 
        ans=p;
}
void dfs(int x,int d)
{
    if (x>ans) check(x-1);//比赛时这一句和下一句打反了QwQ.
    if (x>m) return;
    fo(i,d,m-1)
    if (!bz[i]){
        now[x]=s2[i];
        bz[i]=true;
        dfs(x+1,i+1);
        bz[i]=false;
        now[x]=0;
    }
}
int main()
{
    freopen("lcps.in","r",stdin);
    freopen("lcps.out","w",stdout);
    scanf("%s",&s1);
    scanf("%s",&s2);
    memset(go,255,sizeof(go)); 
    n=strlen(s1); m=strlen(s2);
    fo(i,0,strlen(s1)-1) add(s1[i],i);
    go[0]=1e5+3;
    dfs(1,0);
    printf("%d",ans);
}