P3426 [POI2005]SZA-Template

洛谷博客链接

P3426 [POI2005]SZA-Template - 传送门

题意

给你一字符串，需要你制作一个印章，能用它盖出该字符串（当然不能多盖）。

求最小印章长度。

盖印章：

1. 同一位置，一样字符可以重复盖；

2. 同一位置，不同字符不能重复盖。

例如：

字符串为 ababbababbabababbabababbababbaba，印章为 ababbaba。

盖法如下：

ababbababbabababbabababbababbaba
ababbaba
     ababbaba
            ababbaba
                   ababbaba
                        ababbaba

引理

我们用 $f_i$ 表示能覆盖长度为 i 的前缀的最小印章长度。

用数组 nxt 表示 KMP 算法中的 kmp 数组。

具体思路、后续求法将在思路部分细说，此处先放引理及证明。

Part 1 - $f_i$ 的取值情况

一、$f_i=i$

这个不难理解，当 $f_i=i$ 时当然覆盖完前缀 i 了。

二、$f_i=f_{nx{t}_i}$

1、不难得出 $f_i\le i$。

接下来我们只开始讨论 $f_i<i$ 的情况。

（$f_i=i$ 见引理一。）

---

2、易得 $f_i\le nx{t}_i$。

证明如下：

若 $f_i>nx{t}_i$，为满足题意，我们盖印章的时候必须在最前面和最后面分别盖一个。

例如：

ababbababbabababbabababbababbaba
ababbaba
                        ababbaba

此时，$f_i$ 即为 i 的相同前后缀，且比 $nx{t}_i$ 还长（已假设 $f_i>nx{t}_i$）。

这与 $nx{t}_i$ 为 i 最长相同前后缀的定义相矛盾，

故，$f_i$ 只能小于等于 $nx{t}_i$。

有了这个引理，我们就可以分成两种情况，分别见 3 和 4。

---

3、当 $f_i<nx{t}_i$ 时，$f_i=f_{nx{t}_i}$。

1. 若 $f_i<f_{nx{t}_i}$，则无法用 $f_i$ 填充前 $nx{t}_i$ 个字符，不然 $f_{nx{t}_i}$ 会更小。

   所以，$f_i\ge f_{nx{t}_i}$。

2. 若 $f_i>f_{nx{t}_i}$，则可以用 $f_{nx{t}_i}$ 去填充 $f_i$（$f_i<nx{t}_i$），

   与 $f_i$ 为最短能填充 i 的定义相矛盾。

   所以，$f_i\le f_{nx{t}_i}$。

综上，得 $f_i=f_{nx{t}_i}$。

---

4、当 $f_i=nx{t}_i$ 时，$f_i=f_{nx{t}_i}$。

因为 $f_i$ 的定义，即不可以再出现比 $nx{t}_i$ 还能短且能填充掉 $nx{t}_i$ 的串。

所以此时，$f_{nx{t}_i}=nx{t}_i=f_i$，

即 $f_i=f_{nx{t}_i}$。

---

综上所述，$f_i$ 的取值只可能是 i 或者 $f_{nx{t}_i}$。

---

Part 2 - $f_i$ 什么情况下可以取 $f_{nx{t}_i}$

只要我们在 $[i-nx{t}_i,i]$ 中发现一个 j，使得 $f_j=f_{nx{t}_i}$，

此时就可以让 $f_i\leftarrow f_{nx{t}_i}$。

ababbaba bbabababbabababb ababbaba
ababbaba
                          ababbaba

例如当前字符串 i 的最长公共前后缀是 ababbaba ——也就是我空格出的部分。

ababbababbabababbabababb ababbaba
ababbaba
                         ababbaba

按照刚才所说，

我们需要在后面的 ababbaba （最长公共前后缀中的后缀）中找到 j。

此时，

$f_{nx{t}_i}$ 可以覆盖 $[nx{t}_i,i]$，$f_j$ 可以覆盖 $[1,j]$；

又因为 $j\ge i-nx{t}_i$；

还 $f_j=f_{nx{t}_i}$；

所以，

我们就可以用 $f_{nx{t}_i}$ 覆盖掉整个字符串 i 了。

思路

我们用了 dp 的思路去求解了最小印章长度。

用 KMP 去找最长相同前后缀。

由引理，我们得知：

循环每次求给出字符串的长度为 i 的前缀的 $f_i$ 值：

初始化 $f_i\leftarrow i$；

若出现了 j（$f_j=f_{nx{t}_i}$，$j\ge i-nx{t}_i$），那么 $f_i\leftarrow f_{nx{t}_i}$。

怎么去找 j？

我们开个桶去记录，每次印章长度为 $f_i$ 时最多能覆盖到的字符的下标即可。

具体见代码。

代码

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

const int maxn = 500005;
int n;
char a[maxn];
int kmp[maxn], b[maxn];
int f[maxn];

int main ()
{
	scanf ("%s", a + 1);
	n = strlen (a + 1);
	for (int i = 2, j = 0; i <= n; i++)
	{
		while (j and a[j + 1] != a[i]) j = kmp[j];
		if (a[j + 1] == a[i]) j++;
		kmp[i] = j;
	}
	int j = 0;
	for (int i = 1; i <= n; i++)
	{
		f[i] = i;
		if (b[f[kmp[i]]] >= (i - kmp[i])) f[i] = f[kmp[i]];
		b[f[i]] = i;
	}
	printf ("%d\n", f[n]);
	return 0;
}

算是 KMP 的 kmp 数组的应用，就是思维难度确实挺大。

kmp 数组理解深度 ++;

——$End$——

By 531