#include <bits/stdc++.h>
int max(int a, int b) {
  if (a > b) {
    return a;
  }
  return b;
}
int min(int a, int b) {
  if (a < b) {
    return a;
  }
  return b;
}
int zt(int a, int b) { return max(a, b) - min(a, b); }
int round(int a, int b) {
  if ((a % b) * 2 >= b) {
    return (a / b) + 1;
  }
  return a / b;
}
int ceil(int a, int b) {
  if (a % b == 0) {
    return a / b;
  }
  return (a / b) + 1;
}
int gcd(int a, int b) {
  int c;
  while (b != 0) {
    c = a % b;
    a = b;
    b = c;
  }
  return a;
}
int lcm(int a, int b) {
  int c = gcd(a, b);
  a /= c;
  return a * b;
}
int nCr(int a, int b) {
  int i, r = 1;
  for (i = 1; i <= b; i++) {
    r *= (a + 1 - i);
    r /= i;
  }
  return r;
}
int nHr(int a, int b) { return nCr(a + b - 1, b); }
int fact(int a) {
  int i, r = 1;
  for (i = 1; i <= a; i++) {
    r *= i;
  }
  return r;
}
int pow(int a, int b) {
  int i, r = 1;
  for (i = 1; i <= b; i++) {
    r *= a;
  }
  return r;
}
int dsum(int x) {
  int r = 0;
  while (x) {
    r += (x % 10);
    x /= 10;
  }
  return r;
}
int dsumb(int x, int b) {
  int r = 0;
  while (x) {
    r += (x % b);
    x /= b;
  }
  return r;
}
int sankaku(int x) { return ((1 + x) * x) / 2; }
void swap(int *a, int *b) {
  int c;
  c = (*a);
  (*a) = (*b);
  (*b) = c;
}
long long llmax(long long a, long long b) {
  if (a > b) {
    return a;
  }
  return b;
}
long long llmin(long long a, long long b) {
  if (a < b) {
    return a;
  }
  return b;
}
long long llzt(long long a, long long b) { return llmax(a, b) - llmin(a, b); }
long long llround(long long a, long long b) {
  if ((a % b) * 2 >= b) {
    return (a / b) + 1;
  }
  return a / b;
}
long long llceil(long long a, long long b) {
  if (a % b == 0) {
    return a / b;
  }
  return (a / b) + 1;
}
long long llgcd(long long a, long long b) {
  long long c;
  while (b != 0) {
    c = a % b;
    a = b;
    b = c;
  }
  return a;
}
long long lllcm(long long a, long long b) {
  long long c = llgcd(a, b);
  a /= c;
  return a * b;
}
long long llnCr(long long a, long long b) {
  long long i, r = 1;
  for (i = 1; i <= b; i++) {
    r *= (a + 1 - i);
    r /= i;
  }
  return r;
}
long long llnHr(long long a, long long b) { return llnCr(a + b - 1, b); }
long long llfact(long long a) {
  long long i, r = 1;
  for (i = 1; i <= a; i++) {
    r *= i;
  }
  return r;
}
long long llpow(long long a, long long b) {
  long long i, r = 1;
  for (i = 1; i <= b; i++) {
    r *= a;
  }
  return r;
}
long long lldsum(long long x) {
  long long r = 0;
  while (x) {
    r += (x % 10);
    x /= 10;
  }
  return r;
}
long long lldsumb(long long x, long long b) {
  long long r = 0;
  while (x) {
    r += (x % b);
    x /= b;
  }
  return r;
}
long long llsankaku(long long x) { return ((1 + x) * x) / 2; }
void llswap(long long *a, long long *b) {
  long long c;
  c = (*a);
  (*a) = (*b);
  (*b) = c;
}
double dbmax(double a, double b) {
  if (a > b) {
    return a;
  }
  return b;
}
double dbmin(double a, double b) {
  if (a < b) {
    return a;
  }
  return b;
}
double dbzt(double a, double b) { return dbmax(a, b) - dbmin(a, b); }
void dbswap(double *a, double *b) {
  double c;
  c = (*a);
  (*a) = (*b);
  (*b) = c;
}
void chswap(char *a, char *b) {
  char c;
  c = (*a);
  (*a) = (*b);
  (*b) = c;
}
int sortfncsj(const void *a, const void *b) {
  if (*(int *)a > *(int *)b) {
    return 1;
  }
  if (*(int *)a == *(int *)b) {
    return 0;
  }
  return -1;
}
int sortfnckj(const void *a, const void *b) {
  if (*(int *)a < *(int *)b) {
    return 1;
  }
  if (*(int *)a == *(int *)b) {
    return 0;
  }
  return -1;
}
int llsortfncsj(const void *a, const void *b) {
  if (*(long long *)a > *(long long *)b) {
    return 1;
  }
  if (*(long long *)a == *(long long *)b) {
    return 0;
  }
  return -1;
}
int llsortfnckj(const void *a, const void *b) {
  if (*(long long *)a < *(long long *)b) {
    return 1;
  }
  if (*(long long *)a == *(long long *)b) {
    return 0;
  }
  return -1;
}
int dbsortfncsj(const void *a, const void *b) {
  if (*(double *)a > *(double *)b) {
    return 1;
  }
  if (*(double *)a == *(double *)b) {
    return 0;
  }
  return -1;
}
int dbsortfnckj(const void *a, const void *b) {
  if (*(double *)a < *(double *)b) {
    return 1;
  }
  if (*(double *)a == *(double *)b) {
    return 0;
  }
  return -1;
}
int strsortfncsj(const void *a, const void *b) {
  return strcmp((char *)a, (char *)b);
}
int strsortfnckj(const void *a, const void *b) {
  return strcmp((char *)b, (char *)a);
}
int chsortfncsj(const void *a, const void *b) {
  if (*(char *)a > *(char *)b) {
    return 1;
  }
  if (*(char *)a == *(char *)b) {
    return 0;
  }
  return -1;
}
int chsortfnckj(const void *a, const void *b) {
  if (*(char *)a < *(char *)b) {
    return 1;
  }
  if (*(char *)a == *(char *)b) {
    return 0;
  }
  return -1;
}
void shuffledget(int x[], int n) {
  int i, b[524288], p, c;
  for (i = 0; i < n; i++) {
    b[i] = i;
  }
  for (i = n; i >= 1; i--) {
    p = rand() % i;
    c = b[i - 1];
    b[i - 1] = b[p];
    b[p] = c;
  }
  for (i = 0; i < n; i++) {
    scanf("%d", &x[b[i]]);
  }
}
int dx4[4] = {1, -1, 0, 0};
int dy4[4] = {0, 0, 1, -1};
int dx8[8] = {-1, -1, -1, 0, 0, 1, 1, 1};
int dy8[8] = {-1, 0, 1, -1, 1, -1, 0, 1};
int search(int x, int a[], int n) {
  int st = 0, fi = n - 1, te;
  while (st <= fi) {
    te = (st + fi) / 2;
    if (a[te] < x) {
      st = te + 1;
    } else {
      fi = te - 1;
    }
  }
  return st;
}
void prarr(int arr[], int n) {
  int i;
  for (i = 0; i < n; i++) {
    if (i) {
      printf(" ");
    }
    printf("%d", arr[i]);
  }
  printf("\n");
  return;
}
void getperm(int a[], int n) {
  int i, p;
  for (i = 0; i < n; i++) {
    a[i] = i;
  }
  for (i = n - 1; i >= 1; i--) {
    p = rand() % (i + 1);
    swap(&a[p], &a[i]);
  }
}
typedef struct {
  int val;
  int node;
} sd;
int sdsortfnc(const void *a, const void *b) {
  if (((sd *)a)->val < ((sd *)b)->val) {
    return -1;
  }
  if (((sd *)a)->val > ((sd *)b)->val) {
    return 1;
  }
  return 0;
}
void coordinate_comp(int a[], int n) {
  int i, c = 0;
  sd dat[524288];
  for (i = 0; i < n; i++) {
    dat[i].val = a[i];
    dat[i].node = i;
  }
  qsort(dat, n, sizeof(dat[0]), sdsortfnc);
  a[dat[0].node] = c;
  for (i = 1; i < n; i++) {
    if (dat[i - 1].val != dat[i].val) {
      c++;
    }
    a[dat[i].node] = c;
  }
}
typedef struct {
  long long st;
  long long fi;
  long long cost;
} data;
data heap[524288];
long long hsize = 0;
void resheap() {
  long long i;
  heap[0].cost = -4154118101919364364;
  for (i = 1; i < 524288; i++) {
    heap[i].cost = 4154118101919364364;
  }
}
long long swjud(data high, data low) {
  if (high.cost > low.cost) {
    return 1;
  }
  return -1;
}
void pqpush(data x) {
  long long a, b;
  data m;
  hsize++;
  heap[hsize] = x;
  a = hsize;
  while (1) {
    b = a / 2;
    if (swjud(heap[b], heap[a]) == 1) {
      m = heap[a];
      heap[a] = heap[b];
      heap[b] = m;
    } else {
      break;
    }
    a = b;
  }
}
data pqpop() {
  long long a, b, c;
  data m, sw;
  m = heap[1];
  heap[1] = heap[hsize];
  heap[hsize].cost = 4154118101919364364;
  hsize--;
  a = 1;
  while (1) {
    b = a * 2;
    c = a * 2 + 1;
    if (swjud(heap[a], heap[b]) == -1 && swjud(heap[a], heap[c]) == -1) {
      break;
    }
    if (swjud(heap[b], heap[c]) == -1) {
      sw = heap[a];
      heap[a] = heap[b];
      heap[b] = sw;
      a = b;
    } else {
      sw = heap[a];
      heap[a] = heap[c];
      heap[c] = sw;
      a = c;
    }
  }
  return m;
}
typedef struct {
  long long par;
  long long dep;
  long long size;
} node;
node uft[524288];
void resuf() {
  long long i;
  for (i = 0; i < 524288; i++) {
    uft[i].par = i;
    uft[i].dep = 0;
    uft[i].size = 1;
  }
  return;
}
long long find(long long x) {
  if (uft[x].par == x) {
    return x;
  } else {
    uft[x].par = find(uft[x].par);
    return uft[x].par;
  }
}
void uni(long long x, long long y) {
  long long xp, yp;
  xp = find(x);
  yp = find(y);
  if (xp == yp) {
    return;
  }
  if (uft[xp].dep > uft[yp].dep) {
    uft[yp].par = xp;
    uft[xp].size += uft[yp].size;
  } else {
    uft[xp].par = yp;
    uft[yp].size += uft[xp].size;
    if (uft[xp].dep == uft[yp].dep) {
      uft[yp].dep++;
    }
  }
  return;
}
long long size(long long x) {
  x = find(x);
  return uft[x].size;
}
long long a[524288], d;
long long ccost(long long p, long long q) {
  if (p > q) {
    llswap(&p, &q);
  }
  return a[p] + a[q] + d * (q - p);
}
int main(void) {
  resuf();
  resheap();
  long long i, j, n, m, k, b, c, h, w, r = 0, l, t;
  data pd, od;
  scanf("%lld%lld", &n, &d);
  for (i = 0; i < n; i++) {
    scanf("%lld", &a[i]);
  }
  for (i = 0; i < n - 1; i++) {
    pd.st = i;
    pd.fi = i + 1;
    pd.cost = ccost(i, i + 1);
    pqpush(pd);
  }
  c = 0;
  while (c < (n - 1)) {
    od = pqpop();
    if (find(od.st) == find(od.fi)) {
      continue;
    }
    c++;
    r += od.cost;
    uni(od.st, od.fi);
    if (od.st != 0) {
      if (find(od.st - 1) != find(od.fi)) {
        pd.st = od.st - 1;
        pd.fi = od.fi;
        pd.cost = ccost(od.st - 1, od.fi);
        pqpush(pd);
      }
    }
    if (od.fi != n - 1) {
      if (find(od.st) != find(od.fi + 1)) {
        pd.st = od.st;
        pd.fi = od.fi + 1;
        pd.cost = ccost(od.st, od.fi + 1);
        pqpush(pd);
      }
    }
  }
  printf("%lld\n", r);
  return 0;
}
