結果

問題 No.1074 増殖
ユーザー convexineqconvexineq
提出日時 2020-06-07 03:37:46
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 3,347 bytes
コンパイル時間 242 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 102,640 KB
最終ジャッジ日時 2024-06-06 03:53:34
合計ジャッジ時間 5,189 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
54,144 KB
testcase_01 AC 43 ms
54,400 KB
testcase_02 AC 40 ms
54,016 KB
testcase_03 AC 41 ms
53,504 KB
testcase_04 AC 44 ms
59,520 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 AC 419 ms
98,928 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 RE -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

class BIT: #0-indexed
    def __init__(self, n):
        self.size = n
        self.tree = [0]*(n+1)
        self.depth = n.bit_length()
        self.n0 = 1<<self.depth
#        self.element = [0]*(n+1)
    def get_sum(self, i): #a_0 + ... + a_{i} #閉区間
        s = 0; i += 1
        while i > 0:
            s += self.tree[i]
            i -= i & -i
        return s
    def query(self,l,r): #a_l + ... + a_r 閉区間
        return self.get_sum(r) - self.get_sum(l-1) 
    def add(self, i, x):
        i += 1
        while i <= self.size:
            self.tree[i] += x
            i += i & -i
        # self.element[i] += x
    #def get(self,i): return element[i]        
    def bisect_left(self,w):
        #和が w 以上になる最小の index
        #w が存在しない場合 self.size を返す
        if w <= 0: return 0
        x,k = 0,self.n0
        for _ in range(self.depth):
            k >>= 1
            if x+k <= self.size and self.tree[x+k] < w:
                w -= self.tree[x+k]
                x += k
        return x
    
"""
C++ の std::map ぽいデータ構造
座標圧縮は前提
n: 取りうる最大値（最小値は 0）

self.dic を dict に取り換えることでメモリ節約できる
"""
class stdmap:
    def __init__(self, n):
        self.size = n+1
        self.B = BIT(n+1) #存在すれば 1、しないなら 0
        self.dic = [0]*(n+1) # 値域
    
    def insert(self,a,b): # 値 a に b を上書き
        self.B.add(a,1)
        self.dic[a] = b

    def remove(self,a): # a を取り除く
        self.B.add(a,-1)

    def lower_bound(self,k): # k 以上の最小のkeyを求める
        return self.B.bisect_left(self.B.get_sum(k))
        
    def kth_key(self,k): # k 番目に小さい元のkeyを求める
        return self.B.bisect_left(k)

    def kth_value(self,k): # k 番目に小さい元のmap先を求める
        return self.dic[self.B.bisect_left(k)]

    def prev_key(self,k): #一個前の元のkeyを求める
        idx = self.B.get_sum(k)
        assert idx != 0
        return self.B.bisect_left(idx-1)

    def next_key(self,k):
        idx = self.B.get_sum(k)
        assert idx != self.size
        return self.B.bisect_left(idx+1)

    def __getitem__(self,item):
        return self.dic[item]


# coding: utf-8
# Your code here!
import sys
read = sys.stdin.read
readline = sys.stdin.readline

n = int(input())
M = map(int,read().split())
xyxy = [*zip(M,M,M,M)]


def solve(xy):
    MAX = 2*10**4+10
    mp = stdmap(MAX)
    mp.insert(MAX,0)
    mp.insert(0,1<<31)

    res = [0]*n
    for i in range(n):
        x,y = xy[i]
        k = mp.lower_bound(x)
        v1 = mp[k]
        v2 = mp[mp.next_key(k)]
        if y <= v2:
            continue
        val = (y-v2)*(x-k)
        mp.insert(x,y)
        while y > v1:
            x = k
            v2 = v1
            k1 = mp.prev_key(k)
            mp.remove(k)
            k = k1
            v1 = mp[k]
            val += (y-v2)*(x-k)

        res[i] = val        
    return res


ans = [0]*n
res1 = solve([(x2,y2) for x1,y1,x2,y2 in xyxy])
res2 = solve([(-x1,y2) for x1,y1,x2,y2 in xyxy])
res3 = solve([(x2,-y1) for x1,y1,x2,y2 in xyxy])
res4 = solve([(-x1,-y1) for x1,y1,x2,y2 in xyxy])
for i in range(n):
    ans[i] += res1[i]+res2[i]+res3[i]+res4[i]
print("\n".join(map(str,ans)))
0