// #define _GLIBCXX_DEBUG // for STL debug (optional)
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <cstring>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <utility>
#include <algorithm>
#include <map>
#include <set>
#include <complex>
#include <cmath>
#include <limits>
#include <cfloat>
#include <climits>
#include <ctime>
#include <cassert>
#include <numeric>
#include <fstream>
#include <functional>
#include <bitset>
using namespace std;

#define debug(...) fprintf(stderr, __VA_ARGS__)
#define int long long int
 
template<typename T> void chmax(T &a, T b) {a = max(a, b);}
template<typename T> void chmin(T &a, T b) {a = min(a, b);}
template<typename T> void chadd(T &a, T b) {a = a + b;}
 
typedef pair<int, int> pii;
typedef long long ll;
 
int dx[] = {0, 0, 1, -1};
int dy[] = {1, -1, 0, 0};
const ll INF = 1001001001001001LL;
const ll MOD = 1000000007LL;

// Abstract Segment Tree (抽象・汎化セグメント木)
// Verified: AOJ DSL_2_A: Range Minimum Query (int型のみ)

const int SIZE = 1 << 18;
template<typename T>
struct segtree {
    // ノード、単位元
    T node[2*SIZE], I;

    // オペレーション (update, query の 2 つが必要？)
    T (*upd_f)(T, T), (*qry_f)(T, T);

    // 演算子と単位元をセットし、全ての node を単位元で初期化
    segtree(T (*op1)(T, T), T (*op2)(T, T), T X) {
        upd_f = op1;
        qry_f = op2;
        I = X;
        fill(node, node+2*SIZE, I);
    }

    // 配列 vec の値で初期化
    void init(vector<T> vec) {
        int N = (int)vec.size();
        for(int i=0; i<N; i++) {
            node[SIZE-1+i] = vec[i];
        }
        for(int i=SIZE-2; i>=0; i--) {
            node[i] = qry_f(node[2*i+1], node[2*i+2]);
        }
    }

    // 場所 idx に値 val を反映 (upd_f を用いて update)
    void update(int idx, T val) {
        idx += SIZE - 1;
        node[idx] = upd_f(node[idx], val);
        while(idx > 0) {
            // 登るときは qry_f を使うことに注意！！！
            idx = (idx - 1) / 2;
            node[idx] = qry_f(node[2*idx+1], node[2*idx+2]);
        }
    }

    // 半開区間 [a, b) に対してクエリを投げる
    // (qry_f を用いて処理)
    T query(int a, int b, int l=0, int r=SIZE, int k=0) {
        if(b <= l || r <= a) return I;
        if(a <= l && r <= b) return node[k];
        int mid = (l + r) / 2;
        T vl = query(a, b, l, mid, 2*k+1);
        T vr = query(a, b, mid, r, 2*k+2);
        return qry_f(vl, vr);
    }
};

// [[つかいかた]]
// update 用と query 用の関数と単位元 (関数を噛ませても結果が変わらないもの) を用意して宣言する
// 下の例だと、update 用の関数が upd で、query 用の関数が fnd で、単位元が INT_MAX

using PAIR = pair<int, int>;
PAIR upd(PAIR a, PAIR b) {return b;}
PAIR fnd(PAIR a, PAIR b) {return min(a, b);}
const PAIR INIT(INF, -1);
segtree<PAIR> seg_l(upd, fnd, INIT), seg_r(upd, fnd, INIT);

struct UnionFind {
    int node[200010];
    UnionFind() {
        memset(node, -1, sizeof(node));
    }
    int find(int x) {
        return node[x] < 0 ? x : node[x] = find(node[x]);
    }
    bool unite(int x, int y) {
        x = find(x), y = find(y);
        if(x == y) return false;

        node[x] += node[y];
        node[y] = x;
        return true;
    }
};
UnionFind uf;

signed main() {
    int N, D; cin >> N >> D;

    vector< tuple<int, int> > info;
    for(int i=0; i<N; i++) {
        int val; cin >> val;
        info.emplace_back(val, i);
    }

    sort(info.begin(), info.end());
    vector< tuple<int, int, int> > edges;
    for(int i=0; i<N; i++) {
        int val, idx; tie(val, idx) = info[i];

        // 左を見る
        auto elem_l = seg_l.query(0, idx);
        if(elem_l.first != INF) {
            elem_l.first -= D * (N - idx);
            elem_l.first += val;

            int u = idx, v = elem_l.second, cost = elem_l.first;
            edges.emplace_back(cost, u, v);
        }
        int valL = D * (N - idx) + val;
        seg_l.update(idx, make_pair(valL, idx));

        // 右を見る
        auto elem_r = seg_r.query(idx, N);
        if(elem_r.first != INF) {
            elem_r.first -= D * (idx + 1);
            elem_r.first += val;

            int u = idx, v = elem_r.second, cost = elem_r.first;
            edges.emplace_back(cost, u, v);
        }
        int valR = D * (idx + 1) + val;
        seg_r.update(idx, make_pair(valR, idx));
    }

    sort(edges.begin(), edges.end());
    int ans = 0;
    for(auto e : edges) {
        int u, v, w; tie(w, u, v) = e;
        if(uf.unite(u, v)) ans += w;
    }
    cout << ans << endl;
    return 0;
}


Measured memory_bytes: 38632