【洛谷3759】[TJOI2017] 不勤劳的图书管理员（树套树）

最新推荐文章于 2022-10-12 12:46:22 发布
最新推荐文章于 2022-10-12 12:46:22 发布
阅读量82 收藏
点赞数
原文链接：http://www.cnblogs.com/chenxiaoran666/p/Luogu3759.html
版权

点此看题面

大致题意： 给定一个序列，每个元素有两个属性\(a_i\)和\(v_i\)，每次操作改变两个元素的位置，求每次操作后\(\sum{v_i+v_j}[i<j,a_i>a_j]\)。

关于题意的解读

其实，题目差不多就是让我们求逆序对（只不过每个逆序对有一个权值）。

这会让我们联想起一道题目：【洛谷3157】[CQOI2011] 动态逆序对。

这两道题目似乎只有以下几点不同：

  • 这道题目没有卡树套树。

  • 这道题目貌似不可以用\(CDQ\)分治做。

于是，我们就可以快乐地写 线段树套\(Treap\) 了。

大致思路

这题的思路其实与动态逆序对是差不多的。

如果要将一个数从位置\(x\)换到\(y\)（\(x\le y\)），则增加的贡献值应为\(\sum_{i=x+1}^yv_i+v_x[a_i>a_x]-\sum_{i=x+1}^yv_i+v_x[a_i<a_x]\)，用线段树套\(Treap\)可以在\(O(log^2n)\)的时间内做到这一点。

从\(y\)到\(x\)的情况同理。不过要注意，\(v_x+v_y\)会被重复计算，因此在从\(y\)到\(x\)时，我们就增加贡献\(\sum_{i=x+1}^{y-1}v_i+v_y[a_i<a_y]-\sum_{i=x+1}^{y-1}v_i+v_y[a_i>a_y]\)（注意这里的上界减了\(1\)，避免重复计算）。

具体实现还是见代码吧。

代码

#include<bits/stdc++.h>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define uint unsigned LL
#define LL long long
#define ull unsigned long long
#define swap(x,y) (x^=y,y^=x,x^=y)
#define abs(x) ((x)<0?-(x):(x))
#define INF 1e9
#define Inc(x,y) ((x+=(y))>=MOD&&(x-=MOD))
#define ten(x) (((x)<<3)+((x)<<1)) 
#define MOD 1000000007
#define N 50000 
#define LogN 30
using namespace std;
int n;
struct Book
{
    int pos,val;
}s[N+5];
class FIO
{
    private:
        #define Fsize 100000
        #define tc() (FinNow==FinEnd&&(FinEnd=(FinNow=Fin)+fread(Fin,1,Fsize,stdin),FinNow==FinEnd)?EOF:*FinNow++)
        #define pc(ch) (FoutSize<Fsize?Fout[FoutSize++]=ch:(fwrite(Fout,1,FoutSize,stdout),Fout[(FoutSize=0)++]=ch))
        int f,FoutSize,OutputTop;char ch,Fin[Fsize],*FinNow,*FinEnd,Fout[Fsize],OutputStack[Fsize];
    public:
        FIO() {FinNow=FinEnd=Fin;}
        inline void read(int &x) {x=0,f=1;while(!isdigit(ch=tc())) f=ch^'-'?1:-1;while(x=ten(x)+(ch&15),isdigit(ch=tc()));x*=f;}
        inline void read_char(char &x) {while(isspace(x=tc()));}
        inline void read_string(string &x) {x="";while(isspace(ch=tc()));while(x+=ch,!isspace(ch=tc())) if(!~ch) return;}
        inline void write(int x) {if(!x) return (void)pc('0');if(x<0) pc('-'),x=-x;while(x) OutputStack[++OutputTop]=x%10+48,x/=10;while(OutputTop) pc(OutputStack[OutputTop]),--OutputTop;}
        inline void write_char(char x) {pc(x);}
        inline void write_string(string x) {register int i,len=x.length();for(i=0;i<len;++i) pc(x[i]);}
        inline void end() {fwrite(Fout,1,FoutSize,stdout);}
}F;
class Class_SegmentTreap//线段树套Treap
{
    private:
        int n,Root[N<<2];Book data[N+5];
        class Class_Treap
        {
            private:
                #define Rand() ((r*=233333LL)%=2147483647)
                #define PushUp(x) (node[x].Size=node[node[x].Son[0]].Size+node[node[x].Son[1]].Size+1,node[x].Sum=((node[node[x].Son[0]].Sum+node[node[x].Son[1]].Sum)%MOD+node[x].Val)%MOD)
                #define Rotate(x,d) (k=node[x].Son[d^1],node[x].Son[d^1]=node[k].Son[d],node[k].Son[d]=x,x=k,PushUp(node[x].Son[d]),PushUp(x))
                #define Build(pos,val) ((void)(node[Void[tot]]=Tree(pos,val),node[Void[tot]].Size=1,node[Void[tot]].Data=Rand()),Void[tot--])
                int tot,k,Void[N*LogN+5];ull r;
                struct Tree
                {
                   int Pos,Val,Sum,Size,Data,Son[2];
                   Tree(int x=0,int y=0):Pos(x),Val(y),Sum(y){Son[0]=Son[1]=0;}
                }node[N*LogN+5];
                inline void ins(int &x,Book val)
                {
                    if(!x) return (void)(x=Build(val.pos,val.val));
                    ++node[x].Size;
                    if(node[x].Pos>val.pos) {if(ins(node[x].Son[0],val),node[x].Data<node[node[x].Son[0]].Data) Rotate(x,1);}
                    else {if(ins(node[x].Son[1],val),node[x].Data<node[node[x].Son[1]].Data) Rotate(x,0);}
                    PushUp(x);
                }
                inline void del(int &x,int pos)
                {   
                    if(!x) return;
                    if(node[x].Pos==pos)
                    {
                        if(node[x].Son[0]||node[x].Son[1])
                        {
                            if(!node[x].Son[1]||node[node[x].Son[0]].Data>node[node[x].Son[1]].Data) Rotate(x,1),del(node[x].Son[1],pos);
                            else Rotate(x,0),del(node[x].Son[0],pos);
                        }
                        else return (void)(Void[++tot]=x,x=0);
                    }
                    else node[x].Pos>pos?del(node[x].Son[0],pos):del(node[x].Son[1],pos);
                    PushUp(x);
                }
            public:
                Class_Treap() {r=2333;for(register int i=N*LogN;i;--i) Void[++tot]=i;}
                inline void Insert(int &rt,Book val) {ins(rt,val);}
                inline void Delete(int &rt,int val) {del(rt,val);}
                inline int total_max(int rt,Book val)
                {
                    register int x=rt,res=0;
                    while(x)
                    {
                        if(node[x].Pos>val.pos) Inc(res,((node[node[x].Son[1]].Sum+node[x].Val)%MOD+1LL*(node[node[x].Son[1]].Size+1)*val.val%MOD)%MOD),x=node[x].Son[0];
                        else x=node[x].Son[1]; 
                    }
                    return res;
                }
                inline int total_min(int rt,Book val)
                {
                    register int x=rt,res=0;
                    while(x)
                    {
                        if(node[x].Pos<val.pos) Inc(res,((node[node[x].Son[0]].Sum+node[x].Val)%MOD+1LL*(node[node[x].Son[0]].Size+1)*val.val%MOD)%MOD),x=node[x].Son[1];
                        else x=node[x].Son[0]; 
                    }
                    return res;
                }
                #undef Build
                inline int GetSize(int x) {return node[x].Size;}
        }Treap;
        inline void ins(int l,int r,int rt,int pos,Book val)
        {
            register int i,mid=l+r>>1;Treap.Insert(Root[rt],val);
            if(l^r) pos<=mid?ins(l,mid,rt<<1,pos,val):ins(mid+1,r,rt<<1|1,pos,val);
        }
        inline void Upt(int l,int r,int rt,int pos,Book val)
        {
            register int i,mid=l+r>>1;Treap.Delete(Root[rt],data[pos].pos),Treap.Insert(Root[rt],val);
            if(l^r) pos<=mid?Upt(l,mid,rt<<1,pos,val):Upt(mid+1,r,rt<<1|1,pos,val);
        }
        inline int total_max(int l,int r,int rt,int ql,int qr,Book val)
        {
            if(ql<=l&&r<=qr) return Treap.total_max(Root[rt],val);
            register int mid=l+r>>1,res=0;
            if(ql<=mid) res=total_max(l,mid,rt<<1,ql,qr,val);
            if(qr>mid) Inc(res,total_max(mid+1,r,rt<<1|1,ql,qr,val));
            return res;
        }
        inline int total_min(int l,int r,int rt,int ql,int qr,Book val)
        {
            if(ql<=l&&r<=qr) return Treap.total_min(Root[rt],val);
            register int mid=l+r>>1,res=0;
            if(ql<=mid) res=total_min(l,mid,rt<<1,ql,qr,val);
            if(qr>mid) Inc(res,total_min(mid+1,r,rt<<1|1,ql,qr,val));
            return res;
        }
    public:
        inline void Init(int len) {n=len;}
        inline void Insert(int pos,Book val) {ins(1,n,1,pos,val),data[pos]=val;}
        inline void Update(int pos,Book val) {Upt(1,n,1,pos,val),data[pos]=val;}
        inline int TotalMax(int ql,int qr,Book val) {return ql<=qr?total_max(1,n,1,ql,qr,val):0;}
        inline int TotalMin(int ql,int qr,Book val) {return ql<=qr?total_min(1,n,1,ql,qr,val):0;}
}SegmentTreap;
int main()
{
    register int i,Q,x,y,ans=0;register Book t;
    for(F.read(n),F.read(Q),SegmentTreap.Init(n),i=1;i<=n;++i) 
        F.read(s[i].pos),F.read(s[i].val),Inc(ans,SegmentTreap.TotalMax(1,i-1,s[i])),SegmentTreap.Insert(i,s[i]);
    while(Q--) 
    {
        if(F.read(x),F.read(y),x>y) swap(x,y);//如果x>y，则交换x和y
        Inc(ans,SegmentTreap.TotalMax(x+1,y,s[x])-SegmentTreap.TotalMin(x+1,y,s[x])),//更新ans
        Inc(ans,SegmentTreap.TotalMin(x+1,y-1,s[y])-SegmentTreap.TotalMax(x+1,y-1,s[y])),//更新ans，将上界减1避免重复计算
        SegmentTreap.Update(x,s[y]),SegmentTreap.Update(y,s[x]),t=s[x],s[x]=s[y],s[y]=t,//交换
        F.write(((ans%=MOD)+=MOD)%=MOD),F.write_char('\n');//输出答案，注意取模
    }
    return F.end(),0;
}

转载于:https://www.cnblogs.com/chenxiaoran666/p/Luogu3759.html

weixin_30569001
关注
洛谷P3759 [TJOI2017]不勤劳的图书管理员(二维线段树)
p_b_p_b的博客
06-01 478
裸的动态逆序对,但维护的东西很多写的一不是一般的麻烦采用二维线段树,外层区间,内层权值query的返回值是pair（或用全局变量）代码里有注释// P3759 [TJOI2017]不勤劳的图书管理员 //二维线段树 #include&lt;bits/stdc++.h&gt; #define Pli pair&lt;ll,int&gt; #define mp make_pair #define f...
[TJOI2017][bzoj4889] 不勤劳的图书管理员 [线段树套线段树]
10-20 134
题面 传送门 思路 考虑两本书的位置交换对答案的贡献： （为了方便描述，用“左边那本”和“右边那本”称呼两本我们要交换的书，“中间那本”是我们的讨论对象） 对于位置在两本书中间的书，分情况讨论： 情况1，这本书的权值在两本书中间 如果左边的书比右边的书大（也就是交换的两本书本来是逆序对，交换后变成顺序对），那么这本中间的书对答案的改变为$-(2\ast mid + left + righ...
bzoj4889: [Tjoi2017]不勤劳的图书管理员（树套树）
weixin_34026276的博客
09-01 84
传送门 据说正解线段树套平衡树 然而网上参考（抄）了一个树状数组套动态开点线段树的 思路比较清楚，看代码应该就明白了 1 //minamoto 2 #include<iostream> 3 #include<cstdio> 4 #include<cstring> 5 #define ll long long 6...
4889: [Tjoi2017]不勤劳的图书管理员 树套树
weixin_30840573的博客
04-25 102
国际惯例的题面(Bzoj没有，洛谷找的):动态加权逆序对，一眼树套树。256MB内存，5e4范围，不虚不虚。首先把交换改成两个插入和两个删除。考虑插入和删除的贡献，就是统计前面比这个值大的数的数值和，数量和，后面比这个值小的数的数值和，数量和。然后特判一下当前两个值构成逆序对的情况即可(因为这种情况会被计算两遍)。考虑树状数组套动态开点线段树维护这个东西，线段树只需要单点修改区间求和即可，十分简单...
洛谷P3758 [TJOI2017]可乐
qq_35827379的博客
02-22 241
原题链接 题目描述 加里敦星球的人们特别喜欢喝可乐。因而，他们的敌对星球研发出了一个可乐机器人，并且放在了加里敦星球的 11 号城市上。这个可乐机器人有三种行为： 停在原地，去下一个相邻的城市，自爆。它每一秒都会随机触发一种行为。现在给加里敦星球城市图，在第 00 秒时可乐机器人在 11 号城市，问经过了 tt 秒，可乐机器人的行为方案数是多少？ 输入格式 第一行输入两个正整数 NN，MM。NN 表示城市个数，MM 表示道路个数。 接下来 MM 行每行两个整数 uu，vv，表示 uu，vv 之间有一条道路。
树套树之线段树套线段树（BZOJ3110、洛谷P3332）
forezxl的博客
01-03 1140
线段树可以维护一个序列。当需要维护一个矩阵（即二维平面）内的数值，或者有什么奇怪的区间操作时，就需要用到二维线段树，也就是线段树“套上”线段树。
洛谷P3955 [NOIP2017 普及组] 图书管理员
Nico_Anzio的博客
10-08 576
题目链接 纯粹的暴力枚举,复杂度O(n2) 因为题目数据很小,所以很安全,不会超时. #include<bits/stdc++.h> using namespace std; #define N 1005 #define LL long long int book[N];//放书 int flag;//放需求码长度 int user;//放需求码 int n,m;//放书的数量和读者的数量 int main() { int i,j;//老工具人 cin>>n>&g
P3955 [NOIP2017 普及组] 图书管理员C++题解
xiaoxing06的博客
10-12 1073
P3955 [NOIP2017 普及组] 图书管理员
【bzoj4889】[Tjoi2017]不勤劳的图书管理员 树状数组+分块+二分
weixin_30632089的博客
07-06 83
题目描述（转自洛谷） 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道，如果他前i天不去整理，第i天他的厌烦度是多少，这样他好选...
BZOJ4889 & 洛谷3759：[TJOI2017]不勤劳的图书管理员——题解
weixin_34293246的博客
04-30 105
https://www.lydsy.com/JudgeOnline/problem.php?id=4889 https://www.luogu.org/problemnew/show/P3759 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下...
BZOJ4889：[TJOI2017]不勤劳的图书管理员
aga28832的博客
12-26 89
浅谈树状数组与线段树：https://www.cnblogs.com/AKMer/p/9946944.html 题目传送门：https://www.lydsy.com/JudgeOnline/problem.php?id=4889 树状数组套线段树搞一搞就好了。 时间复杂度：\(O(nlog^2n)\) 空间复杂度：\(O(nlog^2n)\) 代码如下： #include <cs...
[bzoj4889] [Tjoi2017]不勤劳的图书管理员
weixin_30782331的博客
12-20 84
Description 加里敦大学有个帝国图书馆，小豆是图书馆阅览室的一个书籍管理员。他的任务是把书排成有序的，所以无序的书让他产生厌烦，两本乱序的书会让小豆产生这两本书页数的和的厌烦度。现在有n本被打乱顺序的书，在接下来m天中每天都会因为读者的阅览导致书籍顺序改变位置。因为小豆被要求在接下来的m天中至少要整理一次图书。小豆想知道，如果他前i天不去整理，第i天他的厌烦度是多少，这样他好选择厌烦度最...
[TJOI2017]不勤劳的图书管理员 数列分块+二分
Exception
04-26 152
原题:https://loj.ac/problem/2639 题解:就是求动态逆序对的和。考虑暴力的方法，对于区间[L,R] 把[L+1.R-1]的元素和L,R比，其实就是要求[L+1,R-1]中所有大于L的和，所有小于L的和，很容易想到数列分块，把每一个块排序，二分查找求和。要注意不要用STL,注意常数优化。 #include<bits/stdc++.h> #define ll...
树套树-线段树套平衡树
ZigZagK的博客
04-15 5067
树套树-线段树套平衡树的总结。
[学习笔记]树套树
weixin_34068198的博客
11-23 548
说白了，就是在一个树形数据结构上，每个点不再是一个节点，而是另外一个树形数据结构。 空间时间复杂度大多数都是O(nlogn) 线段树套平衡树 许多树套树都可以用线段树套平衡树解决。 空间O(nlogn)是很可观的。 各种区间找值的问题，可以游刃有余解决。 （虽然常数比较大） 例如模板： 【模板】二逼平衡树（树套树） 下标线段树里套权值平衡树 线段树节点O(4*N) 平衡...
基于yolov3结合卡尔曼滤波匈牙利算法虚拟线圈的车流量检测python源码+项目说明.zip
09-21
本项目专注于视频基的车辆跟踪与流量统计，依托深度学习技术，实时追踪路面车辆动态并精准统计各车道车流量。该项目通过以下核心步骤高效处理输入视频数据： 高精度目标检测：采用YOLOv3模型，凭借其卓越的实时性与准确性，快速识别视频中的车辆目标，为后续追踪奠定坚实基础。 稳定目标追踪：集成SORT（Simple Online and Realtime Tracking）算法，结合卡尔曼滤波器预测车辆未来位置，同时运用匈牙利算法优化目标间的匹配，实现连续、稳定的车辆追踪，有效应对遮挡、交叉等复杂场景。 智能车流量统计：借鉴虚拟线圈概念，巧妙设计算法以监测穿越预设区域的车辆，自动计数并准确统计各车道车流量，为交通管理与分析提供可靠数据支持。 整体而言，本项目通过融合先进的目标检测与追踪技术，结合创新的流量统计方法，构建了一个高效、精准的实时车辆跟踪与流量统计系统，为智慧交通领域的发展贡献力量。
本科毕业设计--基于深度学习的人脸面部表情识别.zip
最新发布
09-21
深度学习是机器学习的一个子领域，它基于人工神经网络的研究，特别是利用多层次的神经网络来进行学习和模式识别。深度学习模型能够学习数据的高层次特征，这些特征对于图像和语音识别、自然语言处理、医学图像分析等应用至关重要。以下是深度学习的一些关键概念和组成部分： 1. **神经网络（Neural Networks）**：深度学习的基础是人工神经网络，它是由多个层组成的网络结构，包括输入层、隐藏层和输出层。每个层由多个神经元组成，神经元之间通过权重连接。 2. **前馈神经网络（Feedforward Neural Networks）**：这是最常见的神经网络类型，信息从输入层流向隐藏层，最终到达输出层。 3. **卷积神经网络（Convolutional Neural Networks, CNNs）**：这种网络特别适合处理具有网格结构的数据，如图像。它们使用卷积层来提取图像的特征。 4. **循环神经网络（Recurrent Neural Networks, RNNs）**：这种网络能够处理序列数据，如时间序列或自然语言，因为它们具有记忆功能，能够捕捉数据中的时间依赖性。 5. **长短期记忆网络（Long Short-Term Memory, LSTM）**：LSTM 是一种特殊的 RNN，它能够学习长期依赖关系，非常适合复杂的序列预测任务。 6. **生成对抗网络（Generative Adversarial Networks, GANs）**：由两个网络组成，一个生成器和一个判别器，它们相互竞争，生成器生成数据，判别器评估数据的真实性。 7. **深度学习框架**：如 TensorFlow、Keras、PyTorch 等，这些框架提供了构建、训练和部署深度学习模型的工具和库。 8. **激活函数（Activation Functions）**：如 ReLU、Sigmoid、Tanh 等，它们在神经网络中用于添加非线性，使得网络能够学习复杂的函数。 9. **损失函数（Loss Functions）**：用于评估模型的预测与真实值之间的差异，常见的损失函数包括均方误差（MSE）、交叉熵（Cross-Entropy）等。 10. **优化算法（Optimization Algorithms）**：如梯度下降（Gradient Descent）、随机梯度下降（SGD）、Adam 等，用于更新网络权重，以最小化损失函数。 11. **正则化（Regularization）**：技术如 Dropout、L1/L2 正则化等，用于防止模型过拟合。 12. **迁移学习（Transfer Learning）**：利用在一个任务上训练好的模型来提高另一个相关任务的性能。 深度学习在许多领域都取得了显著的成就，但它也面临着一些挑战，如对大量数据的依赖、模型的解释性差、计算资源消耗大等。研究人员正在不断探索新的方法来解决这些问题。
雷电模拟器改成真机环境
09-21
在Android系统中，模拟器通常会被识别为非真实的设备，这可能导致某些应用或服务无法正常运行，尤其是那些需要硬件信息或root权限的应用。改真机环境意味着通过技术手段使模拟器的行为更接近真实的Android设备，包括模拟设备硬件信息、IMEI、MAC地址等，以绕过检测并提升兼容性。LSPSED（Linux System Property Setter and Editor）模块是Magisk的一个扩展，专门用于修改Android系统的系统属性。通过LSPSED，我们可以更改设备的各种系统属性，如设备型号、制造商信息等，这些修改对于使模拟器看起来更像真实设备至关重要。在使用过程中，还需要添加几个机型模块，这些模块通常包含了特定设备的系统配置信息，如设备代号、硬件信息等。通过安装这些模块，雷电模拟器可以模仿各种流行的Android设备，进一步提高应用的兼容性和真实性。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值