結果
問題 |
No.474 色塗り2
|
ユーザー |
![]() |
提出日時 | 2016-12-21 23:35:02 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,116 bytes |
コンパイル時間 | 460 ms |
コンパイル使用メモリ | 82,456 KB |
実行使用メモリ | 126,236 KB |
最終ジャッジ日時 | 2024-12-14 14:29:31 |
合計ジャッジ時間 | 9,715 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 3 WA * 1 |
ソースコード
def ext_gcd(a, b): if b == 0: return [1, 0] else: y, x = ext_gcd(b, a % b) return [x, y - (a // b) * x] def minv(n, m): return ext_gcd(n, m)[0] mod = 2 ** 21 fact_factor2 = [0] * mod fact_mod = [0] * mod fact_inv_mod = [0] * mod fact_mod[0] = 1 fact_inv_mod[0] = 1 for i in range(1, mod): j = i f2 = 0 while j % 2 == 0: f2 += 1 j = j // 2 fact_mod[i] = fact_mod[i-1] * j % mod fact_inv_mod[i] = minv(fact_mod[i], mod) fact_factor2[i] = fact_factor2[i-1] + f2 def solve(A, B, C): if C == 2: return 0 else: f2 = fact_factor2[C+B-1] - fact_factor2[B] - fact_factor2[C-1] # val = fact_mod[C+B-1] * minv(fact_mod[B], mod) % mod * minv(fact_mod[C-1], mod) % mod val = fact_mod[C+B-1] * fact_inv_mod[B] % mod * fact_inv_mod[C-1] % mod val = ((C * val % mod) << min(32, f2)) % mod + A - 1 if (val & A) == A: return 1 else: return 0 testcase = int(input()) for _ in range(testcase): A, B, C = list(map(int, input().split())) print(solve(A, B, C))