結果
問題 | No.2324 Two Countries within UEC |
ユーザー | moharan627 |
提出日時 | 2023-05-28 15:29:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 281 ms / 2,000 ms |
コード長 | 2,142 bytes |
コンパイル時間 | 600 ms |
コンパイル使用メモリ | 87,212 KB |
実行使用メモリ | 81,308 KB |
最終ジャッジ日時 | 2023-09-09 06:16:06 |
合計ジャッジ時間 | 10,760 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 267 ms
81,276 KB |
testcase_01 | AC | 275 ms
80,900 KB |
testcase_02 | AC | 161 ms
78,632 KB |
testcase_03 | AC | 242 ms
80,232 KB |
testcase_04 | AC | 255 ms
80,692 KB |
testcase_05 | AC | 160 ms
78,348 KB |
testcase_06 | AC | 183 ms
78,884 KB |
testcase_07 | AC | 154 ms
78,608 KB |
testcase_08 | AC | 239 ms
79,740 KB |
testcase_09 | AC | 246 ms
80,980 KB |
testcase_10 | AC | 281 ms
81,124 KB |
testcase_11 | AC | 278 ms
81,132 KB |
testcase_12 | AC | 177 ms
78,856 KB |
testcase_13 | AC | 175 ms
79,176 KB |
testcase_14 | AC | 160 ms
78,516 KB |
testcase_15 | AC | 161 ms
78,460 KB |
testcase_16 | AC | 234 ms
79,848 KB |
testcase_17 | AC | 195 ms
79,132 KB |
testcase_18 | AC | 161 ms
78,604 KB |
testcase_19 | AC | 280 ms
80,896 KB |
testcase_20 | AC | 153 ms
78,364 KB |
testcase_21 | AC | 223 ms
78,496 KB |
testcase_22 | AC | 240 ms
78,984 KB |
testcase_23 | AC | 169 ms
78,004 KB |
testcase_24 | AC | 236 ms
78,596 KB |
testcase_25 | AC | 281 ms
81,308 KB |
testcase_26 | AC | 152 ms
77,560 KB |
testcase_27 | AC | 155 ms
78,104 KB |
testcase_28 | AC | 156 ms
77,660 KB |
testcase_29 | AC | 156 ms
77,708 KB |
testcase_30 | AC | 155 ms
77,952 KB |
testcase_31 | AC | 95 ms
71,484 KB |
testcase_32 | AC | 96 ms
71,388 KB |
testcase_33 | AC | 97 ms
71,580 KB |
testcase_34 | AC | 99 ms
71,452 KB |
testcase_35 | AC | 95 ms
71,572 KB |
testcase_36 | AC | 95 ms
71,580 KB |
testcase_37 | AC | 96 ms
71,640 KB |
testcase_38 | AC | 95 ms
71,332 KB |
testcase_39 | AC | 95 ms
71,452 KB |
testcase_40 | AC | 95 ms
71,620 KB |
testcase_41 | AC | 96 ms
71,440 KB |
testcase_42 | AC | 94 ms
71,636 KB |
ソースコード
import sys #sys.setrecursionlimit(10 ** 6) INF = float('inf') #10**20,2**63に変えるのもあり MOD = 10**9 + 7 MOD2 = 998244353 from collections import defaultdict def ceil(A,B): return -(-A//B) def inv_gcd(a, b): a = a % b if a == 0: return (b, 0) s = b; t = a m0 = 0; m1 = 1 while (t): u = s // t s -= t * u m0 -= m1 * u s, t = t, s m0, m1 = m1, m0 if m0 < 0: m0 += b // s return (s, m0) def inv_mod(x, m): assert 1 <= m z = inv_gcd(x, m) assert z[0] == 1 return z[1] def crt(r, m): assert len(r) == len(m) n = len(r) r0 = 0; m0 = 1 for i in range(n): assert 1 <= m[i] r1 = r[i] % m[i] m1 = m[i] if m0 < m1: r0, r1 = r1, r0 m0, m1 = m1, m0 if (m0 % m1 == 0): if (r0 % m1 != r1): return (0, -1) continue g, im = inv_gcd(m0, m1) u1 = m1 // g if ((r1 - r0) % g): return (0, -1) x = (r1 - r0) // g % u1 * im % u1 r0 += x * m0 m0 *= u1 if r0 < 0: r0 += m0 return (r0, m0) def solve(): def II(): return int(sys.stdin.readline()) def LI(): return list(map(int, sys.stdin.readline().split())) def LC(): return list(input()) def IC(): return [int(c) for c in input()] def MI(): return map(int, sys.stdin.readline().split()) N,M,P,Q = MI() import math for q in range(Q): X,F = MI() X%=P if(X==0): if(F==0): print(M) else: print(0) continue G = math.gcd(X,P) #G<=P #print(G) r,m=crt([F,0],[P,X]) if (m == -1): print(0) continue period = P//G # G<=P #周期 ans = M//period rest_m = M%period #print("R:",r,"M:",m) if(r==0): r+=m #print(r//X,rest_m) if(r//X<=rest_m): print(ans+1) else: print(ans) return solve()