結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
|
提出日時 | 2021-02-15 21:14:31 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 45 ms / 2,000 ms |
コード長 | 1,110 bytes |
コンパイル時間 | 798 ms |
コンパイル使用メモリ | 93,368 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-23 12:56:26 |
合計ジャッジ時間 | 2,645 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include <iostream>#include <string>#include <algorithm>#include <vector>#include <iomanip>#include <cmath>#include <stdio.h>#include <queue>#include <deque>#include <cstdio>#include <set>#include <map>#include <bitset>#include <stack>#include <cctype>using namespace std;long long fac[100010], finv[100010], inv[100010];long long mod = 1000000007;void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < 100010; i++) {fac[i] = fac[i - 1] * i % mod;inv[i] = mod - inv[mod % i] * (mod / i) % mod;finv[i] = finv[i - 1] * inv[i] % mod;}}long long COM(long long n, long long k) {if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % mod) % mod;}int main() {COMinit();long long n;cin >> n;long long ans = 0;for (long long i = 0; i < n; i++) {long long x = n - i;long long y = i;long long ans1 = 1;while (y > 0) {if ((y & 1) == 1) {ans1 = ans1 * x % mod;}x = x * x % mod;y >>= 1;}ans += COM(n, i) * ans1 % mod;ans %= mod;}cout << ans << endl;}