luogu P4145 上帝造题的七分钟2 / 花神游历各国 维护区间和&&区间开根号

因为开根号能使数字减小得非常快

所以开不了几次（6次?）很大的数就会变成1.....

所以我们可以维护区间最大值，若最大值>1,则继续递归子树，暴力修改叶节点，否则直接return

（好像也可以维护区间被开方的次数，但我不会。。。QAQ）

#include<cstdio>
#include<iostream>
#include<cmath>
#define int long long
#define R register int
#define ls (tr<<1)
#define rs (tr<<1|1)
using namespace std;
inline int g() {
    R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix; 
    do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
}
int n,m;
int mx[400010],sum[400010];
inline void build(int tr,int l,int r){
    if(l==r) {mx[tr]=sum[tr]=g(); return ;} 
    R md=(l+r)>>1;
    build(ls,l,md),build(rs,md+1,r);
    mx[tr]=max(mx[ls],mx[rs]);
    sum[tr]=sum[ls]+sum[rs];
}
inline void calc(int tr,int l,int r,int ll,int rr) {
    if(mx[tr]<=1) return ;
    if(l==r) {mx[tr]=sqrt(mx[tr]),sum[tr]=sqrt(sum[tr]); return ;}
    R md=(l+r)>>1;
    if(ll>md) calc(rs,md+1,r,ll,rr);
    else if(rr<md+1) calc(ls,l,md,ll,rr);
    else calc(ls,l,md,ll,md),calc(rs,md+1,r,md+1,rr);
    mx[tr]=max(mx[ls],mx[rs]);
    sum[tr]=sum[ls]+sum[rs];
}
inline int query(int tr,int l,int r,int ll,int rr) {
    if(l==ll&&r==rr) return sum[tr];
    R md=(l+r)>>1;
    if(ll>md) return query(rs,md+1,r,ll,rr);
    else if(rr<md+1) return query(ls,l,md,ll,rr);
    else return query(ls,l,md,ll,md)+query(rs,md+1,r,md+1,rr);
}
signed main() {
    n=g(); build(1,1,n); m=g();
    for(R i=1,k,l,r;i<=m;++i) {
        k=g(),l=g(),r=g();
        if(l>r) swap(l,r);
        if(k&1) printf("%lld
",query(1,1,n,l,r));
        else calc(1,1,n,l,r);
    }
}

upd 2019.06.15

可以用树状数组做，如果这个数已经为$1$就用并查集合并到上一个位置，对于没有被开成$1$的直接暴力开根。。

跑的飞快

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ull unsigned long long
#define ll long long
#define R register ll
#define pause (for(R i=1;i<=10000000000;++i))
#define OUT freopen("out.out","w",stdout);
using namespace std;
namespace Fread {
    static char B[1<<15],*S=B,*D=B;
    #define getchar() (S==D&&(D=(S=B)+fread(B,1,1<<15,stdin),S==D)?EOF:*S++)
    inline ll g() {
        R ret=0,fix=1; register char ch; whaile(!isdigit(ch=getchar())) fix=ch=='-'?-1:fix;
        do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
    } inline bool isempty(const char& ch) {return ch<=36||ch>=127;}
    inline void gs(char* s) {register char ch; while(isempty(ch=getchar())); do *s++=ch; while(!isempty(ch=getchar()));}
}using Fread::g; using Fread::gs;
const int N=100010; int n,m,fa[N];
ll c[N],h[N]; bool flg[N];
inline int lbt(int x) {return x&-x;}
inline void change(int pos) { R tmp=h[pos];
    h[pos]=sqrt(h[pos]); R d=tmp-h[pos];
    if(h[pos]==1) flg[pos]=true;
    for(;pos<=n;pos+=lbt(pos)) c[pos]-=d;
}
inline ll query(int pos) { R ret=0;
    for(;pos;pos-=lbt(pos)) ret+=c[pos]; return ret;
}
inline int getf(int x) {return (!flg[fa[x]])?fa[x]:fa[x]=getf(fa[x]);}
signed main() {
    n=g(); for(R i=1;i<=n;++i) h[i]=c[i]=g(); for(R i=1;i<=n;++i) fa[i]=i-1;
    for(R i=1;i<=n;++i) if(i+lbt(i)<=n) c[i+lbt(i)]+=c[i];
    m=g(); for(R i=1,x,l,r;i<=m;++i) { x=g(),l=g(),r=g(); l>r?swap(l,r):void(0);
        if(x&1) printf("%lld
",query(r)-query(l-1));
        else for(flg[r]?r=getf(r):0;r>=l;r=getf(r)) change(r);
    }
}

---

2019.04.11&&2019.06.15