結果
| 問題 |
No.575 n! / m / m / m...
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-11-02 00:28:05 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 881 bytes |
| コンパイル時間 | 1,527 ms |
| コンパイル使用メモリ | 165,964 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-22 14:54:42 |
| 合計ジャッジ時間 | 4,407 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 10 WA * 13 |
ソースコード
#ifdef __APPLE__
#include "bits:stdc++.h"
#else
#include <bits/stdc++.h>
#endif
using namespace std;
typedef long long LL;
typedef pair<LL, LL> PLL;
auto scan = []{ static LL a; scanf("%lld", &a); return a;};
const double PI = acos(-1);
int main() {
LL n = scan(), m = scan();
if (n <= 10) {
LL ans = 1;
for (LL i = 1; i <= n; ++i) ans *= i;
while (ans % m == 0) ans /= m;
double lgans = log10(ans);
LL flo = (LL) lgans;
lgans -= flo;
printf("%.2fe%lld\n",pow(10,lgans),flo);
}
else {
double lgans = log10(sqrt(2 * PI * n)) + n * log10(n/exp(1)), check = m;
LL div = 0;
//printf("%.4E\n",lgans);
for (LL cur = m; cur <= n && check <= n + 1; cur = cur * m, check = check * m)
div += n / cur;
lgans -= log10(pow(m,div));
LL flo = (LL) lgans;
lgans -= flo;
printf("%.2fe%lld\n",pow(10,lgans),flo);
}
}