結果

問題 No.575 n! / m / m / m...
ユーザー 37zigen37zigen
提出日時 2017-10-07 05:54:02
言語 Java
(openjdk 23)
結果
AC  
実行時間 176 ms / 2,000 ms
コード長 1,540 bytes
コンパイル時間 2,518 ms
コンパイル使用メモリ 82,716 KB
実行使用メモリ 42,972 KB
最終ジャッジ日時 2024-11-17 03:45:51
合計ジャッジ時間 7,724 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

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 = Long.MAX_VALUE;
		ArrayList<long[]> lis = flis((long) (m + 1e-9));
		for (long[] f : lis) {
			long tmp = 0;
			long n_ = (long) (n + 1e-9);
			while (n_ > 1) {
				tmp += n_ / f[0];
				n_ /= f[0];
			}
			c = Math.min(c, tmp / f[1]);
		}
		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 =n * Math.log10((n + 1 / (12 * n - 1 / (10 * n))) / Math.E)
				+ Math.log10(Math.sqrt(2 * Math.PI / n)) - c * Math.log10(m);
		System.out.println(Math.pow(10, v - (long) v) + "e" + (long) v);
	}

	ArrayList<long[]> flis(long n) {
		ArrayList<long[]> ret = new ArrayList<>();
		for (long i = 2; i * i <= n; ++i) {
			if (n % i == 0) {
				int c = 0;
				while (n % i == 0) {
					n /= i;
					++c;
				}
				ret.add(new long[] { i, c });
			}
		}
		if (n > 1)
			ret.add(new long[] { n, 1 });
		return ret;
	}

	long pow(long a, long n) {
		long ret = 1;
		for (; n > 0; n >>= 1, a *= a) {
			if (n % 2 == 1) {
				ret *= a;
			}
		}
		return ret;
	}

	void tr(Object... objects) {
		System.out.println(Arrays.deepToString(objects));
	}

}
0