結果
問題 | No.5010 Better Mo's Algorithm is Needed!! (Unweighted) |
ユーザー | t33f |
提出日時 | 2022-12-17 18:16:54 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 4,974 ms / 5,000 ms |
コード長 | 6,376 bytes |
コンパイル時間 | 2,319 ms |
実行使用メモリ | 7,684 KB |
スコア | 19,044,217,421 |
最終ジャッジ日時 | 2022-12-17 18:27:50 |
合計ジャッジ時間 | 654,940 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge14 |
純コード判定しない問題か言語 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 4,972 ms
7,420 KB |
testcase_01 | AC | 4,974 ms
7,548 KB |
testcase_02 | AC | 4,973 ms
7,500 KB |
testcase_03 | AC | 4,973 ms
7,368 KB |
testcase_04 | AC | 4,972 ms
7,516 KB |
testcase_05 | AC | 4,972 ms
7,584 KB |
testcase_06 | AC | 4,972 ms
7,424 KB |
testcase_07 | AC | 4,973 ms
7,512 KB |
testcase_08 | AC | 4,972 ms
7,324 KB |
testcase_09 | AC | 4,973 ms
7,424 KB |
testcase_10 | AC | 4,972 ms
7,600 KB |
testcase_11 | AC | 4,973 ms
7,336 KB |
testcase_12 | AC | 4,972 ms
7,320 KB |
testcase_13 | AC | 4,973 ms
7,516 KB |
testcase_14 | AC | 4,972 ms
7,388 KB |
testcase_15 | AC | 4,972 ms
7,432 KB |
testcase_16 | AC | 4,973 ms
7,576 KB |
testcase_17 | AC | 4,973 ms
7,392 KB |
testcase_18 | AC | 4,973 ms
7,324 KB |
testcase_19 | AC | 4,972 ms
7,416 KB |
testcase_20 | AC | 4,974 ms
7,324 KB |
testcase_21 | AC | 4,972 ms
7,628 KB |
testcase_22 | AC | 4,973 ms
7,408 KB |
testcase_23 | AC | 4,972 ms
7,572 KB |
testcase_24 | AC | 4,973 ms
7,324 KB |
testcase_25 | AC | 4,972 ms
7,508 KB |
testcase_26 | AC | 4,973 ms
7,372 KB |
testcase_27 | AC | 4,972 ms
7,324 KB |
testcase_28 | AC | 4,973 ms
7,504 KB |
testcase_29 | AC | 4,972 ms
7,428 KB |
testcase_30 | AC | 4,973 ms
7,336 KB |
testcase_31 | AC | 4,972 ms
7,464 KB |
testcase_32 | AC | 4,973 ms
7,540 KB |
testcase_33 | AC | 4,972 ms
7,384 KB |
testcase_34 | AC | 4,973 ms
7,464 KB |
testcase_35 | AC | 4,972 ms
7,492 KB |
testcase_36 | AC | 4,972 ms
7,352 KB |
testcase_37 | AC | 4,973 ms
7,500 KB |
testcase_38 | AC | 4,972 ms
7,464 KB |
testcase_39 | AC | 4,972 ms
7,344 KB |
testcase_40 | AC | 4,972 ms
7,596 KB |
testcase_41 | AC | 4,972 ms
7,324 KB |
testcase_42 | AC | 4,973 ms
7,336 KB |
testcase_43 | AC | 4,972 ms
7,556 KB |
testcase_44 | AC | 4,973 ms
7,444 KB |
testcase_45 | AC | 4,972 ms
7,424 KB |
testcase_46 | AC | 4,974 ms
7,432 KB |
testcase_47 | AC | 4,973 ms
7,524 KB |
testcase_48 | AC | 4,972 ms
7,384 KB |
testcase_49 | AC | 4,972 ms
7,460 KB |
testcase_50 | AC | 4,973 ms
7,548 KB |
testcase_51 | AC | 4,973 ms
7,328 KB |
testcase_52 | AC | 4,973 ms
7,468 KB |
testcase_53 | AC | 4,973 ms
7,612 KB |
testcase_54 | AC | 4,973 ms
7,612 KB |
testcase_55 | AC | 4,973 ms
7,472 KB |
testcase_56 | AC | 4,972 ms
7,320 KB |
testcase_57 | AC | 4,973 ms
7,324 KB |
testcase_58 | AC | 4,972 ms
7,468 KB |
testcase_59 | AC | 4,973 ms
7,384 KB |
testcase_60 | AC | 4,972 ms
7,392 KB |
testcase_61 | AC | 4,972 ms
7,612 KB |
testcase_62 | AC | 4,973 ms
7,328 KB |
testcase_63 | AC | 4,972 ms
7,416 KB |
testcase_64 | AC | 4,973 ms
7,596 KB |
testcase_65 | AC | 4,973 ms
7,460 KB |
testcase_66 | AC | 4,972 ms
7,376 KB |
testcase_67 | AC | 4,973 ms
7,532 KB |
testcase_68 | AC | 4,972 ms
7,536 KB |
testcase_69 | AC | 4,974 ms
7,348 KB |
testcase_70 | AC | 4,972 ms
7,496 KB |
testcase_71 | AC | 4,973 ms
7,512 KB |
testcase_72 | AC | 4,973 ms
7,476 KB |
testcase_73 | AC | 4,973 ms
7,596 KB |
testcase_74 | AC | 4,973 ms
7,480 KB |
testcase_75 | AC | 4,973 ms
7,328 KB |
testcase_76 | AC | 4,971 ms
7,508 KB |
testcase_77 | AC | 4,973 ms
7,524 KB |
testcase_78 | AC | 4,973 ms
7,380 KB |
testcase_79 | AC | 4,973 ms
7,588 KB |
testcase_80 | AC | 4,972 ms
7,336 KB |
testcase_81 | AC | 4,973 ms
7,332 KB |
testcase_82 | AC | 4,972 ms
7,600 KB |
testcase_83 | AC | 4,973 ms
7,628 KB |
testcase_84 | AC | 4,972 ms
7,404 KB |
testcase_85 | AC | 4,973 ms
7,512 KB |
testcase_86 | AC | 4,971 ms
7,684 KB |
testcase_87 | AC | 4,973 ms
7,340 KB |
testcase_88 | AC | 4,972 ms
7,592 KB |
testcase_89 | AC | 4,973 ms
7,540 KB |
testcase_90 | AC | 4,973 ms
7,324 KB |
testcase_91 | AC | 4,973 ms
7,588 KB |
testcase_92 | AC | 4,972 ms
7,552 KB |
testcase_93 | AC | 4,973 ms
7,328 KB |
testcase_94 | AC | 4,973 ms
7,524 KB |
testcase_95 | AC | 4,972 ms
7,504 KB |
testcase_96 | AC | 4,973 ms
7,324 KB |
testcase_97 | AC | 4,973 ms
7,472 KB |
testcase_98 | AC | 4,973 ms
7,436 KB |
testcase_99 | AC | 4,972 ms
7,380 KB |
testcase_100 | AC | 4,973 ms
7,616 KB |
testcase_101 | AC | 4,972 ms
7,324 KB |
testcase_102 | AC | 4,973 ms
7,460 KB |
testcase_103 | AC | 4,972 ms
7,452 KB |
testcase_104 | AC | 4,973 ms
7,540 KB |
testcase_105 | AC | 4,972 ms
7,324 KB |
testcase_106 | AC | 4,973 ms
7,592 KB |
testcase_107 | AC | 4,972 ms
7,348 KB |
testcase_108 | AC | 4,973 ms
7,472 KB |
testcase_109 | AC | 4,972 ms
7,476 KB |
testcase_110 | AC | 4,973 ms
7,492 KB |
testcase_111 | AC | 4,972 ms
7,340 KB |
testcase_112 | AC | 4,973 ms
7,440 KB |
testcase_113 | AC | 4,972 ms
7,452 KB |
testcase_114 | AC | 4,973 ms
7,384 KB |
testcase_115 | AC | 4,973 ms
7,600 KB |
testcase_116 | AC | 4,973 ms
7,472 KB |
testcase_117 | AC | 4,972 ms
7,460 KB |
testcase_118 | AC | 4,973 ms
7,536 KB |
testcase_119 | AC | 4,972 ms
7,540 KB |
ソースコード
#include <cstdlib> #include <cmath> #include <algorithm> #include <chrono> #include <cassert> #include <array> #include <vector> #include <numeric> #include <cstdio> #include <iostream> using namespace std; unsigned int xor128() { static unsigned int x=123456789, y=362436069, z=521288629, w=88675123; unsigned int t; t=(x^(x<<11)); x=y; y=z; z=w; return (w=(w^(w>>19))^(t^(t>>8))); } inline bool rand_bool(double prob) { constexpr double x = 1LL<<32; // uint_max+1 return xor128() < prob * x; } int timelimit = 5 * 1000, elapsed = 0, margin = 50; class Timer { chrono::system_clock::time_point start_time = chrono::system_clock::now(); public: Timer() {} long long get_elapsed_time() { auto diff = chrono::system_clock::now() - start_time; return chrono::duration_cast<chrono::milliseconds>(diff).count(); } } timer; constexpr int N = 200000, Q = 200000; int WT, ST; int L[Q], R[Q]; long long initial_cost(int i) { return 100 * (R[i] - L[i] + 1); } long long adj_cost(int i, int j) { return 100 * (abs(L[i] - L[j]) + abs(R[i] - R[j])); } long long calc_cost(const array<int, Q> &ans) { long long c = initial_cost(ans[0]); for (int i = 1; i < Q; ++i) c += adj_cost(ans[i - 1], ans[i]); return c; } long long score_of_cost(long long cost) { constexpr long long n = (long long)1e18; return n % cost * 2 >= cost ? n / cost + 1 : n / cost; } void read_input() { scanf("%*d %*d %d %d", &WT, &ST); assert(WT == 1); for (int i = 0; i < N; i++) scanf(" %*d"); for (int i = 0; i < Q; i++) scanf(" %d %d", L + i, R + i); } array<int, Q> solve(int B, int S) { array<int, N> a; iota(a.begin(), a.end(), 0); vector<vector<int>> buckets((N + B - 1) / B + 1); for (int i = 0; i < Q; i++) { const int b = L[i] / B; buckets[b].push_back(i); } array<int, Q> ans; int pos = 0; for (int b = 0; b < buckets.size(); ++b) { vector<int> &elems = buckets[b]; sort(elems.begin(), elems.end(), [b](int i, int j) { if (b % 2 == 0) return make_pair(R[i], L[i]) < make_pair(R[j], L[j]); else return make_pair(R[i], L[i]) > make_pair(R[j], L[j]); }); for (int i : elems) ans[pos++] = i; } long long cur_cost = calc_cost(ans); for (;;) { cerr << cur_cost << ' ' << score_of_cost(cur_cost) << endl; for (int i = 0; i < N; ++i) { if (timer.get_elapsed_time() + margin > timelimit) { if (initial_cost(ans.front()) > initial_cost(ans.back())) reverse(ans.begin(), ans.end()); return ans; } const int jmax = min(N, i + S); for (int j = i + 1; j < jmax; ++j) { long long diff = 0; if (i > 0) { diff -= adj_cost(ans[i - 1], ans[i]); diff += adj_cost(ans[i - 1], ans[j]); } else { diff -= initial_cost(ans[0]); diff += initial_cost(ans[j]); } if (j + 1 < N) { diff -= adj_cost(ans[j], ans[j + 1]); diff += adj_cost(ans[i], ans[j + 1]); } if (diff <= 0) { reverse(ans.begin() + i, ans.begin() + j + 1); cur_cost += diff; } } } } } array<int, Q> solve_sa(int B, int S) { array<int, N> a; iota(a.begin(), a.end(), 0); vector<vector<int>> buckets((N + B - 1) / B + 1); for (int i = 0; i < Q; i++) { const int b = L[i] / B; buckets[b].push_back(i); } array<int, Q> ans; int pos = 0; for (int b = 0; b < buckets.size(); ++b) { vector<int> &elems = buckets[b]; sort(elems.begin(), elems.end(), [b](int i, int j) { if (b % 2 == 0) return make_pair(R[i], L[i]) < make_pair(R[j], L[j]); else return make_pair(R[i], L[i]) > make_pair(R[j], L[j]); }); for (int i : elems) ans[pos++] = i; } long long cur_cost = calc_cost(ans); long long best_cost = cur_cost; array<int, N> best_ans = ans; const double start_temp = cur_cost / N * 0.1, end_temp = 0; for (;;) { cerr << cur_cost << ' ' << score_of_cost(cur_cost) << endl; for (int i = 0; i < N; ++i) { if (timer.get_elapsed_time() + margin > timelimit) { if (cur_cost < best_cost) { best_cost = cur_cost; best_ans = ans; } if (initial_cost(best_ans.front()) > initial_cost(best_ans.back())) reverse(best_ans.begin(), best_ans.end()); return best_ans; } double progress = double(timer.get_elapsed_time()) / timelimit, temp = start_temp + (end_temp - start_temp) * progress; const int jmax = min(N, i + S); for (int j = i + 1; j < jmax; ++j) { long long diff = 0; if (i > 0) { diff -= adj_cost(ans[i - 1], ans[i]); diff += adj_cost(ans[i - 1], ans[j]); } else { diff -= initial_cost(ans[0]); diff += initial_cost(ans[j]); } if (j + 1 < N) { diff -= adj_cost(ans[j], ans[j + 1]); diff += adj_cost(ans[i], ans[j + 1]); } bool accept = (diff <= 0); if (!accept && diff < temp * 3) { double prob = exp(-diff / temp); accept = rand_bool(prob); } if (accept) { reverse(ans.begin() + i, ans.begin() + j + 1); cur_cost += diff; } } } if (cur_cost < best_cost) { best_cost = cur_cost; best_ans = ans; // これがボトルネックになる } } } int main(int argc, char **argv) { read_input(); const int B = argc <= 1 ? 3500 : atoi(argv[1]); const int S = argc <= 2 ? 400 : atoi(argv[2]); const array<int, Q> ans = solve(B, S); for (int i : ans) cout << i + 1 << ' '; cout << endl; cerr << score_of_cost(calc_cost(ans)) << endl; }