結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー | 👑 SPD_9X2 |
提出日時 | 2021-01-22 23:19:41 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,407 ms / 2,000 ms |
コード長 | 1,569 bytes |
コンパイル時間 | 280 ms |
コンパイル使用メモリ | 87,280 KB |
実行使用メモリ | 77,196 KB |
最終ジャッジ日時 | 2023-08-27 22:29:45 |
合計ジャッジ時間 | 9,169 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 89 ms
76,240 KB |
testcase_01 | AC | 78 ms
71,580 KB |
testcase_02 | AC | 75 ms
71,392 KB |
testcase_03 | AC | 77 ms
71,380 KB |
testcase_04 | AC | 76 ms
71,444 KB |
testcase_05 | AC | 75 ms
71,504 KB |
testcase_06 | AC | 111 ms
76,832 KB |
testcase_07 | AC | 106 ms
76,588 KB |
testcase_08 | AC | 115 ms
77,080 KB |
testcase_09 | AC | 116 ms
77,108 KB |
testcase_10 | AC | 112 ms
76,968 KB |
testcase_11 | AC | 1,298 ms
76,300 KB |
testcase_12 | AC | 1,107 ms
76,960 KB |
testcase_13 | AC | 1,407 ms
76,956 KB |
testcase_14 | AC | 1,033 ms
76,836 KB |
testcase_15 | AC | 1,174 ms
76,952 KB |
testcase_16 | AC | 833 ms
77,196 KB |
ソースコード
""" """ import sys from sys import stdin import math def extGCD2(a,b): if b: d,y,x = extGCD(b,a % b) y -= a // b * x return d,x,y return a,1,0 def extGCD(a,b): g = math.gcd(a,b) x, y, u, v = 1, 0, 0, 1 while b: k = a // b x -= k * u y -= k * v x, u = u, x y, v = v, y a, b = b, a % b return g ,x, y TT = int(stdin.readline()) mod = 10**9+7 ans = 0 for loop in range(TT): N,K,H,Y = map(int,stdin.readline().split()) A = [N,K,H] A.sort() A.reverse() ans = 0 aa,bb = A[1],A[2] gg,xx,yy = extGCD(aa,bb) for i in range(0,Y // A[0] + 1): rem = Y - A[0] * i #print (rem) a,b,g,x,y = aa,bb,gg,xx,yy if rem % g != 0: continue x *= rem // g y *= rem // g lcm = a * b // g p = lcm // a q = lcm // b #xをp増やし、yをq減らす #xをp減らし、yをq増やす ができる if x >= 0: move = x // p xa,ya = x - move * p , y + move * q else: move = (abs(x)+p-1) // p xa,ya = x + move * p , y - move * q if y >= 0: move = y // q xb,yb = x + move * p, y - move * q else: move = (abs(y)+q-1) // q xb,yb = x - move * p , y + move * q if ya < 0 or xb < 0: continue #print (rem , xa , ya , xb , yb) ans += max(0 , (xb - xa) // p + 1) ans %= mod print (ans)