結果
問題 | No.575 n! / m / m / m... |
ユーザー |
|
提出日時 | 2017-10-07 04:24:34 |
言語 | Java (openjdk 23) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,710 bytes |
コンパイル時間 | 2,280 ms |
コンパイル使用メモリ | 85,552 KB |
実行使用メモリ | 43,116 KB |
最終ジャッジ日時 | 2024-11-17 03:40:16 |
合計ジャッジ時間 | 9,473 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 9 WA * 14 |
ソースコード
import java.math.BigDecimal; import java.util.Arrays; import java.util.Scanner; public class Main { public static void main(String[] args) { new Main().run(); } void run() { Scanner sc = new Scanner(System.in); double n = sc.nextDouble(); double m = sc.nextDouble(); if (n < 100) { double ret = 1; for (int i = 2; i <= n; ++i) { int mul = i; while (mul % m == 0) mul /= m; ret *= mul; } int cnt = 0; while (ret >= 10) { ++cnt; ret /= 10; } System.out.println(ret + "e" + cnt); return; } long c = 0; for (int i = 1; i < 1000; ++i) c += n / Math.pow(m, i); ++n; long d = (long) (n * Math.log10((n + 1 / (12 * n - 1 / (10 * n))) / Math.E) + (Math.log10(2 * Math.PI) - Math.log10(n)) * 0.5 - c * Math.log10(m)); --n; double p = fac(n); for (int i = 0; i < c; ++i) { p /= m; p = regulate(p); } System.out.println(p + "e" + d); } double fac(double n) { if (n < 1e8) { double ret = 1; for (int i = 1; i <= n; ++i) { ret *= i; ret = regulate(ret); } return ret; } double ret = Math.sqrt(2 * Math.PI) * pow(n / Math.E, (long) n) / Math.sqrt(n) * (1 + 1 / (12 * n) + 1 / (288 * n * n) - 139 / (51840 * n * n * n)); ret = regulate(ret); return ret; } double regulate(double a) { while (a >= 10) a /= 10; while (a < 1) a *= 10; return a; } double pow(double a, long n) { double ret = 1; for (; n > 0; n >>= 1, a *= a) { while (a >= 10) a /= 10; if (n % 2 == 1) { ret *= a; while (ret >= 10) ret /= 10; } } return ret; } // n!=lg_10 n!=nlg_10 n void tr(Object... objects) { System.out.println(Arrays.deepToString(objects)); } }