【jzoj 4486】最长公共子串（DP）（字符串）

最长公共子串

题目链接：jzoj 4486

题目大意

给你两个字符串。
其中有一个字符串你可以从给出的几个区间中选一个区间，改变区间内数的相对位置。
问你经过若干次操作后，两个字符串的最长公共子串。

思路

首先我们要发现一个性质：
如果有 $[x_1,y_1],[x_2,y_2]$ 两个区间，而且 $x_2⩽y_1$，那它就相当于 [ min ⁡ { x 1 , x 2 } , max ⁡ { y 1 , y 2 } ] [\min\{x_1,x_2\},\max\{y_1,y_2\}] [min{x1​,x2​},max{y1​,y2​}] 这样的区间。
（因为你操作可以无限次）

那就不难转化到只有字符串长度级别个区间，就可以 $O(n^2)$ 搞。
然后考虑不再区间的，那它其实就相当于一个长度为 $1$ 的区间，你这么搞可以使更好搞。

接着想想要怎么暴力。
枚举开始的段和字符串开始的地方，先往左匹配，然后再往右匹配。
如果可以一整段匹配掉就可以继续匹配下一段，否则就要结束。
（当然你继续匹配只要往左或往右有一个就可以了）

但复杂度是 $O(n^3)$，依然会炸。

那接着不难想到可以预处理，首先预处理你某一段从某个位置开始向左或向右匹配能匹配多长。（就只匹配这一段）
然后你跑的时候如果匹配的长度就是段的长度才继续递归跑。
而且不难看出状态只有 $n^2$ 个，会重复计算，直接上记忆化搜索。
（当然你也可以记搜换成 DP，但记搜好些就写了记搜）

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

struct qj {
	int l, r;
}a[100001], nw[2001];
int tn, sn, k, nl[2001], nr[2001], n, nn;
int f[2001][2001], ans, num[2001][31];
int l[2001][2001], r[2001][2001];
char t[2001], s[2001];

bool cmp(qj x, qj y) {
	return x.l < y.l;
}

int work(int x, int y) {//记忆化搜索（代替DP）
	if (x > n || y > tn) return 0;
	if (f[x][y]) return f[x][y];//记忆化
	if (r[x][y] != nw[x].r - nw[x].l + 1) return f[x][y] = r[x][y];
		else return f[x][y] = r[x][y] + work(x + 1, y + r[x][y]);
		//可以整段匹配就可以继续匹配下一段，否则就不可以直接结束
}

int main() {
//	freopen("lcs.in", "r", stdin);
//	freopen("lcs.out", "w", stdout); 
	
	scanf("%s", t + 1);
	tn = strlen(t + 1); 
	scanf("%s", s + 1);
	sn = strlen(s + 1);
	scanf("%d", &k);
	for (int i = 1; i <= k; i++) {
		scanf("%d %d", &a[i].l, &a[i].r);
		a[i].l++; a[i].r++;
	}
	
	sort(a + 1, a + k + 1, cmp);
	
	nn = 1;
	nl[nn] = a[1].l;
	nr[nn] = a[1].r;
	for (int i = 2; i <= k; i++)
		if (a[i].l <= nr[nn]) nr[nn] = max(nr[nn], a[i].r);
			else nl[++nn] = a[i].l, nr[nn] = a[i].r;
	
	//把每一段弄出来（每个段里面可以互换位置）
	if (nl[1] > 1)
		for (int i = 1; i <= nl[1] - 1; i++)
			nw[++n] = (qj){i, i};
	nl[nn + 1] = sn + 1;
	for (int i = 1; i <= nn; i++) {
		nw[++n] = (qj){nl[i], nr[i]};
		for (int j = nr[i] + 1; j <= nl[i + 1] - 1; j++)
			nw[++n] = (qj){j, j};
	}
	
	for (int i = 1; i <= n; i++) {
		for (int j = nw[i].l; j <= nw[i].r; j++)
			num[i][s[j] - 'a']++;
	}
	
	//预处理某个段与字符串某个地方开始向左 / 向右匹配的最长长度
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= tn; j++) {
			int now = 0;
			while (j + now <= tn && num[i][t[j + now] - 'a']) {
				num[i][t[j + now] - 'a']--;
				now++;
			}
			r[i][j] = now;
			for (int l = 0; l < now; l++)
				num[i][t[j + l] - 'a']++;
			
			now = 0;
			while (j - now >= 1 && num[i][t[j - now] - 'a']) {
				 num[i][t[j - now] - 'a']--;
				 now++;
			}
			l[i][j] = now;
			for (int l = 0; l < now; l++)
				num[i][t[j - l] - 'a']++;
		}
	
	for (int i = 1; i <= n + 1; i++)
		for (int j = 1; j <= tn + 1; j++)
			ans = max(ans, l[i - 1][j - 1] + work(i, j));//记得左边也可以搞一些
	
	printf("%d", ans);
	
	fclose(stdin);
	fclose(stdout);
	
	return 0;
}