[jzoj5635][CF578D]【NOI2018模拟4.7】LCS

三倍经验题当然要优先做啦~~~

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做法

考虑设 $F_{i,j}​$ 表示T串前i项和S串前j项最长公共子串长度

可以发现 状态有效当且仅当 |i-j|<=1或 $min(i,j)-1<=F_{i,j}<=max(i,j)+1$

于是状态数就很少很少

另设$G_{i,st}$ 其中st表示 ${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;
}