백준
백준 14504 수열과 쿼리 18
콩순이냉장고
2025. 2. 10. 21:46
문제 URL : https://www.acmicpc.net/problem/14504
문제 접근법 :
제곱근 분할법을 공부하다가 문제를 풀게됐네요
제곱근 분할이란게 데이터를 루트n개로 분할해서 관리하는형태이더군요
그안에서 데이터를 정렬된 형태로 저장하여
upper_bound를 이용한 데이터개수 찾는 알고리즘으로
해결한다는게 참 인상깊었습니다.
이번건의 문제는
여러 블로그를
참고했고 코드는:
https://david0506.tistory.com/57
[ Query ] Sqrt Decomposition (제곱근 분할법)
Sqrt Decomposition 구간 쿼리를 Segment Tree를 이용해서 처리하면 시간복잡도가 $O(logN)$이다. 이는 트리의 깊이에 비례하는데, 각 노드의 자식 노드의 수를 밑으로 가지는 로그의 시간복잡도를 가지
david0506.tistory.com
이분의 블로그를 참고했습니다.
소스코드:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
//제곱근 분할법 sqrt Decomposition
int n,m;
vector<ll> arr;
vector<vector<ll>>bucket;
vector<tuple<ll,ll,ll,ll>> iv;
int sqrtN;
void input() {
cin >> n;
arr=vector<ll>(n+1);
for(int i =0;i<n;i++)cin>>arr[i];
cin>>m;
int a,b,c,d=0;
for(int i =0;i<m;i++){
cin>>a>>b>>c;
if(a==1)cin>>d;
iv.push_back({a,b,c,d});
}
}
void init(){
sqrtN = sqrt(n);
bucket = vector<vector<ll>>(n+1);
for(int i=0;i<n;i++){
bucket[i/sqrtN].push_back(arr[i]);
}
for(int i =0;i<=sqrtN;i++)sort(bucket[i].begin(),bucket[i].end());
}
void update(int x,ll val){
int idx = x/sqrtN;
int bidx = lower_bound(bucket[idx].begin(),bucket[idx].end(),arr[x])-bucket[idx].begin();
bucket[idx].erase(bucket[idx].begin()+bidx);
arr[x]=val;
bidx = lower_bound(bucket[idx].begin(),bucket[idx].end(),val)-bucket[idx].begin();
bucket[idx].insert(bucket[idx].begin()+bidx,val);
}
ll query(int l,int r,ll k){
ll ret = 0;
while(l%sqrtN!=0 && l<=r){
ret+=arr[l++]>k;
}
while((r+1)%sqrtN!=0&&l<=r){
ret+=arr[r--]>k;
}
while(l<=r){
int idx = l/sqrtN;
ret+=bucket[idx].end()-upper_bound(bucket[idx].begin(),bucket[idx].end(),k);
l+=sqrtN;
}
return ret;
}
void solve() {
init();
for(auto[a,b,c,d]:iv){
if(a==1){
cout<<query(b-1,c-1,d)<<"\n";
}
else{
update(b-1,c);
}
}
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
//freopen("input.txt", "r", stdin);
input();
solve();
}
궁금한점 혹은 모르는점
어떤질문이든 댓글은 언제나 환영입니다.