結果
問題 | No.797 Noelちゃんとピラミッド |
ユーザー |
|
提出日時 | 2019-04-07 13:08:59 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 94 ms / 2,000 ms |
コード長 | 3,543 bytes |
コンパイル時間 | 1,195 ms |
コンパイル使用メモリ | 141,048 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-11-30 14:02:28 |
合計ジャッジ時間 | 5,887 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 60 |
ソースコード
#include <iostream> #include <set> #include <queue> #include <stack> #include <vector> #include <algorithm> #include <map> #include <chrono> #include <math.h> using namespace std; using lli = long long int; using Vint = std::vector<int>; using Vlli = std::vector<lli>; using Wint = std::vector<Vint>; using Wlli = std::vector<Vlli>; using Vbool = std::vector<bool>; using pii = std::pair<int, int>; using pll = std::pair<lli, lli>; constexpr int MOD = 1e9 + 7; constexpr int INFi = 2e9 + 1; constexpr lli INFl = (lli)(9e18) + 1; const vector<pii> DXDY = {std::make_pair(1, 0), std::make_pair(-1, 0), std::make_pair(0, 1), std::make_pair(0, -1)}; constexpr char BR = '\n'; #define FOR(i, a, b) for(int (i) = (a); (i) < (b); (i)++) #define FOReq(i, a, b) for(int (i) = (a); (i) <= (b); (i)++) #define rFOR(i, a, b) for(int (i) = (b); (i) >= (a); i--) #define FORstep(i, a, b, step) for(int (i) = (a); i < (b); i += (step)) #define REP(i, n) FOR(i, 0, n) #define rREP(i, n) rFOR(i, 0, (n-1)) #define vREP(ele, vec) for(auto &(ele) : (vec)) #define vREPcopy(ele, vec) for(auto (ele) : (vec)) #define SORT(A) std::sort((A).begin(), (A).end()) template <class T> inline int argmin(vector<T> vec){return min_element(vec.begin(), vec.end()) - vec.begin();} template <class T> inline int argmax(vector<T> vec){return max_element(vec.begin(), vec.end()) - vec.begin();} template <class T> inline void chmax(T &a, T b){a = max(a, b);} template <class T> inline void chmin(T &a, T b){a = min(a, b);} inline int toInt(string &s){int res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;} inline int toInt(const string s){int res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;} inline long long int toLong(string &s){lli res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;} inline long long int toLong(const string s){lli res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;} template <class T> inline std::string toString(T n){ if(n == 0) return "0"; std::string res = ""; if(n < 0){n = -n;while(n != 0){res += (char)(n % 10 + '0'); n /= 10;} std::reverse(res.begin(), res.end()); return '-' + res;} while(n != 0){res += (char)(n % 10 + '0'); n /= 10;} std::reverse(res.begin(), res.end()); return res; } // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ class Combination_Mod{ private: int n, mod; vector<long long int> S, T; long long int power(long long int a, int index){ if(index == 0) return 1; if(index % 2 == 0) return power((a * a) % mod, index / 2); return (a * power((a * a) % mod, index / 2)) % mod; } long long int calculate_inv(int a){ return power(a, mod - 2); } void init_S(void){ S.at(0) = 1; for(int i = 1; i <= n; i++) S.at(i) = (S.at(i - 1) * i) % mod; } void init_T(void){ T.at(0) = 1; for(int i = 1; i <= n; i++) T.at(i) = calculate_inv(S.at(i)); } public: Combination_Mod(int nn, int modmod){ n = nn; mod = modmod; S.resize(nn + 1); T.resize(nn + 1); init_S(); init_T(); } ~Combination_Mod(void){}; long long int nCr(int nn, int rr){ return ((S.at(nn) * T.at(rr) % mod) * T.at(nn - rr)) % mod; } }; int main(void){ // a1, a2, ..., a_n // ans = ∑ (n-1)_C_k * a_k int n; scanf("%d\n", &n); lli res = 0; Vint A(n); vREP(ele, A) scanf("%d", &ele); Combination_Mod C(n, MOD); REP(i, n) res = (res + (C.nCr(n-1, i) * A.at(i) % MOD)) % MOD; printf("%lld\n", res); return 0; }