Problem: Devu wants to decorate his garden with flowers. He has purchased n boxes, where the i-th box contains fi flowers. All flowers in a single box are of the same color (hence they are indistinguishable). Also, no two boxes have flowers of the same color.

Now Devu wants to select exactly s flowers from the boxes to decorate his garden. Devu would like to know, in how many different ways can he select the flowers from each box? Since this number may be very large, he asks you to find the number modulo (109 + 7). 

Devu considers two ways different if there is at least one box from which different number of flowers are selected in these two ways.

Input

The first line of input contains two space-separated integers n and s (1 ≤ n ≤ 20, 0 ≤ s ≤ 1014).

The second line contains n space-separated integers f1, f2, ... fn (0 ≤ fi ≤ 1012).

Output

Output a single integer — the number of ways in which Devu can select the flowers modulo (109 + 7).

Examples

Input

2 3
1 3


Output

2


Input

2 4
2 2


Output

1


Input

3 5
1 3 2


Output

3

Note

Sample 1. There are two ways of selecting 3 flowers: {1, 2} and {0, 3}.

Sample 2. There is only one way of selecting 4 flowers: {2, 2}.

Sample 3. There are three ways of selecting 5 flowers: {1, 2, 2}, {0, 3, 2}, and {1, 3, 1}.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int MAX = 25;
const int MOD = 1e9 + 7;
int add(int a, int b, int c) {
  int res = a + b;
  return (res >= c ? res - c : res);
}
int mod_neg(int a, int b, int c) {
  int res;
  if (abs(a - b) < c)
    res = a - b;
  else
    res = (a - b) % c;
  return (res < 0 ? res + c : res);
}
int mul(int a, int b, int c) {
  long long res = (long long)a * b;
  return (res >= c ? res % c : res);
}
template <typename T>
T power(T e, T n, T m) {
  T x = 1, p = e;
  while (n) {
    if (n & 1) x = mul(x, p, m);
    p = mul(p, p, m);
    n >>= 1;
  }
  return x;
}
template <typename T>
T extended_euclid(T a, T b, T &x, T &y) {
  T xx = 0, yy = 1;
  y = 0;
  x = 1;
  while (b) {
    T q = a / b, t = b;
    b = a % b;
    a = t;
    t = xx;
    xx = x - q * xx;
    x = t;
    t = yy;
    yy = y - q * yy;
    y = t;
  }
  return a;
}
template <typename T>
T mod_inverse(T a, T n) {
  T x, y, z = 0;
  T d = extended_euclid(a, n, x, y);
  return (d > 1 ? -1 : mod_neg(x, z, n));
}
vector<long long> cnt;
int invp[MAX];
unordered_map<long long, int> polynomial;
void init() {
  invp[0] = invp[1] = 1;
  for (int i = 2; i < MAX; ++i) {
    invp[i] = mod_inverse(i, MOD);
  }
  for (int i = 1; i < MAX; ++i) {
    invp[i] = mul(invp[i], invp[i - 1], MOD);
  }
}
int ncr(long long n, int r) {
  if (r > n) {
    return 0;
  }
  int ans = invp[r], num = n % MOD;
  for (int i = 1; i <= r; ++i) {
    ans = mul(ans, num, MOD);
    num = mod_neg(num, 1, MOD);
  }
  return ans;
}
void gen(int pos, long long expo, int coeff) {
  if (pos == cnt.size()) {
    polynomial[expo] = add(polynomial[expo], coeff, MOD);
    return;
  }
  gen(pos + 1, expo, coeff);
  int new_coeff = (coeff * (-1)) % MOD;
  if (new_coeff < 0) new_coeff += MOD;
  gen(pos + 1, expo + cnt[pos] + 1, new_coeff);
}
int main() {
  init();
  int n;
  long long s;
  scanf("%d %lld", &n, &s);
  cnt.resize(n);
  for (int i = 0; i < n; ++i) {
    scanf("%lld", &cnt[i]);
  }
  int ans = 0;
  gen(0, 0, 1);
  for (auto it : polynomial) {
    long long expo = it.first;
    int coeff = it.second;
    long long rem = s - expo;
    if (rem >= 0) {
      ans = add(ans, mul(coeff, ncr(n + rem - 1, n - 1), MOD), MOD);
    }
  }
  printf("%d\n", ans);
  return 0;
}