結果
問題 | No.2326 Factorial to the Power of Factorial to the... |
ユーザー |
![]() |
提出日時 | 2023-05-28 14:20:20 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 1,208 bytes |
コンパイル時間 | 2,166 ms |
コンパイル使用メモリ | 169,160 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-27 00:35:52 |
合計ジャッジ時間 | 2,397 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 20 |
ソースコード
#include<bits/stdc++.h>#define chmin(x,y) (x) = min((x),(y))#define chmax(x,y) (x) = max((x),(y))using namespace std;using ll = long long;const ll mod = 1000000007;const vector<int> dx = {1,0,-1,0}, dy = {0,1,0,-1};using Graph = vector<vector<int>>;// https://qiita.com/drken/items/3b4fdf0a78e7a138cd9a// a^n mod を計算するlong long modpow(long long a, long long n, long long mod) {long long res = 1;while (n > 0) {if (n & 1) res = res * a % mod;a = a * a % mod;n >>= 1;}return res;}// a^{-1} mod を計算するlong long modinv(long long a, long long mod) {return modpow(a, mod - 2, mod);}int main(){// inputll N,P; cin >> N >> P;// solvell fac = 1, fer = 0;for (ll i = 1; i <= N; i++){// update fac/ferfer *= i;fer += (fac * i) / mod;fac *= i;// take remainder of fac/ferfac %= mod;fer = (fer / mod) + fer % mod;}ll pwr = modpow(fac,fac+fer,mod);ll leg = 0; // N!がPで割り切れる回数を調べるwhile(N > 0){leg += N / P;N /= P;}// output// cout << fac << " " << pwr << " " << leg << endl;cout << (leg * pwr) % mod << endl;}