Problem p02709 (C++) — buggy submission, verdict: Runtime Error

/*
 * じょえチャンネル
 * 高評価・チャンネル登録よろしくおねがいします！
 * https://www.youtube.com/channel/UCRXsI3FL_kvaVL9zoolBfbQ
 */

#include <bits/stdc++.h>

// #pragma GCC target("avx2")
// #pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")

#define f(i, n) for (int i = 0; i < (n); i++)

// here!!!
//  define int long long !!!!!

#define int long long

// define int long long !!!!!

#define mod (int)((1e9) + 7)
// constexpr int mod = 998244353ll;

#ifdef int
#define inf (int)(3e18)
#else
#define inf (int)(5e8)
#endif

#define intt long long
#define itn long long
#define P pair<long long, long long>

#define rep(i, n) for (int i = 0; i < n; i++)
#define REP(i, n) for (int i = 1; i <= n; i++)

#define ALL(v) v.begin(), v.end()

using namespace std;

// Library
// モッドパウ
inline int modpow(int x, int y, int m = mod) {
  int res = 1;
  while (y) {
    if (y & 1) {
      res *= x;
      res %= m;
    }
    x = x * x % m;
    y /= 2;
  }
  return res;
}

int mypow(int x, int y) {
  int res = 1;
  while (y) {
    if (y % 2) {
      res *= x;
    }
    x = x * x;
    y /= 2;
  }
  return res;
}
// is the number (x) a prime number?
bool prime(int x) {
  for (int i = 2; i * i <= x; i++) {
    if (!(x % i)) {
      return false;
    }
  }
  return true;
}

// saidai-kouyakusuu
inline int gcd(int x, int y) {
  if (!y) {
    return x;
  }
  return gcd(y, x % y);
}

// number of keta
int keta(int x) {
  int ans = 0;
  while (x) {
    x /= 10;
    ans++;
  }
  return ans;
}

// sum of keta
int ketasum(int x) {
  int ans = 0;
  while (x) {
    ans += x % 10;
    x /= 10;
  }
  return ans;
}

inline int lcm(int x, int y) {
  int ans = x / gcd(x, y) * y;
  return ans;
}
int twobeki(int x) {
  int ans = 0;
  while (1) {
    if (!(x & 1)) {
      ans++;
      x /= 2;
    } else {
      break;
    }
  }
  return ans;
}

template <class T, class U> inline bool chmax(T &lhs, const U &rhs) {
  if (lhs < rhs) {
    lhs = rhs;
    return 1;
  }
  return 0;
}
template <class T, class U> inline bool chmin(T &lhs, const U &rhs) {
  if (lhs > rhs) {
    lhs = rhs;
    return 1;
  }
  return 0;
}
void Yes() { cout << "Yes" << endl; }
void No() { cout << "No" << endl; }
void YES() { cout << "YES" << endl; }
void NO() { cout << "NO" << endl; }

#define fin(i) scanf("%lld", &i)
#define fout(i) printf("%lld", i)
#define fendl printf("\\n")

int kai(int x, int y) {
  int res = 1;
  for (int i = x - y + 1; i <= x; i++) {
    res *= i;
    res %= mod;
  }
  return res;
}

int comb(int x, int y) {
  if (y > x)
    return 0;
  //    cout<<kai(x, y)<<' '<<modpow(kai(y, y), mod - 2)<<endl;
  return kai(x, y) * modpow(kai(y, y), mod - 2) % mod;
}
// Library-End

#define vecin(v)                                                               \\
  for (int i = 0; i < v.size(); i++)                                           \\
    scanf("%lld", &v[i]);
#define vecout(v)                                                              \\
  {                                                                            \\
    for (int i = 0; i < (int)v.size(); i++)                                    \\
      printf("%lld ", v[i]);                                                   \\
    printf("\\n");                                                              \\
  }

template <typename T, typename U> class SegTree {
protected:
  unsigned int n = 1, rank = 0;
  std::vector<T> node;
  T nodee;
  virtual T nodef(const T &, const T &) const = 0;

public:
  SegTree(unsigned int m, T init, T nodee) : nodee(nodee) {
    while (n < m) {
      n *= 2;
      rank++;
    }
    node.resize(2 * n);
    for (unsigned int i = n; i < 2 * n; i++)
      node[i] = init;
  }
  SegTree(const std::vector<T> &initvec, T nodee) : nodee(nodee) {
    unsigned int m = initvec.size();
    while (n < m) {
      n *= 2;
      rank++;
    }
    node.resize(2 * n);
    for (unsigned int i = n; i < 2 * n; i++) {
      if (i - n < m)
        node[i] = initvec[i - n];
    }
  }
  virtual void update(int i, U x) {
    i += n;
    node[i] = x;
    while (i != 1) {
      i >>= 1;
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
    }
  }
  virtual T query(int l, int r) {
    l += n;
    r += n;
    int ls = nodee, rs = nodee;
    while (l < r) {
      if (l & 1) {
        ls = nodef(ls, node[l]);
        l++;
      }
      if (r & 1) {
        r--;
        rs = nodef(node[r], rs);
      }
      l >>= 1;
      r >>= 1;
    }
    return nodef(ls, rs);
  }
  virtual T operator[](const int &x) { return node[n + x]; }
  void fill(T x) { std::fill(all(node), x); }
  void print() {
    rep(i, n) std::cout << operator[](i) << " ";
    std::cout << std::endl;
  }
};
class RSQ : public SegTree<int, int> {
  int nodef(const int &lhs, const int &rhs) { return lhs + rhs; }

public:
  RSQ(int size, const int &def = 0) : SegTree<int, int>(size, def, 0) {}
  RSQ(const std::vector<int> &initvec) : SegTree<int, int>(initvec, 0) {}
};
class RMiQ : public SegTree<int, int> {
  int nodef(const int &lhs, const int &rhs) { return std::min(lhs, rhs); }

public:
  RMiQ(int size, const int &def = 0) : SegTree<int, int>(size, def, 0) {}
  RMiQ(const std::vector<int> &initvec) : SegTree<int, int>(initvec, 0) {}
};
class RMaQ : public SegTree<int, int> {
  int nodef(const int &lhs, const int &rhs) { return std::max(lhs, rhs); }

public:
  RMaQ(int size, const int &def = 0) : SegTree<int, int>(size, def, 0) {}
  RMaQ(const std::vector<int> &initvec) : SegTree<int, int>(initvec, 0) {}
};
template <typename T, typename U> class IntervalSegTree : public SegTree<T, U> {
protected:
  using SegTree<T, U>::n;
  using SegTree<T, U>::rank;
  using SegTree<T, U>::node;
  using SegTree<T, U>::nodef;
  using SegTree<T, U>::nodee;
  std::vector<U> lazy;
  std::vector<bool> lazyflag;
  std::vector<int> width;
  virtual void lazyf(U &, const U &) = 0;
  virtual void updf(T &, const U &, const unsigned int &) = 0;
  void eval(int k) {
    for (int i = rank; i > 0; i--) {
      int nk = k >> i;
      if (lazyflag[nk]) {
        updf(node[2 * nk], lazy[nk], width[2 * nk]);
        updf(node[2 * nk + 1], lazy[nk], width[2 * nk + 1]);
        if (lazyflag[2 * nk])
          lazyf(lazy[2 * nk], lazy[nk]);
        else
          lazy[2 * nk] = lazy[nk];
        if (lazyflag[2 * nk + 1])
          lazyf(lazy[2 * nk + 1], lazy[nk]);
        else
          lazy[2 * nk + 1] = lazy[nk];
        lazyflag[2 * nk] = lazyflag[2 * nk + 1] = true;
        lazyflag[nk] = false;
      }
    }
  }

public:
  IntervalSegTree(unsigned int m, T init, T nodee)
      : SegTree<T, U>(m, init, nodee) {
    lazy.resize(2 * n);
    lazyflag.resize(2 * n);
    width.resize(2 * n);
    width[1] = n;
    for (unsigned int i = 2; i < 2 * n; i++) {
      width[i] = width[i >> 1] >> 1;
    }
  }
  IntervalSegTree(T nodee, const std::vector<T> &initvec)
      : SegTree<T, U>(initvec, nodee) {
    lazy.resize(2 * n);
    lazyflag.resize(2 * n);
    width.resize(2 * n);
    width[1] = n;
    for (unsigned int i = 2; i < 2 * n; i++) {
      width[i] = width[i >> 1] >> 1;
    }
  }
  void update(int i, U x) {
    i += n;
    eval(i);
    updf(node[i], x, width[i]);
    if (lazyflag[i])
      lazyf(lazy[i], x);
    else {
      lazyflag[i] = true;
      lazy[i] = x;
    }
    while (i /= 2)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  void update(int l, int r, U x) {
    l += n;
    r += n;
    int nl = l, nr = r;
    while (!(nl & 1))
      nl >>= 1;
    while (!(nr & 1))
      nr >>= 1;
    nr--;
    eval(nl);
    eval(nr);
    while (l < r) {
      if (l & 1) {
        updf(node[l], x, width[l]);
        if (lazyflag[l])
          lazyf(lazy[l], x);
        else {
          lazyflag[l] = true;
          lazy[l] = x;
        }
        l++;
      }
      if (r & 1) {
        r--;
        updf(node[r], x, width[r]);
        if (lazyflag[r])
          lazyf(lazy[r], x);
        else {
          lazyflag[r] = true;
          lazy[r] = x;
        }
      }
      l >>= 1;
      r >>= 1;
    }
    while (nl >>= 1)
      node[nl] = nodef(node[2 * nl], node[2 * nl + 1]);
    while (nr >>= 1)
      node[nr] = nodef(node[2 * nr], node[2 * nr + 1]);
  }
  T query(int l, int r) {
    l += n;
    r += n;
    eval(l);
    eval(r - 1);
    int ls = nodee, rs = nodee;
    while (l < r) {
      if (l & 1) {
        ls = nodef(ls, node[l]);
        l++;
      }
      if (r & 1) {
        r--;
        rs = nodef(node[r], rs);
      }
      l >>= 1;
      r >>= 1;
    }
    return nodef(ls, rs);
  }
  T operator[](const int &x) {
    eval(n + x);
    return node[n + x];
  }
  T queryForAll() { return node[1]; }
};
class RAQRSQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return a + b; }
  void lazyf(int &a, const int &b) { a += b; }
  void updf(int &a, const int &b, const unsigned int &width) { a += width * b; }

public:
  RAQRSQ(int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, 0) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RAQRSQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>((int)0, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};
class RAQRMiQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return std::min(a, b); }
  void lazyf(int &a, const int &b) { a += b; }
  void updf(int &a, const int &b, const unsigned int &width) { a += b; }

public:
  RAQRMiQ(int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, inf) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RAQRMiQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>(inf, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};
class RAQRMaQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return std::max(a, b); }
  void lazyf(int &a, const int &b) { a += b; }
  void updf(int &a, const int &b, const unsigned int &width) { a += b; }

public:
  RAQRMaQ(unsigned int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, -inf) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RAQRMaQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>(-inf, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};
class RUQRSQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return a + b; }
  void lazyf(int &a, const int &b) { a = b; }
  void updf(int &a, const int &b, const unsigned int &width) { a = width * b; }

public:
  RUQRSQ(int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, 0) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RUQRSQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>((int)0, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};
class RUQRMiQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return std::min(a, b); }
  void lazyf(int &a, const int &b) { a = b; }
  void updf(int &a, const int &b, const unsigned int &width) { a = b; }

public:
  RUQRMiQ(int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, inf) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RUQRMiQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>(inf, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};
class RUQRMaQ : public IntervalSegTree<int, int> {
  int nodef(const int &a, const int &b) const { return std::max(a, b); }
  void lazyf(int &a, const int &b) { a = b; }
  void updf(int &a, const int &b, const unsigned int &width) { a = b; }

public:
  RUQRMaQ(int size, const int &def = 0)
      : IntervalSegTree<int, int>(size, def, -inf) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
  RUQRMaQ(const std::vector<int> &initvec)
      : IntervalSegTree<int, int>(-inf, initvec) {
    for (int i = n - 1; i > 0; i--)
      node[i] = nodef(node[2 * i], node[2 * i + 1]);
  }
};

#define R_MIN ([](long long a, long long b) { return min(a, b); })
#define R_MAX ([](long long a, long long b) { return max(a, b); })
#define R_SUM ([](long long a, long long b) { return a + b; })

#define NORMAL_UPDATE ([](long long a, long long b) { return b; })
#define ADD_UPDATE ([](long long a, long long b) { return a + b; })
#define MINUS_UPDATE ([](long long a, long long b) { return a - b; }

class Union_Find {
  vector<int> par;
  vector<int> rankmy;
  vector<int> ookisa;

public:
  Union_Find(int size) {
    par = vector<int>(size);
    rankmy = vector<int>(size);
    ookisa = vector<int>(size);
    for (int i = 0; i < size; i++) {
      par[i] = i;
      ookisa[i] = 1;
    }
  }

  int find(int x) {
    if (par[x] == x) {
      return x;
    }
    return par[x] = find(par[x]);
  }

  void unite(int x, int y) {
    x = find(x);
    y = find(y);
    if (x == y) {
      return;
    }
    if (rankmy[x] < rankmy[y]) {
      par[x] = y;
      ookisa[y] += ookisa[x];
      ookisa[x] = 0;
    } else {
      par[y] = x;
      ookisa[x] += ookisa[y];
      ookisa[y] = 0;
      if (rankmy[x] == rankmy[y]) {
        rankmy[x]++;
      }
    }
  }
  int size(int i) {
    i = find(i);
    return ookisa[i];
  }
  bool same(int x, int y) { return find(x) == find(y); }
};

// Union-Find-End

int perm[1000005];
void init_perm() {
  perm[0] = 1;
  REP(i, 1000000) { perm[i] = perm[i - 1] * i % mod; }
}

int nCk(int x, int y) {
  return perm[x] * modpow(perm[x - y], mod - 2) % mod *
         modpow(perm[y], mod - 2) % mod;
}

double kyori(pair<int, int> f, pair<int, int> s) {
  double ans = 0;
  double t = fabs(f.first - s.first);
  double y = fabs(f.second - s.second);
  ans = sqrt(t * t + y * y);
  return ans;
}

inline string stringmax(string &x, string &y) {
  if (x.size() > y.size()) {
    return x;
  }
  if (x.size() < y.size()) {
    return y;
  }
  rep(i, x.size()) {
    if (x[i] > y[i]) {
      return x;
    }
    if (x[i] < y[i]) {
      return y;
    }
  }
  return x;
}

#define endl "\\n" // interactive の時に注意！！！
#define printd cout << fixed << setprecision(10)

itn n, a[2004], ans, dp[2004][2004];
vector<P> v;
signed main() {
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  cin >> n;
  rep(i, n) {
    cin >> a[i];
    v.emplace_back(a[i], i);
  }
  sort(ALL((v)), greater<P>());
  rep(i, n + 1) {
    rep(j, n - i) {
      chmax(dp[i + 1][j],
            dp[i][j] + v[i + j].first * ((n - 1 - (i)) - v[i + j].second));
      chmax(dp[i][j + 1], dp[i][j] + v[i + j].first * (v[i + j].second - (j)));
      chmax(ans, dp[i][j]);
      chmax(ans, dp[i + 1][j]);
      chmax(ans, dp[i][j + 1]);
    }
  }
  cout << ans << endl;
}


Error: 0

Fixed submission:

/*
 * じょえチャンネル
 * 高評価・チャンネル登録よろしくおねがいします！
 * https://www.youtube.com/channel/UCRXsI3FL_kvaVL9zoolBfbQ
 */

#include <bits/stdc++.h>

// #pragma GCC target("avx2")
// #pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")

#define f(i, n) for (int i = 0; i < (n); i++)

#define inf (int)(3e18)

// here!!!
//  define int long long !!!!!

#define int long long

// define int long long !!!!!

#define mod (int)((1e9) + 7)

#define intt long long
#define itn long long
#define P pair<int, int>

#define rep(i, n) for (int i = 0; i < n; i++)
#define REP(i, n) for (int i = 1; i <= n; i++)

#define ALL(v) v.begin(), v.end()

using namespace std;

// Library
// モッドパウ
inline int modpow(int x, int y, int m = mod) {
  int res = 1;
  while (y) {
    if (y % 2) {
      res *= x;
      res %= m;
    }
    x = x * x % m;
    y /= 2;
  }
  return res;
}

int mypow(int x, int y) {
  int res = 1;
  while (y) {
    if (y % 2) {
      res *= x;
    }
    x = x * x;
    y /= 2;
  }
  return res;
}
// is the number (x) a prime number?
bool prime(int x) {
  for (int i = 2; i * i <= x; i++) {
    if (!(x % i)) {
      return false;
    }
  }
  return true;
}

// saidai-kouyakusuu
inline int gcd(int x, int y) {
  if (!y) {
    return x;
  }
  return gcd(y, x % y);
}

// number of keta
int keta(int x) {
  int ans = 0;
  while (x) {
    x /= 10;
    ans++;
  }
  return ans;
}

// sum of keta
int ketasum(int x) {
  int ans = 0;
  while (x) {
    ans += x % 10;
    x /= 10;
  }
  return ans;
}

inline int lcm(int x, int y) {
  int ans = x / gcd(x, y) * y;
  return ans;
}
int twobeki(int x) {
  int ans = 0;
  while (1) {
    if (!(x & 1)) {
      ans++;
      x /= 2;
    } else {
      break;
    }
  }
  return ans;
}

template <class T, class U> inline bool chmax(T &lhs, const U &rhs) {
  if (lhs < rhs) {
    lhs = rhs;
    return 1;
  }
  return 0;
}
template <class T, class U> inline bool chmin(T &lhs, const U &rhs) {
  if (lhs > rhs) {
    lhs = rhs;
    return 1;
  }
  return 0;
}
void Yes() { cout << "Yes" << endl; }
void No() { cout << "No" << endl; }
void YES() { cout << "YES" << endl; }
void NO() { cout << "NO" << endl; }

#define fin(i) scanf("%lld", &i)
#define fout(i) printf("%lld", i)
#define fendl printf("\\n")

int kai(int x, int y) {
  int res = 1;
  for (int i = x - y + 1; i <= x; i++) {
    res *= i;
    res %= mod;
  }
  return res;
}

int comb(int x, int y) {
  if (y > x)
    return 0;
  //    cout<<kai(x, y)<<' '<<modpow(kai(y, y), mod - 2)<<endl;
  return kai(x, y) * modpow(kai(y, y), mod - 2) % mod;
}
// Library-End

#define vecin(v)                                                               \\
  for (int i = 0; i < v.size(); i++)                                           \\
    scanf("%lld", &v[i]);
#define vecout(v)                                                              \\
  {                                                                            \\
    for (int i = 0; i < v.size(); i++)                                         \\
      printf("%lld\\n", v[i]);                                                  \\
  }

// SegTree
template <class T> class SegTree {
  int n;                       // 葉の数
  vector<T> data;              // データを格納するvector
  T def;                       // 初期値かつ単位元
  function<T(T, T)> operation; // 区間クエリで使う処理
  function<T(T, T)> update;    // 点更新で使う処理

  // 区間[a,b)の総和。ノードk=[l,r)に着目している。
  T _query(int a, int b, int k, int l, int r) {
    if (r <= a || b <= l)
      return def; // 交差しない
    if (a <= l && r <= b)
      return data[k]; // a,l,r,bの順で完全に含まれる
    else {
      T c1 = _query(a, b, 2 * k + 1, l, (l + r) / 2); // 左の子
      T c2 = _query(a, b, 2 * k + 2, (l + r) / 2, r); // 右の子
      return operation(c1, c2);
    }
  }

public:
  // _n:必要サイズ, _def:初期値かつ単位元, _operation:クエリ関数,
  // _update:更新関数
  SegTree(size_t _n, T _def, function<T(T, T)> _operation,
          function<T(T, T)> _update)
      : def(_def), operation(_operation), update(_update) {
    n = 1;
    while (n < _n) {
      n *= 2;
    }
    data = vector<T>(2 * n - 1, def);
  }

  // 場所i(0-indexed)の値をxで更新
  void change(int i, T x) {
    i += n - 1;
    data[i] = update(data[i], x);
    while (i > 0) {
      i = (i - 1) / 2;
      data[i] = operation(data[i * 2 + 1], data[i * 2 + 2]);
    }
  }

  // [a, b)の区間クエリを実行
  T query(int a, int b) { return _query(a, b, 0, 0, n); }

  // 添字でアクセス
  T operator[](int i) { return data[i + n - 1]; }
};

#define R_MIN ([](long long a, long long b) { return min(a, b); })
#define R_MAX ([](long long a, long long b) { return max(a, b); })
#define R_SUM ([](long long a, long long b) { return a + b; })

#define NORMAL_UPDATE ([](long long a, long long b) { return b; })
#define ADD_UPDATE ([](long long a, long long b) { return a + b; })
#define MINUS_UPDATE ([](long long a, long long b) { return a - b; }

class Union_Find {
  vector<int> par;
  vector<int> rankmy;
  vector<int> ookisa;

public:
  Union_Find(int size) {
    par = vector<int>(size);
    rankmy = vector<int>(size);
    ookisa = vector<int>(size);
    for (int i = 0; i < size; i++) {
      par[i] = i;
      ookisa[i] = 1;
    }
  }

  int find(int x) {
    if (par[x] == x) {
      return x;
    }
    return par[x] = find(par[x]);
  }

  void unite(int x, int y) {
    x = find(x);
    y = find(y);
    if (x == y) {
      return;
    }
    if (rankmy[x] < rankmy[y]) {
      par[x] = y;
      ookisa[y] += ookisa[x];
      ookisa[x] = 0;
    } else {
      par[y] = x;
      ookisa[x] += ookisa[y];
      ookisa[y] = 0;
      if (rankmy[x] == rankmy[y]) {
        rankmy[x]++;
      }
    }
  }
  int size(int i) {
    i = find(i);
    return ookisa[i];
  }
  bool same(int x, int y) { return find(x) == find(y); }
};

// Union-Find-End
#define endl "\\n" // interactive の時に注意！！！

int n, a[14], mae, ushiro, dp[2004][2004], ans, cnt;
vector<pair<int, P>> v;
bool yoko[2004], tate[2004];
signed main() {
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  cin >> n;
  REP(i, n) { cin >> a[i]; }
  REP(i, n) {
    REP(j, n) {
      dp[i][j] = abs(i - j) * a[j];
      v.push_back({dp[i][j], {i, j}});
    }
  }
  sort(ALL(v));

  cnt = n;
  while (cnt) {
    if (yoko[v.back().second.first] || tate[v.back().second.second]) {
      v.pop_back();
    } else {
      ans += v.back().first;
      yoko[v.back().second.first] = 1;
      tate[v.back().second.second] = 1;
      v.pop_back();
      cnt--;
    }
  }
  cout << ans << endl;
}