結果

問題 No.2474 Empty Quartz
ユーザー 👑 SPD_9X2SPD_9X2
提出日時 2023-07-22 03:39:56
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 167 ms / 2,000 ms
コード長 1,525 bytes
コンパイル時間 316 ms
コンパイル使用メモリ 82,528 KB
実行使用メモリ 92,292 KB
最終ジャッジ日時 2024-06-11 22:48:47
合計ジャッジ時間 2,256 ms
ジャッジサーバーID
(参考情報)
judge5 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 65 ms
77,364 KB
testcase_01 AC 65 ms
77,208 KB
testcase_02 AC 133 ms
92,292 KB
testcase_03 AC 120 ms
91,744 KB
testcase_04 AC 145 ms
91,796 KB
testcase_05 AC 167 ms
91,676 KB
testcase_06 AC 154 ms
91,604 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

"""

Empty Quartz (7) 想定解

"""

import sys
from sys import stdin

mod = 998244353

from collections import deque

def modfac(n, MOD):
 
    f = 1
    factorials = [1]
    for m in range(1, n + 1):
        f *= m
        f %= MOD
        factorials.append(f)
    inv = pow(f, MOD - 2, MOD)
    invs = [1] * (n + 1)
    invs[n] = inv
    for m in range(n, 1, -1):
        inv *= m
        inv %= MOD
        invs[m - 1] = inv
    return factorials, invs

def modnCr(n,r):
    if n < 0 or r < 0 or n < r:
        return 0
    return fac[n] * inv[n-r] * inv[r] % mod

fac,inv = modfac(200000,mod)

T = int(stdin.readline())

assert 1 <= T <= 10**5

ANS = []

for loop in range(T):

    N,K = map(int,stdin.readline().split())
    
    assert 1 <= N <= 10**5
    assert 0 <= K <= N*(N+1)//2

    # x^2 - (N+1)x + K = 0 を解く
    in_root = (N+1)**2 - 4 * K

    if in_root < 0:
        ANS.append(0)
        continue

    # 10^6 > root の時、pythonでこの処理で大丈夫なことはチェック済(下部)
    root = int(in_root**0.5)
    if root**2 != in_root:
        ANS.append(0)
        continue

    topa = (N+1) + root
    topb = (N+1) - root

    ans = 0
    if topa % 2 == 0 and 1 <= topa // 2 <= N+1:
        ans += modnCr(N,topa//2-1)
    if topb % 2 == 0 and topa != topb and 1 <= topb <= N+1:
        ans += modnCr(N,topb//2-1)
    ANS.append(ans % mod)

print (*ANS,sep="\n")

"""    
# 精度チェック用
for i in range(1,10**6):
	
	x = i**2
	
	y = int(x**0.5)
	
	if y**2 != x:
		print (x,y)
"""
0