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;
}

猜你喜欢

转载自blog.csdn.net/YYHS_WSF/article/details/82497286