結果
問題 | No.551 夏休みの思い出(2) |
ユーザー | pekempey |
提出日時 | 2017-07-29 01:12:28 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,873 ms / 4,000 ms |
コード長 | 1,759 bytes |
コンパイル時間 | 770 ms |
コンパイル使用メモリ | 85,712 KB |
実行使用メモリ | 46,200 KB |
最終ジャッジ日時 | 2024-04-18 18:02:51 |
合計ジャッジ時間 | 21,022 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 23 ms
5,248 KB |
testcase_01 | AC | 23 ms
5,376 KB |
testcase_02 | AC | 21 ms
5,376 KB |
testcase_03 | AC | 23 ms
5,376 KB |
testcase_04 | AC | 24 ms
5,376 KB |
testcase_05 | AC | 26 ms
5,376 KB |
testcase_06 | AC | 27 ms
5,376 KB |
testcase_07 | AC | 23 ms
5,376 KB |
testcase_08 | AC | 23 ms
5,376 KB |
testcase_09 | AC | 24 ms
5,376 KB |
testcase_10 | AC | 23 ms
5,376 KB |
testcase_11 | AC | 23 ms
5,376 KB |
testcase_12 | AC | 23 ms
5,376 KB |
testcase_13 | AC | 23 ms
5,376 KB |
testcase_14 | AC | 22 ms
5,376 KB |
testcase_15 | AC | 23 ms
5,376 KB |
testcase_16 | AC | 21 ms
5,376 KB |
testcase_17 | AC | 73 ms
14,104 KB |
testcase_18 | AC | 180 ms
30,948 KB |
testcase_19 | AC | 174 ms
29,412 KB |
testcase_20 | AC | 193 ms
31,076 KB |
testcase_21 | AC | 294 ms
44,788 KB |
testcase_22 | AC | 39 ms
6,792 KB |
testcase_23 | AC | 53 ms
10,104 KB |
testcase_24 | AC | 174 ms
24,036 KB |
testcase_25 | AC | 84 ms
13,720 KB |
testcase_26 | AC | 146 ms
23,144 KB |
testcase_27 | AC | 607 ms
46,068 KB |
testcase_28 | AC | 585 ms
45,936 KB |
testcase_29 | AC | 474 ms
46,068 KB |
testcase_30 | AC | 596 ms
46,068 KB |
testcase_31 | AC | 491 ms
46,072 KB |
testcase_32 | AC | 561 ms
45,940 KB |
testcase_33 | AC | 629 ms
46,072 KB |
testcase_34 | AC | 611 ms
45,944 KB |
testcase_35 | AC | 539 ms
46,200 KB |
testcase_36 | AC | 550 ms
46,196 KB |
testcase_37 | AC | 1,873 ms
46,068 KB |
testcase_38 | AC | 1,456 ms
46,068 KB |
testcase_39 | AC | 692 ms
45,896 KB |
testcase_40 | AC | 1,110 ms
45,940 KB |
testcase_41 | AC | 635 ms
46,072 KB |
testcase_42 | AC | 1,505 ms
46,068 KB |
testcase_43 | AC | 693 ms
45,936 KB |
testcase_44 | AC | 889 ms
46,068 KB |
testcase_45 | AC | 771 ms
46,068 KB |
testcase_46 | AC | 1,725 ms
45,944 KB |
testcase_47 | AC | 21 ms
5,376 KB |
testcase_48 | AC | 22 ms
5,376 KB |
ソースコード
#include <iostream> #include <algorithm> #include <vector> #include <unordered_map> int64_t modpow(int64_t a, int64_t b, int64_t m) { int64_t ret = 1; for (; b > 0; b >>= 1) { if (b & 1) ret = ret * a % m; a = a * a % m; } return ret; } int64_t modinv(int64_t a, int64_t p) { return modpow(a, p - 2, p); } int main() { int64_t p, r; std::cin >> p >> r; auto modulo = [&](int64_t a) { a %= p; if (a < 0) a += p; return a; }; constexpr int S = 1000000; constexpr int T = (1000000000 + S - 1) / S; int64_t rs = 1; std::unordered_map<int64_t, int> baby; for (int i = 0; i < S; i++) { if (!baby.count(rs)) baby[rs] = i; rs = rs * r % p; } rs = modinv(rs, p); auto modlog = [&](int64_t a) -> int64_t { if (a == 0) return 1000000000; int64_t giant = 1; for (int i = 0; i < T; i++) { if (baby.count(a * giant % p)) return baby[a * giant % p] + i * S; giant = giant * rs % p; } return -1; }; int q; std::cin >> q; while (q--) { int64_t a, b, c; scanf("%lld %lld %lld", &a, &b, &c); int64_t d = modulo(b * b - 4 * a * c); int64_t e = modlog(d); if (e == -1 || e % 2 == 1) { puts("-1"); } else { d = e == 1000000000 ? 0LL : modpow(r, e / 2, p); int64_t x1 = modulo((-b + d) * modinv(2 * a % p, p)); int64_t x2 = modulo((-b - d) * modinv(2 * a % p, p)); if (x1 > x2) std::swap(x1, x2); if (x1 == x2) { printf("%lld\n", x1); } else { printf("%lld %lld\n", x1, x2); } } } }