結果

問題 No.575 n! / m / m / m...
ユーザー koba-e964koba-e964
提出日時 2017-10-06 23:22:46
言語 C++11
(gcc 8.4.1)
結果
AC  
実行時間 7 ms / 2,000 ms
コード長 1,929 Byte
コンパイル時間 605 ms
使用メモリ 4,184 KB
最終ジャッジ日時 2021-08-25 19:28:31
合計ジャッジ時間 2,094 ms
ジャッジサーバーID
(参考情報)
judge10 / judge12
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 2 ms
4,132 KB
testcase_01 AC 2 ms
3,968 KB
testcase_02 AC 1 ms
4,032 KB
testcase_03 AC 1 ms
4,112 KB
testcase_04 AC 1 ms
4,036 KB
testcase_05 AC 1 ms
3,904 KB
testcase_06 AC 1 ms
3,856 KB
testcase_07 AC 1 ms
4,068 KB
testcase_08 AC 1 ms
4,016 KB
testcase_09 AC 1 ms
3,828 KB
testcase_10 AC 1 ms
4,100 KB
testcase_11 AC 1 ms
4,068 KB
testcase_12 AC 6 ms
4,012 KB
testcase_13 AC 2 ms
4,020 KB
testcase_14 AC 3 ms
4,032 KB
testcase_15 AC 2 ms
4,020 KB
testcase_16 AC 4 ms
4,092 KB
testcase_17 AC 2 ms
4,092 KB
testcase_18 AC 2 ms
3,904 KB
testcase_19 AC 1 ms
4,016 KB
testcase_20 AC 2 ms
4,184 KB
testcase_21 AC 1 ms
4,048 KB
testcase_22 AC 1 ms
4,016 KB
testcase_23 AC 6 ms
3,904 KB
testcase_24 AC 1 ms
4,020 KB
testcase_25 AC 7 ms
4,020 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <algorithm>
#include <cassert>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <ctime>
#include <deque>
#include <functional>
#include <iomanip>
#include <iostream>
#include <list>
#include <map>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <utility>
#include <vector>

#define REP(i,s,n) for(int i=(int)(s);i<(int)(n);i++)

using namespace std;
typedef long long int ll;
typedef vector<int> VI;
typedef vector<ll> VL;
typedef pair<int, int> PI;

/**
 * For every prime factor p, perform result[p] += e. (where n = \prod p^e)
 * Verified by: Codeforces #400 E
 *              (http://codeforces.com/contest/776/submission/24956471)
 */
void factorize(long long v, std::map<long long, int> &result) {
  long long p = 2;
  while (v > 1 && p * p <= v) {
    int cnt = 0;
    while (v % p == 0) {
      cnt++;
      v /= p;
    }
    if (cnt > 0) {
      if (result.count(p) == 0) {
	result[p] = 0;
      }
      result[p] += cnt;
    }
    p += p == 2 ? 1 : 2;
  }
  if (v > 1) {
    if (result.count(v) == 0) {
      result[v] = 0;
    }
    result[v] += 1;
  }
}

const ll THRESHOLD = 100000;

int main(void) {
  ios::sync_with_stdio(false);
  cin.tie(0);
  ll n, m;
  cin >> n >> m;
  map<ll, int> fact;
  factorize(m, fact);
  ll mi = 1e16;
  for (auto v: fact) {
    ll p = v.first;
    int e = v.second;
    ll einn = 0;
    ll k = n;
    while (k > 0) {
      k /= p;
      einn += k;
    }
    mi = min(mi, einn / e);
  }
  double ret = 0;
  if (n < THRESHOLD) {
    REP(i, 1, n + 1) {
      ret += log10(i);
    }
  } else {
    // Stirling's approximation
    ret = n * (log10(n / exp(1)) + log10(n * sinh(1.0 / n) + 1.0 / (810.0 * pow(n, 6))) / 2.0);
    ret += log10(2 * acos(-1) * n) / 2.0;
  }
  ret -= mi * log10(m);
  ll tenexp = floor(ret);
  ret -= tenexp;
  ret = pow(10, ret);
  cout << ret << "e" << tenexp << "\n";
}
0