結果
| 問題 |
No.2326 Factorial to the Power of Factorial to the...
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-05-28 15:52:12 |
| 言語 | D (dmd 2.109.1) |
| 結果 |
AC
|
| 実行時間 | 46 ms / 2,000 ms |
| コード長 | 819 bytes |
| コンパイル時間 | 1,497 ms |
| コンパイル使用メモリ | 174,516 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-06-22 18:03:29 |
| 合計ジャッジ時間 | 2,728 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 20 |
ソースコード
import std;
void main () {
auto buf = readln.split.to!(int[]);
int N = buf[0], P = buf[1];
solve(N, P);
}
void solve (int N, int P) {
const int MOD = 1_000_000_007;
long factorialN = 1;
foreach (i; 2..N+1) {
factorialN *= i;
factorialN %= MOD;
}
foreach (i; 2..N+1) {
factorialN = modpow(factorialN, i, MOD);
}
long X = 0;
{
long base = P;
while (N/base != 0) {
X += N/P;
base *= P;
}
}
writeln((X*factorialN) % MOD);
}
long modpow (long x, long a, const long MOD) {
if (a == 0) {
return 1;
}
if (a == 1) {
return x;
}
if (a % 2 == 1) {
return (x * modpow(x, a-1, MOD)) % MOD;
}
long t = modpow(x, a/2, MOD);
return (t * t) % MOD;
}