結果
問題 | No.1099 Range Square Sum |
ユーザー |
![]() |
提出日時 | 2020-09-25 11:19:18 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 4,498 bytes |
コンパイル時間 | 150 ms |
コンパイル使用メモリ | 82,388 KB |
実行使用メモリ | 229,484 KB |
最終ジャッジ日時 | 2024-06-28 05:47:18 |
合計ジャッジ時間 | 8,083 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 2 WA * 18 TLE * 1 -- * 9 |
ソースコード
mod = 1000000007eps = 10**-9mask = (1 << 60) - 1def main():import sysinput = sys.stdin.buffer.readlinedef op(a, b):sqsum_a, sum_a = a >> 60, a & maskif sum_a > 10**12:sum_a -= 1 << 60sqsum_a += 1sqsum_b, sum_b = b >> 60, b & maskif sum_b > 10**12:sum_b -= 1 << 60sqsum_b += 1return ((sqsum_a + sqsum_b) << 60) + sum_a + sum_be = 0def mapping(a, x, seg_len):sqsum, sum_ = a >> 60, a & maskif sum_ > 10**12:sum_ -= 1 << 60sqsum += 1sqsum_new = sqsum + (seg_len * x) ** 2 + 2 * x * sum_sum_new = sum_ + seg_len * xreturn (sqsum_new << 60) + sum_newdef composition(x, y):return x + yid = 0class LazySegTree:# Range update querydef __init__(self, A, op=op, e=e, mapping=mapping, composition=composition, id=id, initialize=True):self.N = len(A)self.LV = (self.N - 1).bit_length()self.N0 = 1 << self.LVself.op = opself.e = eself.mapping = mappingself.composition = compositionself.id = idif initialize:self.data = [self.e] * self.N0 + A + [self.e] * (self.N0 - self.N)for i in range(self.N0 - 1, 0, -1):self.data[i] = op(self.data[i * 2], self.data[i * 2 + 1])else:self.data = [self.e] * (self.N0 * 2)self.lazy = [id] * (self.N0 * 2)def _ascend(self, i):for _ in range(i.bit_length() - 1):i >>= 1self.data[i] = self.op(self.data[i * 2], self.data[i * 2 + 1])def _descend(self, idx):lv = idx.bit_length()seg_len = 1 << self.LVfor j in range(lv - 1, 0, -1):seg_len >>= 1i = idx >> jx = self.lazy[i]if x == self.id:continueself.lazy[i * 2] = self.composition(self.lazy[i * 2], x)self.lazy[i * 2 + 1] = self.composition(self.lazy[i * 2 + 1], x)self.lazy[i] = self.idself.data[i * 2] = self.mapping(self.data[i * 2], x, seg_len)self.data[i * 2 + 1] = self.mapping(self.data[i * 2 + 1], x, seg_len)# open interval [l, r)def apply(self, l, r, x):l += self.N0 - 1r += self.N0 - 1self._descend(l // (l & -l))self._descend(r // (r & -r) - 1)l_ori = lr_ori = rseg_len = 1while l < r:if l & 1:self.data[l] = self.mapping(self.data[l], x, seg_len)self.lazy[l] = self.composition(self.lazy[l], x)l += 1if r & 1:r -= 1self.data[r] = self.mapping(self.data[r], x, seg_len)self.lazy[r] = self.composition(self.lazy[r], x)l >>= 1r >>= 1seg_len <<= 1self._ascend(l_ori // (l_ori & -l_ori))self._ascend(r_ori // (r_ori & -r_ori) - 1)# open interval [l, r)def query(self, l, r):l += self.N0 - 1r += self.N0 - 1self._descend(l // (l & -l))self._descend(r // (r & -r) - 1)ret_l = self.eret_r = self.ewhile l < r:if l & 1:ret_l = self.op(ret_l, self.data[l])l += 1if r & 1:ret_r = self.op(self.data[r - 1], ret_r)r -= 1l >>= 1r >>= 1return self.op(ret_l, ret_r)N = int(input())A = list(map(int, input().split()))B = [((a*a) << 60) + a for a in A]ST = LazySegTree(B)ans = []for _ in range(int(input())):query = list(map(int, input().split()))if query[0] == 1:l, r, x = query[1:]ST.apply(l, r+1, x)else:l, r = query[1:]val = ST.query(l, r + 1)sqsum, sum_ = val >> 60, val & maskif sum_ > 10**12:sqsum += 1ans.append(sqsum)print(*ans, sep="\n")if __name__ == '__main__':main()