結果
| 問題 |
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));
}
}