题目描述

最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分。超级计算机中的任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si秒开始,在第Ei秒后结束(第Si秒和Ei秒任务也在运行),其优先级为Pi。同一时间可能有多个任务同时执行,它们的优先级可能相同,也可能不同。调度系统会经常向查询系统询问,第Xi秒正在运行的任务中,优先级最小的Ki个任务(即将任务按照优先级从小到大排序后取前Ki个)的优先级之和是多少。特别的,如果Ki大于第Xi秒正在运行的任务总数,则直接回答第Xi秒正在运行的任务优先
级之和。上述所有参数均为整数,时间的范围在1到n之间(包含1和n)。

输入

输入文件第一行包含两个空格分开的正整数m和n,分别表示任务总数和时间范围。接下来m行,每行包含三个空格分开的正整数Si、Ei和Pi(Si≤Ei),描述一个任务。接下来n行,每行包含四个空格分开的整数Xi、Ai、Bi和Ci,描述一次查询。查询的参数Ki需要由公式 Ki=1+(Ai*Pre+Bi) mod Ci计算得到。其中Pre表示上一次查询的结果,对于第一次查询,Pre=1。

输出

输出共n行,每行一个整数,表示查询结果。

样例输入

4 3
1 2 6
2 3 3
1 3 2
3 3 4
3 1 3 2
1 1 3 4
2 2 4 3

样例输出

2
8
11


题解

主席树(不需要修改)

先将所有任务读进来,将优先级离散化,再按时间排序加入主席树中。

开始时间加进来,结束时间+1减去。

查询时,相当于查询[1,x]这一段区间,所以不用变为两树相减的形式。

先二分查找出最后一个小于等于x的(即第一个大于x的数-1),然后查询。

注意:需要开long long;sum数组修改时要加减原值,不能是离散化后的值;查询时如果l和r相等,需要返回l的原值乘p

upperbound码错调了半天。。。

#include <cstdio>
#include <algorithm>
#define N 200010
using namespace std;
typedef long long lint;
struct data
{
lint ti , p , mk;
}a[N];
lint v[N] , top , root[N] , lp[N << 6] , rp[N << 6] , si[N << 6] , sum[N << 6] , tot;
lint vti[N];
bool cmp1(data a , data b)
{
return a.p < b.p;
}
bool cmp2(data a , data b)
{
return a.ti < b.ti;
}
void pushup(lint x)
{
si[x] = si[lp[x]] + si[rp[x]];
sum[x] = sum[lp[x]] + sum[rp[x]];
}
void update(lint x , lint &y , lint l , lint r , lint p , lint a)
{
y = ++tot;
if(l == r)
{
si[y] = si[x] + a;
sum[y] = sum[x] + v[p] * a;
return;
}
lint mid = (l + r) >> 1;
if(p <= mid) rp[y] = rp[x] , update(lp[x] , lp[y] , l , mid , p , a);
else lp[y] = lp[x] , update(rp[x] , rp[y] , mid + 1 , r , p , a);
pushup(y);
}
lint query(lint y , lint l , lint r , lint p)
{
if(l == r)
return v[l] * p;
lint mid = (l + r) >> 1;
if(p <= si[lp[y]]) return query(lp[y] , l , mid , p);
else return query(rp[y] , mid + 1 , r , p - si[lp[y]]) + sum[lp[y]];
}
int main()
{
lint n , m , i , x , y , z , pre = 1 , c , k , p;
scanf("%lld%lld" , &n , &m);
for(i = 1 ; i <= n ; i ++ )
{
scanf("%lld%lld%lld" , &x , &y , &z);
a[2 * i - 1].ti = x;
a[2 * i].ti = y + 1;
a[2 * i - 1].p = a[2 * i].p = z;
a[2 * i - 1].mk = 1;
a[2 * i].mk = -1;
}
sort(a + 1 , a + 2 * n + 1 , cmp1);
for(i = 1 ; i <= 2 * n ; i ++ )
{
if(a[i].p != v[top]) v[++top] = a[i].p;
a[i].p = top;
}
sort(a + 1 , a + 2 * n + 1 , cmp2);
for(i = 1 ; i <= 2 * n ; i ++ )
{
vti[i] = a[i].ti;
update(root[i - 1] , root[i] , 1 , top , a[i].p , a[i].mk);
}
while(m -- )
{
scanf("%lld%lld%lld%lld" , &c , &x , &y , &z);
k = 1 + (x * pre + y) % z;
p = upper_bound(vti + 1 , vti + 2 * n + 1 , c) - vti - 1;
if(si[root[p]] < k) pre = sum[root[p]];
else pre = query(root[p] , 1 , top , k);
printf("%lld\n" , pre);
}
return 0;
}

【bzoj3932】[CQOI2015]任务查询系统 离散化+主席树的更多相关文章

  1. BZOJ3932 CQOI2015 任务查询系统 【主席树】

    BZOJ3932 CQOI2015 任务查询系统 Description 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的任务用三元组(Si,Ei, ...

  2. [BZOJ3932] [CQOI2015]任务查询系统(主席树 || 树状数组 套 主席树 + 差分 + 离散化)

    传送门 看到这个题有个很暴力的想法, 可以每一个时间点都建一颗主席树,主席树上叶子节点 i 表示优先级为 i 的任务有多少个. 当 x 到 y 有个优先级为 k 的任务时,循环 x 到 y 的每个点, ...

  3. [bzoj3932][CQOI2015]任务查询系统_主席树

    任务查询系统 bzoj-3932 CQOI-2015 题目大意:最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的任务用三元组(Si,Ei,Pi)描述 ...

  4. [bzoj3932][CQOI2015][任务查询系统] (主席树)

    Description 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si ...

  5. [BZOJ3932][CQOI2015]任务查询系统(差分+主席树)

    题面 分析 对于一个区间修改(s,e,v),我们可以将它差分,这样就变成了单点修改s和e+1(s插入,t+1删除) 我们用主席树维护差分数组的前缀和,第i棵主席树维护区间[1,i]之间的所有差分值 那 ...

  6. BZOJ_3932_[CQOI2015]任务查询系统_主席树

    BZOJ_3932_[CQOI2015]任务查询系统_主席树 题意: 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei,P ...

  7. 【BZOJ3932】任务查询系统(主席树)

    [BZOJ3923]任务查询系统(主席树) 题面 Description 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei ...

  8. bzoj3932 / P3168 [CQOI2015]任务查询系统(主席树+差分)

    P3168 [CQOI2015]任务查询系统 看到第k小,就是主席树辣 对于每一段任务(a,b,k),在版本a的主席树+k,版本b+1的主席树-k 同一时间可能有多次修改,所以开个vector存操作, ...

  9. 2018.06.30 BZOJ 3932: [CQOI2015]任务查询系统(主席树)

    3932: [CQOI2015]任务查询系统 Time Limit: 20 Sec Memory Limit: 512 MB Description 最近实验室正在为其管理的超级计算机编制一套任务管理 ...

随机推荐

  1. sas简单使用

    1 数据存取: 逻辑库: libname  自己起的名字 ‘文件所在的路径’,若无这步数据则存在默认的work中. 另一个方法在sas里自己建立一个逻辑库,但是关闭后就消失了. 新建数据:data  ...

  2. Homebrew安装Redis找不到redis.conf文件

    使用Homebrew安装redis完成后,使用命令 redis-server 启动redis,如下图所示: 启动信息中存在一条警告信息:没有指定的配置文件 然而在安装目录中并没有发现redis.con ...

  3. Intellij打包jar文件,“java.lang.SecurityException: Invalid signature file digest for Manifest main attrib

    下面是使用Intellij 打包jar文件的步骤,之后会有运行jar文件时遇到的错误. 打包完成. ================================================== ...

  4. 转:asp.net mvc ef 性能监控调试工具 MiniProfiler

    MiniProfiler官网:http://miniprofiler.com/ MiniProfiler的一个特别有用的功能是它与数据库框架的集成.除了.NET原生的 DbConnection类,Mi ...

  5. 用python读取配置文件config.ini

    还在学习中...写的有点凌乱 感觉还是应该先学会读取配置文件才行,把一些经常需要修改的但是又经常需要用到的参数放到配置文件中方便使用(我是这么觉得的) 首先是config.ini的存放位置,我们把它放 ...

  6. 【system.date】使用说明

    对象:system.date 说明:提供一系列针对日期类型的操作 目录: 方法 返回 说明  system.date.isDate( date_string )  [True | False]  判断 ...

  7. HADOOP docker(二):HDFS 高可用原理

        1.环境简述2.QJM HA简述2.1为什么要做HDFS HA?2.2 HDFS HA的方式2.2 HSFS HA的结构2.3 机器要求3.部署HDFS HA3.1 详细配置3.2 部署HDF ...

  8. CentOS6 安装VNCserver

    1.下载vncserver yum install tigervnc tigervnc-server -y 2.配置 vncserver vi /etc/sysconfig/vncserver 在文件 ...

  9. FZU 1844 Earthquake Damage(最大流最小割)

    Problem Description Open Source Tools help earthquake researchers stay a step ahead. Many geological ...

  10. PHP 5.6.32 增加pdo_dblib.so拓展

    首先说明,php增加pdo_dblib.so拓展不需要重新编译php源文件,只需要增加dblib源包即可. 1.下载安装所需包 1.#下载 wget http://mirrors.ibiblio.or ...