結果
問題 | 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();}