結果
問題 | No.1693 Invasion |
ユーザー | shotoyoo |
提出日時 | 2021-10-01 21:58:37 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,329 bytes |
コンパイル時間 | 147 ms |
コンパイル使用メモリ | 82,364 KB |
実行使用メモリ | 317,132 KB |
最終ジャッジ日時 | 2024-07-19 11:13:42 |
合計ジャッジ時間 | 7,237 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 38 ms
57,344 KB |
testcase_01 | AC | 36 ms
52,224 KB |
testcase_02 | AC | 37 ms
52,992 KB |
testcase_03 | AC | 42 ms
59,520 KB |
testcase_04 | AC | 38 ms
52,736 KB |
testcase_05 | AC | 43 ms
60,800 KB |
testcase_06 | AC | 44 ms
61,696 KB |
testcase_07 | AC | 44 ms
60,416 KB |
testcase_08 | AC | 42 ms
59,392 KB |
testcase_09 | AC | 456 ms
141,240 KB |
testcase_10 | AC | 238 ms
108,536 KB |
testcase_11 | AC | 131 ms
85,504 KB |
testcase_12 | AC | 138 ms
78,928 KB |
testcase_13 | AC | 97 ms
77,156 KB |
testcase_14 | AC | 114 ms
77,320 KB |
testcase_15 | AC | 89 ms
76,416 KB |
testcase_16 | AC | 471 ms
159,872 KB |
testcase_17 | WA | - |
testcase_18 | AC | 261 ms
79,164 KB |
testcase_19 | AC | 51 ms
65,920 KB |
testcase_20 | AC | 42 ms
59,776 KB |
testcase_21 | AC | 760 ms
232,072 KB |
testcase_22 | TLE | - |
testcase_23 | -- | - |
ソースコード
import sys input = lambda : sys.stdin.readline().rstrip() write = lambda x: sys.stdout.write(x+"\n") debug = lambda x: sys.stderr.write(x+"\n") writef = lambda x: print("{:.12f}".format(x)) n,m = list(map(int, input().split())) a = list(map(int, input().split())) M = 998244353 ### 素数の逆元とCombination N = n+m+10 # 必要なテーブルサイズ g1 = [0] * (N+1) # 元テーブル g2 = [0] * (N+1) #逆元テーブル inverse = [0] * (N+1) #逆元テーブル計算用テーブル g1[0] = g1[1] = g2[0] = g2[1] = 1 inverse[0], inverse[1] = [0, 1] for i in range( 2, N + 1 ): g1[i] = ( g1[i-1] * i ) % M inverse[i] = ( -inverse[M % i] * (M//i) ) % M # ai+b==0 mod M <=> i==-b*a^(-1) <=> i^(-1)==-b^(-1)*aより g2[i] = (g2[i-1] * inverse[i]) % M def cmb(n, r, M=M): if ( r<0 or r>n ): return 0 r = min(r, n-r) return ((g1[n] * g2[r] % M) * g2[n-r]) % M def perm(n, r, M=M): if (r<0 or r>n): return 0 return (g1[n] * g2[n-r]) % M dp = [0]*(m+1) dp[0] = 1 mask = (1<<m) - 1 for i in range(n): v = a[i] for j in range(m-v+1): dp[j+v] |= dp[j]<<1 dp[j+v] &= mask ans = 0 for l in range(m+1): val = dp[l] if val==0: continue k = (val&(-val)).bit_length()-1 ans += cmb(m-k, l-k) ans %= M # print(l,k) print(ans%M)