#include <bits/stdc++.h>
using namespace std;
const int block_size = 320;
const long long mod = 1e9 + 7;
const long long inf = 1e9 + 7;
const long double eps = 1e-9;
const double PI = atan(1) * 4;
template <typename T>
inline int sign(const T &a) {
  if (a < 0) return -1;
  if (a > 0) return 1;
  return 0;
}
string to_string(string s) { return '"' + s + '"'; }
string to_string(const char *s) { return to_string((string)s); }
string to_string(bool b) { return (b ? "true" : "false"); }
template <typename A, typename B>
string to_string(pair<A, B> p) {
  return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template <typename A>
string to_string(A v) {
  bool first = true;
  string res = "{";
  for (const auto &x : v) {
    if (!first) {
      res += ", ";
    }
    first = false;
    res += to_string(x);
  }
  res += "}";
  return res;
}
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
  cerr << " " << to_string(H);
  debug_out(T...);
}
template <typename T, typename S>
inline bool upmin(T &a, const S &b) {
  return a > b ? a = b, 1 : 0;
}
template <typename T, typename S>
inline bool upmax(T &a, const S &b) {
  return a < b ? a = b, 1 : 0;
}
template <typename T>
inline void in(T &x) {
  x = 0;
  T f = 1;
  char ch = getchar();
  while (!isdigit(ch)) {
    if (ch == '-') f = -1;
    ch = getchar();
  }
  while (isdigit(ch)) {
    x = x * 10 + ch - '0';
    ch = getchar();
  }
  x *= f;
}
long long twop(int x) { return 1LL << x; }
template <typename A, typename B>
inline void in(A &x, B &y) {
  in(x);
  in(y);
}
template <typename A, typename B, typename C>
inline void in(A &x, B &y, C &z) {
  in(x);
  in(y);
  in(z);
}
template <typename A, typename B, typename C, typename D>
inline void in(A &x, B &y, C &z, D &d) {
  in(x);
  in(y);
  in(z);
  in(d);
}
long long gcd(long long a, long long b) { return b == 0 ? a : gcd(b, a % b); }
namespace SOLVE {
long long a[200010];
template <int maxn>
struct dsu {
  int f[maxn];
  dsu() { init(); }
  void init() { memset(f, -1, sizeof(f)); }
  int fa(int c) {
    if (f[c] < 0) return c;
    return f[c] = fa(f[c]);
  }
  void merge(int a, int b) {
    a = fa(a);
    b = fa(b);
    if (a == b) return;
    if (f[a] < f[b]) swap(a, b);
    f[b] += f[a];
    f[a] = b;
  }
  bool same(int a, int b) { return fa(a) == fa(b); }
};
dsu<200010> bcj;
void main() {
  long long n, d;
  in(n, d);
  for (long long i = 0; i < n; i++) in(a[i]);
  pair<long long, long long> cur = make_pair(a[0], 0);
  vector<pair<long long, pair<long long, long long>>> edges;
  for (long long i = 1; i < n; i++) {
    edges.push_back(
        make_pair(cur.first + i * d + a[i], make_pair(cur.second, i)));
    upmin(cur, make_pair(a[i] - i * d, i));
  }
  cur = make_pair((n - 1) * d + a[n - 1], n - 1);
  for (long long i = n - 2; i >= 0; i--) {
    edges.push_back(
        make_pair(cur.first - i * d + a[i], make_pair(cur.second, i)));
    upmin(cur, make_pair(a[i] + i * d, i));
  }
  sort(edges.begin(), edges.end());
  long long ans = 0;
  for (auto e : edges) {
    if (bcj.same(e.second.first, e.second.second))
      continue;
    else {
      bcj.merge(e.second.first, e.second.second);
      ans += e.first;
    }
  }
  cout << ans << endl;
}
}  // namespace SOLVE
signed main() {
  int t;
  t = 1;
  while (t--) {
    SOLVE::main();
  }
  return 0;
}
