Problem: While discussing a proper problem A for a Codeforces Round, Kostya created a cyclic array of positive integers a_1, a_2, …, a_n. Since the talk was long and not promising, Kostya created a new cyclic array b_1, b_2, …, b_{n} so that b_i = (a_i mod a_{i + 1}), where we take a_{n+1} = a_1. Here mod is the [modulo operation](https://en.wikipedia.org/wiki/Modulo_operation). When the talk became interesting, Kostya completely forgot how array a had looked like. Suddenly, he thought that restoring array a from array b would be an interesting problem (unfortunately, not A).

Input

The first line contains a single integer n (2 ≤ n ≤ 140582) — the length of the array a.

The second line contains n integers b_1, b_2, …, b_{n} (0 ≤ b_i ≤ 187126).

Output

If it is possible to restore some array a of length n so that b_i = a_i mod a_{(i mod n) + 1} holds for all i = 1, 2, …, n, print «YES» in the first line and the integers a_1, a_2, …, a_n in the second line. All a_i should satisfy 1 ≤ a_i ≤ 10^{18}. We can show that if an answer exists, then an answer with such constraint exists as well.

It it impossible to restore any valid a, print «NO» in one line.

You can print each letter in any case (upper or lower).

Examples

Input

4
1 3 1 0


Output

YES
1 3 5 2


Input

2
4 4


Output

NO

Note

In the first example:

  * 1 mod 3 = 1 
  * 3 mod 5 = 3 
  * 5 mod 2 = 1 
  * 2 mod 1 = 0 

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <class TH>
void _dbg(const char *sdbg, TH h) {
  cerr << sdbg << " = " << h << endl;
}
template <class TH, class... TA>
void _dbg(const char *sdbg, TH h, TA... a) {
  while (*sdbg != ',') cerr << *sdbg++;
  cerr << " = " << h << ", ";
  _dbg(sdbg + 1, a...);
}
const int maxn = (1e6) + 7;
const int maxk = 20;
const int inf = (1e9) + 7;
const long long LLinf = (1e18) + 7;
const long double eps = 1e-9;
const long long mod = 1e9 + 7;
long long odp[maxn], tab[maxn], b[maxn];
int n;
int nast(int x) {
  if (x == n) return 1;
  return x + 1;
}
long long naj(long long a, long long b, long long c) {
  long long d = (a / b) * b;
  d += c;
  while (d < a) d += b;
  return d;
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  cin >> n;
  long long ma = 0;
  for (int i = 1; i <= n; i++) {
    tab[i] = max(tab[i], 1LL);
    cin >> b[i];
    ma = max(ma, b[i]);
    tab[i] = max(tab[i], b[i]);
    tab[i + 1] = max(tab[i + 1], b[i] + 1);
  }
  if (ma == 0) {
    cout << "YES\n";
    for (int i = 1; i <= n; i++) cout << "1 ";
    return 0;
  }
  tab[1] = max(tab[1], tab[n + 1]);
  bool f = 0;
  for (int i = 1; i <= n; i++)
    if (tab[i] == b[i]) {
      f = 1;
      for (int j = i - 1; j > 0; j--) {
        tab[j] = naj(tab[j], tab[nast(j)], b[j]);
      }
      for (int j = n; j > i; j--) {
        tab[j] = naj(tab[j], tab[nast(j)], b[j]);
      }
      while ((tab[i] % tab[nast(i)] != b[i]))
        tab[nast(i)] += tab[nast(nast(i))];
      break;
    }
  for (int i = 1; f && i <= n; i++)
    if (tab[i] % tab[nast(i)] != b[i]) assert(0);
  if (f) {
    cout << "YES\n";
    for (int i = 1; i <= n; i++) cout << tab[i] << " ";
    return 0;
  }
  cout << "NO\n";
  return 0;
}