結果
問題 | No.1388 Less than K |
ユーザー | vwxyz |
提出日時 | 2023-11-19 02:24:04 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,356 bytes |
コンパイル時間 | 261 ms |
コンパイル使用メモリ | 81,700 KB |
実行使用メモリ | 98,228 KB |
最終ジャッジ日時 | 2023-11-19 02:24:19 |
合計ジャッジ時間 | 14,560 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 37 ms
53,460 KB |
testcase_01 | AC | 36 ms
53,460 KB |
testcase_02 | AC | 60 ms
68,624 KB |
testcase_03 | AC | 36 ms
53,460 KB |
testcase_04 | AC | 51 ms
64,040 KB |
testcase_05 | AC | 79 ms
75,080 KB |
testcase_06 | AC | 72 ms
73,292 KB |
testcase_07 | AC | 44 ms
59,756 KB |
testcase_08 | AC | 102 ms
76,240 KB |
testcase_09 | AC | 67 ms
71,236 KB |
testcase_10 | AC | 64 ms
71,228 KB |
testcase_11 | AC | 61 ms
68,624 KB |
testcase_12 | AC | 96 ms
76,776 KB |
testcase_13 | AC | 120 ms
77,452 KB |
testcase_14 | AC | 147 ms
77,676 KB |
testcase_15 | AC | 299 ms
79,200 KB |
testcase_16 | AC | 395 ms
80,260 KB |
testcase_17 | AC | 584 ms
81,720 KB |
testcase_18 | AC | 954 ms
84,332 KB |
testcase_19 | AC | 1,389 ms
86,680 KB |
testcase_20 | AC | 1,452 ms
87,240 KB |
testcase_21 | AC | 121 ms
87,388 KB |
testcase_22 | AC | 69 ms
74,136 KB |
testcase_23 | AC | 79 ms
78,196 KB |
testcase_24 | AC | 102 ms
81,880 KB |
testcase_25 | AC | 112 ms
83,200 KB |
testcase_26 | AC | 99 ms
83,900 KB |
testcase_27 | AC | 83 ms
80,272 KB |
testcase_28 | AC | 63 ms
72,380 KB |
testcase_29 | AC | 63 ms
72,180 KB |
testcase_30 | AC | 52 ms
64,296 KB |
testcase_31 | AC | 56 ms
66,840 KB |
testcase_32 | AC | 44 ms
59,748 KB |
testcase_33 | AC | 71 ms
75,476 KB |
testcase_34 | AC | 70 ms
74,696 KB |
testcase_35 | AC | 68 ms
73,920 KB |
testcase_36 | AC | 86 ms
83,188 KB |
testcase_37 | AC | 72 ms
77,096 KB |
testcase_38 | AC | 155 ms
87,464 KB |
testcase_39 | AC | 106 ms
85,764 KB |
testcase_40 | AC | 128 ms
84,972 KB |
testcase_41 | AC | 251 ms
91,324 KB |
testcase_42 | AC | 96 ms
87,448 KB |
testcase_43 | TLE | - |
testcase_44 | -- | - |
testcase_45 | -- | - |
testcase_46 | -- | - |
testcase_47 | -- | - |
testcase_48 | -- | - |
testcase_49 | -- | - |
testcase_50 | -- | - |
testcase_51 | -- | - |
testcase_52 | -- | - |
testcase_53 | -- | - |
testcase_54 | -- | - |
testcase_55 | -- | - |
testcase_56 | -- | - |
testcase_57 | -- | - |
testcase_58 | -- | - |
testcase_59 | -- | - |
testcase_60 | -- | - |
testcase_61 | -- | - |
testcase_62 | -- | - |
testcase_63 | -- | - |
testcase_64 | -- | - |
testcase_65 | -- | - |
testcase_66 | -- | - |
testcase_67 | -- | - |
testcase_68 | -- | - |
testcase_69 | -- | - |
testcase_70 | -- | - |
testcase_71 | -- | - |
testcase_72 | -- | - |
testcase_73 | -- | - |
testcase_74 | -- | - |
testcase_75 | -- | - |
testcase_76 | -- | - |
ソースコード
import sys readline=sys.stdin.readline H,W,K=map(int,readline().split()) H-=1;W-=1 K//=2 mod=998244353 fact=[1]+[i for i in range(1,H+W+1)] for i in range(1,H+W+1): fact[i]*=fact[i-1] fact[i]%=mod fact_inve=[i for i in range(1,H+W+1)]+[pow(fact[H+W],mod-2,mod)] for i in range(H+W-1,-1,-1): fact_inve[i]*=fact_inve[i+1] fact_inve[i]%=mod N=min(H,W) ans=0 if K<=600: for h in range(N+1): if h: prev=dp dp=[0]*(2*K+1) else: dp=[0]*(2*K+1) dp[K]=1 for w in range(max(h-K,0),min(h+K,N)+1): if h and abs((h-1)-w)<=K: dp[w-h+K]+=prev[w-(h-1)+K] if w and abs(h-(w-1))<=K: dp[w-h+K]+=dp[(w-1)-h+K] dp[w-h+K]%=mod ans+=fact[H+W]*fact_inve[H-h]%mod*fact_inve[W-h]%mod*fact_inve[2*h]%mod*dp[K]%mod ans%=mod else: for cnt in range(N+1): s=fact[2*cnt]*fact_inve[cnt]%mod*fact_inve[cnt]%mod for i in range(1,cnt//(K+1)+1): if i%2: s-=2*fact[2*cnt]*fact_inve[cnt-(K+1)*i]%mod*fact_inve[cnt+(K+1)*i]%mod else: s+=2*fact[2*cnt]*fact_inve[cnt-(K+1)*i]%mod*fact_inve[cnt+(K+1)*i]%mod s%=mod ans+=fact[H+W]*fact_inve[H-cnt]%mod*fact_inve[W-cnt]%mod*fact_inve[2*cnt]%mod*s%mod ans%=mod print(ans)