結果

問題 No.1625 三角形の質問
ユーザー chocoruskchocorusk
提出日時 2021-07-23 23:25:27
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2,602 ms / 6,000 ms
コード長 5,356 bytes
コンパイル時間 4,023 ms
コンパイル使用メモリ 204,392 KB
実行使用メモリ 194,100 KB
最終ジャッジ日時 2023-09-26 02:00:34
合計ジャッジ時間 40,506 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 122 ms
19,608 KB
testcase_02 AC 904 ms
82,464 KB
testcase_03 AC 991 ms
85,668 KB
testcase_04 AC 701 ms
72,828 KB
testcase_05 AC 1,499 ms
107,500 KB
testcase_06 AC 2,542 ms
182,540 KB
testcase_07 AC 2,521 ms
182,748 KB
testcase_08 AC 2,546 ms
183,232 KB
testcase_09 AC 2,546 ms
182,736 KB
testcase_10 AC 2,532 ms
182,144 KB
testcase_11 AC 2,520 ms
183,060 KB
testcase_12 AC 2,602 ms
183,436 KB
testcase_13 AC 2,586 ms
183,872 KB
testcase_14 AC 2,587 ms
182,536 KB
testcase_15 AC 2,551 ms
183,332 KB
testcase_16 AC 582 ms
107,880 KB
testcase_17 AC 808 ms
170,496 KB
testcase_18 AC 1,331 ms
105,476 KB
testcase_19 AC 999 ms
194,100 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#include <list>
#include <atcoder/all>
#define popcount __builtin_popcount
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef pair<ll, ll> P;
template<typename Monoid>
struct SegmentTree{
	using F=function<Monoid(Monoid, Monoid)>;
	int sz;
	vector<Monoid> seg;
	F f;
	Monoid e;

    SegmentTree(){}

	SegmentTree(int n, const F f, const Monoid &e): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
	}

	SegmentTree(int n, const F f, const Monoid &e, vector<Monoid> v): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
		for(int i=0; i<n; i++) seg[i+sz]=v[i];
		for(int i=sz-1; i>=1; i--){
			seg[i]=f(seg[2*i], seg[2*i+1]);
		}
	}

    void init(int n, const F f1, const Monoid &e1){
        f=f1, e=e1;
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
	}

	void update(int k, const Monoid &x){
		k+=sz;
		seg[k]=x;
		while(k>1){
			k>>=1;
			seg[k]=f(seg[2*k], seg[2*k+1]);
		}
	}

	Monoid query(int a, int b){
		a+=sz, b+=sz;
		Monoid ret=e;
		for(;a<b; a>>=1, b>>=1){
			if(b&1) ret=f(ret, seg[--b]);
			if(a&1) ret=f(ret, seg[a++]);
		}
		return ret;
	}

	/*Monoid query(int a, int b){ //演算が非可換のとき
		a+=sz, b+=sz;
		Monoid retl=e, retr=e;
		for(;a<b; a>>=1, b>>=1){
			if(b&1) retr=f(seg[--b], retr);
			if(a&1) retl=f(retl, seg[a++]);
		}
		return f(retl, retr);
	}*/

	Monoid operator[](const int &k) const{
		return seg[k+sz];
	}
};
int l[100010], r[100010]; ll s[100010];
int t[100010], lq[100010], rq[100010]; ll sq[100010];
int main()
{
    int n, q; cin>>n>>q;
    

    for(int i=0; i<n; i++){
        ll a, b, c, d, e, f;
        cin>>a>>b>>c>>d>>e>>f;
        l[i]=min({a, c, e}), r[i]=max({a, c, e});
        c-=a, e-=a, d-=b, f-=b;
        s[i]=abs(c*f-d*e);
    }
    for(int i=0; i<q; i++){
        cin>>t[i];
        if(t[i]==1){
            ll a, b, c, d, e, f;
            cin>>a>>b>>c>>d>>e>>f;
            lq[i]=min({a, c, e}), rq[i]=max({a, c, e});
            c-=a, e-=a, d-=b, f-=b;
            sq[i]=abs(c*f-d*e);
        }else{
            cin>>lq[i]>>rq[i];
            rq[i]++;
        }
    }
    vector<int> vl, vr;
    for(int i=0; i<n; i++){
        vl.push_back(l[i]);
        //vr.push_back(r[i]);
    }
    for(int i=0; i<q; i++){
        vl.push_back(lq[i]);
        //vr.push_back(rq[i]);
    }
    sort(vl.begin(), vl.end());
    vl.erase(unique(vl.begin(), vl.end()), vl.end());
    // sort(vr.begin(), vr.end());
    // vr.erase(unique(vr.begin(), vr.end()), vr.end());
    int m=vl.size();
    int sz=1;
    while(sz<m) sz<<=1;
    vector<vector<P>> seg(2*sz);
    
    for(int i=0; i<n; i++){
        int k=lower_bound(vl.begin(), vl.end(), l[i])-vl.begin();
        k+=sz;
        seg[k].push_back({r[i], s[i]});
        while(k>1){
            k>>=1;
            seg[k].push_back({r[i], s[i]});
        }
    }
    for(int i=0; i<q; i++){
        if(t[i]==2) continue;
        int k=lower_bound(vl.begin(), vl.end(), lq[i])-vl.begin();
        k+=sz;
        seg[k].push_back({rq[i], sq[i]});
        while(k>1){
            k>>=1;
            seg[k].push_back({rq[i], sq[i]});
        }
    }
    for(int i=2*sz-1; i>=0; i--) sort(seg[i].begin(), seg[i].end());

    
    vector<SegmentTree<ll>> segm(2*sz);
    for(int i=1; i<2*sz; i++){
        segm[i].init((int)seg[i].size(), [](ll x, ll y){ return max(x, y);}, -1ll);
    }
    for(int i=0; i<n; i++){
        int k=lower_bound(vl.begin(), vl.end(), l[i])-vl.begin();
        k+=sz;
        int u=lower_bound(seg[k].begin(), seg[k].end(), P(r[i], s[i]))-seg[k].begin();
        segm[k].update(u, s[i]);
        while(k>1){
            k>>=1;
            int u=lower_bound(seg[k].begin(), seg[k].end(), P(r[i], s[i]))-seg[k].begin();
            segm[k].update(u, s[i]);
        }
    }
    for(int i=0; i<q; i++){
        if(t[i]==1){
            int k=lower_bound(vl.begin(), vl.end(), lq[i])-vl.begin();
            k+=sz;
            int u=lower_bound(seg[k].begin(), seg[k].end(), P(rq[i], sq[i]))-seg[k].begin();
            segm[k].update(u, sq[i]);
            while(k>1){
                k>>=1;
                int u=lower_bound(seg[k].begin(), seg[k].end(), P(rq[i], sq[i]))-seg[k].begin();
                segm[k].update(u, sq[i]);
            }
        }else{
            int l1=lower_bound(vl.begin(), vl.end(), lq[i])-vl.begin();
            int a=l1+sz, b=2*sz;
            ll mx=-1;
            for(;a<b; a>>=1, b>>=1){
                if(b&1){
                    b--;
                    int u=lower_bound(seg[b].begin(), seg[b].end(), P(rq[i], 0))-seg[b].begin();
                    mx=max(mx, segm[b].query(0, u));
                }
                if(a&1){
                    int u=lower_bound(seg[a].begin(), seg[a].end(), P(rq[i], 0))-seg[a].begin();
                    mx=max(mx, segm[a].query(0, u));
                    a++;
                }
            }
            cout<<mx<<endl;
        }
    }
    return 0;
}
0