hdu6774 String Distance

hdu6774
因为只能插删,所以最优的肯定是留下最长公共子序列。
dp i j 表示到b串的第i个位置,长度为j的子序列,在a串的最小位置。
序列自动机预处理一下a串,然后每次询问做一遍m^2的dp即可。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e5 + 10;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
char a[N],b[N];
int nx[N][26],t,f[22][22];
int main() {
  ios::sync_with_stdio(false);
  for(cin>>t;t;t--){
    cin>>a+1>>b+1;
    int n=strlen(a+1);
    int m=strlen(b+1);
    for(int i=0;i<26;i++)nx[n+1][i]=0;
    for(int i=n;i>=1;i--){
      for(int j=0;j<26;j++)nx[i][j]=nx[i+1][j];
      nx[i][a[i]-'a']=i;
    }
    int Q;
    for(cin>>Q;Q;Q--){
      int l,r;
      cin>>l>>r;
      for(int i=0;i<=m;i++){
        for(int j=0;j<=m;j++){
          f[i][j]=-1;
        }
      }
      for(int i=0;i<=m;i++)f[i][0]=l-1;// 前i个字符 匹配了j位子序列 的最短 a串位置
      int y=0;
      for(int i=1;i<=m;i++){
        for(int j=i-1;j>=0;j--){
          int x=f[i-1][j];
          f[i][j+1]=f[i-1][j+1];
          if(x==-1)continue;
          if(nx[x+1][b[i]-'a']&&nx[x+1][b[i]-'a']<=r){
            if(f[i][j+1]==-1)f[i][j+1]=nx[x+1][b[i]-'a'];
            else f[i][j+1]=min(f[i][j+1],nx[x+1][b[i]-'a']);
            y=max(y,j+1);
          }
        }
      }
      cout<<r-l+1+m-2*y<<'\n';
    }
  }
     return 0;
}

猜你喜欢

转载自blog.csdn.net/qq_40655981/article/details/107570361
今日推荐