BZOJ 3038: 上帝造题的七分钟2

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，谢谢！

（第一遍没看见，wrong。。。。：（）GJ太强

 

ORZ lyz大佬！

查询是树状数组的基本操作，所以不再谈

因为sqrt可以是数据衰减特别快，6次就可以把1e12变成1

所以最多每个数可以被修改6次，总共最多修改MAXN*6次

所以可以暴力修改

修改时用并查集找出向右第一个大于1的数进行修改

O(n lg n)

细节详见代码

#include<cstdio>
#include<iostream>
#include<cmath>
#define ll long long
using namespace std;

inline ll read()
{
	ll ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n;
const int N=1e5+5;
int f[N];
ll a[N],c[N];

inline void add(int x,ll k)
{
	for(int i=x;i<=n;i+=i&-i) c[i]+=k;
}

inline ll sum(int x)
{
	ll ret=0;
	for(int i=x;i;i-=i&-i) ret+=c[i];
	return ret; 
}

inline int find(int x)
{
	if(f[x]!=x) f[x]=find(f[x]);
	return f[x];
}

int main()
{
	n=read();
	for(int i=1;i<=n;i++)
		a[i]=read(),add(i,a[i]),f[i]=i;
	int m=read();
	while(m--)
	{
		int t=read(),l=read(),r=read();
		if (l>r) swap(l,r);
		if(!t)
			for(int i=find(l);(i<=r)&&(i!=0);i=find(i+1))
			 //没有i！=,当i==find（n+1）时 i=0;就会死循环，样例都过不了 
			{
				
				ll k=(int)sqrt(a[i]);
				add(i,k-a[i]);
				a[i]=k;
				if(a[i]<=1) f[i]=find(i+1);
			} else printf("%lld\n",sum(r)-sum(l-1));
	}
	return 0;
}