結果
問題 | No.551 夏休みの思い出(2) |
ユーザー |
![]() |
提出日時 | 2019-02-02 04:01:51 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,826 ms / 4,000 ms |
コード長 | 2,464 bytes |
コンパイル時間 | 1,634 ms |
コンパイル使用メモリ | 176,924 KB |
実行使用メモリ | 24,572 KB |
最終ジャッジ日時 | 2024-11-24 20:27:46 |
合計ジャッジ時間 | 15,466 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 11 ms
6,820 KB |
testcase_01 | AC | 11 ms
6,820 KB |
testcase_02 | AC | 12 ms
6,820 KB |
testcase_03 | AC | 12 ms
6,816 KB |
testcase_04 | AC | 13 ms
6,816 KB |
testcase_05 | AC | 16 ms
6,816 KB |
testcase_06 | AC | 18 ms
6,816 KB |
testcase_07 | AC | 12 ms
6,816 KB |
testcase_08 | AC | 12 ms
6,816 KB |
testcase_09 | AC | 12 ms
6,820 KB |
testcase_10 | AC | 11 ms
6,816 KB |
testcase_11 | AC | 12 ms
6,816 KB |
testcase_12 | AC | 12 ms
6,820 KB |
testcase_13 | AC | 11 ms
6,816 KB |
testcase_14 | AC | 12 ms
6,820 KB |
testcase_15 | AC | 11 ms
6,820 KB |
testcase_16 | AC | 12 ms
6,816 KB |
testcase_17 | AC | 52 ms
13,876 KB |
testcase_18 | AC | 109 ms
24,448 KB |
testcase_19 | AC | 104 ms
24,448 KB |
testcase_20 | AC | 110 ms
24,444 KB |
testcase_21 | AC | 126 ms
24,328 KB |
testcase_22 | AC | 26 ms
6,820 KB |
testcase_23 | AC | 37 ms
9,872 KB |
testcase_24 | AC | 99 ms
23,676 KB |
testcase_25 | AC | 50 ms
13,492 KB |
testcase_26 | AC | 95 ms
23,036 KB |
testcase_27 | AC | 272 ms
24,444 KB |
testcase_28 | AC | 271 ms
24,480 KB |
testcase_29 | AC | 165 ms
24,448 KB |
testcase_30 | AC | 250 ms
24,452 KB |
testcase_31 | AC | 169 ms
24,320 KB |
testcase_32 | AC | 252 ms
24,444 KB |
testcase_33 | AC | 286 ms
24,448 KB |
testcase_34 | AC | 248 ms
24,448 KB |
testcase_35 | AC | 164 ms
24,448 KB |
testcase_36 | AC | 208 ms
24,448 KB |
testcase_37 | AC | 1,826 ms
24,448 KB |
testcase_38 | AC | 1,306 ms
24,448 KB |
testcase_39 | AC | 419 ms
24,572 KB |
testcase_40 | AC | 975 ms
24,448 KB |
testcase_41 | AC | 324 ms
24,448 KB |
testcase_42 | AC | 1,443 ms
24,440 KB |
testcase_43 | AC | 405 ms
24,320 KB |
testcase_44 | AC | 584 ms
24,448 KB |
testcase_45 | AC | 422 ms
24,420 KB |
testcase_46 | AC | 1,444 ms
24,444 KB |
testcase_47 | AC | 12 ms
6,816 KB |
testcase_48 | AC | 11 ms
6,816 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; using uint = unsigned int; using lint = long long int; using ulint = unsigned long long int; template<class T = int> using V = vector<T>; template<class T = int> using VV = V< V<T> >; template<class T, class U> void assign(V<T>& v, int n, const U& a) { v.assign(n, a); } template<class T, class... Args> void assign(V<T>& v, int n, const Args&... args) { v.resize(n); for (auto&& e : v) assign(e, args...); } lint tmod(lint a, lint p) { assert(p > 0); return (a %= p) < 0 ? a + p : a; } int mod_pow(lint a, lint n, int p) { assert(n >= 0); a = tmod(a, p); lint res = 1; while (n > 0) { if (n & 1) { (res *= a) %= p; } (a *= a) %= p; n >>= 1; } return res; } lint mod_inv(lint a, lint p) { a = tmod(a, p); lint b = p, x = 1, u = 0; while (b) { lint q = a / b; swap(a -= q * b, b); swap(x -= q * u, u); } return a == 1 ? tmod(x, p) : -1; } int mod_log(lint g, lint x, int p) { g = tmod(g, p); x = tmod(x, p); int m = 5e5; static unordered_map<int, int> mp; lint a = 1; if (mp.empty()) { for (int j = 0; j < m; ++j) { mp[a] = j; (a *= g) %= p; } } g = mod_inv(mod_pow(g, m, p), p), a = 1; for (int i = 0; i < (p - 1 + m - 1) / m; ++i) { if (mp.count(a * x % p)) { return i * m + mp[a * x % p]; } (a *= g) %= p; } return -1; } template<class Int> Int ext_gcd(Int a, Int b, Int& x, Int& y) { Int u = y = 0, v = x = 1; while (b) { Int q = a / b; swap(a -= q * b, b); swap(x -= q * u, u); swap(y -= q * v, v); } return a; } int mod_root(lint a, lint n, int p, int g) { a = tmod(a, p); if (a == 0) { return n > 0 ? 0 : -1; } n = tmod(n, p - 1); int e = mod_log(g, a, p); assert(e != -1); lint x, y; int d = ext_gcd<lint>(n, p - 1, x, y); if (e % d) return -1; return mod_pow(g, e / d * x, p); } int main() { cin.tie(nullptr); ios_base::sync_with_stdio(false); int p, g, q; cin >> p >> g >> q; while (q--) { lint a, b, c; cin >> a >> b >> c; lint D = tmod(b * b - 4 * a * c, p); lint sqrt_D = mod_root(D, 2, p, g); if (sqrt_D == -1) { cout << -1 << '\n'; continue; } lint x = -b + sqrt_D; lint y = -b - sqrt_D; x = tmod(mod_inv(2 * a, p) * x, p); y = tmod(mod_inv(2 * a, p) * y, p); if (x > y) swap(x, y); cout << x; if (y != x) { cout << ' ' << y; } cout << '\n'; } }