結果

問題 No.2950 Max Min Product
ユーザー sasa8uyauya
提出日時 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other TLE * 1 -- * 36
権限があれば一括ダウンロードができます

ソースコード

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