#include using namespace std; using LL = long long int; #define incID(i, l, r) for(int i = (l) ; i < (r); ++i) #define decID(i, l, r) for(int i = (r) - 1; i >= (l); --i) #define incII(i, l, r) for(int i = (l) ; i <= (r); ++i) #define decII(i, l, r) for(int i = (r) ; i >= (l); --i) #define inc(i, n) incID(i, 0, n) #define dec(i, n) decID(i, 0, n) #define inc1(i, n) incII(i, 1, n) #define dec1(i, n) decII(i, 1, n) #define inID(v, l, r) ((l) <= (v) && (v) < (r)) #define inII(v, l, r) ((l) <= (v) && (v) <= (r)) #define PB push_back #define EB emplace_back #define MP make_pair #define FI first #define SE second #define FR front() #define BA back() #define ALL(v) v.begin(), v.end() #define RALL(v) v.rbegin(), v.rend() auto setmin = [](auto & a, auto b) { return (b < a ? a = b, true : false); }; auto setmax = [](auto & a, auto b) { return (b > a ? a = b, true : false); }; auto setmineq = [](auto & a, auto b) { return (b <= a ? a = b, true : false); }; auto setmaxeq = [](auto & a, auto b) { return (b >= a ? a = b, true : false); }; #define SI(v) static_cast(v.size()) #define RF(e, v) for(auto & e: v) #define until(e) while(! (e)) #define if_not(e) if(! (e)) #define ef else if #define UR assert(false) #define IN(T, ...) T __VA_ARGS__; IN_(__VA_ARGS__); void IN_() { }; template void IN_(T & a, U & ... b) { cin >> a; IN_(b ...); }; template void OUT(T && a) { cout << a << endl; } template void OUT(T && a, U && ... b) { cout << a << " "; OUT(b ...); } // ---- ---- template istream & operator>>(istream & s, vector & v) { RF(e, v) { s >> e; } return s; } template ostream & operator<<(ostream & s, vector const & v) { inc(i, SI(v)) { s << (i == 0 ? "" : " ") << v[i]; } return s; } auto mo = [](auto a, auto b) { assert(b != 0); a %= b; if(a < 0) { a += abs(b); } return a; }; vector> prime_factorization(LL x) { assert(x > 0); vector> f; for(LL i = 2; i <= x; i++) { if(i * i > x) { i = x; } if(x % i == 0) { f.EB(i, 0); while(x % i == 0) { f.back().SE++; x /= i; } } } return f; } vector divisors(LL x) { auto pf = prime_factorization(x); vector d = { 1 }; for(auto e: pf) { int ds = d.size(); inc(i, ds) { LL v = d[i]; inc(j, e.SE) { v *= e.FI; d.PB(v); } } } sort(ALL(d)); return d; } template T gcd(T a, T b) { return (b == 0 ? a : gcd(b, a % b)); } template T lcm(T a, T b) { return a / gcd(a, b) * b; } int main() { IN(LL, n, m, k); IN(string, s); vector a(n), b(m); cin >> b >> a; sort(ALL(a)); sort(ALL(b)); LL ans = 0; if(s == "+") { map mp; inc(j, m) { mp[b[j] % k]++; } inc(i, n) { ans += mp[mo(k - a[i], k)]; } } else { map mp; auto D = divisors(k); inc(j, m) { b[j] = gcd(b[j], k); mp[b[j]]++; } RF(d1, D) { RF(d2, D) { if(d1 > d2 && d1 % d2 == 0) { mp[d2] += mp[d1]; } } } inc(i, n) { a[i] = gcd(a[i], k); ans += mp[k / a[i]]; } } OUT(ans); }