結果
問題 | No.849 yuki国の分割統治 |
ユーザー | titia |
提出日時 | 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 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 26 |
ソースコード
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)))