#include //#include //#include using namespace std; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(n);i++) #define rep1(i,n) for(int i=1;i<=(n);i++) #define rrep(i,n) for(int i=(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; #define Sp(p) cout< a) { swap(a, b); } ll r = a%b; while (r != 0) { a = b; b = r; r = a%b; } return b; } ll lcm(ll a, ll b) { return (a / gcd(a, b))*b; } map mp; void dfs(vector &yaku, int i, ll now) { if (i == yaku.size()) { //DEBUG(now); mp[now]++; return; } for (int j = 0; j <= yaku[i].second; j++) { dfs(yaku, i + 1, now); now *= yaku[i].first; } } void solve2(int n, int m, int k, vi a, vi b) { vi primes; vector is_prime(33000, true); is_prime[0] = is_prime[1] = false; for (int x = 2; x < 33000; x++) { if (is_prime[x]) { primes.push_back(x); for (int j = 2; x * j < 33000; j++) { is_prime[x * j] = false; } } } //DEBUG(primes.size()); rep (i, m) { int x = b[i]; //DEBUG(x); vector yaku; for (int p: primes) { if (x < p) break; int cnt = 0; while (x % p == 0) { x /= p; cnt++; } if (cnt > 0) { yaku.push_back(pii(p, cnt)); } } if (x > 1) yaku.push_back(pii(x, 1)); dfs(yaku, 0, 1); //DEBUG_VEC(yaku); } /* for (auto itr = mp.begin(); itr != mp.end(); itr++) { DEBUG(*itr); } //*/ ll ans = 0; rep (i, n) { ll g = gcd(k, a[i]); if (g == 1) { ans += m; continue; } else { ans += mp[k / g]; } } cout << ans << endl; } signed main() { fio(); int n, m, k; cin >> n >> m >> k; char op; cin >> op; vi a(n), b(m); rep (i, m) cin >> b[i]; rep (i, n) cin >> a[i]; if (op == '+') solve1(n, m, k, a, b); else solve2(n, m, k, a, b); }