-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathd_kmp.cpp
54 lines (49 loc) · 1.38 KB
/
d_kmp.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
#include <iostream>
#include <fstream>
#include <list>
#include <vector>
int main() {
std::ifstream fin("temp.in");
std::string str, main;
fin >> str >> main;
int zf[str.size()];
std::fill(zf, zf + str.size(), 0);
int left = 0;
int right = 0;
for (int i = 1; i < str.size(); ++i) {
zf[i] = std::max(0, std::min(right - i, zf[i - left]));
while (i + zf[i] < str.size() && str[zf[i]] == str[i + zf[i]]) {
zf[i]++;
}
if (i + zf[i] > right) {
left = i;
right = i + zf[i];
}
}
const int size = str.size();
left = 0;
right = 0;
std::vector<int> z_cache(size, 0);
std::list<int> coordinates;
for (int i = 0; i < main.size(); ++i) {
int cur = i % size;
z_cache[cur] = std::max(0, std::min(right - i, zf[i - left]));
while (i < main.size() && z_cache[cur] < size &&
i + z_cache[cur] < main.size() &&
main[i + z_cache[cur]] ==
str[z_cache[cur]]) {
z_cache[cur]++;
}
if (i + z_cache[cur] > right) {
left = i;
right = i + z_cache[cur];
}
if (z_cache[cur] == size) {
coordinates.push_back(i + 1);
}
}
std::cout << coordinates.size() << '\n';
for (auto i: coordinates) {
// std::cout << i << ' ';
}
}