結果

問題 No.502 階乗を計算するだけ
ユーザー Min_25Min_25
提出日時 2017-04-08 11:49:55
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 709 ms / 1,000 ms
コード長 1,338 bytes
コンパイル時間 520 ms
コンパイル使用メモリ 56,132 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-24 23:19:36
合計ジャッジ時間 5,911 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 1 ms
4,380 KB
testcase_14 AC 1 ms
4,380 KB
testcase_15 AC 1 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 1 ms
4,380 KB
testcase_18 AC 2 ms
4,380 KB
testcase_19 AC 1 ms
4,380 KB
testcase_20 AC 1 ms
4,380 KB
testcase_21 AC 1 ms
4,380 KB
testcase_22 AC 3 ms
4,376 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 3 ms
4,376 KB
testcase_25 AC 2 ms
4,376 KB
testcase_26 AC 3 ms
4,380 KB
testcase_27 AC 2 ms
4,376 KB
testcase_28 AC 2 ms
4,380 KB
testcase_29 AC 1 ms
4,380 KB
testcase_30 AC 3 ms
4,380 KB
testcase_31 AC 2 ms
4,376 KB
testcase_32 AC 551 ms
4,376 KB
testcase_33 AC 153 ms
4,376 KB
testcase_34 AC 709 ms
4,376 KB
testcase_35 AC 347 ms
4,376 KB
testcase_36 AC 357 ms
4,376 KB
testcase_37 AC 637 ms
4,380 KB
testcase_38 AC 56 ms
4,380 KB
testcase_39 AC 306 ms
4,376 KB
testcase_40 AC 154 ms
4,376 KB
testcase_41 AC 1 ms
4,376 KB
testcase_42 AC 1 ms
4,380 KB
testcase_43 AC 1 ms
4,376 KB
testcase_44 AC 2 ms
4,376 KB
testcase_45 AC 1 ms
4,376 KB
testcase_46 AC 1 ms
4,380 KB
testcase_47 AC 2 ms
4,376 KB
testcase_48 AC 1 ms
4,376 KB
testcase_49 AC 1 ms
4,380 KB
testcase_50 AC 1 ms
4,376 KB
testcase_51 AC 1 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <vector>
#include <algorithm>
#include <cassert>

using namespace std;

using i64 = long long;

int pow_mod(int b, int e, int mod) {
  int ret = 1;
  for (; e; e >>= 1, b = i64(b) * b % mod) {
    if (e & 1) ret = i64(ret) * b % mod;
  }
  return ret;
}

int fact_valuation(int n, int p) {
  int ret = n /= p;
  while (n >= p) {
    n /= p;
    ret += n;
  }
  return ret;
}

template <int p>
int fact_mod(i64 n) {
  if (n >= p) return 0;
  if (n * 2 > p) {
    int ret = fact_mod<p>(p - n - 1);
    ret = (p - n - 1) & 1 ? ret : p - ret;
    return pow_mod(ret, p - 2, p);
  }
  vector<int> ns;
  for (int i = n; i > 1; i /= 5) 
    for (int j = i; j > 1; j /= 3) 
      for (int k = j; k > 1; k >>= 1) 
        ns.push_back(k);

  std::sort(ns.begin(), ns.end());
  const int d[8] = {6, 4, 2, 4, 2, 4, 6, 2};
  int ret = 1, prod = 1, i = 1, r = 0;
  for (auto end : ns) {
    for (; i <= end; i += d[r & 7], ++r) {
      prod = i64(prod) * i % p;
    }
    ret = i64(ret) * prod % p;
  }
  ret = i64(ret) * pow_mod(2, fact_valuation(n, 2), p) % p;
  ret = i64(ret) * pow_mod(3, fact_valuation(n, 3), p) % p;
  ret = i64(ret) * pow_mod(5, fact_valuation(n, 5), p) % p;
  return ret;
}

int main() {
  const int mod = 1e9 + 7;
  i64 N;
  while (~scanf("%lld", &N)) {
    printf("%d\n", fact_mod<mod>(N));
  }
}
0