loj#6281. 数列分块入门 5

Posted CHerish_OI

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了loj#6281. 数列分块入门 5相关的知识,希望对你有一定的参考价值。

#6281. 数列分块入门 5

题目:传送门 

简要题意:

   给出一个长为 n 的数列,以及 n 个操作,操作涉及区间开方,区间求和。

 


 

 

题解:

   怎么说...这道题loj的数据有点水。

   和bzoj花神游历各国是一样的...但是loj没有卡掉不完全优化的代码。

   基础操作就不说了(同分块4),主要讲优化吧:

   不难发现,因为开方后我们是向下取整的,那么如果一直开方下去,整个序列一定会变成0或1。那么如果当前区间都是0或1就没有操作的必要了嘛。

   讲到这里思路就已经很清晰了,但是一开始的时候本蒟蒻并没有优化头尾的块,依然在loj AC了...bzoj就挂了,事实证明:暴力加优化才是分块真正的难点!

 


 

 

代码:

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cstdlib>
 4 #include<cmath>
 5 #include<algorithm>
 6 using namespace std;
 7 typedef long long LL;
 8 LL n,m,a[110000],sum[110000];
 9 LL block,pos[110000];
10 bool v[110000];
11 void reset(LL x)
12 {
13     if(v[x])return ;
14     v[x]=1;sum[x]=0;
15     for(int i=(x-1)*block+1;i<=x*block;i++)
16     {
17         a[i]=sqrt(a[i]);sum[x]+=a[i];
18         if(a[i]>1)v[x]=0;
19     }
20 }
21 void update(LL l,LL r)
22 {
23     if(!v[pos[l]])
24     {
25         for(int i=l;i<=min(pos[l]*block,r);i++)
26         {
27             sum[pos[l]]-=a[i];
28             a[i]=sqrt(a[i]);
29             sum[pos[l]]+=a[i];
30         }
31         v[pos[l]]=1;for(int i=(pos[l]-1)*block+1;i<=pos[l]*block;i++)if(a[i]>1)v[pos[l]]=0;
32     }
33     if(pos[l]!=pos[r])
34     {
35         if(!v[pos[r]])
36         {
37             for(int i=(pos[r]-1)*block+1;i<=r;i++)
38             {
39                 sum[pos[r]]-=a[i];
40                 a[i]=sqrt(a[i]);
41                 sum[pos[r]]+=a[i];
42             }
43             v[pos[r]]=1;for(int i=(pos[r]-1)*block;i<=pos[r]*block;i++)if(a[i]>1)v[pos[r]]=0;
44         }
45     }
46     for(int i=pos[l]+1;i<=pos[r]-1;i++)reset(i);
47 }
48 void sol(LL l,LL r)
49 {
50     LL ans=0;
51     for(int i=l;i<=min(pos[l]*block,r);i++)ans+=a[i];
52     if(pos[l]!=pos[r])
53         for(int i=(pos[r]-1)*block+1;i<=r;i++)
54             ans+=a[i];
55     for(int i=pos[l]+1;i<=pos[r]-1;i++)ans+=sum[i];
56     printf("%lld\n",ans);
57 }
58 int main()
59 {
60     memset(v,false,sizeof(v));
61     scanf("%lld",&n);
62     for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
63     block=sqrt(n);
64     for(int i=1;i<=n;i++)
65     {
66         pos[i]=(i-1)/block+1;
67         sum[pos[i]]+=a[i];
68     }
69     for(int i=1;i<=n;i++)
70     {
71         LL opt,l,r,c;scanf("%lld%lld%lld%lld",&opt,&l,&r,&c);
72         if(opt==0)update(l,r);
73         else sol(l,r);
74     }
75     return 0;
76 }

 

以上是关于loj#6281. 数列分块入门 5的主要内容,如果未能解决你的问题,请参考以下文章

LOJ.6281.数列分块入门5(分块 区间开方)

LOJ#6281. 数列分块入门 5

#6281. 数列分块入门 5

#6280. 数列分块入门 4 #6281. 数列分块入门 5

Loj 6285. 数列分块入门 9

Loj 6282. 数列分块入门 6