BZOJ 3211 弗洛拉前往国家 树阵+并检查集合

标题效果:给定一个序列,它提供了以下操作:

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的数

此外此题增加读入优化能够快一倍 令人非常费解0.0

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 100100
using namespace std;
typedef long long ll;
int n,m,a[M],fa[M];
ll c[M];
inline int getc() {
static const int L = 1 << 15;
static char buf[L], *S = buf, *T = buf;
if (S == T) {
T = (S = buf) + fread(buf, 1, L, stdin);
if (S == T)
return EOF;
}
return *S++;
}
inline int getint() {
int c;
while(!isdigit(c = getc()) && c != '-');
bool sign = c == '-';
int tmp = sign ? 0 : c - '0';
while(isdigit(c = getc()))
tmp = (tmp << 1) + (tmp << 3) + c - '0';
return sign ? -tmp : tmp;
}
int Find(int x)
{
if(!fa[x]||fa[x]==x)
return fa[x]=x;
return fa[x]=Find(fa[x]);
}
void Update(int x,int y)
{
for(;x<=n;x+=x&-x)
c[x]+=y;
}
ll Get_Ans(int x)
{
ll re=0;
for(;x;x-=x&-x)
re+=c[x];
return re;
}
void Modify(int x,int y)
{
int i=x;
for( i=x ; i<=y ; i=Find(i+1) )
{
int temp=(int)sqrt(a[i]);
Update(i,temp-a[i]);
a[i]=temp;
if(a[i]<=1)
fa[i]=Find(i+1);
}
}
int main()
{
int i,p,x,y;
cin>>n;
for(i=1;i<=n;i++)
{
a[i]=getint();
if(a[i]==1||!a[i])
fa[i]=i+1;
Update(i,a[i]);
}
m=getint();
for(i=1;i<=m;i++)
{
p=getint();
x=getint();
y=getint();
if(p==1)
printf("%lld\n", Get_Ans(y)-Get_Ans(x-1) );
else
Modify(x,y);
}
}
上一篇:原生JS实现图片放大镜插件


下一篇:ACM/ICPC 之 Floyd练习六道(ZOJ2027-POJ2253-POJ2472-POJ1125-POJ1603-POJ2607)