// template version 1.2
using namespace std;
#include <iostream>
#include <bits/stdc++.h>

#define infile "../test/sample-1.in"
#define int long long
// #define INF 2147483647
#define INF 1000000000000000000LL
#define MOD 1000000007LL

// {{{ define basic macro
#define _overload3(_1,_2,_3,name,...) name
#define _rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(int i=(int)(a);i<(int)(b);++i)
#define rep(...) _overload3(__VA_ARGS__,repi,_rep,)(__VA_ARGS__)
#define _rrep(i,n) rrepi(i,0,n)
#define rrepi(i,a,b) for(int i=(int)(b-1);i>=(int)(a);--i)
#define rrep(...) _overload3(__VA_ARGS__,rrepi,_rrep,)(__VA_ARGS__)

#define foreach(x,a) for(auto& (x) : (a) )
#define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++)

#define all(x) (x).begin(),(x).end()
#define sum(v) accumulate(all(v), 0LL)
#define sz(x) ((int)(x).size())
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define uni(x) sort(all(x));x.erase(unique(all(x)),x.end())
#define ten(n) ((int)1e##n)

template<class T> inline void chmax(T &a, const T &b) { if(a < b) a = b; }
template<class T> inline void chmin(T &a, const T &b) { if(a > b) a = b; }

template<typename A, size_t N, typename T> // n次元配列の初期化。第２引数の型のサイズごとに初期化していく。
void Fill(A (&array)[N], const T &val){
    std::fill( (T*)array, (T*)(array+N), val );
}

template <class T = int>
T in() {T x; cin>>x; return (x);}

struct Fast {
  Fast(){
    std::cin.tie(0);
    ios::sync_with_stdio(false);
  }
} fast;
// }}}
//{{{ dump macro
#ifdef PCM
  #include "dump.hpp"
#else
  #define dump(...) 42
  #define dump_1d(...) 42
  #define dump_2d(...) 42
#endif
//}}}
//{{{ others
typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef long double ld;
typedef pair<int,int> pii;
typedef tuple<int,int,int> iii;
template<typename T>
using PQ = priority_queue<T, vector<T>, greater<T>>;

int dx[]={1, -1, 0, 0};
int dy[]={0, 0, 1, -1};
//}}}

// mod lib {{{
#define MAX_P 200005
ll fact[MAX_P]; // fact[n]: n!
ll inv[MAX_P];  // inv[n]: n^(-1)
ll finv[MAX_P]; // finv[n]: (n!)^(-1)

void init(ll mod){  // 階乗関連テーブルの初期化 O(n)
  fact[0]=1;
  for(ll i=1;i<MAX_P;i++)
    fact[i]=(fact[i-1]*i)%mod;

  inv[1]=1;
  for(ll i=2;i<MAX_P;i++)
    inv[i]=inv[mod%i]*(mod-mod/i)%mod;

  finv[0]=1;  //　高速(10倍程度)な逆元階乗テーブルの作成。
  for(ll i=1;i<MAX_P;i++)
    finv[i]=finv[i-1]*inv[i]%mod;
}

ll extgcd(ll a,ll b,ll& x,ll& y){ //O(log(max(a,b)))
  ll d=a;
  if(b!=0){
    d=extgcd(b,a%b,y,x);
    y-=(a/b)*x;
  }else{
    x=1;y=0;
  }
  return d;
}
ll mod_inverse(ll a,ll mod){  // ユークリッドの拡張互除法による逆元
  ll x,y;
  extgcd(a,mod,x,y);
  return (mod+x%mod)%mod;
}

ll mod_pow(ll x,ll n,ll mod){ // 二分累乗法による階乗
  ll res=1;
  while(n){
    if(n&1) (res*=x)%=mod;
    (x*=x)%=mod;
    n>>=1;
  }
  return res;
}

ll mod_inverse2(ll a,ll mod){  // フェルマーの小定理による逆元
  return mod_pow(a,mod-2,mod);
}


ll mod_fact(ll n,ll mod,ll& e){ // n!=a*p^eとしたときのa mod pを返す。O(p)
  e=0;
  if(n==0) return 1;
  ll res=mod_fact(n/mod,mod,e);
  e+=n/mod;
  if(n/mod%2!=0)return res*(mod-fact[n%mod]) %mod;
  return res*fact[n%mod]%mod;
}

ll mod_comb(ll n,ll k,ll mod){
  if(n==k||k==0) return 1;
  ll e1,e2,e3;
  ll a1=mod_fact(n,mod,e1),a2=mod_fact(k,mod,e2),a3=mod_fact(n-k,mod,e3);
  if(e1>e2+e3) return 0;
  return a1*mod_inverse(a2*a3%mod,mod)%mod;
}

ll mod_comb2(ll n,ll k,ll mod){
  ll res=1;
  for(ll i=0;i<k;i++){
    res*=(n-i)%mod;
    res%=mod;
    res*=mod_inverse(i+1,mod);
    res%=mod;
  }
  return res;
}

//only for prime mod
ll mod_comb_pm(ll n,ll k,ll mod){  // テーブルがあればO(1)
  if(k<0||k>n) return 0;
  return fact[n]*finv[k]%mod*finv[n-k]%mod;
}

ll montmort(ll n,ll mod){
  ll res=0,inv=1;
  for(ll k=2;k<=n;k++){
    (inv*=mod_inverse(k,mod))%=mod;
    if(k%2) (res+=mod-inv)%=mod;
    else (res+=inv)%=mod;
  }
  for(ll i=1;i<=n;i++)
    (res*=i)%=mod;
  return res;
}

// calculate P(t) from given points in [0,N]
ll LagrangePolynomial(vector<ll> &y,ll t,const ll mod){
  init(mod);
  ll n=y.size()-1;
  ll num=1;
  for(ll i=0;i<=n;i++)
    num=num*((t-i)%mod)%mod;
  ll res=0;
  for(ll i=0;i<=n;i++){
    ll tmp=(y[i]*num%mod)*mod_inverse((t-i)%mod,mod)%mod;
    tmp=tmp*finv[i]%mod;
    tmp=tmp*finv[n-i]%mod;
    if((n-i)&1) tmp=mod-tmp;
    res=(res+tmp)%mod;
  }
  return res;
}
// }}}

int solve(){
  int n;cin>>n;
  int a[n]; rep(i, n) {cin>>a[i];}

  // prepare n!
  int fac[n+1];
  int ps[n+1];
  fac[0]=1;
  ps[0]=1;
  rep(i, 1, n+1){
    fac[i] = fac[i-1]*i;
    fac[i] %= MOD;
    ps[i] = (int)(M_E*fac[i]);
    ps[i] %= MOD;
  }

  int g[n];
  int G[n];
  rep(d, n){
    // assert(fac[n]%(1+d)==0);
    g[d] = fac[n] * mod_inverse(1+d, MOD);
    g[d]%=MOD;
  }
  G[0] = g[0];
  rep(d, 1, n){
    G[d] = G[d-1]+g[d];
    G[d] %= MOD;
  }
  int res = 0;
  rep(i, n){
    int tmp=0;
    tmp += G[i] + G[n-i-1];
    tmp -= g[0];
    tmp %= MOD;
    res += a[i]*tmp;
    res %= MOD;
    dump(res);
  }

  cout << res << endl;
  return 0;
}

int test(){ //{{{
  return 0;
}
//}}}

signed main() { //{{{
#ifdef INPUT_FROM_FILE
  std::ifstream in(infile);
  std::cin.rdbuf(in.rdbuf());
#endif
#ifdef PCM
  test();
#endif
  solve();
  return 0;
} //}}}


Measured memory_bytes: 4224