#pragma GCC optimize("Ofast") #include using namespace std; #define rep(i, n) for(int i = 0; i < (int)(n); ++i) #define all(x) x.begin(),x.end() #define ln '\n' const long double PI = acos(-1.0L); const long long MOD = 1000000007LL; // const long long MOD = 998244353LL; typedef long long ll; typedef pair pii; typedef pair pll; template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true;} return false; } template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true;} return false; } //////////////////////////////////////////////////////////////////////////////////////////////////////// ll gcd(ll x, ll y) { return x ? gcd(y%x,x) : y; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); ll N,M,K; cin >> N >> M >> K; char c; cin >> c; vector B(M),A(N); rep(i,M) cin >> B[i]; rep(i,N) cin >> A[i]; if (c == '+') { ll ans = 0; map dic; rep(i,N) dic[A[i]%K]++; rep(i,M) ans += dic[(K-B[i]%K)%K]; cout << ans << ln; } else { /* vector C; for (int i = 1; i*i <= K; ++i) { if (K%i==0) { C.emplace_back(i); C.emplace_back(K/i); } } ll ans = 0; sort(all(C)); int L = C.size(); vector cntA(L),cntB(L); rep(i,N) { rep(j,L) { if (A[i]%C[j]==0) { cntA[j]++; } } } rep(i,M) { rep(j,L) { if (B[i]%C[j]==0) { cntB[j]++; } } } rep(i,L) { ans += cntA[i]*cntB[L-1-i]; } cout << ans << ln; */ map dicA,dicB; rep(i,N) { ll num = gcd(K,A[i]); dicA[num]++; } rep(i,M) { ll num = gcd(K,B[i]); dicB[num]++; } ll ans = 0; for (auto i : dicA) { ll num = K / i.first; vector C; for (int j = 1; j*j <= i.first; ++j) { if (i.first%j==0) { C.emplace_back(j); if (j*j < i.first) C.emplace_back(i.first/j); } } for (auto j : C) { ans += i.second*dicB[j*num]; } } cout << ans << ln; } }