#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; class range {private: struct I{int x;int operator*(){return x;}bool operator!=(I& lhs){return x ostream& operator<<(ostream& os, const pair& p){ return os << "{" << p.first << ", " << p.second << "}"; } template ostream& operator<<(ostream& os, const vector& obj) { os << "{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const set& obj) { os << "set{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const map& obj) { os << "map{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } #ifdef ONLINE_JUDGE #define dump(expr) ; #else #define dump(expr) { cerr << "\033[33m#L" << __LINE__ << ": " << expr << "\033[39m" << endl; } #endif ll extgcd(ll a, ll b, ll &x, ll &y) { ll g = a; x = 1; y = 0; if (b != 0) g = extgcd(b, a % b, y, x), y -= (a / b) * x; return g; } namespace solver { ll a, b, c, d; ll n; vector xs, ys; void read() { cin >> a >> b >> c >> d; cin >> n; xs.resize(n); ys.resize(n); for (int i : range(n)) cin >> xs[i] >> ys[i]; } void experiment() { set> st; dump(a) dump(b) dump(c) dump(d) for (ll s : range(-100, 100)) for (ll t : range(-100, 100)) { ll x = a * s + c * t; ll y = b * s + d * t; if (x >= 0 && y >= 0 && x <= 15 && y <= 15) st.insert({x, y}); } cout << st << endl; } using RetType = ll; RetType run() { ll A = abs(a * d - b * c); ll gx = gcd(a, c); ll gy = gcd(b, d); set> st; dump("A " << A); if (A == 0) { if (gx) a /= gx; if (gy) b /= gy; if (a == 0) { for (int i : range(n)) st.insert({xs[i], ys[i] % b}); } else { for (int i : range(n)) { ll x = xs[i], y = ys[i]; ll m = x / a; x -= m * a; y -= m * b; st.insert({x, y}); } } return st.size(); } ll h = gy; ll w = A / gy; ll s0, t0; extgcd(b, d, s0, t0); ll u = a * s0 + c * t0; for (int i : range(n)) { ll x = xs[i], y = ys[i]; ll m = y / gy; x -= m * u; y -= m * h; x = ((x % w) + w) % w; st.insert({x, y}); } return st.size(); } } // namespace template void run(F f) { if constexpr (std::is_same_v) f(); else cout << f() << endl; } int main(int argc, char** argv) { cerr << fixed << setprecision(12); cout << fixed << setprecision(12); int testcase = 1; if (argc > 1) testcase = atoi(argv[1]); while (testcase--) { solver::read(); } // solver::experiment(); run(solver::run); }