結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
|
提出日時 | 2017-10-06 23:15:17 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 25 ms / 2,000 ms |
コード長 | 2,205 bytes |
コンパイル時間 | 1,160 ms |
コンパイル使用メモリ | 106,520 KB |
実行使用メモリ | 7,936 KB |
最終ジャッジ日時 | 2024-11-17 01:56:59 |
合計ジャッジ時間 | 2,768 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#define _USE_MATH_DEFINES#include <cstdio>#include <iostream>#include <sstream>#include <fstream>#include <iomanip>#include <algorithm>#include <cmath>#include <complex>#include <string>#include <vector>#include <list>#include <queue>#include <stack>#include <set>#include <map>#include <bitset>#include <numeric>#include <limits>#include <climits>#include <cfloat>#include <functional>#include <iterator>using namespace std;void mod_inverse(int n, vector<long long>& inv, int mod){inv.assign(n+1, -1);inv[1] = 1;for(int i=2; i<=n; ++i)inv[i] = inv[mod % i] * (mod - mod / i) % mod;}class FactorialCalculation{private:const int mod;vector<long long> factorial;vector<long long> invFactorial;public:FactorialCalculation(int n, int mod) : mod(mod){factorial.resize(n+1, 1);invFactorial.resize(n+1, 1);vector<long long> inv;mod_inverse(n, inv, mod);for(int i=1; i<=n; ++i){factorial[i] = factorial[i-1] * i % mod;invFactorial[i] = invFactorial[i-1] * inv[i] % mod;}}long long getFactorial(int n){return factorial[n];}long long getInvFactorial(int n){return invFactorial[n];}long long getPermutation(int n, int r){if(n < r)return 0;return factorial[n] * invFactorial[n-r] % mod;}long long getCombination(int n, int r){if(n < r)return 0;return getPermutation(n, r) * invFactorial[r] % mod;}long long getHomogeneous(int n, int r){return getCombination(n+r-1, r);}};const int MOD = 1000000007;long long power(int a, int b){long long ret = 1;long long tmp = a;while(b > 0){if(b & 1){ret *= tmp;ret %= MOD;}tmp *= tmp;tmp %= MOD;b >>= 1;}return ret;}int main(){int n;cin >> n;FactorialCalculation fc(2*n, MOD);long long ans = 0;for(int i=1; i<=n; ++i){ans += fc.getCombination(n, i) * power(i, n-i);ans %= MOD;}cout << ans << endl;return 0;}