HDU - 6774 String Distance

传送门：HDU - 6774

思路分析

实际上就是求询问串和$t$串的$lcs$
如果按照普通的求$lcs$的方法，肯定是不行的
$dp[i][j]$表示在t串的长度为$i$的前缀里，匹配到长度为$j$的$lcs$的最小右端点
可以通过序列自动机预处理来优化匹配的过程，这样匹配的复杂度是$O(m*m)$的

AC代码

#include <bits/stdc++.h>
#define fi first
#define se second
#define ll long long
#define pb push_back
#define mp make_pair
#define fun function
#define sz(x) (x).size()
#define lowbit(x) (x)&(-x)
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))

namespace FastIO {
#define BUF_SIZE 100000
#define OUT_SIZE 100000
    bool IOerror=0;
    inline char nc() {
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE;
        if(p1==pend) {
            p1=buf;
            pend=buf+fread(buf,1,BUF_SIZE,stdin);
            if(pend==p1) {
                IOerror=1;
                return -1;
            }
        }
        return *p1++;
    }
    inline bool blank(char ch) {
        return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';
    }
    template<class T> inline bool read(T &x) {
        bool sign=0;
        char ch=nc();
        x=0;
        for(; blank(ch); ch=nc());
        if(IOerror)return false;
        if(ch=='-')sign=1,ch=nc();
        for(; ch>='0'&&ch<='9'; ch=nc())x=x*10+ch-'0';
        if(sign)x=-x;
        return true;
    }
    template<class T,class... U>bool read(T& h,U&... t) {
        return read(h)&&read(t...);
    }
#undef OUT_SIZE
#undef BUF_SIZE
};
using namespace std;
using namespace FastIO;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

const int INF = 0x3f3f3f3f;
const int N = 1e5+10;

int nex[N][30],dp[30][30];
char s[N],t[N];

int main() {

#ifdef xiaofan
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
#endif

    int T;
    cin>>T;
    while(T--) {
        cin>>s+1;
        cin>>t+1;
        int n=strlen(s+1);
        int m=strlen(t+1);

        for(int j=0; j<26; j++) {
            nex[n][j]=n+1;
        }
        nex[n][s[n]-'a']=n;
        for(int i=n-1; i>=1; i--) {
            for(int j=0; j<26; j++) {
                nex[i][j]=nex[i+1][j];
            }
            nex[i][s[i]-'a'] = i;
        }
        
        int q;
        cin>>q;
        while(q--){
            int l,r;
            cin>>l>>r;
            mem(dp,INF);
            for(int i=1;i<=m;i++) dp[i][1]=nex[l][t[i]-'a'];
            for(int i=1;i<=m;i++){
                for(int j=1;j<=i;j++){
                    if(dp[i-1][j]<=r) dp[i][j]=min(dp[i][j],dp[i-1][j]);
                    if(dp[i-1][j-1]<r) dp[i][j]=min(dp[i][j],nex[dp[i-1][j-1]+1][t[i]-'a']);
                    
                }
            }
            int ma=0;
            dp[m][0]=0;
            for(int i=m;i>=0;i--){
                if(dp[m][i]<=r) {
                    ma=i;
                    break;
                    
                }
            }
            int del=r-l+1-ma;
            int add=m-ma;
            cout<<del+add<<endl;
        }
    }




    return 0;
}