LOJ数列分块入门九题(中)
#6281. 数列分块入门 5 - 题目 - LibreOJ (loj.ac)
区间开方,区间求和题。
显然,针对区间维护开方操作很难做到,于是考虑其值的性质,显然,int范围内的值最多开方6次就会变为1,之后再开方依然为1,于是考虑暴力维护一个区间内的值是否全部为1即可,全部为1标记上之后再遇到O(1)处理即可。复杂度大概就为O(n*6+n*sqrt(n))
代码如下
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=(1e5)+5;
int a[maxn];
int sum[maxn];
int bel[maxn];
int L[maxn],R[maxn];
bool tag[maxn];
void build(int n){
int block=sqrt(n);
int num=n/block;
if(n%block)num++;
for(int i=1;i<=num;i++){
L[i]=(i-1)*block+1;
R[i]=i*block;
for(int j=L[i];j<=R[i];j++){
bel[j]=i;
sum[i]+=a[j];
}
}
}
void update(int l,int r){
if(bel[l]==bel[r]){
for(int i=l;i<=r;i++){
if(tag[bel[i]])break;
sum[bel[i]]=sum[bel[i]]-a[i]+int(sqrt(a[i]));
a[i]=int(sqrt(a[i]));
}
}
else{
for(int i=l;i<=R[bel[l]];i++){
if(tag[bel[i]])break;
sum[bel[i]]=sum[bel[i]]-a[i]+int(sqrt(a[i]));
a[i]=int(sqrt(a[i]));
}
for(int i=L[bel[r]];i<=r;i++){
if(tag[bel[i]])break;
sum[bel[i]]=sum[bel[i]]-a[i]+int(sqrt(a[i]));
a[i]=int(sqrt(a[i]));
}
for(int i=bel[l]+1;i<bel[r];i++){
bool fg=1;
for(int j=L[i];j<=R[i];j++){
if(tag[i])break;
sum[i]=sum[i]-a[j]+int(sqrt(a[j]));
a[j]=int(sqrt(a[j]));
if(a[j]>1)fg=0;
}
if(fg)tag[i]=1;
}
}
}
void calc(int l,int r){
int ans=0;
if(bel[l]==bel[r]){
for(int i=l;i<=r;i++){
ans+=a[i];
}
}
else{
for(int i=l;i<=R[bel[l]];i++){
ans+=a[i];
}
for(int i=L[bel[r]];i<=r;i++){
ans+=a[i];
}
for(int i=bel[l]+1;i<bel[r];i++){
ans+=sum[i];
}
}
cout<<ans<<endl;
}
int main(){
int n;
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i];
build(n);
for(int i=1;i<=n;i++){
int opt,l,r,c;
cin>>opt>>l>>r>>c;
if(opt){
calc(l,r);
}
else{
update(l,r);
}
}
return 0;
}
#6283. 数列分块入门 7 - 题目 - LibreOJ (loj.ac)
区间乘法,区间加法,单点询问
类似线段树的实现操作,注意先乘后加。当时我调bug花了很久,主要原因是在处理离散块的时候我直接对a数组进行了修改操作,导致标记数组出错,解决方法为直接对离散块所属的完整块进行下放标记的操作(暴力但是可行)
代码如下
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=(1e5)+5;
const int modd=(1e4)+7;
int a[maxn];
int bel[maxn];
int L[maxn],R[maxn];
int tag1[maxn];
int tag2[maxn];
void build(int n){
int block=sqrt(n);
int num=n/block;
if(n%block)num++;
for(int i=1;i<=num;i++){
L[i]=(i-1)*block+1;
R[i]=min(n,i*block);
tag1[i]=0;
tag2[i]=1;
for(int j=L[i];j<=R[i];j++){
bel[j]=i;
}
}
}
void pushdown(int p){
int l=L[bel[p]],r=R[bel[p]];
for(int i=l;i<=r;i++){
a[i]=(a[i]*tag2[bel[p]]+tag1[bel[p]])%modd;
}
tag1[bel[p]]=0;
tag2[bel[p]]=1;
}
void add(int l,int r,int c){
if(bel[l]==bel[r]){
pushdown(l);
for(int i=l;i<=r;i++){
a[i]=(a[i]+c)%modd;
}
}
else{
pushdown(l);
for(int i=l;i<=R[bel[l]];i++){
a[i]=(a[i]+c)%modd;
}
pushdown(r);
for(int i=L[bel[r]];i<=r;i++){
a[i]=(a[i]+c)%modd;
}
for(int i=bel[l]+1;i<bel[r];i++){
tag1[i]=(tag1[i]+c)%modd;
}
}
}
void multiply(int l,int r,int c){
if(bel[l]==bel[r]){
pushdown(l);
for(int i=l;i<=r;i++){
a[i]=(a[i]*c)%modd;
}
}
else{
pushdown(l);
for(int i=l;i<=R[bel[l]];i++){
a[i]=(a[i]*c)%modd;
}
pushdown(r);
for(int i=L[bel[r]];i<=r;i++){
a[i]=(a[i]*c)%modd;
}
for(int i=bel[l]+1;i<bel[r];i++){
tag1[i]=(tag1[i]*c)%modd;
tag2[i]=(tag2[i]*c)%modd;
}
}
}
void calc(int i){
int ans=(a[i]*tag2[bel[i]]+tag1[bel[i]])%modd;
cout<<ans<<endl;
}
int main(){
int n;
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i];
build(n);
for(int i=1;i<=n;i++){
int opt,l,r,c;
cin>>opt>>l>>r>>c;
if(opt==0){
add(l,r,c);
}
if(opt==1){
multiply(l,r,c);
}
if(opt==2){
calc(r);
}
}
return 0;
}
#6284. 数列分块入门 8 - 题目 - LibreOJ (loj.ac)
区间询问等于一个数 c 的元素,并将这个区间的所有元素改为 c。
这题我折腾了挺久。思路为标记每个区间是否为同一个值,若是则标记为该值,否则标记为1e11。之后在分块的模板中再分三种情况进行分析。1.该区间全部为c,2.该区间全部为另外一个值,3.该区间的值不统一。注意点:在离散块情况2进行修改操作时需要先进行标记下放操作。
代码如下
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=(1e5)+5;
ll a[maxn];
ll tag[maxn];
ll bel[maxn];
ll L[maxn],R[maxn];
void build(ll n){
ll block=sqrt(n);
ll num=n/block;
if(n%block)num++;
for(int i=1;i<=num;i++){
L[i]=(i-1)*block+1;
R[i]=min(n,i*block);
tag[i]=1e11;
for(int j=L[i];j<=R[i];j++){
bel[j]=i;
}
}
}
void pushdown(int p){
for(int i=L[p];i<=R[p];i++){
a[i]=tag[p];
}
}
void ope(ll l,ll r,ll c){
ll ans=0;
if(bel[l]==bel[r]){
if(tag[bel[l]]==c){
ans+=r-l+1;
}
else if(tag[bel[l]]!=1e11){
pushdown(bel[l]);
for(int i=l;i<=r;i++){
a[i]=c;
}
tag[bel[l]]=1e11;
}
else{
for(int i=l;i<=r;i++){
if(a[i]==c)ans++;
else a[i]=c;
}
}
}
else{
if(tag[bel[l]]==c){
ans+=R[bel[l]]-l+1;
}
else if(tag[bel[l]]!=1e11){
pushdown(bel[l]);
for(int i=l;i<=R[bel[l]];i++){
a[i]=c;
}
tag[bel[l]]=1e11;
}
else{
for(int i=l;i<=R[bel[l]];i++){
if(a[i]==c)ans++;
else a[i]=c;
}
}
if(tag[bel[r]]==c){
ans+=r-L[bel[r]]+1;
}
else if(tag[bel[r]]!=1e11){
pushdown(bel[r]);
for(int i=L[bel[r]];i<=r;i++){
a[i]=c;
}
tag[bel[r]]=1e11;
}
else{
for(int i=L[bel[r]];i<=r;i++){
if(a[i]==c)ans++;
else a[i]=c;
}
}
for(int i=bel[l]+1;i<bel[r];i++){
if(tag[i]==c){
ans+=(R[i]-L[i]+1);
}
else if(tag[i]==1e11){
for(int j=L[i];j<=R[i];j++){
if(a[j]==c)ans++;
else a[j]=c;
}
tag[i]=c;
}
else{
tag[i]=c;
}
}
}
cout<<ans<<endl;
}
int main(){
int n;
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i];
build(n);
for(int i=1;i<=n;i++){
ll l,r,c;
cin>>l>>r>>c;
ope(l,r,c);
}
return 0;
}
LOJ数列分块入门九题(中)的更多相关文章
- LOJ 数列分块入门 9 题解题报告
LOJ 数列分块入门 9 题解题报告 \(\text{By DaiRuiChen007}\) I. 数列分块入门 1 题目大意 \(\text{Link}\) 维护一个长度为 \(n\) 的序列,支持 ...
- 数列分块入门九题(三):LOJ6283~6285
Preface 最后一题我一直觉得用莫队是最好的. 数列分块入门 7--区间乘法,区间加法,单点询问 还是很简单的吧,比起数列分块入门 7就多了个区间乘. 类似于线段树,由于乘法的优先级高于加法,因此 ...
- 数列分块入门九题(二):LOJ6280~6282
Preface 个人感觉这中间的三题是最水的没有之一 数列分块入门 4--区间加法,区间求和 这个也是很多数据结构完爆的题目线段树入门题,但是练分块我们就要写吗 修改还是与之前类似,只不过我们要维护每 ...
- 数列分块入门九题(一):LOJ6277~6279
Preface 分块,一个神奇的暴力算法.可以把很多\(O(n^2)\)的数据结构题的暴力优化到常数极小的\(O(n\sqrt n)\).当一些毒瘤题无法用线段树,主席树,平衡树,树状数组...... ...
- [Loj] 数列分块入门 1 - 9
数列分块入门 1 https://loj.ac/problem/6277 区间加 + 单点查询 #include <iostream> #include <cstdio> #i ...
- loj 数列分块入门 6 9(区间众数)
6 题意 给出一个长为\(n\)的数列,以及\(n\)个操作,操作涉及单点插入,单点询问,数据随机生成. 题解 参考:http://hzwer.com/8053.html 每个块内用一个\(vecto ...
- loj 数列分块入门 5 7 8
5 题意 给出一个长为\(n\)的数列,以及\(n\)个操作,操作涉及区间开方,区间求和. 思路 用\(tag\)记录这一块是否已全为\(1\). 除分块外,还可用 树状数组+并查集(链表) 或者 线 ...
- LOJ 数列分块入门系列
目录 1.区间加+单点查 每个块维护tag,散的暴力改. code: #include<bits/stdc++.h> using namespace std; const int maxn ...
- LOJ6285 数列分块入门9(分块)
昨天对着代码看了一晚上 然后今天终于在loj上过了 数列分块入门9题撒花★,°:.☆( ̄▽ ̄)/$:.°★ . 然后相当玄学 块的大小调成\(\sqrt{n}\)会TLE,改成150就过了 啧 然后就 ...
- LOJ——#6277. 数列分块入门 1
~~推荐播客~~ 「分块」数列分块入门1 – 9 by hzwer 浅谈基础根号算法——分块 博主蒟蒻,有缘人可直接观摩以上大佬的博客... #6277. 数列分块入门 1 题目大意: 给出一个长为 ...
随机推荐
- Jmeter 命令执行脚本并输出报告
Jmeter 可以通过cmd 命令执行脚本,并输出报告 执行顺序如下: 第一 进入 安装的jmeter bin 目录下:(我安装的是E盘)
- uniapp改变icon
<!DOCTYPE html><html lang="zh-CN"> <head> <meta charset="utf-8&q ...
- AVL tree rotate
AVL tree single rotate /** * Rotate binary tree node with left child. * For AVL trees, this is a sin ...
- chklist
1. 重复检查的必要性.一段代码如果在测试期间没问题,也要间隔几天再去观察是否有问题.2. 如果是集群式的服务,使用定时任务要采用分布式锁,或使用工具随机发送一台都可以.3. 如果任务跑失败,需要支持 ...
- lxml库和BeautifulSoup库常用点小结
算是本人的学习笔记吧,仅供个人学习使用. 以下内容摘自<Python3网络爬虫开发实战--崔庆才著> 1.lxml库 XPath 常用规则: 表达式 描述 nodename 选取此节点的所 ...
- 大规模并行处理器编程实战_原书第2版_pdf
链接:https://pan.baidu.com/s/1c8ez8dCTz5bUQchwhXAF7w 提取码:tc1f
- 删除v-for方法生成的组件的方法
一.思路 在vue:data中的数据对象添加布尔类型的属性,用来表明其是否被渲染,如果需要删除组件,就把这个属性设置为false: 在计算属性vue:computed中,根据该布尔属性过滤,生成一个新 ...
- Linux工作中最常用命令整理
ls 命令:显示指定工作目录下之内容 ls -a # 显示所有文件夹,包含隐藏的. 和.. ls -l # 显示文件的详细信息,包含文件形态,权限,所属,大小,其实就是平常用的 ll ll -h # ...
- Python学习笔记文件读写之遍历目录树
随笔记录方便自己和同路人查阅. #------------------------------------------------我是可耻的分割线--------------------------- ...
- ant Vue select 多条数据重复的问题
1.问题描述:select实现下拉的功能,并且带前端搜索,但是后端给的数据没去重,对象中所有的参数值有保持一致,试了下,发现只要key值和value值不一致就不会报错: 2.解决:key值用index ...