結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
![]() |
提出日時 | 2019-09-17 12:27:33 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 45 ms / 2,000 ms |
コード長 | 2,232 bytes |
コンパイル時間 | 739 ms |
コンパイル使用メモリ | 107,060 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-07 10:28:32 |
合計ジャッジ時間 | 1,847 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
#include<iostream> #include<string> #include<cstdio> #include<vector> #include<cmath> #include<algorithm> #include<functional> #include<iomanip> #include<queue> #include<ciso646> #include<random> #include<map> #include<set> #include<complex> #include<bitset> #include<stack> #include<unordered_map> #include<utility> #include<tuple> using namespace std; typedef long long ll; typedef unsigned int ui; const ll mod = 1000000007; const ll INF = (ll)1000000007 * 1000000007; typedef pair<int, int> P; #define stop char nyaa;cin>>nyaa; #define rep(i,n) for(int i=0;i<n;i++) #define per(i,n) for(int i=n-1;i>=0;i--) #define Rep(i,sta,n) for(int i=sta;i<n;i++) #define Per(i,sta,n) for(int i=n-1;i>=sta;i--) #define rep1(i,n) for(int i=1;i<=n;i++) #define per1(i,n) for(int i=n;i>=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef long double ld; typedef complex<ld> Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair<ll, ll> LP; int n; struct mod_fac{ private: vector<ll> f,i_f; int n;ll MOD; public: void init(int n,ll m){ MOD=m; f.resize(n+1,0); i_f.resize(n+1,0); f[0]=1; rep(i,n){ f[i+1]=f[i]*(i+1)%MOD; } i_f[n]=power(f[n],MOD-2); per(i,n){ i_f[i]=(i_f[i+1]*(i+1))%MOD; } } ll power(ll a ,ll x){ if (x==0) return 1; if (x==1) return a%MOD; if (x==2) return (a*a)%MOD; if (x%2==1) return (a*power(power(a,(x-1)/2),2))%MOD; return (power(power(a,x/2),2))%MOD; } ll fac(int x){ return f[x]; } ll invfac(int x){ return i_f[x]; } ll comb(int m,int k){ if (m<0 or k<0) return 0; if (m<k) return 0; return (((f[m]*i_f[k])%MOD)*i_f[m-k])%MOD; } }; void solve(){ cin >> n; mod_fac m;m.init(n,mod); ll ans=0; rep(i,n+1){ ans+=m.comb(n,i)*m.power(i,n-i)%mod; ans%=mod; } cout << ans << endl; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(50); solve(); }