数列分块入门 7(LibreOj-6283)
生活随笔
收集整理的這篇文章主要介紹了
数列分块入门 7(LibreOj-6283)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
【題目描述】
給出一個長為 n 的數列,以及 n 個操作,操作涉及區間乘法,區間加法,單點詢問。
【輸入格式】
第一行輸入一個數字 n。
第二行輸入 n?個數字,第 i 個數字為 ai,以空格隔開。
接下來輸入 n 行詢問,每行輸入四個數字 opt、l、r、c,以空格隔開。
若 opt=0,表示將位于 [l,r] 的之間的數字都加 c。
若 opt=1,表示將位于 [l,r] 的之間的數字都乘 c。
若 opt=2,表示詢問 ar 的值 mod 10007(l 和 c 忽略)。
【輸出格式】
對于每次詢問,輸出一行一個數字表示答案。
【樣例】
樣例輸入
7
1 2 2 3 9 3 2
0 1 3 1
2 1 3 1
1 1 4 4
0 1 7 2
1 2 6 4
1 1 6 5
2 2 6 4
樣例輸出
3
100
【數據范圍與提示】
對于 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 = 10000+7; const int N = 1000000+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];//pos記錄第i個元素在第幾個塊中 int mulTag[N];//乘法標記 int addTag[N];//加法標記 void init(int n){block=sqrt(n);//塊的長度sum=n/block;//塊個數if(n%block)sum++;for(int i=1;i<=n;i++)//第i個元素在第幾塊中pos[i]=(i-1)/block+1;for(int i=1;i<=sum;i++)//乘法標記初始化為1mulTag[i]=1; } void reset(int x){//更新不完整塊的值for(int i=(pos[x]-1)*block+1;i<=min(pos[x]*block,n);i++)//更新元素a[i]=(mulTag[pos[x]]*a[i]+addTag[pos[i]])%MOD;//更新標記mulTag[pos[x]]=1;addTag[pos[x]]=0; } void update(int L,int R,int x,int op){reset(L);//更新左端不完整塊for(int i=L;i<=min(pos[L]*block,R);i++){//左邊的邊角料if(op==0)a[i]=(a[i]+x)%MOD;elsea[i]=(a[i]*x)%MOD;}if(pos[L]!=pos[R]){//存在右區間才遍歷,防止重復計算reset(R);//更新右端不完整塊for(int i=(pos[R]-1)*block+1;i<=R;i++){//右邊的邊角料if(op==0)a[i]=(a[i]+x)%MOD;elsea[i]=(a[i]*x)%MOD;}}for(int i=pos[L]+1;i<=pos[R]-1;i++){//中間的整塊if(op==0)addTag[i]=(addTag[i]+x)%MOD;else{addTag[i]=(addTag[i]*x)%MOD;mulTag[i]=(mulTag[i]*x)%MOD;}} } int query(int x){return (a[x]*mulTag[pos[x]]+addTag[pos[x]])%MOD; } int main(){scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%d",&a[i]);a[i]%=MOD;}init(n);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||op==1)//加或乘update(left,right,x,op);else if(op==2)printf("%d\n",query(right));}return 0; }?
總結
以上是生活随笔為你收集整理的数列分块入门 7(LibreOj-6283)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 图论 —— 网络流 —— 最大流 ——
- 下一篇: C++语言基础 —— STL —— 算法