結果

問題 No.878 Range High-Element Query
ユーザー furafura
提出日時 2020-08-07 02:33:05
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 251 ms / 2,000 ms
コード長 2,597 bytes
コンパイル時間 2,693 ms
コンパイル使用メモリ 216,616 KB
実行使用メモリ 23,216 KB
最終ジャッジ日時 2023-10-22 06:24:41
合計ジャッジ時間 6,194 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 3 ms
4,348 KB
testcase_02 AC 3 ms
4,348 KB
testcase_03 AC 3 ms
4,348 KB
testcase_04 AC 3 ms
4,348 KB
testcase_05 AC 3 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 3 ms
4,348 KB
testcase_09 AC 3 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 213 ms
20,776 KB
testcase_12 AC 182 ms
21,744 KB
testcase_13 AC 166 ms
14,804 KB
testcase_14 AC 127 ms
13,736 KB
testcase_15 AC 177 ms
21,056 KB
testcase_16 AC 236 ms
22,912 KB
testcase_17 AC 251 ms
23,216 KB
testcase_18 AC 250 ms
23,216 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;

template<class M>
class segment_tree{
	int sz;
	vector<M> seg;
public:
	segment_tree(){}
	segment_tree(int n){ build(n); }
	template<class T>
	segment_tree(const vector<T>& a){ build(a); }
	void build(int n){
		for(sz=1;sz<n;sz<<=1);
		seg.assign(2*sz,M());
	}
	template<class T>
	void build(const vector<T>& a){
		int n=a.size();
		build(n);
		rep(i,n) seg[sz+i]=a[i];
		for(int i=sz-1;i>0;i--){
			seg[i]=seg[2*i]*seg[2*i+1];
		}
	}
	const M& operator[](int i)const{
		return seg[sz+i];
	}
	void update(int i,const M& x){
		i+=sz;
		seg[i]=x;
		for(i>>=1;i>0;i>>=1) seg[i]=seg[2*i]*seg[2*i+1];
	}
	M get(int l,int r)const{
		M lcum,rcum;
		for(int a=l+sz,b=r+sz;a<b;a>>=1,b>>=1){
			if(a&1) lcum=lcum*seg[a++];
			if(b&1) rcum=seg[--b]*rcum;
		}
		return lcum*rcum;
	}
	template<class N,class F>
	N query(int l,int r,const F& f){
		N lcum,rcum;
		for(int a=l+sz,b=r+sz;a<b;a>>=1,b>>=1){
			if(a&1) lcum=lcum*f(seg[a++]);
			if(b&1) rcum=f(seg[--b])*rcum;
		}
		return lcum*rcum;
	}
};

template<class T>
class max_monoid{
	T a;
public:
	max_monoid():a(numeric_limits<T>::min()){}
	max_monoid(const T& val):a(val){}
	max_monoid operator*(const max_monoid& x)const{
		return max(a,x.a);
	}
	T& get(){ return a; }
	const T& get()const{ return a; }
};

template<class T>
class add_monoid{
	T a;
public:
	add_monoid():a(){}
	add_monoid(const T& val):a(val){}
	add_monoid operator*(const add_monoid& x)const{
		return a+x.a;
	}
	T& get(){ return a; }
	const T& get()const{ return a; }
};

template<class T>
class sequence_monoid{
	vector<T> a;
public:
	sequence_monoid(){}
	sequence_monoid(const T& val):a(1,val){}
	sequence_monoid(const vector<T>& x):a(x){
		sort(a.begin(),a.end());
	}
	sequence_monoid operator*(const sequence_monoid& x)const{
		vector<T> res;
		merge(a.begin(),a.end(),x.a.begin(),x.a.end(),back_inserter(res));
		return res;
	}
	int f(int t){
		return lower_bound(a.begin(),a.end(),t)-a.begin();
	}
};

const int INF=1<<29;

int main(){
	int n,q; scanf("%d%d",&n,&q);
	vector<int> a(n+1);
	a[0]=INF;
	rep(i,n) scanf("%d",&a[i+1]);

	vector<int> pos(n);
	segment_tree<max_monoid<int>> S1(a);
	rep(i,n){
		int lo=0,hi=i+1;
		while(hi-lo>1){
			int mi=(lo+hi)/2;
			if(S1.get(mi,i+1).get()>a[i+1]) lo=mi; else hi=mi;
		}
		pos[i]=lo-1;
	}

	segment_tree<sequence_monoid<int>> S2(pos);
	rep(_,q){
		int l,r; scanf("%*d%d%d",&l,&r); l--;
		auto ans=S2.query<add_monoid<int>>(l,r,[&](sequence_monoid<int>& x)->add_monoid<int>{
			return x.f(l);
		});
		printf("%d\n",ans.get());
	}

	return 0;
}
0