結果

問題 No.1074 増殖
ユーザー leaf_1415leaf_1415
提出日時 2020-06-05 22:17:19
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 460 ms / 2,000 ms
コード長 3,125 bytes
コンパイル時間 876 ms
コンパイル使用メモリ 86,540 KB
実行使用メモリ 12,988 KB
最終ジャッジ日時 2024-12-17 15:56:12
合計ジャッジ時間 4,042 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))
#define mod 1000000007

using namespace std;
typedef pair<llint, llint> P;

// range-add, range-min query
 
struct LazySegTree{
	int size;
	vector<llint> seg, delay;
	
	LazySegTree(){}
	LazySegTree(int size){
		this->size = size;
		seg.resize(1<<(size+1));
		delay.resize(1<<(size+1));
	}
	
	void init()
	{
		for(int i = 0; i < (1<<(size+1)); i++){
			seg[i] = 0;
			delay[i] = -1;
		}
	}
	
	void eval(int l, int r, int k)
	{
		if(delay[k] != -1){
			seg[k] = delay[k] * (r-l+1);  //総和クエリのときは区間幅を乗じる必要がある
			if(l < r){
				delay[k*2] = delay[k];
				delay[k*2+1] = delay[k];
			}
			delay[k] = -1;
		}
	}
	
	void update(int i, llint val)
	{
		int l = 0, r = (1<<size)-1, k = 1;
		eval(l, r, k);
		for(int j = size-1; j >= 0; j--){
			k <<= 1;
			if(i & (1<<j)){
				k++;
				r = (l+r)/2;
			}
			else l = (l+r)/2+1;
			eval(l, r, k);
		}
		
		i += (1 << size);
		seg[i] = val;
		while(i > 1){
			i /= 2;
			seg[i] = (seg[i*2] + seg[i*2+1]);
		}
	}
	
	void add(int a, int b, int k, int l, int r, llint val)
	{
		eval(l, r, k);
		
		if(b < l || r < a) return;
		if(a <= l && r <= b){
			delay[k] = val;
			eval(l, r, k);
			return;
		}
		add(a, b, k*2, l, (l+r)/2, val);
		add(a, b, k*2+1, (l+r)/2+1, r, val);
		seg[k] = (seg[k*2] + seg[k*2+1]);
	}
	void add(int a, int b, llint val){
		if(a > b) return;
		add(a, b, 1, 0, (1<<size)-1, val);
	}
 
	llint query(int a, int b, int k, int l, int r)
	{
		eval(l, r, k);
		
		if(b < l || r < a) return 0;
		if(a <= l && r <= b) return seg[k];
		llint lval = query(a, b, k*2, l, (l+r)/2);
		llint rval = query(a, b, k*2+1, (l+r)/2+1, r);
		return (lval + rval);
	}
	llint query(int a, int b)
	{
		return query(a, b, 1, 0, (1<<size)-1);
	}
};

llint n;
llint lx[40005], ly[40005], rx[40005], ry[40005];
llint ans[40005];
LazySegTree seg(18);

void calc(llint x[], llint y[])
{
	seg.init();
	seg.update(0, inf);
	
	for(int i = 1; i <= n; i++){
		llint X = x[i], Y = y[i];
		if(seg.query(X, X) >= Y) continue;
		
		llint ub = 20005, lb = 0, mid;
		while(ub-lb>1){
			mid = (ub+lb)/2;
			if(seg.query(mid, mid) >= Y) lb = mid;
			else ub = mid;
		}
		llint sum = Y * lb;
		sum += seg.query(lb+1, X);
		ans[i] += X*Y-sum;
		seg.add(lb+1, X, Y);
		//cout << X << " " << Y << " " << sum << endl;
	}
}

int main(void)
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	
	cin >> n;
	for(int i = 1; i <= n; i++){
		cin >> lx[i] >> ly[i] >> rx[i] >> ry[i];
		lx[i] *= -1, ly[i] *= -1;
	}
	
	calc(rx, ry);
	calc(rx, ly);
	calc(lx, ly);
	calc(lx, ry);
	
	for(int i = 1; i <= n; i++) cout << ans[i] << endl;
	
	return 0;
}
0