結果
| 問題 | 
                            No.2271 平方根の13桁精度近似計算
                             | 
                    
| コンテスト | |
| ユーザー | 
                             risujiroh
                         | 
                    
| 提出日時 | 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";
}
            
            
            
        
            
risujiroh