結果
問題 | No.5010 Better Mo's Algorithm is Needed!! (Unweighted) |
ユーザー | t33f |
提出日時 | 2022-12-17 16:40:24 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 4,974 ms / 5,000 ms |
コード長 | 3,291 bytes |
コンパイル時間 | 1,106 ms |
実行使用メモリ | 7,692 KB |
スコア | 19,046,722,799 |
最終ジャッジ日時 | 2022-12-17 16:51:13 |
合計ジャッジ時間 | 644,045 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge15 |
純コード判定しない問題か言語 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 4,972 ms
7,452 KB |
testcase_01 | AC | 4,973 ms
7,628 KB |
testcase_02 | AC | 4,973 ms
7,636 KB |
testcase_03 | AC | 4,972 ms
7,444 KB |
testcase_04 | AC | 4,973 ms
7,452 KB |
testcase_05 | AC | 4,972 ms
7,680 KB |
testcase_06 | AC | 4,973 ms
7,344 KB |
testcase_07 | AC | 4,972 ms
7,392 KB |
testcase_08 | AC | 4,973 ms
7,356 KB |
testcase_09 | AC | 4,972 ms
7,328 KB |
testcase_10 | AC | 4,973 ms
7,552 KB |
testcase_11 | AC | 4,972 ms
7,396 KB |
testcase_12 | AC | 4,972 ms
7,352 KB |
testcase_13 | AC | 4,972 ms
7,476 KB |
testcase_14 | AC | 4,973 ms
7,352 KB |
testcase_15 | AC | 4,972 ms
7,280 KB |
testcase_16 | AC | 4,972 ms
7,580 KB |
testcase_17 | AC | 4,972 ms
7,492 KB |
testcase_18 | AC | 4,972 ms
7,456 KB |
testcase_19 | AC | 4,972 ms
7,432 KB |
testcase_20 | AC | 4,972 ms
7,280 KB |
testcase_21 | AC | 4,972 ms
7,468 KB |
testcase_22 | AC | 4,972 ms
7,380 KB |
testcase_23 | AC | 4,972 ms
7,692 KB |
testcase_24 | AC | 4,973 ms
7,464 KB |
testcase_25 | AC | 4,973 ms
7,500 KB |
testcase_26 | AC | 4,972 ms
7,288 KB |
testcase_27 | AC | 4,973 ms
7,404 KB |
testcase_28 | AC | 4,972 ms
7,456 KB |
testcase_29 | AC | 4,973 ms
7,376 KB |
testcase_30 | AC | 4,972 ms
7,424 KB |
testcase_31 | AC | 4,972 ms
7,420 KB |
testcase_32 | AC | 4,972 ms
7,532 KB |
testcase_33 | AC | 4,973 ms
7,300 KB |
testcase_34 | AC | 4,973 ms
7,476 KB |
testcase_35 | AC | 4,972 ms
7,512 KB |
testcase_36 | AC | 4,971 ms
7,476 KB |
testcase_37 | AC | 4,972 ms
7,592 KB |
testcase_38 | AC | 4,972 ms
7,352 KB |
testcase_39 | AC | 4,972 ms
7,300 KB |
testcase_40 | AC | 4,972 ms
7,612 KB |
testcase_41 | AC | 4,972 ms
7,352 KB |
testcase_42 | AC | 4,973 ms
7,356 KB |
testcase_43 | AC | 4,972 ms
7,420 KB |
testcase_44 | AC | 4,972 ms
7,468 KB |
testcase_45 | AC | 4,971 ms
7,416 KB |
testcase_46 | AC | 4,972 ms
7,436 KB |
testcase_47 | AC | 4,972 ms
7,676 KB |
testcase_48 | AC | 4,973 ms
7,284 KB |
testcase_49 | AC | 4,971 ms
7,456 KB |
testcase_50 | AC | 4,973 ms
7,496 KB |
testcase_51 | AC | 4,972 ms
7,356 KB |
testcase_52 | AC | 4,972 ms
7,424 KB |
testcase_53 | AC | 4,972 ms
7,480 KB |
testcase_54 | AC | 4,971 ms
7,684 KB |
testcase_55 | AC | 4,972 ms
7,492 KB |
testcase_56 | AC | 4,972 ms
7,400 KB |
testcase_57 | AC | 4,972 ms
7,452 KB |
testcase_58 | AC | 4,972 ms
7,448 KB |
testcase_59 | AC | 4,972 ms
7,336 KB |
testcase_60 | AC | 4,972 ms
7,400 KB |
testcase_61 | AC | 4,972 ms
7,540 KB |
testcase_62 | AC | 4,971 ms
7,336 KB |
testcase_63 | AC | 4,973 ms
7,276 KB |
testcase_64 | AC | 4,972 ms
7,596 KB |
testcase_65 | AC | 4,972 ms
7,336 KB |
testcase_66 | AC | 4,972 ms
7,284 KB |
testcase_67 | AC | 4,971 ms
7,576 KB |
testcase_68 | AC | 4,972 ms
7,540 KB |
testcase_69 | AC | 4,973 ms
7,416 KB |
testcase_70 | AC | 4,972 ms
7,448 KB |
testcase_71 | AC | 4,973 ms
7,556 KB |
testcase_72 | AC | 4,972 ms
7,396 KB |
testcase_73 | AC | 4,971 ms
7,552 KB |
testcase_74 | AC | 4,972 ms
7,488 KB |
testcase_75 | AC | 4,972 ms
7,372 KB |
testcase_76 | AC | 4,972 ms
7,512 KB |
testcase_77 | AC | 4,972 ms
7,552 KB |
testcase_78 | AC | 4,972 ms
7,444 KB |
testcase_79 | AC | 4,972 ms
7,452 KB |
testcase_80 | AC | 4,972 ms
7,332 KB |
testcase_81 | AC | 4,972 ms
7,280 KB |
testcase_82 | AC | 4,974 ms
7,640 KB |
testcase_83 | AC | 4,972 ms
7,580 KB |
testcase_84 | AC | 4,972 ms
7,324 KB |
testcase_85 | AC | 4,973 ms
7,428 KB |
testcase_86 | AC | 4,973 ms
7,504 KB |
testcase_87 | AC | 4,972 ms
7,332 KB |
testcase_88 | AC | 4,972 ms
7,564 KB |
testcase_89 | AC | 4,972 ms
7,604 KB |
testcase_90 | AC | 4,973 ms
7,460 KB |
testcase_91 | AC | 4,972 ms
7,500 KB |
testcase_92 | AC | 4,973 ms
7,492 KB |
testcase_93 | AC | 4,972 ms
7,372 KB |
testcase_94 | AC | 4,973 ms
7,532 KB |
testcase_95 | AC | 4,972 ms
7,324 KB |
testcase_96 | AC | 4,972 ms
7,280 KB |
testcase_97 | AC | 4,972 ms
7,480 KB |
testcase_98 | AC | 4,972 ms
7,352 KB |
testcase_99 | AC | 4,972 ms
7,416 KB |
testcase_100 | AC | 4,971 ms
7,660 KB |
testcase_101 | AC | 4,972 ms
7,264 KB |
testcase_102 | AC | 4,973 ms
7,300 KB |
testcase_103 | AC | 4,972 ms
7,444 KB |
testcase_104 | AC | 4,972 ms
7,496 KB |
testcase_105 | AC | 4,972 ms
7,264 KB |
testcase_106 | AC | 4,972 ms
7,560 KB |
testcase_107 | AC | 4,972 ms
7,268 KB |
testcase_108 | AC | 4,972 ms
7,284 KB |
testcase_109 | AC | 4,972 ms
7,480 KB |
testcase_110 | AC | 4,973 ms
7,496 KB |
testcase_111 | AC | 4,973 ms
7,268 KB |
testcase_112 | AC | 4,972 ms
7,396 KB |
testcase_113 | AC | 4,973 ms
7,460 KB |
testcase_114 | AC | 4,972 ms
7,472 KB |
testcase_115 | AC | 4,972 ms
7,540 KB |
testcase_116 | AC | 4,972 ms
7,388 KB |
testcase_117 | AC | 4,972 ms
7,456 KB |
testcase_118 | AC | 4,972 ms
7,476 KB |
testcase_119 | AC | 4,972 ms
7,644 KB |
ソースコード
#include <cstdlib> #include <cmath> #include <algorithm> #include <chrono> #include <cassert> #include <array> #include <vector> #include <numeric> #include <cstdio> #include <iostream> using namespace std; 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, N> &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, N> 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) 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; } } } } return 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, N> ans = solve(B, S); for (int i : ans) cout << i + 1 << ' '; cout << endl; cerr << score_of_cost(calc_cost(ans)) << endl; }