7-2 破译报文（PTA 字符串实验）

OneGrave

已于 2023-11-05 14:57:06 修改

阅读量1.9k

点赞数 4

文章标签：哈希算法散列表算法

于 2023-09-29 00:04:32 首次发布

本文链接：https://blog.csdn.net/OneGrave/article/details/133396346

版权

7-2 破译报文

分数 80

作者朱允刚

单位吉林大学

小明接到一个破解报文的任务：该报文是一串文本，破解出的密文应是在报文串中出现次数大于1的所有子串中的最长者。规定报文本身不能称为自己的子串。请编写效率尽可能高的程序帮小明完成这个棘手的任务。

输入格式:

输入为一个字符串，表示报文，包含不超过10000个字母。

输出格式:

输出为一个整数，表示破解出的密文串的长度。

输入样例1:
xabceabcf
输出样例1:
3

具体算法如下：

首先定义一个哈希函数Hash，用于计算字符串的哈希值。哈希值的计算采用了前缀哈希的方式，使用了一个预先设定的质数k和一个模数mod。对于字符串中的每个位置，计算其前缀和的哈希值。

字符串前缀哈希法，把字符串变成一个p进制数字（哈希值），实现不同的字符串映射到
如 X1X2X3⋯Xn−1Xn 的字符串,采用字符的ascii 码乘上 P 的次方来计算哈希值。相当于将字符串看做一个 P 进制的数字，每个位上的数字值为 ascil码的值。其中 $p_{i}$ 代表 $k$ 的 $i$ 次方

   $\left( X_1 \cdot P_{n-1} + X_2 \cdot P_{n-2} + \ldots + X_{n-1} \cdot P_1 + X_n \cdot P_0 \right) \mod Q$
s为前缀和数组，str为字符串数组
前缀和公式   $s[i+1]=s[i]\times k+ str_{i+1}$
区间和公式 $s[l,r]=s[r]-s[l-1]\times k^{r-l+1}$
定义一个辅助函数check，用于检查长度为x的子串是否在报文串中出现了多次。该函数使用一个哈希表st来记录已经出现过的子串的哈希值。遍历报文串的所有长度为x的子串，计算其哈希值并判断是否在哈希表中已经存在。如果存在，则说明该子串出现了多次，返回true；否则，将该子串的哈希值添加到哈希表中，并继续遍历下一个子串。最后，如果没有找到重复出现的子串，则返回false。使用前缀和，我们可将查询一个子串的哈希值的时间缩减为 $O(1)$ 。
主函数solve中，首先读入报文串并计算其长度。然后，使用前缀哈希的方式计算字符串的哈希值，并保存在数组s中。同时，计算每个位置的权值p。接下来，使用二分查找的方式找到报文串中出现次数大于1的所有子串中的最长者。初始时，将左边界l设为0，右边界r设为报文串的长度。我们不难看出，答案是具有单调性质的。如果长度为 $x$ 的子串能被找出两次以上，那么长度为 $x-1$ 的子串也能被找出两次以上。

时间复杂度 $O(nlogn)$ ，空间复杂度 $O(n)$

#define int long long
#define endl '\n'
#define mem(a,b) memset(a,b,sizeof(a))
#define rep(a, b, c) for (a = b; a <= c; a++)
using namespace std;

typedef long long ll;
typedef double db;
typedef pair<int,int> pii;
const int N = 10010, mod = 1e9 + 7;

int k = 100007;
char str[N];
int s[N],p[N];
int n;
unordered_map <int,bool> st;

int Hash(int l,int r){
    int t = (s[r] - (s[l-1] * p[r-l+1] % mod) + mod) % mod;
    return t;
}

bool check(int x){
    int i;
    st.clear();
    
    rep(i,1,n-x+1){
        int t = Hash(i,i+x-1);
        // cout << ' ' << i << ' ' << i+x-1 <<" HASH is " << t << endl;

        if (st[t]) return 1;
        st[t] = 1;
    }
    return 0;
}

void solve(){
    cin >> str+1;
    n = strlen(str+1);

    int i;
    
    p[0] = 1;
    rep(i,1,n){
        s[i] = (s[i-1] * k + str[i]) % mod;
        p[i] = p[i-1] * k % mod;
    }   

    int l = 0,r = n;

    while (l + 1 != r){
        int mid = l + r >> 1;
        if (check(mid)){
            l = mid;
        }else{
            r = mid;
        }
    }

    cout << l << endl;
}