bzoj3211 花神游历各国 线段树,势能分析

时间:2023-05-15 17:09:08

【bzoj3211】花神游历各国

2014年3月17日2,7230

Description

 bzoj3211 花神游历各国 线段树,势能分析

Input

 bzoj3211 花神游历各国 线段树,势能分析

Output

每次x=1时,每行一个整数,表示这次旅行的开心度

Sample Input

4

1 100 5 5

5

1 1 2

2 1 2

1 1 2

2 2 3

1 1 4

Sample Output

101

11

11

HINT

对于100%的数据, n ≤ 100000,m≤200000 ,data[i]非负且小于10^9

题解

  这道题目记录,因为一个数经过较少的次数就会被开根号到1;

  所以势能分析即可,

  最坏情况是O(log * 10^9 m)

 #include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<cstdio> #define N 100007
#define ll long long
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch>''||ch<''){if (ch=='-') f=-;ch=getchar();}
while(ch<=''&&ch>=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
return x*f;
}
int n,m;
ll a[N];
struct data
{
int l,r;
ll sum;
bool flag;
}tr[N<<]; inline void update(int p)
{
tr[p].sum=tr[p<<].sum+tr[p<<|].sum;
tr[p].flag=tr[p<<].flag&tr[p<<|].flag;
}
void build(int p,int l,int r)
{
tr[p].l=l;tr[p].r=r;
if(l==r)
{
tr[p].sum=a[l];
if(a[l]==||a[l]==)tr[p].flag=;
return;
}
int mid=(l+r)>>;
build(p<<,l,mid),build(p<<|,mid+,r);
update(p);
}
void modify(int p,int x,int y)
{
if(tr[p].flag)return;
int l=tr[p].l,r=tr[p].r;
if(l==r)
{
tr[p].sum=(ll)sqrt(tr[p].sum);
if(tr[p].sum==||tr[p].sum==)tr[p].flag=;
return;
}
int mid=(l+r)>>;
if(mid>=y) modify(p<<,x,y);
else if(mid<x) modify(p<<|,x,y);
else modify(p<<,x,mid),modify(p<<|,mid+,y);
update(p);
}
ll query(int p,int x,int y)
{
int l=tr[p].l,r=tr[p].r;
if(l==x&&r==y)return tr[p].sum;
int mid=(l+r)>>;
if(mid>=y)return query(p<<,x,y);
else if(mid<x)return query(p<<|,x,y);
else return query(p<<,x,mid)+query(p<<|,mid+,y);
}
int main()
{
n=read();
for(int i=;i<=n;i++)
a[i]=read();
build(,,n);
m=read();
for(int i=;i<=m;i++)
{
int k=read(),x=read(),y=read();
if(x>y)swap(x,y);
if(k== )modify(,x,y);
else printf("%lld\n",query(,x,y));
}
}