Problem: Due to the recent popularity of the Deep learning new countries are starting to look like Neural Networks. That is, the countries are being built deep with many layers, each layer possibly having many cities. They also have one entry, and one exit point.

There are exactly L layers, each having N cities. Let us look at the two adjacent layers L1 and L2. Each city from the layer L1 is connected to each city from the layer L2 with the traveling cost cij for <image>, and each pair of adjacent layers has the same cost in between their cities as any other pair (they just stacked the same layers, as usual). Also, the traveling costs to each city from the layer L2 are same for all cities in the L1, that is cij is the same for <image>, and fixed j.

Doctor G. needs to speed up his computations for this country so he asks you to find the number of paths he can take from entry to exit point such that his traveling cost is divisible by given number M.

Input

The first line of input contains N (1 ≤ N ≤ 106), L (2 ≤ L ≤ 105) and M (2 ≤ M ≤ 100), the number of cities in each layer, the number of layers and the number that travelling cost should be divisible by, respectively.

Second, third and fourth line contain N integers each denoting costs 0 ≤ cost ≤ M from entry point to the first layer, costs between adjacent layers as described above, and costs from the last layer to the exit point.

Output

Output a single integer, the number of paths Doctor G. can take which have total cost divisible by M, modulo 109 + 7.

Example

Input

2 3 13
4 6
2 1
3 4


Output

2

Note

<image>

This is a country with 3 layers, each layer having 2 cities. Paths <image>, and <image> are the only paths having total cost divisible by 13. Notice that input edges for layer cities have the same cost, and that they are same for all layers.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const long double epsylon = 1e-9;
const std::string PROGRAM_NAME = "ivo";
long long mod = 1000000000 + 7;
struct matrix {
  matrix(const vector<vector<long long> >& a_) {
    a = a_;
    n = (int)a.size();
    if (a_[0].size() != a_.size()) {
      wrong_input("The matrix is not square!\n");
    }
    for (int i = 0; i < n; ++i) {
      for (int j = 0; j < n; ++j) {
        a[i][j] %= mod;
      }
    }
  }
  matrix(const vector<vector<int> >& a_) {
    n = (int)a_.size();
    a.resize(n);
    if (a_[0].size() != a_.size()) {
      wrong_input("The matrix is not square!\n");
    }
    for (int i = 0; i < n; ++i) {
      a[i].resize(n);
      for (int j = 0; j < n; ++j) {
        a[i][j] = a_[i][j] % mod;
      }
    }
  }
  matrix& operator+=(const matrix& o) {
    if (a.size() != o.a.size()) {
      wrong_input("The sizes of the matrices is not the same!\n");
    }
    for (int i = 0; i < n; ++i) {
      for (int j = 0; j < n; ++j) {
        a[i][j] = (a[i][j] + o.a[i][j]) % mod;
      }
    }
    return *this;
  }
  matrix operator+(const matrix& o) {
    matrix res = *this;
    res += o;
    return res;
  }
  matrix& operator-=(const matrix& o) {
    if (a.size() != o.a.size()) {
      wrong_input("The sizes of the matrices is not the same!\n");
    }
    for (int i = 0; i < n; ++i) {
      for (int j = 0; j < n; ++j) {
        a[i][j] = (a[i][j] - o.a[i][j] + mod) % mod;
      }
    }
    return *this;
  }
  matrix operator-(const matrix& o) {
    matrix res = *this;
    res -= o;
    return res;
  }
  matrix& operator*=(const matrix& o) {
    if (a.size() != o.a.size()) {
      wrong_input("The sizes of the matrices is not the same!\n");
    }
    vector<vector<long long> > res(n);
    for (int i = 0; i < n; ++i) {
      res[i].resize(n, 0);
    }
    for (int i = 0; i < n; ++i) {
      for (int j = 0; j < n; ++j) {
        for (int k = 0; k < n; ++k) {
          res[i][j] = (res[i][j] + (a[i][k] * o.a[k][j]) % mod) % mod;
        }
      }
    }
    a = res;
    return *this;
  }
  matrix operator*(const matrix& o) {
    matrix res = *this;
    res *= o;
    return res;
  }
  long long stepen(long long x, long long y) {
    long long res = 1;
    long long p = x;
    while (y) {
      if (y % 2) {
        res = (res * p) % mod;
      }
      p = (p * p) % mod;
      y /= 2;
    }
    return res;
  }
  matrix stepen(long long power) {
    matrix res = unary(n);
    matrix temp = *this;
    while (power) {
      if (power % 2LL) {
        res *= temp;
      }
      temp *= temp;
      power /= 2LL;
    }
    return res;
  }
  void wrong_input(const string& str) { cerr << str; }
  static matrix unary(int n) {
    vector<vector<long long> > res(n);
    for (int i = 0; i < n; ++i) {
      res[i].resize(n, 0);
      res[i][i] = 1;
    }
    return matrix(res);
  }

 public:
  int n;
  vector<vector<long long> > a;
};
int main() {
  int n, l, m;
  cin >> n >> l >> m;
  vector<long long> cur(m, 0), mid(m, 0), last(n, 0);
  for (int i = 0; i < n; ++i) {
    int temp;
    scanf("%d", &temp);
    cur[temp % m] += 1;
  }
  for (int i = 0; i < n; ++i) {
    int temp;
    scanf("%d", &temp);
    mid[temp % m] += 1;
    last[i] = temp;
  }
  for (int i = 0; i < n; ++i) {
    int temp;
    scanf("%d", &temp);
    last[i] = (last[i] + temp) % m;
  }
  vector<vector<int> > ma(m, vector<int>(m, 0));
  for (int i = 0; i < m; ++i) {
    for (int j = 0; j < m; ++j) {
      ma[i][j] = mid[(j - i + m) % m];
    }
  }
  matrix a(ma);
  matrix s = a.stepen(l - 2);
  vector<int> nxt(m, 0);
  for (int i = 0; i < m; ++i) {
    for (int j = 0; j < m; ++j) {
      nxt[i] = (int)((nxt[i] + (cur[j] * s.a[j][i]) % mod) % mod);
    }
  }
  for (int i = 0; i < m; ++i) {
    cur[i] = nxt[i];
  }
  long long answer = 0;
  for (int i = 0; i < n; ++i) {
    int other = (m - last[i]) % m;
    answer = (answer + cur[other]) % mod;
  }
  cout << answer << endl;
  return 0;
}