// This code is very suspicious #include using namespace std; using VI = vector; using VVI = vector; using PII = pair; using LL = long long; using VL = vector; using VVL = vector; using PLL = pair; using VS = vector; #define ALL(a) begin((a)),end((a)) #define RALL(a) (a).rbegin(), (a).rend() #define PB push_back #define EB emplace_back #define MP make_pair #define SZ(a) int((a).size()) #define SORT(c) sort(ALL((c))) #define RSORT(c) sort(RALL((c))) #define UNIQ(c) (c).erase(unique(ALL((c))), end((c))) #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define REP(i,n) FOR(i,0,n) #define FF first #define SS second #define DUMP(x) cout<<#x<<":"<<(x)< istream& operator>>(istream& is, pair& p){ return is >> p.FF >> p.SS; } template istream& operator>>(istream& is, vector& xs){ for(auto& x: xs) is >> x; return is; } template ostream& operator<<(ostream& os, const pair& p){ return os << p.FF << " " << p.SS; } template ostream& operator<<(ostream& os, const vector& xs){ for(unsigned int i=0;i void maxi(T& x, T y){ if(x < y) x = y; } template void mini(T& x, T y){ if(x > y) x = y; } const double EPS = 1e-10; const double PI = acos(-1.0); const LL MOD = 1e9+7; const LL INF = 4e18; template> struct Line{ static const T NaN = numeric_limits::max(); T a, b; Line(T a_=0, T b_=0) : a(a_), b(b_) {} mutable function next; bool operator<(const Line& rhs) const { if(rhs.b != NaN) return Comp()(a, rhs.a); const Line* nxt = next(); return Comp()(a * rhs.a + b, nxt->a * rhs.a + nxt->b); } }; template> struct ConvexHullTrick{ using LineT = Line; multiset lines; ConvexHullTrick(){ } bool bad(typename multiset::iterator it){ auto nxt = next(it); if(it == begin(lines)){ if(nxt == end(lines)) return false; return it->a == nxt->a && Comp()(nxt->b, it->b); } auto prv = prev(it); if(nxt == end(lines)) return prv->a == it->a && Comp()(nxt->b, it->b); return Comp()(-(double)(prv->a - it->a)*(nxt->b - it->b), -(double)(nxt->a - it->a)*(prv->b - it->b)); } // y = ax + b void add(T a, T b){ auto it = lines.emplace(a, b); if(bad(it)){ lines.erase(it); return; } it->next = [=](){ return &*next(it); }; auto tmp = next(it); while(tmp != end(lines) && bad(tmp)) tmp = lines.erase(tmp); --tmp; while(tmp != begin(lines)){ --tmp; if(!bad(tmp)) break; tmp = lines.erase(tmp); } } T f(const LineT& l, T x){ return l.a * x + l.b; } T getmin(T x){ auto it = lines.lower_bound(LineT(x, LineT::NaN)); return f(*it, x); } }; int main(){ cin.tie(0); ios_base::sync_with_stdio(false); LL N, A, B, W; cin >> N >> A >> B >> W; VL ds(N+1); REP(i,N) cin >> ds[i+1]; VL dp(N+1, INF); ConvexHullTrick cht; LL ans = N*(N+1)/2*B - N*A; for(LL i=1;i<=N;++i){ dp[i] = i * (i-1) / 2 * B - (i-1)*A + ds[i]; if(i>1) mini(dp[i], ds[i] - (i-1)*A + (i*i-i)/2*B + cht.getmin(i)); cht.add(-B * i, dp[i] + i*A + i*(i+1)/2*B); mini(ans, dp[i] + (N-i)*(N-i+1)/2*B - (N-i)*A); } cout << ans+W << endl; return 0; }