結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
![]() |
提出日時 | 2017-11-22 20:47:09 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 14 ms / 2,000 ms |
コード長 | 1,889 bytes |
コンパイル時間 | 783 ms |
コンパイル使用メモリ | 92,248 KB |
実行使用メモリ | 5,888 KB |
最終ジャッジ日時 | 2024-11-26 11:22:57 |
合計ジャッジ時間 | 2,385 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include <algorithm>#include <climits>#include <cmath>#include <cstdio>#include <cstdlib>#include <ctime>#include <iostream>#include <sstream>#include <functional>#include <map>#include <string>#include <cstring>#include <vector>#include <queue>#include <stack>#include <deque>#include <set>#include <list>#include <numeric>using namespace std;typedef long long ll;typedef unsigned long long ull;typedef pair<ll,ll> P;const double PI = 3.14159265358979323846;const double EPS = 1e-12;const ll INF = 1LL<<29;const ll mod = 1e9+7;#define rep(i,n) for(int (i)=0;(i)<(ll)(n);++(i))#define repd(i,n,d) for(ll (i)=0;(i)<(ll)(n);(i)+=(d))#define all(v) (v).begin(), (v).end()#define pb(x) push_back(x)#define mp(x,y) make_pair((x),(y))#define mset(m,v) memset((m),(v),sizeof(m))#define chmin(X,Y) ((X)>(Y)?X=(Y),true:false)#define chmax(X,Y) ((X)<(Y)?X=(Y),true:false)#define fst first#define snd second#define UNIQUE(x) (x).erase(unique(all(x)),(x).end())template<class T> ostream &operator<<(ostream &os, const vector<T> &v){int n=v.size();rep(i,n)os<<v[i]<<(i==n-1?"":" ");return os;}#define N 100010ll inv[N], fact[N], ifact[N];void init_fact(ll n = N){inv[1] = 1;for(int i = 2; i < n;i++) inv[i] = inv[mod%i] * (mod - mod/i) % mod;fact[0] = ifact[0] = 1;for(int i = 1; i < n; i++){fact[i] = (fact[i-1]*i)%mod;ifact[i]=(ifact[i-1]*inv[i])%mod;}}ll comb(ll a, ll b){return fact[a+b]*ifact[a]%mod*ifact[b]%mod;}ll comb_nk(ll n, ll k){return comb(n-k, k);}ll pow_mod(ll a, ll r, ll m){ll x = 1;while(r){if(r&1) (x*=a)%=m;(a*=a)%=m;r>>=1;}return x;}int main(){ll x;cin>>x;init_fact();ll res = 0;rep(i, x){res += comb(x-i, i)*pow_mod(x-i, i, mod)%mod;res %= mod;}cout<<res<<endl;return 0;}