結果

問題 No.575 n! / m / m / m...
ユーザー 37zigen37zigen
提出日時 2017-10-07 03:57:33
言語 Java21
(openjdk 21)
結果
WA  
実行時間 -
コード長 1,661 bytes
コンパイル時間 2,921 ms
コンパイル使用メモリ 81,820 KB
実行使用メモリ 43,088 KB
最終ジャッジ日時 2024-11-17 03:29:32
合計ジャッジ時間 9,507 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 163 ms
42,836 KB
testcase_01 WA -
testcase_02 AC 161 ms
42,808 KB
testcase_03 AC 166 ms
42,768 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 AC 155 ms
42,652 KB
testcase_07 AC 167 ms
42,576 KB
testcase_08 AC 150 ms
42,388 KB
testcase_09 AC 161 ms
42,728 KB
testcase_10 AC 155 ms
42,760 KB
testcase_11 AC 153 ms
42,644 KB
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
権限があれば一括ダウンロードができます

ソースコード

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);
		long d = (long) (n * Math.log10((n + 1 / (12 * n - 1 / (10 * n))) / Math.E)
				+ Math.log10(Math.sqrt(2 * Math.PI / n)) - c * Math.log10(m));
		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;
	}

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

}
0