結果

問題 No.575 n! / m / m / m...
ユーザー e869120e869120
提出日時 2018-06-10 12:21:33
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,115 bytes
コンパイル時間 675 ms
コンパイル使用メモリ 79,024 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-30 13:14:23
合計ジャッジ時間 3,148 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 WA * 13
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;

long double factorial(long long x) {
	if (x <= 10000) {
		long double ret = 0;
		for (int i = 1; i <= x; i++) ret += log10(1.0L*i);
		return ret;
	}
	long double ans = 1.0L*x*(log(1.0L*x) - 1) / log(10.0L);
	long double dig = log10(1.0L * x); dig /= 2;
	ans += dig + 0.399L;
	return ans;
}

long long n, m;

void calc(long double G) {
	long long T = (long long)(1.0L*G);
	G -= 1.0L*T;
	printf("%.10Lf", pow(10.0L, G));
	cout << "e" << T << endl;
}

int main() {
	cin >> n >> m;
	long double E = factorial(n);

	long long I = m;
	vector<pair<long long, long long>>F;
	for (long long i = 2; i*i <= m; i++) {
		long long cnt = 0;
		while (I%i == 0) { I /= i; cnt++; }
		if (cnt >= 1) F.push_back(make_pair(i, cnt));
	}
	if (I >= 2) F.push_back(make_pair(I, 1));

	long long maxn = 0;
	for (int i = 0; i < F.size(); i++) {
		long long J = n, sum = 0;
		while (J >= 1) {
			J /= F[i].first; sum += J;
		}
		maxn = max(maxn, sum / F[i].second);
	}

	long double T = log10(1.0L*m);
	E -= 1.0L*T*maxn;

	calc(E);
	return 0;
}
0