結果
問題 |
No.1935 Water Simulation
|
ユーザー |
|
提出日時 | 2022-05-13 22:58:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 35 ms / 2,000 ms |
コード長 | 1,657 bytes |
コンパイル時間 | 2,245 ms |
コンパイル使用メモリ | 209,532 KB |
最終ジャッジ日時 | 2025-01-29 07:27:30 |
ジャッジサーバーID (参考情報) |
judge5 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 29 |
ソースコード
#include <bits/stdc++.h> using namespace std; #ifdef _RUTHEN #include "debug.hpp" #else #define show(...) true #endif int kth_vertex(vector<int> &a, long long K, int s) { // K <= 10^18 int LOG = 60; int N = (int)a.size(); vector<vector<int>> db(LOG, vector<int>(N)); db[0] = a; for (int k = 0; k < LOG - 1; k++) { for (int i = 0; i < N; i++) { db[k + 1][i] = db[k][db[k][i]]; } } for (int k = 0; k < LOG; k++) { if (K >> k & 1) { s = db[k][s]; } } return s; } using ll = long long; #define rep(i, n) for (int i = 0; i < (n); i++) template <class T> using V = vector<T>; int main() { ios::sync_with_stdio(false); cin.tie(0); V<int> A(4); rep(i, 4) cin >> A[i]; ll N; cin >> N; map<V<int>, int> id; V<int> B = {A[0], 0, 0, 0, 0}; int Nc = 0; id[B] = Nc++; while (1) { int pre = B[4] % 4, nex = (B[4] + 1) % 4; int t = min(B[pre], A[nex] - B[nex]); B[pre] -= t; B[nex] += t; B[4] = (B[4] + 1) % 4; if (id.count(B)) break; id[B] = Nc; Nc++; } V<V<int>> db(60, V<int>(Nc, -1)); V<int> C(Nc); for (auto &[Bi, c] : id) { auto B = Bi; int pre = B[4] % 4, nex = (B[4] + 1) % 4; int t = min(B[pre], A[nex] - B[nex]); B[pre] -= t; B[nex] += t; B[4] = (B[4] + 1) % 4; C[c] = id[B]; } auto ans = kth_vertex(C, N, 0); for (auto &[B, c] : id) { if (c == ans) { cout << B[0] << ' ' << B[1] << ' ' << B[2] << ' ' << B[3] << '\n'; } } return 0; }