「6月雅礼集训 2017 Day1」说无可说

【题目大意】

给出n个字符串，求有多少组字符串之间编辑距离为1~8。

n<=200，∑|S| <= 10^6

【题解】

首先找编辑距离有一个n^2的dp，由于发现只找小于等于8的，所以搜旁边16个状态即可。

复杂度O(n^2|S| * 16)

# include <vector>
# include <stdio.h>
# include <iostream>
# include <string.h>
# include <algorithm>

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;

const int N =  + , M = 1e6 + , AN = ;
const int mod = 1e9 + ;

int n, len[M], ans[AN];
vector<char> a[N];
char str[M];

int u, v;
int dp[][M];
inline int f(int i, int j) {
    if(i == ) return j;
    if(j == ) return i;
    if(j < max(, i-) || j > min(len[v], i+)) return 1e9;
    return dp[i&][j];
}

inline void fin(int i, int j, int s) {
    dp[i&][j] = s;
}

inline int gans() {
    for (int i=; i<=len[v] || i<=len[u]; ++i) fin(, i, ), fin(, i, );
    for (int i=; i<=len[u]; ++i)
        for (int j=max(, i-), jto = min(i+, len[v]); j<=jto; ++j) {
            fin(i, j, min(f(i-, j), f(i, j-)) + );
            if(a[u][i-] != a[v][j-]) fin(i, j, min(f(i, j), f(i-, j-) + ));
            else fin(i, j, min(f(i, j), f(i-, j-)));
        }
    return f(len[u], len[v]);
}

int main() {
//    freopen("say.in", "r", stdin);
//    freopen("say.out", "w", stdout);
    cin >> n;
    for (int i=; i<=n; ++i) {
        scanf("%s", str);
        len[i] = strlen(str);
        for (int j=; str[j]; ++j) a[i].push_back(str[j]);
    }
    for (int i=, tem; i<=n; ++i) {
        for (int j=i+; j<=n; ++j) {
            u = i, v = j;
            tem = gans();
            if(tem >=  && tem <= ) ++ ans[tem];
        }
    }

    for (int i=; i<=; ++i) cout << ans[i] << ' ';
    cout << endl;
    return ;
}

当然，过不了，有60分。

所以考虑dp有很多空余状态，改成dfs加剪枝就可以过了。。

# include <vector>
# include <stdio.h>
# include <iostream>
# include <string.h>
# include <algorithm>

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;

const int N =  + , M = 1e6 + , AN = ;
const int mod = 1e9 + ;

int n, len[M], ans[AN];
vector<char> a[N];
char str[M];

# define ABS(x) ((x) >  ? (x) : -(x))

int tem;
inline void solve(int u, int v, int cur_u, int cur_v, int cnt) {
    if(cnt + ABS(len[v] - len[u] - (cur_v - cur_u)) >= tem) return;
    while(cur_u < len[u] && cur_v < len[v]) {
        if(a[u][cur_u] != a[v][cur_v]) {
            solve(u, v, cur_u+, cur_v, cnt+);
            solve(u, v, cur_u, cur_v+, cnt+);
            solve(u, v, cur_u+, cur_v+, cnt+);
            return ;
        }
        ++cur_u, ++cur_v;
    }
    if(cur_u == len[u]) tem = min(tem, cnt + len[v] - cur_v);
    if(cur_v == len[v]) tem = min(tem, cnt + len[u] - cur_u);
}

int main() {
//    freopen("say.in", "r", stdin);
//    freopen("say.out", "w", stdout);
    cin >> n;
    for (int i=; i<=n; ++i) {
        scanf("%s", str);
        len[i] = strlen(str);
        for (int j=; str[j]; ++j) a[i].push_back(str[j]);
    }
    for (int i=; i<=n; ++i) {
        for (int j=i+; j<=n; ++j) {
            tem = ;
            solve(i, j, , , );
            ++ ans[tem];
        }
    }

    for (int i=; i<=; ++i) cout << ans[i] << ' ';
    cout << endl;

    return ;
}

这里dfs的时候需要注意一点，要一段一段跳，不能一格一格跳，会T。。