結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
![]() |
提出日時 | 2021-01-24 00:26:11 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 634 ms / 2,000 ms |
コード長 | 1,114 bytes |
コンパイル時間 | 318 ms |
コンパイル使用メモリ | 82,716 KB |
実行使用メモリ | 64,728 KB |
最終ジャッジ日時 | 2025-01-01 22:20:10 |
合計ジャッジ時間 | 4,623 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
import syssys.setrecursionlimit(10**7)def I(): return int(sys.stdin.readline().rstrip())def MI(): return map(int,sys.stdin.readline().rstrip().split())def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))def LI2(): return list(map(int,sys.stdin.readline().rstrip()))def S(): return sys.stdin.readline().rstrip()def LS(): return list(sys.stdin.readline().rstrip().split())def LS2(): return list(sys.stdin.readline().rstrip())def ext_gcd(a, b): # a*x+b*y == gcd(a,b) たる gcd(a,b),x,yif b > 0:d,x,y = ext_gcd(b,a % b)return d,y,x-(a//b)*yreturn a,1,0T = I()mod = 10**9+7for _ in range(T):N,K,H,Y = MI()# Nを最大としてよいM = max(N,K,H)if M == K:N,K = K,Nelif M == H:N,H = H,Ng,a,b = ext_gcd(K,H)K //= gH //= gans = 0for i in range(Y//N+1):Z = Y-N*iif Z % g != 0:continueZ //= g# K*j+H*k == Z の自然数解の個数c = a*Zd = b*Zans += max(1+d//K-(-c+H-1)//H,0)ans %= modprint(ans)