# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
657400 | bicsi | Necklace (Subtask 4) (BOI19_necklace4) | C++14 | 353 ms | 632 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
using ld = long double;
const ld EPS = 6e-11;
int main() {
string s, t; cin >> s >> t;
mt19937_64 rng(22);
vector<complex<ld>> rnd(256);
for (int i = 0; i < 256; ++i) {
ld real = 1.0L * rng() / rng.max();
ld imag = 1.0L * rng() / rng.max();
rnd[i] = {real, imag};
}
vector<pair<ld, int>> vs, vt;
auto check = [&](int len) {
complex<ld> eps = polar(1.0L, 2.0L * M_PI / len);
pair<int, int> ans = {-1, -1};
for (int rev = 0; rev < 2; ++rev) {
vs.clear(); vt.clear();
for (int it = 0; it < 2; ++it) {
complex<ld> h = 0;
for (int i = 0; i < (int)s.size(); ++i) {
h = (h + rnd[s[i]]) * eps;
if (i >= len - 1) {
vs.emplace_back(norm(h), i - len + 1);
h -= rnd[s[i - len + 1]];
}
}
swap(s, t); swap(vs, vt);
}
sort(vs.begin(), vs.end());
sort(vt.begin(), vt.end());
for (int i = 0, j = 0; i < (int)vs.size(); ++i) {
while (j < vt.size() && vs[i].first + EPS > vt[j].first) {
if (abs(vt[j].first - vs[i].first) < EPS) {
ans.first = vs[i].second;
ans.second = vt[j].second;
if (rev) ans.first = s.size() - ans.first - len;
}
++j;
}
}
reverse(s.begin(), s.end());
}
return ans;
};
int ans = 0, ps = -1, pt = -1;
int bound = 1;
while (check(bound).first >= 0) bound *= 2;
for (int i = bound; i >= bound/2; --i) {
auto [a, b] = check(i);
if (a >= 0) {
ans = i, ps = a, pt = b;
break;
}
}
cout << ans << endl << ps << " " << pt << endl;
return 0;
}
컴파일 시 표준 에러 (stderr) 메시지
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |