字符串哈希-k-substrings-CodeForces-961F

字符串哈希-k-substrings-CodeForces-961F

题意：

$$\begin{gathered} 给定一个长度为n的字符串s。 \\ 定义k串为s中，第k个字符到倒数第k个字符构成的子串。 \\ 抽象地，即sub{s}_k=s[k,n+1-k]。 \end{gathered}$$

$$\begin{gathered} 对长度为n的字符串而言，共有\lceil \frac{n}{2}\rceil 个k串，分别是sub{s}_1,sub{s}_2,...,sub{s}_{\frac{n}{2}}。 \\ 显然，sub{s}_1=s[1,n]=s。 \end{gathered}$$

$求每一个k串中，最长相同前后缀的长度。$

Sample Input:
15
bcabcabcabcabca

Sample Output:
9 7 5 3 1 -1 -1 -1


Sample Input:
24
abaaabaaaabaaabaaaabaaab

Sample Output:
15 13 11 9 7 5 3 1 1 -1 -1 1


Sample Input:
19
cabcabbcabcabbcabca

Sample Output:
5 3 1 -1 -1 1 1 -1 -1 -1

Note:
The answer for first sample test is folowing:

1-substring: bcabcabcabcabca
2-substring: cabcabcabcabc
3-substring: abcabcabcab
4-substring: bcabcabca
5-substring: cabcabc
6-substring: abcab
7-substring: bca
8-substring: c

数据范围：

$s的长度n\in [2,1000000]。$

题解:

$从短到长，枚举每一个k串，再从大到小，枚举最大相同前后缀的长度，用哈希来判断是否相等。$

$$\begin{gathered} 有一个性质需要注意： \\ 设ans[i]表示第i长的k串的最大相同前后缀的长度，则ans[i]<=ans[i+1]+2。 \\ 这很好理解，枚举k串时每次只在前一个k串的首尾增加2个字符， \\ 最好情况就是增加的2个字符恰好在原来的基础上把相同前后缀拓宽了2个字符，如样例一的解释。 \\ 那么第i个k串的答案范围就在[-1,ans[i+1]+2]之间来取。 \end{gathered}$$

具体落实：

$$\begin{gathered} ①、首先根据n的奇偶性，解决最短的k串的即最后一个k串的相同前后缀的长度。 \\ 若n为奇数，最短的k串长度为1，那么ans[\lceil \frac{n}{2}\rceil ]=-1， \\ 若n为偶数，最短的k串长度为2，判断s[\frac{n}{2}]与s[\frac{n}{2}+1]是否相等，若相等，则ans[\frac{n}{2}]=1,否则ans[\frac{n}{2}]=-1。 \\ ②、以2为步长，拓宽k串的长度，对第i个k串，在第i+1个k串的ans[i+1]+2到-1之间枚举答案ans[i]。 \end{gathered}$$

---

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<map>
#define ll long long
using namespace std;
const int N=2e6+10;
const int base=131;
const int mod=1e9+7;
char s[N];
int n,ans[N/2];
ll h[N],p[N];

ll get(int l,int r)
{
    return (h[r]-h[l-1]*p[r-l+1]%mod+mod)%mod;
}

int main()
{
    scanf("%d",&n);
    scanf("%s",s+1);

    p[0]=1;
    for(int i=1;i<=n;i++) p[i]=p[i-1]*base%mod;

    for(int i=1;i<=n;i++) h[i]=(h[i-1]*base%mod+s[i]-'a')%mod;

    ///最中间的串，考虑奇偶性
    int l,r;
    int mid=n+1>>1;
    if(n%2==0)
    {
        l=n>>1,r=l+1;
        if(s[l]==s[r]) ans[mid]=1;
        else ans[mid]=-1;
    }
    else ans[mid]=-1,l=r=mid;

    for(int i=mid-1;i>=1;i--)
    {
        l--,r++;  ///k串的左右端点
        for(int j=ans[i+1]+2;j>=-1;j-=2)  ///枚举最大匹配的长度
        {
            if(get(l,l+j-1)==get(r-j+1,r)||j==-1)
            {
                ans[i]=j;
                break;
            }
        }
    }

    for(int i=1;i<=mid;i++) printf("%d ",ans[i]);
    puts("");

    return 0;
}