結果

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

ソースコード

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 minx = (1LL << 60);
    for (int i = 0; i < F.size(); i++) {
        long long J = n, sum = 0;
        while (J >= 1) {
            J /= F[i].first; sum += J;
        }
        minx = min(minx, sum / F[i].second);
    }
    long double T = log10(1.0L*m);
    E -= 1.0L*T*minx;
    calc(E);
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0