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

所以开不了几次(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=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
}
int n,m;
int mx[],sum[];
inline void build(int tr,int l,int r){
if(l==r) {mx[tr]=sum[tr]=g(); return ;}
R md=(l+r)>>;
build(ls,l,md),build(rs,md+,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]<=) return ;
if(l==r) {mx[tr]=sqrt(mx[tr]),sum[tr]=sqrt(sum[tr]); return ;}
R md=(l+r)>>;
if(ll>md) calc(rs,md+,r,ll,rr);
else if(rr<md+) calc(ls,l,md,ll,rr);
else calc(ls,l,md,ll,md),calc(rs,md+,r,md+,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)>>;
if(ll>md) return query(rs,md+,r,ll,rr);
else if(rr<md+) return query(ls,l,md,ll,rr);
else return query(ls,l,md,ll,md)+query(rs,md+,r,md+,rr);
}
signed main() {
n=g(); build(,,n); m=g();
for(R i=,k,l,r;i<=m;++i) {
k=g(),l=g(),r=g();
if(l>r) swap(l,r);
if(k&) printf("%lld\n",query(,,n,l,r));
else calc(,,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[<<],*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=,fix=; register char ch; whaile(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
} inline bool isempty(const char& ch) {return ch<=||ch>=;}
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=; 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]==) flg[pos]=true;
for(;pos<=n;pos+=lbt(pos)) c[pos]-=d;
}
inline ll query(int pos) { R ret=;
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=;i<=n;++i) h[i]=c[i]=g(); for(R i=;i<=n;++i) fa[i]=i-;
for(R i=;i<=n;++i) if(i+lbt(i)<=n) c[i+lbt(i)]+=c[i];
m=g(); for(R i=,x,l,r;i<=m;++i) { x=g(),l=g(),r=g(); l>r?swap(l,r):void();
if(x&) printf("%lld\n",query(r)-query(l-));
else for(flg[r]?r=getf(r):;r>=l;r=getf(r)) change(r);
}
}

2019.04.11&&2019.06.15

最新文章

  1. java中的代码块执行顺序
  2. GCD学习之dispatch_barrier_async
  3. unix 常用命令
  4. Android Studio中新建项目时Your android sdk is out of date or is missing templates的解决办法
  5. 数据库分库分表(sharding)系列(四) 多数据源的事务处理
  6. python标准库 sysconfig模块
  7. 【filezilla】 ubuntu下安装filezilla
  8. git 初步
  9. 【NOIP2003提高组】加分二叉树
  10. FreeSWITCH 内线拨号 总是使用 dialplan/public 拨号计划,而对 dialplan/default 视而不见
  11. [译]Ocelot - Headers Transformation
  12. react 路由导航栏 withRouter
  13. Glide4 高效加载图片的配置【转】
  14. pip 安装pandas报UnicodeDecodeError: &#39;ascii&#39; codec can&#39;t decode byte 0xd5错
  15. vue组件 订单支付15分钟倒计时
  16. 10、java初始化顺序
  17. USB 之传输编码格式 NRZI 介绍
  18. FFmpeg 入门(6):音频同步
  19. 小白学CMD下运行MySQL
  20. 【转载】Direct3D纹理映射

热门文章

  1. python处理txt文件的一种情况
  2. matlab写txt文件
  3. kvm初体验之六:克隆
  4. POSTGRESQL 导入导出
  5. 关于STM32中GPIO的8种工作模式
  6. 递归------python实现列表创建二叉树
  7. 基于STM32的uCGUI移植和优化
  8. Mysql常用命令行大全(二)
  9. SpringBoot @RequestBody 中文乱码
  10. 《精通Spring4.X企业应用开发实战》读后感第四章