結果
問題 |
No.573 a^2[i] = a[i]
|
ユーザー |
|
提出日時 | 2019-02-03 14:51:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 61 ms / 2,000 ms |
コード長 | 904 bytes |
コンパイル時間 | 764 ms |
コンパイル使用メモリ | 71,456 KB |
実行使用メモリ | 6,656 KB |
最終ジャッジ日時 | 2024-12-16 11:00:50 |
合計ジャッジ時間 | 2,549 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include<iostream> #include<vector> using namespace std; typedef long long ll; const ll mod = 1e9 + 7; ll modpow(ll a, ll b, ll p = 1e9+7){ if(b == 0) return 1; if(b % 2 == 0){ ll d = modpow(a, b/2, p); return (d*d) % p; }else{ return (a%p * modpow(a, b-1, p)) % p; } } int main(){ vector<ll> fact(100001), inv(100001); fact[0] = 1; for(ll i = 1; i <= 100000; i++){ fact[i] = (fact[i-1] * i)%mod; } inv[100000] = modpow(fact[100000], mod-2, mod); for(ll i = 99999; i >= 0; i--){ inv[i] = (inv[i+1] * (i+1))%mod; } int n; cin >> n; auto ncm = [=](ll n, ll m) -> ll{ return ((fact[n]*inv[m])%mod*inv[n-m])%mod; }; ll ans = 0; for(int i = 1; i <= n; i++){ ans += (ncm(n,i)*modpow(i, n-i, mod))%mod; ans %= mod; } cout << ans << endl; return 0; }