/* -*- coding: utf-8 -*- * * 1681.cc: No.1681 +-* - yukicoder */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; /* constant */ const int MAX_N = 200000; const int MOD = 1000000007; /* typedef */ typedef long long ll; template struct MI { int v; MI(): v() {} MI(int _v): v(_v) {} MI operator+(const MI m) const { return MI((v + m.v) % MOD); } MI operator-(const MI m) const { return MI((v + MOD - m.v) % MOD); } MI operator*(const MI m) const { return MI((ll)v * m.v % MOD); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; typedef MI mi; /* global variables */ int as[MAX_N]; mi es[MAX_N + 1]; /* subroutines */ /* main */ int main() { int n; scanf("%d", &n); for (int i = 0; i < n; i++) scanf("%d", as + i); es[0] = 1; for (int i = 0; i < n; i++) es[i + 1] = es[i] * 3; mi sum = 0, p = as[0]; for (int i = 1; i < n; i++) { sum += p * 2 * es[n - 1 - i]; p *= as[i]; } sum += p; printf("%d\n", sum.v); return 0; }