Jzoj4889 最长公共回文子序列

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

以后写题解都开大字体

嗯，看到数据范围发现可以做呢

枚举b每一个子序列并在A中匹配，复杂度2^m*n

明显T，怎么加速匹配呢？

一个SAM不就行了嘛

这个可是子序列，搞一个数组就可以了，存下每个字符的位置让后二分就好了

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
char a[100010],b[30];
int f[26][100010]={0},t[26]={0},n,m;
int w[30],c=0,ans=0;
bool ok(int x){
	for(int i=c=0;x;i++,x>>=1) if(x&1) w[c++]=i;
	for(int i=0;(i<<1)<c;++i) if(b[w[i]]!=b[w[c-i-1]]) return 0;
}
int cal(){
	int pos=-1;
	for(int i=0;i<c;++i){
		int j=b[w[i]]-'a',l=0,r=t[j];
		for(int M;l<r;){
			M=l+r>>1;
			if(f[j][M]<=pos) l=M+1;
			else r=M;
		}
		if(l==t[j]) return 0;
		pos=f[j][l];
	}
	return c;
}
int main(){
	freopen("lcps.in","r",stdin);
	freopen("lcps.out","w",stdout);
	scanf("%s%s",a,b);
	n=strlen(b); m=strlen(a);
	for(int i=0;i<m;++i)
		f[a[i]-'a'][t[a[i]-'a']++]=i;
	for(int i=0;i<26;++i) f[i][t[i]]=m+1;
	for(int i=1;i<(1<<n);++i)
		if(ok(i)) ans=max(ans,cal());
	printf("%d
",ans);
}