結果
問題 |
No.2271 平方根の13桁精度近似計算
|
ユーザー |
![]() |
提出日時 | 2023-04-14 23:49:24 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,010 bytes |
コンパイル時間 | 3,058 ms |
コンパイル使用メモリ | 247,952 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-10-10 14:50:25 |
合計ジャッジ時間 | 4,660 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 38 WA * 2 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using atcoder::modint; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int N, E; cin >> N >> E; if (E == 0) { cout << "0\n"; return 0; } int mod = 1; modint ::set_mod(mod); modint ans = 0; while (E--) { modint::set_mod(mod * 5); int r = ans.val(); bool found = false; for (int i = 0; i < 5; ++i) { auto x = modint(r + i * mod); if (x * x == N) { ans = x; found = true; break; } } if (!found) { cout << "NaN\n"; return 0; } mod *= 5; } auto go = [&](modint x) { if (x * x != N) { return; } if (x.val() <= 1 << 29) { cout << x.val() << '\n'; exit(0); } if (-(1 << 29) <= x.val() - mod) { cout << x.val() - mod << '\n'; exit(0); } }; for (int i = 0; i < 5; ++i) { go(ans + mod / 5 * i); go(-ans + mod / 5 * i); } cout << "NaN\n"; }