結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
![]() |
提出日時 | 2020-02-04 03:15:22 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 28 ms / 2,000 ms |
コード長 | 1,811 bytes |
コンパイル時間 | 1,835 ms |
コンパイル使用メモリ | 158,340 KB |
実行使用メモリ | 15,468 KB |
最終ジャッジ日時 | 2024-09-19 12:11:18 |
合計ジャッジ時間 | 4,449 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include "bits/stdc++.h"using namespace std;#define int long long#define FOR(i, a, b) for(int i=(a);i<(b);i++)#define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--)#define REP(i, n) for(int i=0; i<(n); i++)#define RREP(i, n) for(int i=(n-1); i>=0; i--)#define ALL(a) (a).begin(),(a).end()#define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end());#define CONTAIN(a, b) find(ALL(a), (b)) != (a).end()#define array2(type, x, y) array<array<type, y>, x>#define vector2(type) vector<vector<type> >#define out(...) printf(__VA_ARGS__)int dxy[] = {0, 1, 0, -1, 0};void solve();signed main(){#if DEBUGstd::ifstream in("input.txt");std::cin.rdbuf(in.rdbuf());#endifcin.tie(0);ios::sync_with_stdio(false);solve();return 0;}/*================================*/#if DEBUG#define SIZE 100#else#define SIZE 123450#endifint N,M,Q;const int MAX = 510000;const int MOD = 1000000007;long long fac[MAX], finv[MAX], inv[MAX];void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < MAX; 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(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;}int mod_pow(int b, int e, int m) {int result = 1;while (e > 0) {if ((e & 1) == 1) result = (result * b) % m;e >>= 1;b = (b * b) % m;}return result;}void solve() {cin>>N;COMinit();int ans = 0;REP(i,N) {int k = i+1;ans += COM(N,k) * mod_pow(k, N-k, MOD);ans %= MOD;}cout << ans << endl;}