結果

問題 No.1919 Many Monster Battles
コンテスト
ユーザー googol_S0
提出日時 2022-04-29 22:55:57
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,931 bytes
コンパイル時間 209 ms
コンパイル使用メモリ 82,460 KB
実行使用メモリ 274,556 KB
最終ジャッジ日時 2024-06-29 04:33:39
合計ジャッジ時間 6,895 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 1 -- * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

from copy import *
def init(N,node,A,unit,func):
  n=1
  while n<N:
    n<<=1
  for i in range(n*2-1):
    if len(node)<=i:
      node.append(deepcopy(unit))
    else:
      node[i]=deepcopy(unit)
  for i in range(len(A)):
    node[n-1+i]=deepcopy(A[i])
  for i in range(n-2,-1,-1):
    node[i]=func(node[(i<<1)+1],node[(i<<1)+2])
  node.append(func)
  node.append(unit)
  node.append(n)

def upd(node,x,a):
  y=node[-1]+x
  node[y-1]+=a
  while y>1:
    y=y>>1
    node[y-1]=node[-3](node[(y<<1)-1],node[y<<1])

def query(node,l,r):
  x,y=l,r
  z=node[-1]-1
  rr=deepcopy(node[-2])
  rl=deepcopy(node[-2])
  while True:
    if x==y:
      return node[-3](rl,rr)
    if x&1:
      rl=node[-3](rl,node[x+z])
      x+=1
    if y&1:
      rr=node[-3](node[y+z-1],rr)
    if z==0:
      return node[-3](rl,rr)
    x>>=1
    y>>=1
    z>>=1

mod=10**9+7
def comp(x):
  y=sorted(set(x))
  d=dict()
  for i in range(len(y)):
    d[y[i]]=i
  return d

def solve(A,B):
  N=len(A)
  X=[(A[i],B[i]) for i in range(N)]
  X.sort()
  Y=[X[i][0]-X[i][1] for i in range(N)]
  Z=[X[i][0]+X[i][1] for i in range(N)]
  y,z=comp(Y),comp(Z)
  ANS=0
  segy=[]
  segz=[]
  init(N,segy,[],0,lambda x,y:x+y)
  init(N,segz,[],0,lambda x,y:x+y)
  D=dict()
  for i in range(N):
    a,b=y[Y[i]],z[Z[i]]
    ANS=(ANS+X[i][0]*(query(segy,0,a)+query(segz,0,b)-i+D.get((X[i][0],X[i][1]),0)))%mod
    upd(segy,a,1)
    upd(segz,b,1)
    D[(X[i][0],X[i][1])]=D.get((X[i][0],X[i][1]),0)+1
  segy=[]
  segz=[]
  init(N,segy,[],0,lambda x,y:x+y)
  init(N,segz,[],0,lambda x,y:x+y)
  D=dict()
  for i in range(N-1,-1,-1):
    a,b=y[Y[i]],z[Z[i]]
    ANS=(ANS-X[i][0]*(query(segy,a+1,N)+query(segz,b+1,N)-(N-1-i)+D.get((X[i][0],X[i][1]),0)))%mod
    upd(segy,a,1)
    upd(segz,b,1)
    D[(X[i][0],X[i][1])]=D.get((X[i][0],X[i][1]),0)+1
  return ANS

N=int(input())
A=list(map(int,input().split()))
B=list(map(int,input().split()))
print(solve(A,B)*2%mod,solve(B,A)*2%mod)
0