結果
問題 | No.797 Noelちゃんとピラミッド |
ユーザー |
|
提出日時 | 2019-03-15 22:10:54 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 16 ms / 2,000 ms |
コード長 | 1,417 bytes |
コンパイル時間 | 1,278 ms |
コンパイル使用メモリ | 72,380 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-01 21:04:19 |
合計ジャッジ時間 | 3,394 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 60 |
ソースコード
#include<iostream> #include<string> #include<iomanip> #include<cmath> #include<vector> #include<algorithm> using namespace std; #define int long long #define rep(i,n) for(int i = 0; i < (n); i++) #define INF ((long long)1e18) #define MOD ((int)1e9+7) #define endl "\n" #define yn(f) ((f)?"Yes":"No") #define YN(f) ((f)?"YES":"NO") #define MAX 110000 #define MAX_VAL (int)(1e5+10) long long fac[MAX_VAL], mmi[MAX_VAL]; void factorial_mod(){ fac[0]=fac[1]=1; for(long long i = 2; i < MAX_VAL; fac[i]%=MOD,i++) fac[i] = fac[i-1]*(i%MOD); } long long power_mod(long long x, long long n){ long long ans = 1; for(;n;n>>=1,x*=x,ans%=MOD,x%=MOD) if(n&1)ans*=x; return ans%MOD; } void exgcd(int a, int b, int &x, int &y){ if(b == 0){ x = 1; y = 0; return ; } exgcd(b,a%b,y,x); y -= a/b * x; } void modular_multiplicatibe_inverse(){ int x, y; exgcd(fac[MAX_VAL-1],MOD,x,y); mmi[MAX_VAL-1] = x; for(long long i = MAX_VAL-2; i >= 0; mmi[i]%=MOD,i--) mmi[i] = mmi[i+1]*((i+1)%MOD); } long long combination_mod(long long n, long long r){ return fac[n]*(mmi[r]*mmi[n-r]%MOD)%MOD; } signed main(){ cin.tie(0); ios::sync_with_stdio(false); cout<<fixed<<setprecision(10); int N; int ans = 0, a; cin>>N; factorial_mod(); modular_multiplicatibe_inverse(); for(int i = 0; i < N; i++){ cin>>a; ans += a*combination_mod(N-1,i); ans %= MOD; } cout<<ans<<endl; return 0; }