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>
inline int64_t ceil_div(int64_t a, int64_t b) {
  return a / b + (a % b ? 1 : 0);
}
int main(int argc, char* argv[]) {
  int n;
  scanf("%d", &n);
  std::vector<int64_t> a(n), b(n);
  for (int i = 0; i < n; ++i) scanf("%lld", &b[i]);
  auto map_index = [=](int ind) -> int {
    while (ind >= n) ind -= n;
    while (ind < 0) ind += n;
    return ind;
  };
  if (std::all_of(b.begin() + 1, b.end(), [&](int x) { return x == b[0]; })) {
    if (b[0] == 0) {
      puts("YES");
      for (int i = 1; i <= n; ++i) printf("1 ");
      putchar('\n');
    } else
      puts("NO");
    return 0;
  }
  int64_t maxvalue = *std::max_element(b.begin(), b.end());
  int pos = 0;
  while (pos < n) {
    if (b[pos] == maxvalue && b[map_index(pos - 1)] < b[pos])
      break;
    else
      ++pos;
  }
  a[pos] = b[pos];
  for (int i = map_index(pos - 1); i != pos; i = map_index(i - 1)) {
    int64_t bi = b[i];
    int64_t bim1 = b[map_index(i - 1)];
    int64_t aip1 = a[map_index(i + 1)];
    int64_t& ai = a[i];
    int64_t ki = (bi > bim1) ? 0 : ceil_div(bim1 - bi + 1, aip1);
    ai = bi + ki * aip1;
  }
  puts("YES");
  for (int64_t value : a) printf("%lld ", value);
  putchar('\n');
  return 0;
}