luogu P4145 上帝造题的七分钟2 / 花神游历各国 维护区间和&&区间开根号
2024-10-19 23:38:22
因为开根号能使数字减小得非常快
所以开不了几次(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
最新文章
- java中的代码块执行顺序
- GCD学习之dispatch_barrier_async
- unix 常用命令
- Android Studio中新建项目时Your android sdk is out of date or is missing templates的解决办法
- 数据库分库分表(sharding)系列(四) 多数据源的事务处理
- python标准库 sysconfig模块
- 【filezilla】 ubuntu下安装filezilla
- git 初步
- 【NOIP2003提高组】加分二叉树
- FreeSWITCH 内线拨号 总是使用 dialplan/public 拨号计划,而对 dialplan/default 视而不见
- [译]Ocelot - Headers Transformation
- react 路由导航栏 withRouter
- Glide4 高效加载图片的配置【转】
- pip 安装pandas报UnicodeDecodeError: &#39;ascii&#39; codec can&#39;t decode byte 0xd5错
- vue组件 订单支付15分钟倒计时
- 10、java初始化顺序
- USB 之传输编码格式 NRZI 介绍
- FFmpeg 入门(6):音频同步
- 小白学CMD下运行MySQL
- 【转载】Direct3D纹理映射