题面

\rm HandInDevil

HandInDevil 的头发很油，因此随时有跳蚤跳上

\rm HandInDevil

HandInDevil 的头发。

现在把

\rm HandInDevil

HandInDevil 的头发抽象成一个数轴，一开始上面没有跳蚤。每一个时刻，有三种事件之一会发生：

有一个每次会跳

t

t

t 格的跳蚤跳上

H

a

n

d

I

n

D

e

v

i

l

\rm HandInDevil

HandInDevil 的头发的位置

x

x

x 。
H
a

n

d

I

n

D

e

v

i

l

\rm HandInDevil

HandInDevil 开始挠头发，每一只跳蚤向右跳一次，跳的距离为各自的

t

t

t 。
H
a

n

d

I

n

D

e

v

i

l

\rm HandInDevil

HandInDevil 没有挠到跳蚤，打算找下一个目标，向你询问头发上区间

[

l

,

r

]

[l,r]

[l,r] 内跳蚤的个数。

输入

第一行一个整数

Q，表示事件个数。

接下来

Q 行描述事件，输入若干个整数：

若第一个数为 1 ，则接下来跟着两个整数

x

i

x_i

xi 和

t

i

t_i

ti ，表示一个跳蚤跳了上去。
若第一个数为 2 ，则
H

a

n

d

I

n

D

e

v

i

l

\rm HandInDevil

HandInDevil 挠了一次头发，所有跳蚤往右跳。
若第一个数为 3 ，则接下来跟着两个整数
l

i

l_i

li 和

r

i

r_i

ri ，你需要输出此时区间

[

l

i

,

r

i

]

[l_i,r_i]

[li,ri] 中跳蚤的数量。

输出

对于每一个事件 3 ，输出一行一个整数，表示询问的答案。

Sample Input

Sample Output

数据范围

≤

1\leq Q,x_i,t_i,l_i,r_i\leq 10^5\;,\;l_i\leq r_i

1≤Q,xi,ti,li,ri≤105,li≤ri

2 ms , 512 mb

\text{2 ms , 512 mb}

2 ms , 512 mb

题解

-既然跟 HandInDevil 扯上关系那肯定是分块题啊-

考虑分块，令块大小为

B，头发有效长度为

n，我们把

t_i> B

ti>B 的跳蚤暴力跳，然后用线段树或树状数组等数据结构维护区间和，这部分复杂度为

(

log

⁡

)

O(Q\frac{n}{B}\log n)

O(QBnlogn) 。

然后对于

≤

t_i\leq B

ti≤B 的跳蚤，我们建立

B 棵平衡树，假设跳蚤

i 跳上时已经发生了

cnt_i

cnti 次事件 2，那么我们就把

−

∗

x_i-cnt_i*t_i

xi−cnti∗ti 加进第

t_i

ti 棵平衡树里。

询问的时候，可以先加上区间内

t_i>B

ti>B 的跳蚤数量，然后在每棵平衡树里进行询问。假设此时事件 2 发生了

′

cnt'

cnt′ 次，对于第

i 棵平衡树，先把

[

]

[l_i,r_i]

[li,ri] 变成

[

−

′

∗

−

′

∗

]

[l_i-cnt'*i,r_i-cnt'*i]

[li−cnt′∗i,ri−cnt′∗i] ，然后在该平衡树里询问值在这个区间内的结点数量，加进答案。这部分复杂度是

(

log

⁡

)

O(QB\log Q)

O(QBlogQ) 。

由于

Q 和

n 同阶，所以取

B=\sqrt{n}

B=n

时复杂度最小，总复杂度为

(

log

⁡

)

O(Q\sqrt n\log n)

O(Qn

logn) ，险过。

CODE

#include<set>

#include<map>

#include<cmath>

#include<ctime>

#include<stack>

#include<queue>

#include<vector>

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

using namespace std;

#define MAXN 100005

#define DB double

#define LL long long

#define ENDL putchar('\n')

#define lowbit(x) (-(x) & (x))

#define FI first

#define SE second

#pragma GCC optimize(2)

LL read() {

	LL f = 1,x = 0;char s = getchar();

	while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}

	while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}

	return f * x;

}

const int MOD = 998244353;

const int SQ = 317;

int n = 100000,m,i,j,s,o,k;

struct np{int s[2];np(){s[0]=s[1]=0;}np(int A,int B){s[0]=A;s[1]=B;}};

struct tr{

	int s[2];

	LL ky;int hp,siz;

	tr(){s[0]=s[1]=ky=hp=siz=0;}

}tre[MAXN];

int buc[MAXN],cntb;

int newnode(LL key) {

	int x = buc[cntb --]; tre[x] = tr();

	tre[x].ky = key; tre[x].hp = rand()*1ll*rand()%MOD;

	tre[x].siz = 1; return x;

}

int update(int x) {

	if(!x) return x;

	int ls = tre[x].s[0],rs = tre[x].s[1];

	tre[x].siz = tre[ls].siz + tre[rs].siz + 1;

	return x;

}

np spli(int x,LL key) {

	np as(0,0); if(!x) return as;

	int d = (tre[x].ky <= key);

	as = spli(tre[x].s[d],key);

	tre[x].s[d] = as.s[d^1];

	as.s[d^1] = update(x);

	return as;

}

np spli2(int x,int kth) {

	np as(0,0); if(!x) return as;

	int d = (tre[tre[x].s[0]].siz+1 <= kth);

	if(d) kth -= tre[tre[x].s[0]].siz+1;

	as = spli(tre[x].s[d],kth);

	tre[x].s[d] = as.s[d^1];

	as.s[d^1] = update(x);

	return as;

}

int merg(int p1,int p2) {

	if(!p1 || !p2) return p1+p2;

	if(tre[p1].hp < tre[p2].hp) {

		tre[p1].s[1] = merg(tre[p1].s[1],p2);

		return update(p1);

	}

	tre[p2].s[0] = merg(p1,tre[p2].s[0]);

	return update(p2);

}

int ins(int x,LL key) {

	np p = spli(x,key);

	return merg(p.s[0],merg(newnode(key),p.s[1]));

}

int del(int x,LL key) {

	np pl = spli(x,key-1);

	np pr = spli2(pl.s[1],1);

	if(pr.s[0] && tre[pr.s[0]].ky == key) buc[++ cntb] = pr.s[0];

	else pr.s[1] = merg(pr.s[0],pr.s[1]);

	return merg(pl.s[0],pr.s[1]);

}

int query(int &x,LL l,LL r) {

	if(l > r) return 0;

	np p1 = spli(x,l-1);

	np p2 = spli(p1.s[1],r);

	int as = tre[p2.s[0]].siz;

	x = merg(p1.s[0],merg(p2.s[0],p2.s[1]));

	return as;

}

int bl[SQ+5],li[SQ+5],rt[SQ+5];

int c[MAXN];

void addt(int x,int y){while(x<=n)c[x]+=y,x+=lowbit(x);}

int sum(int x) {int as=0;while(x>0)as+=c[x],x-=lowbit(x);return as;}

pair<int,int> PAIR(int A,int B) {pair<int,int> a;a.FI=A;a.SE=B;return a;}

pair<int,int> q1[MAXN];

int cnq;

int main() {

	freopen("flea.in","r",stdin);

	freopen("flea.out","w",stdout);

	m = 0;

	for(int i = 100000;i > 0;i --) {

		buc[++ cntb] = i;

		int bid = (i+SQ-1)/SQ;

		m = max(m,bid); li[bid] = i;

	}

	li[m+1] = 100001;

	srand(time(0));

	int Q = read(),cnt2 = 0;

	while(Q --) {

		k = read();

		if(k == 1) {

			s = read();o = read();

			if(o > SQ) q1[++ cnq] = PAIR(s,o),addt(s,1);

			else {

				LL ad = (LL)s-cnt2*1ll*o;

				rt[o] = ins(rt[o],ad);

			}

		}

		else if(k == 2) {

			cnt2 ++;

			int leq = cnq;cnq = 0;

			for(int i = 1;i <= leq;i ++) {

				int ss = q1[i].FI,tt = q1[i].SE;

				addt(ss,-1);

				ss += tt;

				if(ss <= 100000) {

					addt(ss,1);

					q1[++ cnq] = PAIR(ss,tt);

				}

			}

		}

		else {

			s = read();o = read();

			int ans = sum(o)-sum(s-1);

			for(int i = 1;i <= SQ;i ++) {

				LL ad1 = (LL)s-cnt2*1ll*i , ad2 = (LL)o-cnt2*1ll*i;

				ans += query(rt[i],ad1,ad2);

			}

			printf("%d\n",ans);

		}

	}

	return 0;

}

Better Solution

有一个更快的做法，我们把

t_i>B

ti>B 的跳蚤暴力跳时，用分块维护区间和。把头发序列分块，每个块大小为

B，修改

(

)

O(1)

O(1) ，询问

(

)

O(\frac{n}{B})

O(Bn)，这部分修改总复杂度

(

)

O(Q\frac{n}{B})

O(QBn)，询问总复杂度

(

)

O(Q\frac{n}{B})

O(QBn) 。

那么取块大小为

log

⁡

B=\sqrt{\frac{n}{\log n}}

B=lognn

，总复杂度就为

(

log

⁡

)

O(Q\sqrt{n\log n})

O(Qnlogn

)，实际块大小还可以根据平衡树的常数微调，比之前的做法快上不少。

CODE

#include<set>

#include<map>

#include<cmath>

#include<ctime>

#include<stack>

#include<queue>

#include<vector>

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

using namespace std;

#define MAXN 100005

#define DB double

#define LL long long

#define ENDL putchar('\n')

#define lowbit(x) (-(x) & (x))

#define FI first

#define SE second

#pragma GCC optimize(2)

LL read() {

	LL f = 1,x = 0;char s = getchar();

	while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}

	while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}

	return f * x;

}

const int MOD = 998244353;

const int SQ = 77;

int n = 100000,m,i,j,s,o,k;

struct np{int s[2];np(){s[0]=s[1]=0;}np(int A,int B){s[0]=A;s[1]=B;}};

struct tr{

	int s[2];

	LL ky;int hp,siz;

	tr(){s[0]=s[1]=ky=hp=siz=0;}

}tre[MAXN];

int buc[MAXN],cntb;

int newnode(LL key) {

	int x = buc[cntb --]; tre[x] = tr();

	tre[x].ky = key; tre[x].hp = rand()*1ll*rand()%MOD;

	tre[x].siz = 1; return x;

}

int update(int x) {

	if(!x) return x;

	int ls = tre[x].s[0],rs = tre[x].s[1];

	tre[x].siz = tre[ls].siz + tre[rs].siz + 1;

	return x;

}

np spli(int x,LL key) {

	np as(0,0); if(!x) return as;

	int d = (tre[x].ky <= key);

	as = spli(tre[x].s[d],key);

	tre[x].s[d] = as.s[d^1];

	as.s[d^1] = update(x);

	return as;

}

np spli2(int x,int kth) {

	np as(0,0); if(!x) return as;

	int d = (tre[tre[x].s[0]].siz+1 <= kth);

	if(d) kth -= tre[tre[x].s[0]].siz+1;

	as = spli(tre[x].s[d],kth);

	tre[x].s[d] = as.s[d^1];

	as.s[d^1] = update(x);

	return as;

}

int merg(int p1,int p2) {

	if(!p1 || !p2) return p1+p2;

	if(tre[p1].hp < tre[p2].hp) {

		tre[p1].s[1] = merg(tre[p1].s[1],p2);

		return update(p1);

	}

	tre[p2].s[0] = merg(p1,tre[p2].s[0]);

	return update(p2);

}

int ins(int x,LL key) {

	np p = spli(x,key);

	return merg(p.s[0],merg(newnode(key),p.s[1]));

}

int del(int x,LL key) {

	np pl = spli(x,key-1);

	np pr = spli2(pl.s[1],1);

	if(pr.s[0] && tre[pr.s[0]].ky == key) buc[++ cntb] = pr.s[0];

	else pr.s[1] = merg(pr.s[0],pr.s[1]);

	return merg(pl.s[0],pr.s[1]);

}

int query(int &x,LL l,LL r) {

	if(l > r) return 0;

	np p1 = spli(x,l-1);

	np p2 = spli(p1.s[1],r);

	int as = tre[p2.s[0]].siz;

	x = merg(p1.s[0],merg(p2.s[0],p2.s[1]));

	return as;

}

int bl[MAXN],li[MAXN],rt[MAXN];

int ct[MAXN];

pair<int,int> PAIR(int A,int B) {pair<int,int> a;a.FI=A;a.SE=B;return a;}

pair<int,int> q1[MAXN];

int cnq;

int main() {

	freopen("flea.in","r",stdin);

	freopen("flea.out","w",stdout);

	m = 0;

	for(int i = 100000;i > 0;i --) {

		buc[++ cntb] = i;

		int bid = (i+SQ-1)/SQ;

		m = max(m,bid); li[bid] = i;

	}

	li[m+1] = 100001;

	srand(time(0));

	int Q = read(),cnt2 = 0;

	while(Q --) {

		k = read();

		if(k == 1) {

			s = read();o = read();

			if(o > SQ) q1[++ cnq] = PAIR(s,o),ct[s] ++,bl[(s+SQ-1)/SQ] ++;

			else {

				LL ad = (LL)s-cnt2*1ll*o;

				rt[o] = ins(rt[o],ad);

			}

		}

		else if(k == 2) {

			cnt2 ++;

			int leq = cnq;cnq = 0;

			for(int i = 1;i <= leq;i ++) {

				int ss = q1[i].FI,tt = q1[i].SE;

				ct[ss] --; bl[(ss+SQ-1)/SQ] --;

				ss += tt;

				if(ss <= 100000) {

					ct[ss] ++; bl[(ss+SQ-1)/SQ] ++;

					q1[++ cnq] = PAIR(ss,tt);

				}

			}

		}

		else {

			s = read();o = read();

			int ans = 0;

			int ll = (s+SQ-1)/SQ,rr = (o+SQ-1)/SQ;

			for(int i = ll+1;i <= rr-1;i ++) ans += bl[i];

			if(ll == rr) for(int i = s;i <= o;i ++) ans += ct[i];

			else {

				for(int i = s;(i+SQ-1)/SQ == ll;i ++) ans += ct[i];

				for(int i = o;(i+SQ-1)/SQ == rr;i --) ans += ct[i];

			}

			for(int i = 1;i <= SQ;i ++) {

				LL ad1 = (LL)s-cnt2*1ll*i , ad2 = (LL)o-cnt2*1ll*i;

				ans += query(rt[i],ad1,ad2);

			}

			printf("%d\n",ans);

		}

	}

	return 0;

}

HandInDevil 的头发（分块）的更多相关文章

列表 ul ol dl 和块级标签和行及标签之间的转换
1. 无序列表有序列表自定义列表 1,无序列表第一你不必须有子标签 <li></li> 第二 ul天生自带内外边距 List-style的属性值 circle(空心圆 ...
js瀑布流原理实现揭秘 javascript 原生实现
web,js瀑布流揭秘瀑布流再很久之前流行,可能如我一样入行晚的 ,可能就没有机会去使用.但是这个技术终究是个挺炫酷的东西,花了一个上午来研究,用原生js实现了一个,下面会附上源码,供大家解读. 说 ...
ASP.NET 文件下载
using System; using System.Web; using System.IO; public partial class _Default : System.Web.UI.Page ...
漫话JavaScript与异步·第一话——异步：何处惹尘埃
自JavaScript诞生之日起,频繁与异步打交道便是这门语言的使命,并为此衍生出了许多设计和理念.因此,深入理解异步的概念对于前端工程师来说极为重要. 什么是异步? 程序是分"块" ...
css小知识 2
效果为为什么还出现出现不同的效果? 浏览器在解析第二个p的时候,因为第二个字母见没有空格,它会认为这是一个单词没有写完,所以不会换行列表 1.无序列表ul 第二,内部必须有子代标签<li&g ...
http状态码以及请求响应头相关
1xx消息[编辑] 这一类型的状态码,代表请求已被接受,需要继续处理.这类响应是临时响应,只包含状态行和某些可选的响应头信息,并以空行结束.由于HTTP/1.0协议中没有定义任何1xx状态码,所以除非 ...
HTML知识梳理（笔记）
HTML常见元素 meta 定义和用法<meta> 元素可提供有关页面的元信息(meta-information),比如针对搜索引擎和更新频度的描述和关键词. <meta> 标 ...
操作系统之IO管理
IO系统结构设备的分类按数据组织分块设备: 信息的存取总是以数据块为单位. 它属于有结构设备,如磁盘等. 磁盘设备的基本特征是传输速率较高,以及可寻址,即对它可随机地读/写任一块. 字符设备: ...
TCP/IP基础总结性学习（6）
HTTP 首部一. HTTP 报文首部 1.HTTP 报文的结构: 2.HTTP 请求报文图示: 举例子: 3.HTTP 响应报文: 下面的示例是访问 http://hackr.jp 时,请求报文 ...

随机推荐

ER图/模型转换为关系模型
ER图中的主要成分是实体类型和联系类型,转换规则就是如何把实体类型.联系类型转换成关系模式. 1. 二元联系转换规则1.1(实体类型的转换):将每个实体类型转换成一个关系模式,实体的属性即为关系模式 ...
camunda流程引擎概念术语
前言本文重点介绍开源流程引擎camunda的核心概念,这些概念同样适用于JBMP.Activiti.Flowable流程引擎,了解这些基本概念和原理,使用流程引擎API将更得心应手. 一.Proce ...
12.web基础与HTTP协议
web基础与HTTP协议目录 web基础与HTTP协议 web基础域名概述 HTML概述 HTML基本标签 HTML语法规则 HTML文件结构头标签中常用标签内容标签中常用标签静态网页与动态 ...
1.windows编程入门MessageBox使用 -windows编程
引言:刚开始入门windows编程的时候,我记得当时我对MFC的给出的一大堆代码感到束手无策.因为历史的缘故,windows编程入门的代码并没有体现出C++语言的简洁性,相反一上来就给了我们一大堆代码 ...
Vue MD5加密你用吗?
安装 npm install --save js-md5 1.按需引入(在你需要的项目中引入) 引入: import md5 from 'js-md5' 使用: md5('加密信息') 2.全局引入( ...
【python基础】第09回数据类型内置方法 01
本章内容概要 1.数据类型的内置方法简介 2.整型相关方法 3.浮点型相关方法 4.字符串相关方法 5.列表相关方法本章内容详情 1.数据类型的内置方法简介数据类型是用来记录事物状态的,而事物的状 ...
NC201605 Bits
NC201605 Bits 题目题目描述 Nancy喜欢做游戏! 汉诺塔是一个神奇的游戏,神奇在哪里呢? 给出 \(3\) 根柱子,最开始时 \(n\) 个盘子按照大小被置于最左的柱子. 如果盘子数 ...
让你的Nginx支持分布式追踪
Background NGINX 是一个通用且流行的应用程序.也是最流行的 Web 服务器,它可用于提供静态文件内容,但也通常与其他服务一起用作分布式系统中的组件,在其中它用作反向代理.负载均衡或 ...
C# / VB.NET 将Html转为Word
本文分享以C#程序代码为例,实现将Html文件转换Word文档的方法(附VB.NET代码).在实际转换场景中可参考本文的方法,转换前,请按照如下方法引用Word API的dll文件到Visual St ...
CodeForce——Deltix Round, Autumn 2021 (open for everyone, rated, Div. 1 + Div. 2)前三道题目题解
目录 A: B: C: 题目链接 A Divide and Multiply standard input/output 1 s, 256 MB 正在上传-重新上传取消 x13036 B Willia ...

HandInDevil 的头发 （分 块）

题面

输入

输出

Sample Input

Sample Output

数据范围

题解

CODE

Better Solution

CODE

HandInDevil 的头发 （分 块）的更多相关文章

随机推荐

热门专题

HandInDevil 的头发（分块）

HandInDevil 的头发（分块）的更多相关文章