Time Limit: 3 Sec Memory Limit: 128 MB

Description

XLk覺得《上帝造題的七分鐘》不太過癮，于是有了第二部。
"第一分鐘，X說，要有數列，于是便給定了一個正整數數列。
第二分鐘，L說，要能修改，于是便有了對一段數中每個數都開平方(下取整)的操作。
第三分鐘，k說，要能查詢，于是便有了求一段數的和的操作。
第四分鐘，彩虹喵說，要是noip難度，于是便有了數據范圍。
第五分鐘，詩人說，要有韻律，于是便有了時間限制和內存限制。
第六分鐘，和雪說，要省點事，于是便有了保證運算過程中及最終結果均不超過64位有符號整數類型的表示范圍的限制。
第七分鐘，這道題終于造完了，然而，造題的神牛們再也不想寫這道題的程序了。"
——《上帝造題的七分鐘·第二部》
所以這個神圣的任務就交給你了。

Input

第一行一個整數n，代表數列中數的個數。
第二行n個正整數，表示初始狀態下數列中的數。
第三行一個整數m，表示有m次操作。
接下來m行每行三個整數k,l,r，k=0表示給[l,r]中的每個數開平方(下取整)，k=1表示詢問[l,r]中各個數的和。

Output

對于詢問操作，每行輸出一個回答。

Sample Input

10 1 2 3 4 5 6 7 8 9 10 5 0 1 10 1 1 10 1 1 5 0 5 8 1 4 8

Sample Output

19 7 6

HINT

1：對于100%的數據，1<=n<=100000，1<=l<=r<=n，數列中的數大于0，且不超過1e12。

2：數據不保證L<=R 若L>R，請自行交換L,R，謝謝！

Solution

顯然的，1e12的數可以在最多10次以內就可以被開為1，所以到了線段樹一個區間，如果這個區間全是1，不管是詢問還是修改都直接return回去。不然直接在線段樹里面暴力修改。

#include<cstdio> #include<algorithm> #include<cmath> #include<cstring> #define lowbit(x) ((x)&(-(x))) #define lc o<<1 #define rc o<<1|1 using namespace std; #define REP(i,a,n) for(register int i=(a);i<=(n);++i) #define PER(i,a,n) for(register int i=(a);i>=(n);--i) #define FEC(i,x) for(register int i=head[x];i;i=g[i].ne) template<typename A>inline void read(A&a){a=0;char c=0;int f=1;while(c<'0'||c>'9')(c=getchar())=='-'?f=-1:0;while(c>='0'&&c<='9')a=(a<<3)+(a<<1)+c-'0',c=getchar();f==-1?a=-a:0;} char buf[30];template<typename A>inline void write(A a){if(a<0)putchar('-'),a=-a;int top=0;if(!a)buf[top=1]='0';while(a)buf[++top]=a%10+'0',a/=10;while(top)putchar(buf[top--]);} typedef long long ll;typedef unsigned long long ull; template<typename A,typename B>inline bool SMAX(A&x,const B&y){return y>x?x=y,1:0;} template<typename A,typename B>inline bool SMIN(A&x,const B&y){return y<x?x=y,1:0;}const int N=100000+7; int n,Q,opt,x,y;ll a[N];//錯誤筆記：似乎a數組就要開ll了 struct Node{ll maxv;}t[N<<2]; inline void Build(int o,int L,int R){if(L==R)t[o].maxv=a[L];else{int M=(L+R)>>1;Build(lc,L,M);Build(rc,M+1,R);t[o].maxv=max(t[lc].maxv,t[rc].maxv);}} inline void Update(int o,int L,int R,int l,int r){if(t[o].maxv<=1)return;if(L==R)return (void)(t[o].maxv=sqrt(t[o].maxv));int M=(L+R)>>1;if(l<=M)Update(lc,L,M,l,r);if(r>M)Update(rc,M+1,R,l,r);t[o].maxv=max(t[lc].maxv,t[rc].maxv); } inline ll Sum(int o,int L,int R,int l,int r){if(t[o].maxv<=1)return (min(R,r)-max(L,l)+1)*t[o].maxv;if(L==R)return t[o].maxv;int M=(L+R)>>1;ll ans=0;if(l<=M)ans+=Sum(lc,L,M,l,r);if(r>M)ans+=Sum(rc,M+1,R,l,r);return ans; }int main(){read(n);for(register int i=1;i<=n;++i)read(a[i]);Build(1,1,n);read(Q);for(register int i=1;i<=Q;++i){read(opt),read(x),read(y);if(x>y)x^=y^=x^=y;if(opt)write(Sum(1,1,n,x,y)),putchar('\n');else Update(1,1,n,x,y);} }

轉載于:https://www.cnblogs.com/hankeke/p/BZOJ3038.html

總結

以上是生活随笔為你收集整理的BZOJ3038 上帝造题的七分钟的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。