結果

問題 No.321 (P,Q)-サンタと街の子供たち
ユーザー NoneNone
提出日時 2021-04-28 11:37:17
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 3,444 bytes
コンパイル時間 1,333 ms
コンパイル使用メモリ 82,256 KB
実行使用メモリ 76,644 KB
最終ジャッジ日時 2024-07-07 05:52:25
合計ジャッジ時間 7,056 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 AC 35 ms
52,864 KB
testcase_02 AC 33 ms
52,480 KB
testcase_03 AC 33 ms
52,992 KB
testcase_04 AC 33 ms
52,864 KB
testcase_05 RE -
testcase_06 WA -
testcase_07 AC 32 ms
52,736 KB
testcase_08 WA -
testcase_09 AC 34 ms
52,992 KB
testcase_10 AC 33 ms
53,120 KB
testcase_11 WA -
testcase_12 AC 34 ms
53,120 KB
testcase_13 AC 33 ms
53,248 KB
testcase_14 AC 89 ms
76,464 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 AC 90 ms
76,288 KB
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 AC 104 ms
75,868 KB
testcase_27 WA -
testcase_28 AC 91 ms
76,416 KB
testcase_29 WA -
testcase_30 AC 58 ms
70,016 KB
testcase_31 WA -
testcase_32 AC 82 ms
76,612 KB
testcase_33 AC 83 ms
76,160 KB
testcase_34 AC 81 ms
75,880 KB
testcase_35 WA -
testcase_36 AC 55 ms
69,248 KB
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 AC 103 ms
76,288 KB
testcase_41 WA -
testcase_42 WA -
testcase_43 WA -
testcase_44 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

def my_gcd(a, b):
    while b != 0:
        a, b = b, a % b
    return a

def gcd2(array):
    temp = array[0]
    for a in array[1:]:
        temp = gcd(temp, a)
    return temp

def lcm(a,b):
    return a*b//gcd(a,b)

def lcm2(array):
    temp = 1
    for a in array:
        temp = lcm(temp, a)
    return temp

def gcd_ext(a0,b0):
    """
    ax + by = gcd(a,b) を満たす|x|が最小の(x,y)
    一般解は (x+k*b//gcd(a,b),y-k*a//gcd(a,b)) と表される

    備考:
    ax + by = G (=m*gcd(a,b)) の場合は、
    単に m=G//gcd を (x,y) に掛ければよい
    """
    a,b=abs(a0),abs(b0)
    sign_a,sing_b=(a0>0)-(a0<0),(b0>0)-(b0<0)
    x0, y0, x, y = 0, 1, 1, 0
    while b!=0:
        q=a//b
        a,b=b,a%b
        x0, y0, x, y = x-q*x0, y-q*y0, x0, y0
    return (x*sign_a,y*sing_b)

def bezout_identity(a,b,G,cut=10**6,positive=True):
    """
    ベズーの等式
    ax + by = G
    を満たす解(x,y)
    G は gcdの倍数でなければならない

    計算量: O(cut)
    """
    x,y = gcd_ext(a,b)
    g=a*x+b*y
    if G%g: return []
    m=G//g
    xmax,xmin=cut,-cut*(positive==False)
    kmax,kmin=(xmax-x*m)//abs(b//g), (xmin-x*m-1)//abs(b//g)+1
    table=[]
    if positive: # 正数解のみ
        for k in range(kmin,kmax+1):
            xx,yy=x*m+k*b//g,y*m-k*a//g
            if xx>cut or xx<0 or yy>cut or yy<0: continue
            table.append((xx,yy))
        return table
    else:
        for k in range(kmin,kmax+1):
            xx,yy=x*m+k*b//g,y*m-k*a//g
            if abs(xx)>cut or abs(yy)>cut: continue
            table.append((xx,yy))
        return table

def bezout_identity_lower_bound(a,b,G,x0=0,equal=True):
    """
    ax + by = G を満たすxが x0 以上（より大きい）のうちで最小の(x,y)
    一般解は (x+k*b//gcd(a,b),y-k*a//gcd(a,b)) と表される
    """
    x,y = gcd_ext(a,b)
    g=a*x+b*y

    if G%g: return None
    m=G//g
    bg=abs(b//g)
    if equal==True:
        if b>=0:
            k=(x0-x*m+bg-1)//bg
        else:
            k=(x*m-x0)//bg
    else:
        if b>=0:
            k=(x0-x*m+bg)//bg
        else:
            k=(x*m-x0-1)//bg
    return x*m+k*b//g,y*m-k*a//g

def bezout_identity_upper_bound(a,b,G,x0=0,equal=True):
    """
    ax + by = G を満たす xが x0 以下(未満)のうちで最大の(x,y)
    一般解は (x+k*b//gcd(a,b),y-k*a//gcd(a,b)) と表される
    """
    x,y = gcd_ext(a,b)
    g=a*x+b*y

    if G%g: return None
    m=G//g
    bg=abs(b//g)
    if equal==True:
        if b>=0:
            k=(x0-x*m)//bg
        else:
            k=(x*m-x0-1)//bg+1
    else:
        if b>=0:
            k=(x0-x*m+bg-1)//bg-1
        else:
            k=(x*m-x0)//bg+1

    return x*m+k*b//g,y*m-k*a//g

def mod_inv(a,MOD):
    """
    gcd(a,MOD)=1 を満たす必要あり
    """
    x,y=gcd_ext(a,MOD)
    return x%MOD

def mod_solve(a,b,MOD):
    """
    ax=b (mod MOD)
    """
    x,y=gcd_ext(a,MOD)
    g=a*x+MOD*y
    a_inv=x%MOD
    if b%g: return None
    a,b,MOD=a//g,b//g,MOD//g
    return b*a_inv%MOD

######################################################

from math import gcd
P,Q=map(int, input().split())
g=gcd(P,Q)
a,b=gcd_ext(P,Q)
N=int(input())
cnt=0
for _ in range(N):
    x,y=map(int, input().split())
    if x%g==0 and y%g==0:
        mx=x//g
        my=y//g
        if (((x//g+y//g)*a%2==0 and (x//g+y//g)*b%2==0)) or ((x//g)*(b//g)+(y//g)*(a//g))%2==1:
            cnt+=1
print(cnt)
0