首页 > 其他好文 > 详细

bzoj3211: 花神游历各国

时间：2016-09-17 23:34:30 阅读：188 评论：0 收藏：0 [点我收藏+]

标签：

/*向下取整smg！
Popoqqq：题目大意：给定一个序列，提供下列操作：
1.将[l.r]区间内每个数a[i]变为sqrt(a[i])
2.查询[l,r]区间的和
根号是不支持区间修改的，于是我们选择单点修改区间查询的树状数组，但是这样是O(n^2)的，怎么办？
我们发现一个数x最多开loglogx次根号就会变为1 也就是一个int范围内的数只要开6次根号就会变为1 于是修改的总时间复杂度为O(nloglogn)
但是单次修改怎么办？我们维护一个并查集，一旦一个数为1或0，我们就把这个位置的father设为它右面的那个位置即可这样可以均摊O(n)时间找到一个数后面第一个>1的数
*/

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#include<cmath>
using namespace std;
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define clr(x,c) memset(x,c,sizeof(x))
#define ll long long
int read(){
	int x=0,f=1;char c=getchar();
	while(!isdigit(c)){
		if(c==‘-‘) f=-1;c=getchar();
	}
	while(isdigit(c)) x=x*10+c-‘0‘,c=getchar();
	return x*f;
}
const int nmax=1e5+5;
ll sm[nmax];int fa[nmax],a[nmax];
int find(int x){
	return fa[x]==x?x:fa[x]=find(fa[x])
;}
ll query(int x){
	ll ans=0;
	for(int i=x;i;i-=(i&-i)) ans+=sm[i];
	return ans;
}
void update(int x,int n,ll add){
	for(int i=x;i<=n;i+=(i&-i)) sm[i]+=add;
}
int main(){
	int n=read(),u,v,d;
	rep(i,1,n) a[i]=read(),update(i,n,a[i]),fa[i]=i;fa[n+1]=n+1;
	int m=read();
	rep(i,1,m){
		u=read(),v=read(),d=read();
		if(u==1) printf("%lld\n",query(d)-query(v-1));
		else{
			int j=v;ll t;
			while(j<=d){
				j=find(j);if(j>d) break;t=a[j];a[j]=(ll)sqrt(t);update(j,n,a[j]-t);
				if(a[j]<=1) fa[j]=find(j+1);++j;
			}
		}
	}
	return 0;
}