Problem: There are n types of coins in Byteland. Conveniently, the denomination of the coin type k divides the denomination of the coin type k + 1, the denomination of the coin type 1 equals 1 tugrick. The ratio of the denominations of coin types k + 1 and k equals ak. It is known that for each x there are at most 20 coin types of denomination x.

Byteasar has bk coins of type k with him, and he needs to pay exactly m tugricks. It is known that Byteasar never has more than 3·105 coins with him. Byteasar want to know how many ways there are to pay exactly m tugricks. Two ways are different if there is an integer k such that the amount of coins of type k differs in these two ways. As all Byteland citizens, Byteasar wants to know the number of ways modulo 109 + 7.

Input

The first line contains single integer n (1 ≤ n ≤ 3·105) — the number of coin types.

The second line contains n - 1 integers a1, a2, ..., an - 1 (1 ≤ ak ≤ 109) — the ratios between the coin types denominations. It is guaranteed that for each x there are at most 20 coin types of denomination x.

The third line contains n non-negative integers b1, b2, ..., bn — the number of coins of each type Byteasar has. It is guaranteed that the sum of these integers doesn't exceed 3·105.

The fourth line contains single integer m (0 ≤ m < 1010000) — the amount in tugricks Byteasar needs to pay.

Output

Print single integer — the number of ways to pay exactly m tugricks modulo 109 + 7.

Examples

Input

1

4
2


Output

1


Input

2
1
4 4
2


Output

3


Input

3
3 3
10 10 10
17


Output

6

Note

In the first example Byteasar has 4 coins of denomination 1, and he has to pay 2 tugricks. There is only one way.

In the second example Byteasar has 4 coins of each of two different types of denomination 1, he has to pay 2 tugricks. There are 3 ways: pay one coin of the first type and one coin of the other, pay two coins of the first type, and pay two coins of the second type.

In the third example the denominations are equal to 1, 3, 9.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int N = 600005, P = 1000000007, B = 1000000000;
using ll = long long;
int n, b[N], c[N];
ll a[N], f[N], g[N], sum[N];
string s;
struct BigInt {
  int pos, val[1205];
  BigInt() {
    memset(val, 0, sizeof(val));
    val[0] = 1;
  }
  void get(string s) {
    int sz = s.size();
    val[0] = 0;
    for (int i = 0; i < sz; i++) {
      for (int j = 1200; j >= 0; --j) {
        ll tmp = val[j] * 10ll;
        val[j] = tmp % B;
        val[j + 1] += tmp / B;
      }
      val[0] += s[i] - '0';
      for (int j = 0; j <= 1200; j++)
        if (val[j] >= B) val[j] -= B, ++val[j + 1];
    }
    pos = 1200;
  }
  bool empty() {
    for (int i = 0; i <= 1200; i++)
      if (val[i]) return 0;
    return 1;
  }
  inline int &operator[](const size_t &s) { return val[s]; }
  inline const int &operator[](const size_t &s) const { return val[s]; }
} m, d[10005];
BigInt operator+(const BigInt &a, const BigInt &b) {
  BigInt c;
  c[0] = 0;
  for (int i = 0; i <= 1200; i++) {
    c[i] += a[i] + b[i];
    if (c[i] >= B) c[i] -= B, ++c[i + 1];
  }
  return c;
}
BigInt operator-(const BigInt &a, const BigInt &b) {
  BigInt c;
  c[0] = 0;
  for (int i = 0; i <= 1200; i++) {
    c[i] += a[i] - b[i];
    if (c[i] < 0) c[i] += B, --c[i + 1];
  }
  return c;
}
BigInt operator*(const BigInt &a, ll b) {
  BigInt c;
  c[0] = 0;
  for (int i = 0; i <= 1200; i++) {
    ll tmp = a[i] * b;
    c[i] += tmp % B;
    c[i + 1] += tmp / B;
    if (c[i] >= B) c[i] -= B, ++c[i + 1];
  }
  return c;
}
BigInt operator/(const BigInt &a, ll b) {
  BigInt c;
  c[0] = 0;
  c.pos = a.pos;
  ll tmp = 0;
  for (int i = a.pos; i >= 0; --i) {
    c[i] = (tmp + a[i]) / b;
    tmp = (tmp + a[i]) % b * B;
    if (!c[i] && i == c.pos) --c.pos;
  }
  return c;
}
int main() {
  ios::sync_with_stdio(false);
  cin >> n;
  a[1] = 1;
  for (int i = 2; i <= n; i++) cin >> a[i];
  for (int i = 1; i <= n; i++) cin >> b[i];
  cin >> s;
  m.get(s);
  BigInt cur = m;
  for (int i = 2; i <= n; i++) {
    if (cur.empty()) {
      n = i - 1;
      break;
    }
    if (a[i] != 1) {
      BigInt tmp = cur / a[i];
      c[i] = (cur - tmp * a[i])[0];
      cur = tmp;
    }
  }
  int lim = 0;
  f[0] = 1;
  for (int i = 1; i <= n; i++) {
    memcpy(g, f, sizeof(ll) * (lim + 1));
    memset(f, 0, sizeof(ll) * (lim + 1));
    for (int j = c[i]; j <= lim; j += a[i]) f[j / a[i]] = g[j];
    lim = lim / a[i] + b[i];
    for (int j = 1; j <= lim + 1; j++) sum[j] = (sum[j - 1] + f[j - 1]) % P;
    for (int j = 0; j <= lim; j++)
      f[j] = (sum[j + 1] - sum[max(0, j - b[i])] + P) % P;
  }
  for (int i = 1; i <= 1200; i++)
    if (cur[i]) {
      puts("0");
      return 0;
    }
  if (cur[0] > lim) {
    puts("0");
    return 0;
  }
  cout << f[cur[0]] << endl;
  return 0;
}