結果

問題 No.1554 array_and_me
ユーザー ygd.ygd.
提出日時 2021-06-20 13:57:42
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,883 bytes
コンパイル時間 308 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 95,744 KB
最終ジャッジ日時 2024-06-22 22:28:10
合計ジャッジ時間 5,926 ms
ジャッジサーバーID
(参考情報)
judge4 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 53 ms
63,872 KB
testcase_01 AC 97 ms
80,700 KB
testcase_02 AC 96 ms
80,756 KB
testcase_03 AC 96 ms
80,836 KB
testcase_04 AC 97 ms
80,832 KB
testcase_05 AC 98 ms
80,768 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 98 ms
89,088 KB
testcase_12 AC 100 ms
88,960 KB
testcase_13 AC 99 ms
89,344 KB
testcase_14 AC 99 ms
89,344 KB
testcase_15 AC 99 ms
88,960 KB
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 WA -
testcase_41 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

def cmb(n, r, p):
  if (r < 0) or (n < r):
    return 0
  r = min(r, n - r)
  return fac[n]*finv[r]*finv[n-r]%p

def perm(n,r,p):
  if (r < 0) or (n < r):
    return 0
  return fac[n]*finv[n-r]%p

N = 2*pow(10,5) + 1000
MOD = 998244353

fac = [-1]*(N+1); fac[0] = 1; fac[1] = 1 #階乗
finv = [-1]*(N+1); finv[0] = 1; finv[1] = 1 #階乗の逆元
inv = [-1]*(N+1); inv[0] = 0; inv[1] = 1 #逆元
for i in range(2,N+1):
  fac[i] = fac[i-1]*i%MOD
  inv[i] = MOD - inv[MOD%i]*(MOD//i)%MOD
  finv[i] = finv[i-1]*inv[i]%MOD

#print(cmb(5, 4, MOD))

def main():
    T = int(input())
    for _ in range(T):
        N,K = map(int,input().split())
        A = list(map(int,input().split()))
        S = sum(A)
        A.sort(reverse=True)
        num = [0]*N
        nokori = K
        for i in range(N):
            num[i] += K*A[i]//S
            nokori -= num[i]
        #print(num)
        idx = 0
        while nokori > 0:
            num[idx] += 1
            nokori -= 1
            idx += 1
        #print("UP",num)
        

        def calc(A,B,mod=998244353):
            k = sum(B) #回数
            s = sum(A) #全ての和
            chi = 1
            par = pow(s,k,mod)
            for i in range(len(A)):
                chi *= pow(A[i],B[i],mod)*cmb(k,B[i],mod)
                chi %= mod
                k -= B[i] #残りの回数が減る
            #print(chi,par)
            ret = chi*pow(par,mod-2,mod)
            ret %= mod
            return ret

        #for i in range(10):
        #    for j in range(10):
        #        for k in range(10):
        #            for l in range(10):
        #                num = [i,j,k,l,10-(i+j+k+l)]
        #                ans = calc(A,num)
        #                if ans == 484466660:
        #                    print(num)

        ans = calc(A,num)
        print(ans%MOD)



if __name__ == '__main__':
    main()
0