結果
| 問題 |
No.573 a^2[i] = a[i]
|
| コンテスト | |
| ユーザー |
bal4u
|
| 提出日時 | 2019-07-06 12:25:14 |
| 言語 | C (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 15 ms / 2,000 ms |
| コード長 | 791 bytes |
| コンパイル時間 | 274 ms |
| コンパイル使用メモリ | 30,720 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-25 03:59:28 |
| 合計ジャッジ時間 | 1,378 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 47 |
ソースコード
// yukicoder: No.573 a^2[i] = a[i]
// 2019.7.6 bal4u
#include <stdio.h>
typedef long long ll;
#define M 1000000007
#define MAX 100000
ll fact[MAX+2], inv[MAX+2], factinv[MAX+2];
ll comb(int n, int k) {
return fact[n] * factinv[k] % M * factinv[n-k] % M;
}
ll bigPow(int a, int p) {
ll r = 1;
while (p) {
if (p & 1) r = r * a % M;
a = (ll)a * a % M;
p >>= 1;
}
return r;
}
int main()
{
int i, N, ans;
scanf("%d", &N);
fact[0] = 1, inv[1] = 1, factinv[0] = 1;
for (i = 1; i <= N; i++) fact[i] = fact[i-1] * i % M;
for (i = 2; i <= N; i++) inv[i] = (M + (-(M/i) * inv[M % i]) % M) % M;
for (i = 1; i <= N; i++) factinv[i] = factinv[i-1] * inv[i] % M;
ans = 0;
for (i = 0; i <= N; i++) ans = (ans + comb(N, i) * bigPow(N-i, i)) % M;
printf("%d\n", ans);
}
bal4u