結果
| 問題 |
No.986 Present
|
| コンテスト | |
| ユーザー |
convexineq
|
| 提出日時 | 2020-02-11 15:47:27 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 253 ms / 2,000 ms |
| コード長 | 908 bytes |
| コンパイル時間 | 116 ms |
| コンパイル使用メモリ | 12,544 KB |
| 実行使用メモリ | 23,808 KB |
| 最終ジャッジ日時 | 2024-10-01 07:59:33 |
| 合計ジャッジ時間 | 4,411 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 30 |
ソースコード
def extgcd(x,y):
if y==0: return 1,0 #g=x
r0,r1,s0,s1 = x,y,1,0
while r1 != 0:
r0,r1, s0,s1 = r1,r0%r1, s1,s0-r0//r1*s1
#g = r0
return s0,(r0-s0*x)//y
def modinv(a,MOD):
x,y = extgcd(a,MOD)
return x%MOD
# coding: utf-8
# Your code here!
import sys
readline = sys.stdin.readline
read = sys.stdin.read
n,m = [int(i) if i != "-1" else 0 for i in read().split()]
MOD = 998244353
p2 = [1]*(m+1)
for i in range(m):
p2[i+1] = p2[i]*2%MOD
# ans1: 2^d
# ans2: dim<S> == d なる S の数
# ans3: dim V == d なる V \subset F2_^M の数
d = min(n,m)
ans1 = pow(2,d,MOD)
ans3 = [0]*(m+1)
ans3[0] = vd = fac = 1
for i in range(d):
ans3[i+1] = ans3[i]*(p2[m] - p2[i])%MOD
vd *= p2[d] - p2[i]
vd %= MOD
fac *= (i+1)
fac %= MOD
ans2 = ans3[d]*modinv(fac,MOD)%MOD
vdinv = modinv(vd,MOD)
ans3[d] *= vdinv
ans3[d] %= MOD
print(ans1,ans2,ans3[d])
convexineq