結果

問題 No.321 (P,Q)-サンタと街の子供たち
ユーザー drken1215drken1215
提出日時 2020-10-07 03:17:02
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 3,082 bytes
コンパイル時間 2,251 ms
コンパイル使用メモリ 204,020 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-20 03:13:31
合計ジャッジ時間 6,023 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 RE -
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 1 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 1 ms
5,376 KB
testcase_14 AC 39 ms
5,376 KB
testcase_15 AC 79 ms
5,376 KB
testcase_16 AC 31 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 42 ms
5,376 KB
testcase_19 AC 56 ms
5,376 KB
testcase_20 AC 79 ms
5,376 KB
testcase_21 AC 45 ms
5,376 KB
testcase_22 AC 14 ms
5,376 KB
testcase_23 AC 62 ms
5,376 KB
testcase_24 AC 17 ms
5,376 KB
testcase_25 AC 47 ms
5,376 KB
testcase_26 AC 75 ms
5,376 KB
testcase_27 AC 57 ms
5,376 KB
testcase_28 AC 55 ms
5,376 KB
testcase_29 AC 79 ms
5,376 KB
testcase_30 AC 4 ms
5,376 KB
testcase_31 AC 2 ms
5,376 KB
testcase_32 AC 41 ms
5,376 KB
testcase_33 AC 38 ms
5,376 KB
testcase_34 AC 29 ms
5,376 KB
testcase_35 AC 61 ms
5,376 KB
testcase_36 AC 4 ms
5,376 KB
testcase_37 AC 48 ms
5,376 KB
testcase_38 AC 38 ms
5,376 KB
testcase_39 AC 14 ms
5,376 KB
testcase_40 AC 73 ms
5,376 KB
testcase_41 AC 21 ms
5,376 KB
testcase_42 AC 65 ms
5,376 KB
testcase_43 AC 30 ms
5,376 KB
testcase_44 AC 59 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

// Gauss Integer
struct gint {
    long long x, y;
    constexpr gint() : x(0), y(0) { }
    constexpr gint(long long x) : x(x), y(0) { }
    constexpr gint(long long x, long long y) : x(x), y(y) { }
    friend constexpr long long abs(const gint& r) noexcept {
        return r.x * r.x + r.y * r.y;
    }
    constexpr gint operator - () const noexcept {
        return gint(-x, -y);
    }
    constexpr gint operator + (const gint& r) const noexcept { return gint(*this) += r; }
    constexpr gint operator - (const gint& r) const noexcept { return gint(*this) -= r; }
    constexpr gint operator * (const gint& r) const noexcept { return gint(*this) *= r; }
    constexpr gint operator / (const gint& r) const noexcept { return gint(*this) /= r; }
    constexpr gint operator % (const gint& r) const noexcept { return gint(*this) %= r; }
    constexpr gint& operator += (const gint& r) noexcept {
        x += r.x, y += r.y;
        return *this;
    }
    constexpr gint& operator -= (const gint& r) noexcept {
        x -= r.x, y -= r.y;
        return *this;
    }
    constexpr gint& operator *= (const gint& r) noexcept {
        auto tx = x * r.x - y * r.y;
        auto ty = x * r.y + y * r.x;
        x = tx, y = ty;
        return *this;
    }
    constexpr gint& operator /= (const gint& r) noexcept {
        long long a = x, b = y, c = r.x, d = r.y;        
    	assert(c != 0 || d != 0);
        long long tx = (a * c + b * d) / (c * c + d * d);
        long long ty = (b * c - a * d) / (c * c + d * d);
        gint res;
        long long dmin = -1;
        for (long long nx = tx - 1; nx <= tx + 1; ++nx) {
            for (long long ny = ty - 1; ny <= ty + 1; ++ny) {
                gint q(nx, ny);
                long long d = abs((*this) - q * r);
                if (dmin == -1 || dmin > d) {
                    dmin = d;
                    res = q;
                }
            }
        }
        return *this = res;
    }
    constexpr gint& operator %= (const gint& r) noexcept {
        gint q = (*this) / r;
        return (*this) -= q * r;
    }
    constexpr bool operator == (const gint& r) const noexcept {
        return x == r.x && y == r.y;
    }
    constexpr bool operator != (const gint& r) const noexcept {
        return x != r.x || y != r.y;
    }
    friend ostream& operator << (ostream& os, const gint& r) noexcept {
        if (r.x == 0 && r.y == 0) return os << "0";
        else if (r.x == 0) return os << r.y << "i";
        else os << r.x;
        if (r.y == 0) return os;
        else if (r.y > 0) return os << " + " << r.y << "i";
        else return os << " - " << (-r.y) << "i";
    }
    friend gint gcd(const gint& x, const gint& y) {
        if (y == 0) return x;
        else return gcd(y, x % y);
    }
};

int main() {
    long long N, P, Q;
    cin >> P >> Q;
    gint G = gcd(gint(P, Q), gint(P, -Q));
    cin >> N;
    int res = 0;
    while (N--) {
        long long X, Y;
        cin >> X >> Y;
        if (gint(X, Y) % G == 0) ++res;
    }
    cout << res << endl;
}
0