結果
問題 | No.551 夏休みの思い出(2) |
ユーザー | pekempey |
提出日時 | 2017-07-29 01:30:54 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,214 ms / 4,000 ms |
コード長 | 1,914 bytes |
コンパイル時間 | 741 ms |
コンパイル使用メモリ | 78,512 KB |
実行使用メモリ | 134,736 KB |
最終ジャッジ日時 | 2024-04-19 02:12:01 |
合計ジャッジ時間 | 43,628 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 491 ms
134,480 KB |
testcase_01 | AC | 611 ms
134,476 KB |
testcase_02 | AC | 696 ms
134,604 KB |
testcase_03 | AC | 769 ms
134,604 KB |
testcase_04 | AC | 776 ms
134,396 KB |
testcase_05 | AC | 794 ms
134,476 KB |
testcase_06 | AC | 792 ms
134,548 KB |
testcase_07 | AC | 849 ms
134,604 KB |
testcase_08 | AC | 857 ms
134,476 KB |
testcase_09 | AC | 843 ms
134,732 KB |
testcase_10 | AC | 820 ms
134,604 KB |
testcase_11 | AC | 855 ms
134,632 KB |
testcase_12 | AC | 846 ms
134,736 KB |
testcase_13 | AC | 851 ms
134,604 KB |
testcase_14 | AC | 848 ms
134,604 KB |
testcase_15 | AC | 853 ms
134,604 KB |
testcase_16 | AC | 845 ms
134,728 KB |
testcase_17 | AC | 839 ms
134,512 KB |
testcase_18 | AC | 829 ms
134,548 KB |
testcase_19 | AC | 818 ms
134,544 KB |
testcase_20 | AC | 827 ms
134,600 KB |
testcase_21 | AC | 813 ms
134,476 KB |
testcase_22 | AC | 855 ms
134,604 KB |
testcase_23 | AC | 850 ms
134,604 KB |
testcase_24 | AC | 829 ms
134,608 KB |
testcase_25 | AC | 831 ms
134,476 KB |
testcase_26 | AC | 831 ms
134,604 KB |
testcase_27 | AC | 798 ms
134,732 KB |
testcase_28 | AC | 787 ms
134,500 KB |
testcase_29 | AC | 755 ms
134,732 KB |
testcase_30 | AC | 780 ms
134,492 KB |
testcase_31 | AC | 748 ms
134,608 KB |
testcase_32 | AC | 860 ms
134,696 KB |
testcase_33 | AC | 777 ms
134,480 KB |
testcase_34 | AC | 768 ms
134,476 KB |
testcase_35 | AC | 750 ms
134,732 KB |
testcase_36 | AC | 759 ms
134,536 KB |
testcase_37 | AC | 1,213 ms
134,604 KB |
testcase_38 | AC | 1,102 ms
134,604 KB |
testcase_39 | AC | 813 ms
134,604 KB |
testcase_40 | AC | 1,016 ms
134,476 KB |
testcase_41 | AC | 796 ms
134,604 KB |
testcase_42 | AC | 1,110 ms
134,504 KB |
testcase_43 | AC | 814 ms
134,732 KB |
testcase_44 | AC | 877 ms
134,472 KB |
testcase_45 | AC | 838 ms
134,604 KB |
testcase_46 | AC | 1,214 ms
134,604 KB |
testcase_47 | AC | 630 ms
134,728 KB |
testcase_48 | AC | 630 ms
134,476 KB |
ソースコード
#include <iostream> #include <algorithm> #include <vector> 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 = 5000000; constexpr int T = (1000000000 + S - 1) / S; int64_t rs = 1; std::vector<std::pair<int64_t, int>> baby; for (int i = 0; i < S; i++) { baby.emplace_back(rs, i); rs = rs * r % p; } rs = modinv(rs, p); std::sort(baby.begin(), baby.end()); auto modlog = [&](int64_t a) -> int64_t { if (a == 0) return 1000000000; int64_t giant = 1; for (int i = 0; i < T; i++) { auto it = std::lower_bound(baby.begin(), baby.end(), std::make_pair(a * giant % p, 0)); if (it != baby.end() && a * giant % p == it->first) { return it->second + 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); } } } }