結果

問題 No.1625 三角形の質問
ユーザー chocorusk
提出日時 2021-07-23 23:25:27
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2,872 ms / 6,000 ms
コード長 5,356 bytes
コンパイル時間 4,702 ms
コンパイル使用メモリ 200,344 KB
最終ジャッジ日時 2025-01-23 09:02:58
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

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