結果
問題 |
No.2271 平方根の13桁精度近似計算
|
ユーザー |
![]() |
提出日時 | 2023-04-15 00:26:56 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 221 ms / 2,000 ms |
コード長 | 871 bytes |
コンパイル時間 | 1,046 ms |
コンパイル使用メモリ | 29,824 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-10-10 15:20:03 |
合計ジャッジ時間 | 2,257 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 40 |
ソースコード
#include <stdio.h> int main () { long long n = 0LL; int e = 0; int res = 0; long long ans = -1LL; long long d = 1LL; long long tmp = 0LL; long long g[625][625] = {}; int gcnt[625] = {}; res = scanf("%lld", &n); res = scanf("%d", &e); while (e > 0) { d *= 5LL; e--; } if (n > 0LL) { n %= d; } else { n = (d-(-n)%d)%d; } for (int i = 0; i < 625; i++) { g[(i*i)%625][gcnt[(i*i)%625]] = (long long)i; gcnt[(i*i)%625]++; } while (ans < 0LL && tmp < d && tmp <= (1LL<<29LL)) { for (int i = 0; i < gcnt[(int)(n%625LL)]; i++) { long long tmp2 = tmp+g[(int)(n%625LL)][i]; if ((tmp2*tmp2)%d == n && tmp2 <= (1LL<<29LL)) { ans = tmp2; } } tmp += 625LL; } if (ans < 0LL) { printf("NaN\n"); } else { printf("%lld\n", ans); } return 0; }