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