結果
問題 | No.1758 Lazy Segment Tree...? |
ユーザー | 👑 PCTprobability |
提出日時 | 2021-10-16 01:48:15 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 4,280 ms / 8,000 ms |
コード長 | 1,242 bytes |
コンパイル時間 | 271 ms |
コンパイル使用メモリ | 87,324 KB |
実行使用メモリ | 80,660 KB |
最終ジャッジ日時 | 2023-09-01 23:38:21 |
合計ジャッジ時間 | 62,622 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge16 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 72 ms
71,292 KB |
testcase_01 | AC | 70 ms
71,116 KB |
testcase_02 | AC | 164 ms
78,492 KB |
testcase_03 | AC | 71 ms
71,436 KB |
testcase_04 | AC | 71 ms
71,488 KB |
testcase_05 | AC | 71 ms
71,344 KB |
testcase_06 | AC | 101 ms
77,252 KB |
testcase_07 | AC | 133 ms
77,732 KB |
testcase_08 | AC | 131 ms
77,780 KB |
testcase_09 | AC | 176 ms
79,228 KB |
testcase_10 | AC | 145 ms
78,048 KB |
testcase_11 | AC | 164 ms
78,212 KB |
testcase_12 | AC | 163 ms
78,112 KB |
testcase_13 | AC | 171 ms
78,200 KB |
testcase_14 | AC | 131 ms
77,616 KB |
testcase_15 | AC | 163 ms
78,300 KB |
testcase_16 | AC | 170 ms
78,184 KB |
testcase_17 | AC | 174 ms
77,948 KB |
testcase_18 | AC | 151 ms
77,820 KB |
testcase_19 | AC | 159 ms
78,244 KB |
testcase_20 | AC | 152 ms
77,728 KB |
testcase_21 | AC | 3,306 ms
79,488 KB |
testcase_22 | AC | 2,394 ms
79,224 KB |
testcase_23 | AC | 3,290 ms
80,660 KB |
testcase_24 | AC | 2,795 ms
79,084 KB |
testcase_25 | AC | 2,038 ms
78,968 KB |
testcase_26 | AC | 3,270 ms
79,628 KB |
testcase_27 | AC | 2,346 ms
79,468 KB |
testcase_28 | AC | 3,309 ms
79,476 KB |
testcase_29 | AC | 3,475 ms
79,332 KB |
testcase_30 | AC | 4,001 ms
79,724 KB |
testcase_31 | AC | 2,583 ms
79,288 KB |
testcase_32 | AC | 3,238 ms
80,212 KB |
testcase_33 | AC | 3,344 ms
79,156 KB |
testcase_34 | AC | 2,685 ms
79,284 KB |
testcase_35 | AC | 2,106 ms
79,564 KB |
testcase_36 | AC | 4,280 ms
79,784 KB |
testcase_37 | AC | 4,102 ms
79,520 KB |
testcase_38 | AC | 3,630 ms
79,580 KB |
ソースコード
mod=998244353 def modpow(a,b,m): res=1 while b: if b%2: res*=a res%=m a*=a a%=m b//=2 return res def powermat(a,n): x=len(a) b=[[0 for i in range(3)] for j in range(3)] c=[[0 for i in range(3)] for j in range(3)] for i in range(3): b[i][i]=1 while n: if n%2: for i in range(3): for j in range(3): c[i][j]=0 for i in range(3): for j in range(3): for k in range(3): c[i][j]+=a[i][k]*b[k][j] for i in range(3): for j in range(3): b[i][j]=c[i][j]%mod for i in range(3): for j in range(3): c[i][j]=0 for i in range(3): for j in range(3): for k in range(3): c[i][j]+=a[i][k]*a[k][j] for i in range(3): for j in range(3): a[i][j]=c[i][j]%mod n//=2 return b n,q=map(int,input().split()) ans=0 for i in range(1,n): x=i*(n-i+1) y=n*(n+1)//2-x z=(n-i)*(n-i+1)*(n-i+2)//6 x%=mod y%=mod z%=mod a=[[0 for i in range(3)] for j in range(3)] a[1][2]+=z; a[2][2]+=n*(n+1) a[0][0]+=y a[0][0]+=n*(n+1)//2 a[1][1]+=y; a[1][1]+=n*(n+1)//2 a[1][0]+=x a[0][1]+=x a=powermat(a,q) ans+=a[0][2] print(ans*modpow(2,mod-2,mod)%mod)