数列分块入门 3(LibreOj-6279)
生活随笔
收集整理的這篇文章主要介紹了
数列分块入门 3(LibreOj-6279)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
【題目描述】
給出一個長為 n 的數列,以及 n 個操作,操作涉及區間加法,詢問區間內小于某個值 x 的前驅(比其小的最大元素)。
【輸入格式】
第一行輸入一個數字 n。
第二行輸入 n?個數字,第 i 個數字為 ai,以空格隔開。
接下來輸入 n 行詢問,每行輸入四個數字 opt、l、r、c,以空格隔開。
若 opt=0,表示將位于 [l,r] 的之間的數字都加 c。
若 opt=1,表示詢問位于 [l,r] 中 c 的前驅的值,不存在輸出 -1。
【輸出格式】
對于每次詢問,輸出一行一個數字表示答案。
【樣例】
樣例輸入
4
1 2 2 3
0 1 3 1
1 1 4 4
0 1 2 2
1 1 2 4
樣例輸出
3
-1
【數據范圍與提示】
對于 100%?的數據,1<=n<=100000,-2^31<=other,ans<=2^31-1。
【源代碼】
#include<iostream> #include<cstdio> #include<cstdlib> #include<string> #include<cstring> #include<cmath> #include<ctime> #include<algorithm> #include<utility> #include<stack> #include<queue> #include<vector> #include<set> #include<map> #include<bitset> #define EPS 1e-9 #define PI acos(-1.0) #define INF 0x3f3f3f3f #define LL long long const int MOD = 1E9+7; const int N = 100000+5; const int dx[] = {-1,1,0,0,-1,-1,1,1}; const int dy[] = {0,0,-1,1,-1,1,-1,1}; using namespace std;int n; int block,sum;//block為塊的長度,sum為塊的個數 int a[N];//存放數列元素 int pos[N],tag[N];//pos記錄第i個元素在第幾個塊中,tag為操作標記 int ans[N];//維護整塊和 vector<int> v[N]; void init(){block=sqrt(n);//塊的長度sum=n/block;//塊個數if(n%block)sum++;for(int i=1;i<=n;i++){pos[i]=(i-1)/block+1;//第i個元素在第幾塊中v[pos[i]].push_back(a[i]);//保存每個數分塊的序號}for(int i=1;i<=sum;i++)//對整塊進行排序sort(v[i].begin(),v[i].end()); } void resort(int x){v[pos[x]].clear();for(int i=(pos[x]-1)*block+1;i<=min(pos[x]*block,n);i++)v[pos[x]].push_back(a[i]);sort(v[pos[x]].begin(),v[pos[x]].end()); } void update(int L,int R,int x){for(int i=L;i<=min(pos[L]*block,R);i++)//左邊的邊角料a[i]+=x;resort(L);//對不完整塊排序if(pos[L]!=pos[R]){//存在右區間才遍歷,防止重復計算for(int i=(pos[R]-1)*block+1;i<=R;i++)//右邊的邊角料a[i]+=x;resort(R);//對不完整塊排序}for(int i=pos[L]+1;i<=pos[R]-1;i++)//中間的整塊tag[i]+=x; } void query(int L,int R,int x){int res=-INF;bool flag=false;for(int i=L;i<=min(pos[L]*block,R);i++){//左邊的邊角料if(a[i]+tag[pos[L]]<x){res=max(a[i]+tag[pos[L]],res);flag=true;}}if(pos[L]!=pos[R]){//存在右區間才遍歷,防止重復計算for(int i=(pos[R]-1)*block+1;i<=R;i++){//右邊的邊角料if(a[i]+tag[pos[R]]<x){res=max(a[i]+tag[pos[R]],res);flag=true;}}}for(int i=pos[L]+1;i<=pos[R]-1;i++){//中間的整塊進行二分查找int temp=x-tag[i];int k=lower_bound(v[i].begin(),v[i].end(),temp)-v[i].begin();if(k!=0&&v[i][k-1]+tag[i]<x){res=max(res,v[i][k-1]+tag[i]);flag=true;}}if(!flag)printf("-1\n");elseprintf("%d\n",res); } int main(){scanf("%d",&n);for(int i=1;i<=n;i++)scanf("%d",&a[i]);init();for(int i=1;i<=n;i++){int op;int left,right,x;scanf("%d",&op);scanf("%d%d%d",&left,&right,&x);if(op==0)update(left,right,x);else if(op==1)query(left,right,x);}return 0; }?
總結
以上是生活随笔為你收集整理的数列分块入门 3(LibreOj-6279)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 理论基础 —— 栈
- 下一篇: 基础算法 —— 调度问题 —— 流水调度