結果
問題 | No.1081 和の和 |
ユーザー | ate |
提出日時 | 2020-06-19 23:06:30 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 175 ms / 2,000 ms |
コード長 | 3,496 bytes |
コンパイル時間 | 2,511 ms |
コンパイル使用メモリ | 196,224 KB |
最終ジャッジ日時 | 2025-01-11 07:48:07 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 8 |
ソースコード
#include<bits/stdc++.h> using namespace std; template<int64_t MOD> struct mod_int{ int64_t val; constexpr mod_int():val(0){} constexpr mod_int(int64_t x)noexcept:val(x>=0?x%MOD:(MOD-(-x)%MOD)%MOD){} constexpr int64_t value()const noexcept{return val;} constexpr int64_t get_MOD()const noexcept{return MOD;} constexpr mod_int operator+(mod_int const& rhs){ return mod_int(*this)+=rhs; } constexpr mod_int operator-(mod_int const& rhs){ return mod_int(*this)-=rhs; } constexpr mod_int operator*(mod_int const& rhs){ return mod_int(*this)*=rhs; } constexpr mod_int operator/(mod_int const& rhs){ return mod_int(*this)/=rhs; } constexpr mod_int &operator+=(mod_int const& rhs)noexcept{ val += rhs.val; if(val>=MOD)val-=MOD; return *this; } constexpr mod_int &operator-=(mod_int const& rhs)noexcept{ if(val<rhs.val)val+=MOD; val-=rhs.val; return *this; } constexpr mod_int &operator*=(mod_int const& rhs)noexcept{ (val*=rhs.val)%=MOD; return *this; } constexpr mod_int &operator/=(mod_int rhs)noexcept{ *this*=rhs.inverse(); return *this; } constexpr bool operator==(mod_int const& rhs)const{ return val==rhs.val; } constexpr bool operator!=(mod_int const& rhs)const{ return !(val==rhs.val); } constexpr bool operator<(mod_int const& rhs)const{ return val<rhs.val; } constexpr bool operator>(mod_int const& rhs)const{ return val>rhs.val; } constexpr bool operator<=(mod_int const& rhs)const{ return !(val>rhs.val); } constexpr bool operator>=(mod_int const& rhs)const{ return !(val<rhs.val); } constexpr friend std::ostream &operator<<(std::ostream& os, mod_int const& mi){ return os<<mi.val; } constexpr friend std::istream &operator>>(std::istream& is, mod_int & mi){ int64_t t=0;is>>t; mi = mod_int<MOD>(t); return is; } constexpr mod_int inverse(){ int64_t a=val,b=MOD,u=1,v=0,t=0; while(b>0){ t=a/b; std::swap(a-=t*b,b); std::swap(u-=t*v,v); } return mod_int(u); } constexpr mod_int mpow(int64_t n)const{ mod_int res(1),e(*this); while(n){ if(n&1)res*=e; e*=e; n>>=1; } return res; } }; using mint = mod_int<1000000007>; mint mpow(int64_t p,int64_t n){ mint res=1; mint e = p; while(n){ if(n&1)res*=e; e*=e; n>>=1; } return res; } struct combination{ using i64 = int_fast64_t; static constexpr i64 COM_MAX = 2e6; static constexpr i64 MOD = 1e9+7; vector<i64> fac,inv,finv; combination(){ fac.resize(COM_MAX); inv.resize(COM_MAX); finv.resize(COM_MAX); fac[0] = fac[1] = 1; inv[1] = 1; finv[0] = finv[1] = 1; for(int i=2;i<COM_MAX;++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; } } template<class T> mint nCk(T n,T k){ if(n<k)return 0; if(n<0||k<0)return 0; if(n>=COM_MAX){ mint res1 = mint(1); mint res2 = mint(fac[k]); for(int i=1;i<=k;++i){ res1 *= n-(i-1); } return res1/res2; } return mint(fac[n]*(finv[k]*finv[n-k]%MOD)); } template<class T> mint nHk(T n,T k){ return nCk(n+k-1,k-1); } template<class T> mint operator()(T n,T k){ return nCk(n,k); } }; signed main(){ int n; cin>>n; combination com; mint ans; for(int i=1;i<=n;++i){ mint ai; cin>>ai; ans += ai*com(n-1,i-1); } cout<<(ans)<<endl; }