結果
問題 | No.1099 Range Square Sum |
ユーザー | vwxyz |
提出日時 | 2021-09-01 11:08:40 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,389 ms / 2,000 ms |
コード長 | 6,871 bytes |
コンパイル時間 | 1,650 ms |
コンパイル使用メモリ | 82,192 KB |
実行使用メモリ | 235,552 KB |
最終ジャッジ日時 | 2024-05-05 13:30:38 |
合計ジャッジ時間 | 15,704 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 131 ms
89,128 KB |
testcase_01 | AC | 126 ms
89,104 KB |
testcase_02 | AC | 132 ms
89,336 KB |
testcase_03 | AC | 133 ms
88,900 KB |
testcase_04 | AC | 133 ms
89,156 KB |
testcase_05 | AC | 135 ms
88,868 KB |
testcase_06 | AC | 130 ms
89,112 KB |
testcase_07 | AC | 125 ms
88,972 KB |
testcase_08 | AC | 128 ms
89,052 KB |
testcase_09 | AC | 123 ms
89,080 KB |
testcase_10 | AC | 124 ms
89,052 KB |
testcase_11 | AC | 188 ms
92,168 KB |
testcase_12 | AC | 182 ms
91,776 KB |
testcase_13 | AC | 185 ms
92,296 KB |
testcase_14 | AC | 192 ms
91,596 KB |
testcase_15 | AC | 194 ms
91,580 KB |
testcase_16 | AC | 191 ms
91,860 KB |
testcase_17 | AC | 184 ms
91,732 KB |
testcase_18 | AC | 181 ms
91,688 KB |
testcase_19 | AC | 181 ms
91,744 KB |
testcase_20 | AC | 185 ms
91,660 KB |
testcase_21 | AC | 1,377 ms
235,552 KB |
testcase_22 | AC | 1,389 ms
233,012 KB |
testcase_23 | AC | 1,310 ms
233,604 KB |
testcase_24 | AC | 1,385 ms
234,740 KB |
testcase_25 | AC | 1,363 ms
232,304 KB |
testcase_26 | AC | 714 ms
154,956 KB |
testcase_27 | AC | 707 ms
154,428 KB |
testcase_28 | AC | 729 ms
155,192 KB |
testcase_29 | AC | 758 ms
155,048 KB |
testcase_30 | AC | 718 ms
153,824 KB |
ソースコード
import bisect import copy import decimal import fractions import functools import heapq import itertools import math import random import sys from collections import Counter,deque,defaultdict from functools import lru_cache,reduce from heapq import heappush,heappop,heapify,heappushpop,_heappop_max,_heapify_max def _heappush_max(heap,item): heap.append(item) heapq._siftdown_max(heap, 0, len(heap)-1) def _heappushpop_max(heap, item): if heap and item < heap[0]: item, heap[0] = heap[0], item heapq._siftup_max(heap, 0) return item from math import gcd as GCD read=sys.stdin.read readline=sys.stdin.readline readlines=sys.stdin.readlines class Lazy_Segment_Tree: def __init__(self,N,f,e,f_act,e_act,operate): self.N=N self.f=f self.e=e self.f_act=f_act self.e_act=e_act self.operate=operate self.segment_tree=[self.e]*(self.N+self.N) self.segment_tree_act=[self.e_act]*(self.N+self.N) def __getitem__(self,i): if type(i) is int: if -self.N<=i<0: i+=self.N*2 elif 0<=i<self.N: i+=self.N else: raise IndexError('list index out of range') self.Propagate_Above(i) self.Recalculate_Above(i) return self.Operate_At(i) else: a,b,c=i.start,i.stop,i.step if a==None or a<-self.N: a=self.N elif self.N<=a: a=self.N*2 elif a<0: a+=self.N*2 else: a+=self.N if b==None or self.N<=b: b=self.N*2 elif b<-self.N: b=self.N elif b<0: b+=self.N*2 else: b+=self.N return self.segment_tree[slice(a,b,c)] def __setitem__(self,i,x): if -self.N<=i<0: i+=self.N*2 elif 0<=i<self.N: i+=self.N else: raise IndexError('list index out of range') self.Propagate_Above(i) self.segment_tree[i]=x self.segment_tree_act[i]=self.e_act self.Recalculate_Above(i) def Operate_At(self,i): return self.operate(self.segment_tree[i],self.segment_tree_act[i]) def Propagate_At(self,i): self.segment_tree[i]=self.Operate_At(i) self.segment_tree_act[i<<1]=self.f_act(self.segment_tree_act[i<<1],self.segment_tree_act[i]) self.segment_tree_act[i<<1|1]=self.f_act(self.segment_tree_act[i<<1|1],self.segment_tree_act[i]) self.segment_tree_act[i]=self.e_act def Propagate_Above(self,i): H=i.bit_length()-1 for h in range(H,0,-1): self.Propagate_At(i>>h) def Recalculate_Above(self,i): while i>1: i>>=1 self.segment_tree[i]=self.f(self.Operate_At(i<<1),self.Operate_At(i<<1|1)) def Build(self,lst): for i,x in enumerate(lst,self.N): self.segment_tree[i]=x for i in range(self.N-1,0,-1): self.segment_tree[i]=self.f(self.segment_tree[i<<1],self.segment_tree[i<<1|1]) self.segment_tree_act=[self.e_act]*(self.N+self.N) def Fold(self,L=None,R=None): if L==None or L<-self.N: L=self.N elif self.N<=L: L=self.N*2 elif L<0: L+=self.N*2 else: L+=self.N if R==None or self.N<=R: R=self.N*2 elif R<-self.N: R=self.N elif R<0: R+=self.N*2 else: R+=self.N self.Propagate_Above(L//(L&-L)) self.Propagate_Above(R//(R&-R)-1) vL=self.e vR=self.e while L<R: if L&1: vL=self.f(vL,self.Operate_At(L)) L+=1 if R&1: R-=1 vR=self.f(self.Operate_At(R),vR) L>>=1 R>>=1 return self.f(vL,vR) def Fold_Index(self,L=None,R=None): if L==None or L<-self.N: L=self.N elif self.N<=L: L=self.N*2 elif L<0: L+=self.N*2 else: L+=self.N if R==None or self.N<=R: R=self.N*2 elif R<-self.N: R=self.N elif R<0: R+=self.N*2 else: R+=self.N if L==R: return None x=self.Fold(L-self.N,R-self.N) while L<R: if L&1: if self.segment_tree[L]==x: i=L break L+=1 if R&1: R-=1 if self.segment_tree[R]==x: i=R break L>>=1 R>>=1 while i<self.N: if self.segment_tree[i]==self.segment_tree[i<<1]: i<<=1 else: i<<=1 i|=1 i-=self.N return i def Operate_Range(self,a,L=None,R=None): if L==None or L<-self.N: L=self.N elif self.N<=L: L=self.N*2 elif L<0: L+=self.N*2 else: L+=self.N if R==None or self.N<=R: R=self.N*2 elif R<-self.N: R=self.N elif R<0: R+=self.N*2 else: R+=self.N L0=L//(L&-L) R0=R//(R&-R)-1 self.Propagate_Above(L0) self.Propagate_Above(R0) while L<R: if L&1: self.segment_tree_act[L]=self.f_act(self.segment_tree_act[L],a) L+=1 if R&1: R-=1 self.segment_tree_act[R]=self.f_act(self.segment_tree_act[R],a) L>>=1 R>>=1 self.Recalculate_Above(L0) self.Recalculate_Above(R0) def Update(self): for i in range(1,self.N): self.Propagate_At(i) for i in range(self.N,self.N*2): self.segment_tree[i]=self.Operate_At(i) self.segment_tree_act[i]=self.e_act for i in range(self.N-1,0,-1): self.segment_tree[i]=self.f(self.segment_tree[i<<1],self.segment_tree[i<<1|1]) N=int(readline()) def f(tpl0,tpl1): x0,y0,z0=tpl0 x1,y1,z1=tpl1 return (x0+x1,y0+y1,z0+z1) e=(0,0,1) def f_act(a,b): return a+b e_act=0 def operate(tpl,a): x,y,z=tpl return (x+z*a,y+2*x*a+z*a**2,z) LST=Lazy_Segment_Tree(N,f,e,f_act,e_act,operate) A=list(map(int,readline().split())) LST.Build([(a,a**2,1) for a in A]) Q=int(readline()) for _ in range(Q): query=tuple(map(int,readline().split())) if query[0]==1: _,l,r,x=query l-=1;r-=1 LST.Operate_Range(x,l,r+1) else: _,l,r=query l-=1;r-=1 ans=LST.Fold(l,r+1)[1] print(ans)