結果

問題 No.2950 Max Min Product
ユーザー sasa8uyauyasasa8uyauya
提出日時 2024-10-26 00:12:48
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,633 bytes
コンパイル時間 365 ms
コンパイル使用メモリ 82,480 KB
実行使用メモリ 115,428 KB
最終ジャッジ日時 2024-10-26 00:12:55
合計ジャッジ時間 6,007 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
61,252 KB
testcase_01 AC 40 ms
53,272 KB
testcase_02 AC 41 ms
53,824 KB
testcase_03 TLE -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

n=int(input())
a=list(map(int,input().split()))
M=998244353
q=[(a[i],i) for i in range(n)]
q.sort()
p=[[] for i in range(n+1)]
B=int(n**0.5)+1
st1=[0]*B*B
st2=[0]*B
for v,i in q[::-1]:
  st1[i]+=1
  st2[i//B]+=1
  ok=i
  ng=-1
  while abs(ok-ng)>1:
    m=(ok+ng)//2
    g=0
    l=m
    r=i
    yl=l//B
    yr=r//B
    if yl==yr:
      g+=sum(st1[l:r+1])
    else:
      g+=sum(st1[l:yl*B+B])
      g+=sum(st2[yl+1:yr])
      g+=sum(st1[yr*B:r+1])
    if g==1:
      ok=m
    else:
      ng=m
  L=ok
  ok=i
  ng=n
  while abs(ok-ng)>1:
    m=(ok+ng)//2
    g=0
    l=i
    r=m
    yl=l//B
    yr=r//B
    if yl==yr:
      g+=sum(st1[l:r+1])
    else:
      g+=sum(st1[l:yl*B+B])
      g+=sum(st2[yl+1:yr])
      g+=sum(st1[yr*B:r+1])
    if g==1:
      ok=m
    else:
      ng=m
  R=ok
  p[i]+=[(L,i,v)]
  p[R+1]+=[(L,i,pow(v,M-2,M))]
st1=[0]*B*B
st2=[0]*B
for v,i in q:
  st1[i]+=1
  st2[i//B]+=1
  ok=i
  ng=-1
  while abs(ok-ng)>1:
    m=(ok+ng)//2
    g=0
    l=m
    r=i
    yl=l//B
    yr=r//B
    if yl==yr:
      g+=sum(st1[l:r+1])
    else:
      g+=sum(st1[l:yl*B+B])
      g+=sum(st2[yl+1:yr])
      g+=sum(st1[yr*B:r+1])
    if g==1:
      ok=m
    else:
      ng=m
  L=ok
  ok=i
  ng=n
  while abs(ok-ng)>1:
    m=(ok+ng)//2
    g=0
    l=i
    r=m
    yl=l//B
    yr=r//B
    if yl==yr:
      g+=sum(st1[l:r+1])
    else:
      g+=sum(st1[l:yl*B+B])
      g+=sum(st2[yl+1:yr])
      g+=sum(st1[yr*B:r+1])
    if g==1:
      ok=m
    else:
      ng=m
  R=ok
  p[i]+=[(L,i,v)]
  p[R+1]+=[(L,i,pow(v,M-2,M))]
ans=0
st1=[1]*B*B
st2=[B]*B
lt=[1]*B
for y in range(n):
  for l,r,v in p[y]:
    yl=l//B
    yr=r//B
    for i in range(yl*B,yl*B+B):
      st1[i]*=lt[yl]
      st1[i]%=M
    lt[yl]=1
    for i in range(yr*B,yr*B+B):
      st1[i]*=lt[yr]
      st1[i]%=M
    lt[yr]=1
    if yl==yr:
      for i in range(l,r+1):
        st2[yl]-=st1[i]
        st1[i]*=v
        st1[i]%=M
        st2[yl]+=st1[i]
        st2[yl]%=M
    else:
      for i in range(l,yl*B+B):
        st2[yl]-=st1[i]
        st1[i]*=v
        st1[i]%=M
        st2[yl]+=st1[i]
        st2[yl]%=M
      for i in range(yl+1,yr):
        st2[i]*=v
        st2[i]%=M
        lt[i]=v
      for i in range(yr*B,r+1):
        st2[yr]-=st1[i]
        st1[i]*=v
        st1[i]%=M
        st2[yr]+=st1[i]
        st2[yr]%=M
  g=0
  l=0
  r=y
  yl=l//B
  yr=r//B
  for i in range(yl*B,yl*B+B):
    st1[i]*=lt[yl]
    st1[i]%=M
  lt[yl]=1
  for i in range(yr*B,yr*B+B):
    st1[i]*=lt[yr]
    st1[i]%=M
  lt[yr]=1
  if yl==yr:
    g+=sum(st1[l:r+1])
  else:
    g+=sum(st1[l:yl*B+B])
    g+=sum(st2[yl+1:yr])
    g+=sum(st1[yr*B:r+1])
  ans+=g
  ans%=M
print(ans)
0