LOJ#6281. 數列分塊入門 5

題目描述

給出一個長為 nnn 的數列，以及 nnn 個操作，操作涉及區間開方，區間求和。

輸入格式

第一行輸入一個數字 nnn。

第二行輸入 nnn 個數字，第 i 個數字為 aia_ia​i​​，以空格隔開。

接下來輸入 nnn 行詢問，每行輸入四個數字 opt\mathrm{opt}opt、lll、rrr、ccc，以空格隔開。

若 opt=0\mathrm{opt} = 0opt=0，表示將位於 [l,r][l, r][l,r] 的之間的數字都開方。

若 opt=1\mathrm{opt} = 1opt=1，表示詢問位於 [l,r][l, r][l,r] 的所有數字的和。

輸出格式

對於每次詢問，輸出一行一個數字表示答案。

樣例

樣例輸入

4
1 2 2 3
0 1 3 1
1 1 4 4
0 1 2 2
1 1 2 4

樣例輸出

6
2

數據範圍與提示

對於 100% 100\%100% 的數據，1≤n≤50000,−231≤others 1 \leq n \leq 50000, -2^{31} \leq \mathrm{others}1≤n≤50000,−2​31​​≤others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans≤2​31​​−1。

這道題的難點在於如何維護開根這個神奇的操作

我自己測的是1e7的數差不多開五六次根就會變成1，所以我們直接維護整個塊內的數是否變成了1就可以了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define int long long 
using namespace std;
const int MAXN=1e5+10;
const int INF=1e8+10;
inline char nc()
{
    static char buf[MAXN],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXN,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    char c=nc();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=nc();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=nc();}
    return x*f;
}
int N;
int a[MAXN],block,L[MAXN],R[MAXN],belong[MAXN],sum[MAXN],flag[MAXN];

void Sqrt(int l,int r)
{
	for(int i=l;i<=min(r,R[l]);i++)
	{
		sum[belong[i]]-=a[i];
		a[i]=sqrt(a[i]);
		sum[belong[i]]+=a[i];
	}
	if(belong[l]!=belong[r])
		for(int i=L[r];i<=r;i++)
			sum[belong[i]]-=a[i],a[i]=sqrt(a[i]),sum[belong[i]]+=a[i];
	for(int i=belong[l]+1;i<=belong[r]-1;i++)
	{
		if(flag[i]) {continue;}
		flag[i]=1;
		for(int j=L[i*block];j<=R[i*block];j++)
		{
			sum[i]-=a[j];
			a[j]=sqrt(a[j]);
			sum[i]+=a[j];
			if(a[j]>1) flag[i]=0;
		}
	}
}
int Query(int l,int r)
{
	int ans=0;
	for(int i=l;i<=min(r,R[l]);i++)
		ans+=a[i];
	if(belong[l]!=belong[r])
		for(int i=L[r];i<=r;i++)
			ans+=a[i];
	for(int i=belong[l]+1;i<=belong[r]-1;i++)
		ans+=sum[i];
	return ans;
}
main()
{
    #ifdef WIN32
    freopen("a.in","r",stdin);
   // freopen("b.out","w",stdout);
    #else
    #endif
	N=read();block=sqrt(N);
	for(int i=1;i<=N;i++) a[i]=read();
	for(int i=1;i<=N;i++) belong[i]=(i-1)/block+1,L[i]=(belong[i]-1)*block+1,R[i]=belong[i]*block;
	for(int i=1;i<=N;i++) sum[belong[i]]+=a[i];
	for(int i=1;i<=N;i++)
	{
		int opt=read(),l=read(),r=read(),c=read();
		if(opt==0) 
			Sqrt(l,r); 
		else  
			printf("%d\n",Query(l,r));
	}
    return 0;
}