結果
| 問題 |
No.575 n! / m / m / m...
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-10-07 05:19:36 |
| 言語 | Java (openjdk 23) |
| 結果 |
WA
(最新)
J_TLE
(最初)
|
| 実行時間 | - |
| コード長 | 1,356 bytes |
| コンパイル時間 | 2,330 ms |
| コンパイル使用メモリ | 83,400 KB |
| 実行使用メモリ | 55,272 KB |
| 最終ジャッジ日時 | 2024-11-25 05:39:32 |
| 合計ジャッジ時間 | 8,641 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 21 WA * 2 |
ソースコード
import java.util.ArrayList;
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();
long c = Integer.MAX_VALUE;
ArrayList<Long> lis = plis((long) m);
for (long p : lis) {
long tmp = 0;
for (int i = 1; i < 1000; ++i)
tmp += n / Math.pow(p, i);
c = Math.min(c, tmp);
}
if (n < 170) {
double ret = 1;
for (int i = 2; i <= n; ++i) {
ret *= i;
}
while (ret % m == 0)
ret /= m;
int cnt = 0;
while (ret >= 10) {
++cnt;
ret /= 10;
}
System.out.println(ret + "e" + cnt);
return;
}
++n;
double v = (Math.log10(Math.sqrt(2 * Math.PI)) + n * Math.log10(n / Math.E) - 0.5 * Math.log10(n)
+ Math.log10(1 + 1 / (12 * n) + 1 / (288 * n * n) - 139 / (51840 * n * n * n))) - c * Math.log10(m);
System.out.println(Math.pow(10, v - (long) v) + "e" + (long)v);
}
ArrayList<Long> plis(long n) {
long ori = n;
ArrayList<Long> ret = new ArrayList<>();
for (long i = 2; i * i <= n; ++i) {
if (n % i == 0) {
while (n % i == 0) {
n /= i;
}
ret.add(i);
}
}
ret.add(n);
return ret;
}
void tr(Object... objects) {
System.out.println(Arrays.deepToString(objects));
}
}