#include using namespace std; struct iofast_t { iofast_t() { ios::sync_with_stdio(false); cin.tie(nullptr); } } iofast; struct uns_t {} uns; template auto vec(Element init, Head arg, Args ...args) { if constexpr (sizeof...(Args) == 0) return vector(arg, init); else return vector(arg, vec(init, args...)); } template auto vec(uns_t, Head arg, Args ...args) { return vec(Element(), arg, args...); } template auto distance(const Container &c, decltype(begin(c)) iter) { return distance(begin(c), iter); } template ::value_type>> auto isort(RIter first, RIter last, Compare comp = Compare()) { vector i(distance(first, last)); iota(begin(i), end(i), 0); sort(begin(i), end(i), [&](auto x, auto y) { return comp(*(first + x), *(first + y)); }); return i; } template typename, typename = void_t<>> struct detect : false_type {}; template typename Check> struct detect>> : true_type {}; template typename Check> constexpr inline bool detect_v = detect::value; template using has_member_sort = decltype(declval().sort()); template > auto sorted(Container c, Compare comp = Compare()) { if constexpr (detect_v) { c.sort(comp); return c; } else { sort(begin(c), end(c), comp); return c; } } template > auto uniqued(Container c, Compare comp = Compare()) { c.erase(unique(begin(c), end(c), comp), end(c)); return c; } template > T &chmin(T &l, T r, Compare &&f = less()) { return l = min(l, r, f); } template > T &chmax(T &l, T r, Compare &&f = less()) { return l = max(l, r, f); } template constexpr auto fix(F &&f) noexcept { return [f = std::tuple(std::forward(f))](auto &&...args) mutable { return std::get<0>(f)(fix(std::get<0>(f)), std::forward(args)...); }; } #include int64_t sign(int64_t x) { return clamp(x, -1, 1); } int64_t floor(int64_t x, int64_t y) { if (0 < sign(x) * sign(y)) { return x / y; } return -((abs(x) + abs(y) - 1) / abs(y)); } int64_t ceil(int64_t x, int64_t y) { if (0 < sign(x) * sign(y)) { return (abs(x) + abs(y) - 1) / abs(y); } return x / y; } int main() { int n, m, p, q; cin >> n >> m >> p >> q; while (q--) { int x, f; cin >> x >> f; if (x % p == 0) { cout << (f == 0 ? m : 0) << '\n'; continue; } auto k = (f * atcoder::pow_mod(x, p - 2, p)) % p; auto l = max(0, ceil(1 - k, p)); auto r = floor(m - k, p); cout << max(0, r - l + 1) << '\n'; } cout << flush; }