单词检索(search)

\(Description\)

小可可是学校图书馆的管理员，现在他接手了一个十分棘手的任务。
由于学校需要一些材料，校长需要在文章中检索一些信息。校长一共给了小可可N篇文章，每篇文章为一个字符串。现在，校长需要他找到这样的单词，它至少在这N篇文章中的M篇文章里出现过，且单词长度为L。可是，工作量十分庞大，但校长又急需小可可完成这项任务。
现在他向你求助，需要你编写程序完成这项艰巨的任务。

\(Input\)

第1行3个正整数N，M，L，表示文章的数目，单词至少出现在M篇文章中和每个单词的长度。
接下来N行，每行一个字符串，表示一篇文章。

\(Output\)

仅一行，表示满足检索条件的单词数。

\(Sample Input\)

3 2 2
noip
istudycpp
imacppstudent

\(Sample Output\)

5
【样例解释】
这5个单词分别为：st，tu，ud，pp，cp。

\(Data Constraint\)

对于20%的数据有1≤N,M≤10；
对于60%的数据有1≤N,M≤100；
对于100%的数据有1≤N,M≤2000，L≤1000。每篇文章长度不大于1000，均有小写字母组成。

思路

字符串？
既然限定了答案字符的长度，那很容易想到分离字符串，统计出现个数
那就是 \(O(nl)\) 级别的个数，空间可以承受
然后怎么存分离出的字符
如果可以
对于每一篇文章，我们把分离的字符去重（因为答案要求的字符至少在N篇文章中的M篇文章里出现过）
最后合在一起
此时统计出现次数超过 \(M\) 的字符数量就是答案
想想，没有问题
关键是分离
如果把每个字符截取下来，那么一片文章处理的复杂度就是 \(O(|S|l)\)
\(n\) 篇就是 \(O(n|S|l)\) 的，显然不可以承受
然后呢？

正解

既然有了思路，那就直接打吧（也有 \(60\) 分了）
可是，现在是追求正解的时候
思路中的瓶颈在于分离字符
而我们知道要获取一段区间的信息，如果满足可加性，那就可以前缀和维护
字符串？
哦，哈希！

具体方式

对于一个字符串，我们对其进行哈希时，考虑映射成一个数字
一个很巧妙的方法是：把字符本身当做一串数字，一个字符对应权值
从高到低每一位的位权就参考十进制的方法，记作 \(B^{k-1}\)
其中 \(k\) 为第 \(k\) 位，\(B\) 是几进制
因为题中字符串均由小写字母构成
所以设为 \(26\) 进制（不这样也行，哈希映射的方法自己选择）
那么对于一个字符串S，有：

\[ H(S) = \sum_{i=1}^{|S|} c_iB^{|S|-i} \]

截取 S 中后 \(l\) 个字符组成的字符串设为 C‘，前 \(|S| - l\) 个字符组成的字符串设为 C
则

\[ H(C') = H(S) - H(C) * B^l \]

于是维护 \(s_i =\sum_{k=1}^i H(|1...k|)\)

双模数，以防万一

代码

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL p1 = 1e9 + 7 , p2 = 1e9 + 9 , B = 26;
LL s1[1005] , s2[1005] , b1[1005] , b2[1005];
struct node{
    LL x , y;
}d[2000005] , c[2000005];
int cnt , tot , n , m , l , ans;
char ch;

inline bool cmp(node a , node b){return (a.x < b.x || (a.x == b.x && a.y < b.y));}

int main()
{
//  freopen("单词检索.in" , "r" , stdin);
    scanf("%d%d%d" , &n , &m , &l);
    b1[0] = b2[0] = 1;
    for(register int i = 1; i <= 1001; i++) b1[i] = b1[i - 1] * B % p1;
    for(register int i = 1; i <= 1001; i++) b2[i] = b2[i - 1] * B % p2;
    int x;
    for(register int i = 1; i <= n; i++)
    {
        x = 0 , s1[0] = 0 , s2[0] = 0 , cnt = 0;
        ch = getchar();
        while (ch < 'a' || ch > 'z') ch = getchar();
        for(; ch >= 'a' && ch <= 'z'; ch = getchar())
        {
            x++;
            s1[x] = (s1[x - 1] * B % p1 + (ch - 'a' + 1)) % p1;
            s2[x] = (s2[x - 1] * B % p2 + (ch - 'a' + 1)) % p2;
            if (x >= l)
            {
                d[++cnt].x = (s1[x] + p1 - s1[x - l] * b1[l] % p1) % p1;
                d[cnt].y = (s2[x] + p2 - s2[x - l] * b2[l] % p2) % p2;
            }
        }
        sort(d + 1 , d + cnt + 1 , cmp);
        c[++tot] = d[1];
        for(register int j = 2; j <= cnt; j++)
        if (d[j].x != d[j - 1].x || d[j].y != d[j - 1].y)
            c[++tot] = d[j];
    }
    sort(c + 1 , c + tot + 1 , cmp);
    x = 1;
    while (x <= tot)
    {
        int s = 1;
        while (c[x].x == c[x + 1].x && c[x].y == c[x + 1].y && x < tot) x++ , s++;
        if (s >= m) ans++;
        x++;
    }
    printf("%d" , ans);
}

原文地址：https://www.cnblogs.com/leiyuanze/p/12368625.html