結果

問題 No.1074 増殖
ユーザー convexineqconvexineq
提出日時 2020-06-07 03:54:33
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 3,612 bytes
コンパイル時間 302 ms
コンパイル使用メモリ 87,052 KB
実行使用メモリ 111,712 KB
最終ジャッジ日時 2023-08-25 03:04:26
合計ジャッジ時間 6,156 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 70 ms
72,748 KB
testcase_01 AC 70 ms
72,472 KB
testcase_02 AC 72 ms
72,744 KB
testcase_03 AC 71 ms
72,624 KB
testcase_04 AC 74 ms
77,116 KB
testcase_05 AC 74 ms
77,120 KB
testcase_06 WA -
testcase_07 AC 467 ms
107,812 KB
testcase_08 AC 570 ms
111,712 KB
testcase_09 AC 709 ms
103,580 KB
testcase_10 AC 339 ms
99,652 KB
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.keys = set()
        self.B = BIT(n+1) #存在すれば 1、しないなら 0
        self.dic = [0]*(n+1) # 値域
    
    def insert(self,a,b): # 値 a に b を上書き
        if a not in self.keys:
            self.B.add(a,1)
            self.keys.add(a)
        self.dic[a] = b

    def remove(self,a): # a を取り除く
        self.keys.remove(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,n):
    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]
        x0,y0 =x,y
        k = mp.lower_bound(x)
        v1 = mp[k]
        v2 = mp[mp.next_key(k)]
        if y <= v2:
            continue
        val = (y-v2)*(x-k)
        while y > v1:
            x = k
            v2 = v1
            k1 = mp.prev_key(k)
            #print(k,k1)
            mp.remove(k)
            k = k1
            v1 = mp[k]
            val += (y-v2)*(x-k)
            #print(x,y,k,k1,v2,val)
        mp.insert(x0,y0)

        res[i] = val        
    return res


ans = [0]*n

#a = [(1,10),(3,5),(3,6),(4,10)]
#print(solve(a,4))

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