結果

問題 No.575 n! / m / m / m...
ユーザー 37zigen37zigen
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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));
	}

}
0