結果
| 問題 |
No.1705 Mode of long array
|
| コンテスト | |
| ユーザー |
moririn2528_c
|
| 提出日時 | 2021-10-08 22:57:21 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 175 ms / 3,000 ms |
| コード長 | 5,335 bytes |
| コンパイル時間 | 1,370 ms |
| コンパイル使用メモリ | 114,760 KB |
| 実行使用メモリ | 11,300 KB |
| 最終ジャッジ日時 | 2024-07-23 06:10:45 |
| 合計ジャッジ時間 | 7,048 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 51 |
ソースコード
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<vector>
#include<cmath>
#include<algorithm>
#include<map>
#include<queue>
#include<deque>
#include<iomanip>
#include<tuple>
#include<cassert>
#include<set>
#include<complex>
#include<numeric>
#include<functional>
#include<unordered_map>
#include<unordered_set>
using namespace std;
typedef long long int LL;
typedef pair<int,int> P;
typedef pair<LL,LL> LP;
const int INF=1<<30;
const LL MAX=1e9+7;
void array_show(int *array,int array_n,char middle=' '){
for(int i=0;i<array_n;i++)printf("%d%c",array[i],(i!=array_n-1?middle:'\n'));
}
void array_show(LL *array,int array_n,char middle=' '){
for(int i=0;i<array_n;i++)printf("%lld%c",array[i],(i!=array_n-1?middle:'\n'));
}
void array_show(vector<int> &vec_s,int vec_n=-1,char middle=' '){
if(vec_n==-1)vec_n=vec_s.size();
for(int i=0;i<vec_n;i++)printf("%d%c",vec_s[i],(i!=vec_n-1?middle:'\n'));
}
void array_show(vector<LL> &vec_s,int vec_n=-1,char middle=' '){
if(vec_n==-1)vec_n=vec_s.size();
for(int i=0;i<vec_n;i++)printf("%lld%c",vec_s[i],(i!=vec_n-1?middle:'\n'));
}
template<typename T> ostream& operator<<(ostream& os,const vector<T>& v1){
int n=v1.size();
for(int i=0;i<n;i++){
if(i)os<<" ";
os<<v1[i];
}
return os;
}
template<typename T1,typename T2> ostream& operator<<(ostream& os,const pair<T1,T2>& p){
os<<p.first<<" "<<p.second;
return os;
}
template<typename T> istream& operator>>(istream& is,vector<T>& v1){
int n=v1.size();
for(int i=0;i<n;i++)is>>v1[i];
return is;
}
template<typename T1,typename T2> istream& operator>>(istream& is,pair<T1,T2>& p){
is>>p.first>>p.second;
return is;
}
template <typename T> class seg_tree{
//monoid
private:
function<T(T,T)> func;
T e;
int N,n=-1;
vector<T> seg;
void init(){
assert(n>=0);
int i;
for(i=0;(1<<i)<n;i++);
N=(1<<i)-1;
seg.assign(2*(N+1),e);
}
void init_reload(){
for(int i=N-1;i>=0;i--){
seg[i]=func(seg[2*i+1],seg[2*i+2]);
}
}
void update(int pos){
T a=func(seg[pos*2+1],seg[pos*2+2]);
if(seg[pos]==a)return;
seg[pos]=a;
if(pos==0)return;
update((pos-1)/2);
}
public:
seg_tree(function<T(T,T)> _func,T _e,int _n):func(_func),e(_e),n(_n){
init();
}
seg_tree(function<T(T,T)> _func,T _e,vector<T> vec):func(_func),e(_e){
n=vec.size();
init();
for(int i=0;i<n;i++){
seg[N+i]=vec[i];
}
init_reload();
}
seg_tree(function<T(T,T)> _func,T _e,int _n,T a):func(_func),e(_e),n(_n){
init(e);
for(int i=0;i<n;i++){
seg[N+i]=a;
}
init_reload();
}
void set(int pos,T a){
assert(pos>=0 && pos<=N);
pos+=N;
seg[pos]=a;
update((pos-1)/2);
}
T search(int a,int b,int l,int r,int x){//[a,b) search
if(a<=l && r<=b)return seg[x];
int m=(l+r)/2;
if(b<=m)return search(a,b,l,m,2*x+1);
if(m<=a)return search(a,b,m,r,2*x+2);
return func(search(a,m,l,m,2*x+1),search(m,b,m,r,2*x+2));
}
T search(int a,int b){
assert(a<b);
return search(a,b,0,N+1,0);
}
int max_right(function<bool(T)>& g,int pos,int l,int r,int x,T& y){
//suppose that S is return value, g(func(pos,..,S-1))=true,g(func(pos,..,S))=false
if(pos<=l && g(func(y,seg[x]))){
y=func(y,seg[x]);
return r;
}
if(l+1==r)return l;
int m=(l+r)/2;
if(pos<m){
int s=max_right(g,pos,l,m,2*x+1,y);
if(s<m)return s;
}
return max_right(g,pos,m,r,2*x+2,y);
}
int max_right(function<bool(T)> g,int pos){
T y=e;
int s=max_right(g,pos,0,N+1,0,y);
return min(s,n);
}
int min_left(function<bool(T)>& g,int pos,int l,int r,int x,T& y){
//suppose that S is return value, g(func(S,..,pos-1))=true,g(func(S-1,..,pos-1))=false
int s;
if(r<=pos && g(func(seg[x],y))){
y=func(seg[x],y);
return l;
}
if(l+1==r)return r;
int m=(l+r)/2;
if(m<pos){
s=min_left(g,pos,m,r,2*x+2,y);
if(m<s)return s;
}
return min_left(g,pos,l,m,2*x+1,y);
}
int min_left(function<bool(T)> g,int pos){
assert(pos>=0);
if(pos==0)return 0;
T y=e;
return min_left(g,pos,0,N+1,0,y);
}
};
namespace sol{
LP op(LP a,LP b){
if(a<b)return b;
else return a;
}
void solve(){
LL n,m,q;
int i,j,k;
LL a,b,c;
cin>>n>>m;
vector<LL> v1(m);
cin>>v1;
vector<LP> va(m);
for(i=0;i<m;i++){
va[i]={v1[i],i};
}
seg_tree<LP> seg(op,{-1,-1},va);
cin>>q;
for(i=0;i<q;i++){
cin>>a>>b>>c;
b--;
if(a==1)v1[b]+=c;
if(a==2)v1[b]-=c;
if(a<=2){
seg.set(b,{v1[b],b});
continue;
}
cout<<seg.search(0,m).second+1<<endl;
}
}
}
int main(){
cin.tie(0);
ios::sync_with_stdio(false);
sol::solve();
}
moririn2528_c