結果

問題 No.3122 Median of Medians of Division
ユーザー 沙耶花
提出日時 2025-04-18 23:51:37
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,255 bytes
コンパイル時間 4,569 ms
コンパイル使用メモリ 256,456 KB
実行使用メモリ 442,164 KB
最終ジャッジ日時 2025-04-18 23:51:51
合計ジャッジ時間 13,426 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 1
other TLE * 1 -- * 39
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:18:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   18 |         rep(i,n)scanf("%d",&a[i]);
      |                 ~~~~~^~~~~~~~~~~~
main.cpp:22:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   22 |                 scanf("%d %d %d",&ts[i],&ls[i],&rs[i]);
      |                 ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 4000000000000000001LL
int mm;
int xs[320][400005];

int main(){
	
	int n,q;
	cin>>n>>q;
	vector<int> a(n);
	rep(i,n)scanf("%d",&a[i]);
	vector<int> ts(q),ls(q),rs(q);
	vector<int> vs = a;
	rep(i,q){
		scanf("%d %d %d",&ts[i],&ls[i],&rs[i]);
		ls[i]--;
		if(ts[i]==1)vs.push_back(rs[i]);
	}
	sort(vs.begin(),vs.end());
	vs.erase(unique(vs.begin(),vs.end()),vs.end());
	mm = vs.size();
	rep(i,n){
		a[i] = lower_bound(vs.begin(),vs.end(),a[i])-vs.begin();
	}
	rep(i,q){
		if(ts[i]==1)rs[i] = lower_bound(vs.begin(),vs.end(),rs[i])-vs.begin();
	}
	int B = 635;
	vector<vector<int>> as(n/B+1,vector<int>());
	rep(i,n){
		as[i/B].push_back(a[i]);
	}
	rep(i,as.size()){
		rep(j,as[i].size()){
			xs[i][as[i][j]]+=2;
			if(j+1!=as[i].size()){
				xs[i][max(as[i][j],as[i][j+1])]--;
			}
		}
		rep(j,mm){
			xs[i][j+1] += xs[i][j];
		}
	}
	rep(_,q){
		int t = ts[_];
		int l = ls[_],r = rs[_];
		if(t==1){
			as[l/B][l%B] = r;
			int ii = l/B;
			rep(j,mm){
				xs[ii][j] = 0;
			}
			rep(j,as[ii].size()){
				xs[ii][as[ii][j]]+=2;
				if(j+1!=as[ii].size()){
					xs[ii][max(as[ii][j],as[ii][j+1])]++;
				}
			}
			rep(j,mm){
				xs[ii][j+1] += xs[ii][j];
			}
		}
		else{
			int ok = 0,ng = mm;
			while(ng-ok>1){
				int mid = (ok+ng)/2;
				int tl = l,tr = r;
				int x = 0;
				while(tl%B!=0 && tl<tr){
					if(as[tl/B][tl%B]<mid)x+=2;
					if(tl+1 < tr){
						if(max(as[tl/B][tl%B], as[(tl+1)/B][(tl+1)%B])<mid)x--;
					}
					tl++;
				}
				while(tl + B < tr){
					x += xs[tl/B][mid-1];//distance(xs[tl/B].begin(),lower_bound(xs[tl/B].begin(),xs[tl/B].end(),mid));
					//y += ys[tl/B][mid-1];//distance(ys[tl/B].begin(),lower_bound(ys[tl/B].begin(),ys[tl/B].end(),mid));
					tl += B-1;
					if(max(as[tl/B][tl%B], as[(tl+1)/B][(tl+1)%B])<mid)x--;
					tl++;
				}
				while(tl<tr){
					if (as[tl/B][tl%B]<mid)x+=2;
					if(tl+1 < tr){
						if(max(as[tl/B][tl%B], as[(tl+1)/B][(tl+1)%B])<mid)x--;
					}
					tl++;
				}
				if(r-l -x  > 0)ok = mid;
				else ng= mid;
			}
			printf("%d\n",ok);
		}
		
	}
	return 0;
}
0