您的位置:首页 > 其它

【hdu 3518】Boring counting

2017-10-04 18:44 399 查看

【链接】h在这里写链接


【题意】


给出一个字符串,求出至少不重叠出现2次以上的子串有多少个。

【题解】


枚举要找的子串的长度i; 根据height数组,找出连续>=i的height; 这几个起始的位置的后缀的最长公共前缀都大于等于i; 且它们起始位置开始的长度为i的串ts都是一样的。 (且没有其他和它们一样起始i个位置的字符串为ts的了); 则,找到最小的起始位置和最大的起始位置。 只要这两个位置的差大于等于i; 就说明这个ts,至少出现了两次(不重叠); 否则就全是重叠的。

【错的次数】


0

【反思】


熟练了就很自然了吧。

【代码】

#include<bits/stdc++.h>
using namespace std;

const int N = 2e5;
const int MAX_CHAR = 255;//每个数字的最大值。
char s[N + 10];//如果是数字,就写成int s[N+10]就好,从0开始存
int Sa[N + 10], T1[N + 10], T2[N + 10], C[N + 10];
int Height[N + 10], Rank[N + 10];

void build_Sa(int n, int m) {
int i, *x = T1, *y = T2;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[i] = s[i]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[i]]] = i;
for (int k = 1; k <= n; k <<= 1)
{
int p = 0;
for (i = n - k; i<n; i++) y[p++] = i;
for (i = 0; i<n; i++) if (Sa[i] >= k) y[p++] = Sa[i] - k;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[y[i]]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[y[i]]]] = y[i];
swap(x, y);
p = 1; x[Sa[0]] = 0;
for (i = 1; i<n; i++)
x[Sa[i]] = y[Sa[i - 1]] == y[Sa[i]] && y[Sa[i - 1] + k] == y[Sa[i] + k] ? p - 1 : p++;
if (p >= n) break;
m = p;
}
}

void getHeight(int n)
{
int i, j, k = 0;
for (i = 1; i <= n; i++) Rank[Sa[i]] = i;
for (i = 0; i<n; i++) {
if (k) k--;
j = Sa[Rank[i] - 1];
while (s[i + k] == s[j + k]) k++;
Height[Rank[i]] = k;
}
}

int main() {
int n;
while (~scanf("%s", s)) {
if (s[0] == '#') break;
n = strlen(s);
s
 = 0;
build_Sa(n + 1, MAX_CHAR);//注意调用n+1
getHeight(n);
int ans = 0;
for (int i = 1; i <= n / 2; i++) {
int mi = n+2, ma = 0;
for (int j = 2;j <= n;j++)
if (Height[j] >= i) {
mi = min(min(Sa[j], Sa[j - 1]), mi);
ma = max(max(Sa[j], Sa[j - 1]), ma);
}
else {
if (ma - mi >= i) ans++;
mi = n+2, ma = 0;
}
if (ma - mi >= i) ans++;
}
printf("%d\n", ans);
}
return 0;
}


内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: