結果

問題 No.849 yuki国の分割統治
ユーザー titiatitia
提出日時 2019-07-12 07:17:17
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 278 ms / 2,000 ms
コード長 1,007 bytes
コンパイル時間 308 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 49,992 KB
最終ジャッジ日時 2024-10-07 01:48:24
合計ジャッジ時間 6,622 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 27 ms
11,008 KB
testcase_01 AC 27 ms
10,752 KB
testcase_02 AC 27 ms
10,880 KB
testcase_03 AC 29 ms
10,880 KB
testcase_04 AC 30 ms
10,752 KB
testcase_05 AC 28 ms
10,880 KB
testcase_06 AC 27 ms
10,880 KB
testcase_07 AC 261 ms
44,416 KB
testcase_08 AC 238 ms
40,064 KB
testcase_09 AC 251 ms
42,496 KB
testcase_10 AC 266 ms
42,240 KB
testcase_11 AC 227 ms
38,144 KB
testcase_12 AC 252 ms
37,760 KB
testcase_13 AC 269 ms
43,776 KB
testcase_14 AC 245 ms
38,912 KB
testcase_15 AC 229 ms
38,528 KB
testcase_16 AC 261 ms
46,464 KB
testcase_17 AC 273 ms
43,776 KB
testcase_18 AC 267 ms
46,336 KB
testcase_19 AC 268 ms
44,288 KB
testcase_20 AC 278 ms
46,464 KB
testcase_21 AC 251 ms
37,632 KB
testcase_22 AC 271 ms
46,464 KB
testcase_23 AC 265 ms
46,336 KB
testcase_24 AC 261 ms
44,288 KB
testcase_25 AC 261 ms
49,992 KB
testcase_26 AC 27 ms
10,752 KB
testcase_27 AC 29 ms
10,880 KB
testcase_28 AC 29 ms
10,880 KB
testcase_29 AC 27 ms
10,752 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import math

input = sys.stdin.readline

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

a,b,c,d=map(int,input().split())
N=int(input())
P=[list(map(int,input().split())) for i in range(N)]


if a*d-b*c==0 and a!=0:
    a1=math.gcd(a,c)
    b1=b//(a//a1)

    def rep_point(p,q):
        k=p//a1

        return p-k*a1,q-k*b1

    Q=[rep_point(P[i][0],P[i][1]) for i in range(N)]

    print(len(set(Q)))
    sys.exit()
    
elif a*d-b*c==0:
    b1=math.gcd(b,d)
    a1=a//(b//b1)

    def rep_point(p,q):
        k=q//b1

        return p-k*a1,q-k*b1

    Q=[rep_point(P[i][0],P[i][1]) for i in range(N)]

    print(len(set(Q)))
    sys.exit()

x = lcm(a,c)

if a==0:
    t=abs(b)
    a,b=c,d
    
    
else:

    while c!=0:
        if a>c:
            a,b,c,d=c,d,a,b
            
        r,q=divmod(c,a)
        c,d=q,d-b*r

    t = abs(d)   

def rep_point(p,q):
    k=p//a

    p-=k*a
    q-=k*b

    q=q%t

    return p,q

Q=[rep_point(P[i][0],P[i][1]) for i in range(N)]

print(len(set(Q)))
0