結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
|
提出日時 | 2020-06-19 16:06:34 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 14 ms / 2,000 ms |
コード長 | 1,210 bytes |
コンパイル時間 | 1,417 ms |
コンパイル使用メモリ | 168,980 KB |
実行使用メモリ | 8,280 KB |
最終ジャッジ日時 | 2024-07-03 13:27:22 |
合計ジャッジ時間 | 2,922 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include<bits/stdc++.h>using namespace std;typedef long long ll;typedef long double ld;#define rep(i,a,b) for(int i=a;i<b;i++)#define rrep(i,a,b) for(int i=a;i>=b;i--)template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }const ll INF=1LL<<60;const int inf=(1<<30)-1;const int mod=1e9+7;int dx[8]={1,0,-1,0,-1,-1,1,1};int dy[8]={0,1,0,-1,-1,1,-1,1};const int nmax=200005;ll fac[nmax],finv[nmax],inv[nmax];void COMinit(){fac[0]=fac[1]=1;finv[0]=finv[1]=1;inv[1]=1;for(int i=2;i<nmax;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;}}ll com(int n,int k){if(n<k||n<0||k<0){return 0;}return fac[n]*(finv[k]*finv[n-k]%mod)%mod;}ll pow_mod(ll n,ll k,ll m){ll res=1;for(;k>0;k>>=1){if(k&1){res=(res*n)%m;}n=(n*n)%m;}return res;}int main(){cin.tie(0);ios::sync_with_stdio(false);ll n;cin >> n;ll ans=0;COMinit();rep(i,0,n+1){ans+=com(n,i)*pow_mod(n-i,i,mod)%mod;ans%=mod;}cout << ans << endl;}