#include #include using namespace std; using namespace atcoder; istream &operator>>(istream &is, modint &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint &a) { return os << a.val(); } istream &operator>>(istream &is, modint998244353 &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint998244353 &a) { return os << a.val(); } istream &operator>>(istream &is, modint1000000007 &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint1000000007 &a) { return os << a.val(); } typedef long long ll; typedef vector> Graph; typedef pair pii; typedef pair pll; #define FOR(i,l,r) for (int i = l;i < (int)(r); i++) #define rep(i,n) for (int i = 0;i < (int)(n); i++) #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define my_sort(x) sort(x.begin(), x.end()) #define my_max(x) *max_element(all(x)) #define my_min(x) *min_element(all(x)) template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } const int INF = (1<<30) - 1; const ll LINF = (1LL<<62) - 1; const double PI = acos(-1); vector di = {1,0,-1,0}; vector dj = {0,1,0,-1}; #ifdef LOCAL # include # define debug(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define debug(...) (static_cast(0)) #endif //https://drken1215.hatenablog.com/entry/2018/06/08/210000 //COMinit()を忘れない!!! const ll NMAX = 202020; const ll MOD = 998244353; //const int MOD = 1e9+7; ll fac[NMAX],finv[NMAX],inv[NMAX]; void COMinit(){ fac[0] = fac[1] = 1LL; finv[0] = finv[1] = 1LL; inv[1] = 1LL; for (int i=2;i> N; vector A(N); rep(i,N) cin >> A[i]; mint ans = 0; if(N <= 3){ for(int i = 0; i < N; i++){ for(int l = 0; l <= i; l++){ for(int r = 0; r <= N - i - 1; r++){ mint tmp = A[i]; tmp *= l + 1; tmp *= r + 1; tmp *= nCr(i, l); tmp *= nCr(N - i - 1, r); ans += tmp; debug(i, l, r, tmp); } } } } else{ for(int i = 0; i < N; i++){ mint tmp = A[i]; tmp *= (i + 2) * (-mint(2).pow(N - 3)) * (i - N - 1); ans += tmp; } } cout << ans << endl; }