三倍经验题当然要优先做啦~~~
jzoj5636 CF578D
Solution
T可以看成S去掉一个字母再加一个字母.
首先把S分成一些相同字母组成的段(假设有k段),去掉这些段中的任意一个字母是等价的.
之后我们可以在任何一个位置添加任何字母,除非和原来的S一模一样.
那么方案数就是k(nm-n).
然而,有些T会被计算2次,这些串都有子串长成ababa..aba这样(只要长度大于等于2),它和bab..bab是一致的.(注:对于字符串ab 可以删去a得到 ba 也可以删去b得到 ba 对于字符串qaq 也有类似性质)容易发现 冗余状态只有这一种 尝试统计形如“aba…bab”的字符串个数即可
以上题解做法 这里讲一种DP套DP做法
考虑设 Fi,j F i , j 表示T串前i项和S串前j项最长公共子串长度
可以发现 状态有效当且仅当 |i-j|<=1或 min(i,j)−1<=Fi,j<=max(i,j)+1 m i n ( i , j ) − 1 <= F i , j <= m a x ( i , j ) + 1
于是状态数就很少很少
另设 Gi,st G i , s t 其中st表示 Fi F i 状态
每次把 st 强制转出来进行转移 再转回去即可
Code
#include <cmath>
#include <ctime>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define oo 2139062143
#define fo(i,x,y) for (ll i=(x),asdfgh=(y);i<=asdfgh;++i)
#define fd(i,x,y) for (ll i=(x),asdfgh=(y);i>=asdfgh;--i)
#define fi first
#define se second
using namespace std;
typedef double db;
typedef long long ll;
typedef pair<ll,ll> PI;
ll abs(ll x) {return(((x)>=0)?(x):(-(x)));}
ll max(ll x,ll y) {return(((x)>(y))?(x):(y));}
ll min(ll x,ll y) {return(((x)<(y))?(x):(y));}
ll lowbit(ll x) {return((x)&(-x));}
const ll N=100100,M=26;
ll n,m;
ll f[2][4];//[i][i-1] [i][i] [i][i+1]
ll g[N][8],st[N];
ll mi(ll o)
{
ll rt=0;
fo(i,1,3)
{
if(f[o][i]) rt|=(1<<(i-1))*f[o][i];
}
return rt;
}
void je(ll o,ll cd)
{
ll rt=cd;
fo(i,1,3) f[o][i]=(bool)(rt&(1<<(i-1)));
}
bool check(ll o)
{
fo(i,1,3)
if(f[o][i]<0) return (0);
return 1;
}
int main()
{
freopen("lcs.in","r",stdin);
freopen("lcs.out","w",stdout);
//freopen("D:/LiuYuanHao/b.in","r",stdin);
scanf("%lld%lld\n",&n,&m);
fo(i,1,n)
{
char ch=getchar();while(!(ch>='a'&&ch<='z')) ch=getchar();
st[i]=ch-'a'+1;
}
//1 1/0 1/0
f[0][1]=1;
fo(ch,1,m)
{
if(ch==st[1]) f[0][2]=1,f[0][3]=1;
else if(ch==st[2]) f[0][2]=0,f[0][3]=1;
else f[0][2]=f[0][3]=0;
g[1][mi(0)]++;
}
fo(i,1,n-1)
{
fo(ch,1,m)/*T[1]*/
fo(s,0,7)/*i*/
if(g[i][s])
{
memset(f[1],255,sizeof f[1]);
je(0,s);
if(ch==st[i])f[1][1]=f[0][1]/*+1*/;
f[1][1]=max(f[1][1],f[0][2]);
if(ch==st[i+1])f[1][2]=f[0][2];
f[1][2]=max(f[1][2],max(f[1][1]-1,f[0][3]-1));
if(ch==st[i+2]) f[1][3]=f[0][3];
f[1][3]=max(f[1][3],f[1][2]);
ll tmp=0;
if(check(1))
g[i+1][tmp=mi(1)]+=g[i][s];
}
}
ll ans=0;
fo(s,0,7)
{
je(0,s);
if(f[0][2]!=0) continue;
ans+=g[n][s];
}
printf("%lld",ans);
return 0;
}