题意

https://www.lydsy.com/JudgeOnline/problem.php?id=4399

思路

码农题,需要一定代码功底。方法很暴力,先将权值离散,表示在线段树里储存的位置,每个连通块用一棵动点线段树存储,合并两个连通块直接对两个动点线段树进行合并,查询操作在当前连通块的线段树上进行,只不过有询问乘积大小,直接权值取原权值的 \(\ln​\) ,比较和的大小即可。

现在分析线段树合并的复杂度,举一个最基本的例子:权值为\([1,n]\) ,\(n\) 棵动点线段树,每个线段树插入了一个权值,那么总共有 \(n\log n\) 个点,而每一次合并相当于少掉了一个点,那么合并完这 \(n\) 棵线段树后复杂度就是消失的点的个数,不会超过总共的点数,所以复杂度是 \(n \log n\) 的。

类似的,对最一般的情况,也是一样的分析方法,最后得出线段树合并的复杂度为节点个数的结论。

代码

#include<bits/stdc++.h>
#define FOR(i,x,y) for(int i=(x),i##END=(y);i<=i##END;++i)
#define DOR(i,x,y) for(int i=(x),i##END=(y);i>=i##END;--i)
typedef long long LL;
using namespace std;
const int N=4e5+5;
const int NN=N*12; struct SegmentTree
{
struct node
{
int cnt;double sum;
node operator +(const node &_)const
{
return (node){cnt+_.cnt,sum+_.sum};
}
};
node nd[NN];
int lson[NN],rson[NN],rt[N],tot;
int &operator [](const int x){return rt[x];}
void build()
{
memset(rt,0,sizeof(rt));
nd[tot=0]=(node){0,0};
lson[0]=rson[0]=0;
}
void create(int &k)
{
if(!k)k=++tot,nd[k]=nd[0],lson[k]=rson[k]=0;
}
void update(int &k,int x,int addcnt,double addsum,int l,int r)
{
create(k);
if(l==r)
{
nd[k].cnt+=addcnt;
nd[k].sum+=addsum;
return;
}
int mid=(l+r)>>1;
if(x<=mid)update(lson[k],x,addcnt,addsum,l,mid);
else update(rson[k],x,addcnt,addsum,mid+1,r);
nd[k]=nd[lson[k]]+nd[rson[k]];
}
int sweep_off(int &k,int L,int R,int l,int r)
{
if(!k)return 0;
if(L<=l&&r<=R){int res=nd[k].cnt;k=0;return res;}
int mid=(l+r)>>1,res=0;
if(L<=mid)res+=sweep_off(lson[k],L,R,l,mid);
if(R>mid)res+=sweep_off(rson[k],L,R,mid+1,r);
nd[k]=nd[lson[k]]+nd[rson[k]];
return res;
}
int querycnt(int k){return nd[k].cnt;}
double querysum(int k){return nd[k].sum;}
int queryKth(int k,int K,int l,int r)
{
if(l==r)return l;
int mid=(l+r)>>1;
if(nd[lson[k]].cnt>=K)return queryKth(lson[k],K,l,mid);
else return queryKth(rson[k],K-nd[lson[k]].cnt,mid+1,r);
}
void merge(int &x,int y,int l,int r)
{
if(!x||!y){x=(x|y);return;}
if(l==r){nd[x]=nd[x]+nd[y];return;}
int mid=(l+r)>>1;
merge(lson[x],lson[y],l,mid);
merge(rson[x],rson[y],mid+1,r);
nd[x]=nd[lson[x]]+nd[rson[x]];
}
}ST;
int n,m,fa[N];
int op[N],a[N],b[N];
int disc[N],tot;
double logdisc[N];
int getfa(int k){return k==fa[k]?k:fa[k]=getfa(fa[k]);} int main()
{
scanf("%d",&m);
FOR(i,1,m)
{
scanf("%d%d",&op[i],&a[i]);
if(op[i]!=1&&op[i]!=7)scanf("%d",&b[i]);
} FOR(i,1,m)
{
if(op[i]==1)disc[++tot]=a[i];
else if(op[i]==3||op[i]==4)disc[++tot]=b[i];
}
sort(disc+1,disc+1+tot);
tot=unique(disc+1,disc+1+tot)-disc-1;
FOR(i,1,m)
{
if(op[i]==1)a[i]=lower_bound(disc+1,disc+1+tot,a[i])-disc;
else if(op[i]==3||op[i]==4)b[i]=lower_bound(disc+1,disc+1+tot,b[i])-disc;
}
FOR(i,1,tot)logdisc[i]=log(disc[i]); ST.build();
FOR(i,1,m)
{
if(op[i]==1)
{
n++;
fa[n]=n;
ST.update(ST[n],a[i],1,logdisc[a[i]],1,tot);
}
else if(op[i]==2)
{
a[i]=getfa(a[i]),b[i]=getfa(b[i]);
if(a[i]==b[i])continue;
ST.merge(ST[a[i]],ST[b[i]],1,tot);
fa[b[i]]=a[i];
}
else if(op[i]==3)
{
if(b[i]==1)continue;
a[i]=getfa(a[i]);
int cnt=ST.sweep_off(ST[a[i]],1,b[i]-1,1,tot);
ST.update(ST[a[i]],b[i],cnt,cnt*logdisc[b[i]],1,tot);
}
else if(op[i]==4)
{
if(b[i]==tot)continue;
a[i]=getfa(a[i]);
int cnt=ST.sweep_off(ST[a[i]],b[i]+1,tot,1,tot);
ST.update(ST[a[i]],b[i],cnt,cnt*logdisc[b[i]],1,tot);
}
else if(op[i]==5)
{
a[i]=getfa(a[i]);
printf("%d\n",disc[ST.queryKth(ST[a[i]],b[i],1,tot)]);
}
else if(op[i]==6)
{
a[i]=getfa(a[i]),b[i]=getfa(b[i]);
if(ST.querysum(ST[a[i]])>ST.querysum(ST[b[i]]))
puts("1");
else puts("0");
}
else if(op[i]==7)
{
a[i]=getfa(a[i]);
printf("%d\n",ST.querycnt(ST[a[i]]));
}
} return 0;
}

BZOJ 4399 魔法少女LJJ(线段树合并)的更多相关文章

  1. BZOJ.4399.魔法少女LJJ(线段树合并)

    BZOJ 注意\(c\leq7\)→_→ 然后就是裸的权值线段树+线段树合并了. 对于取\(\max/\min\)操作可以直接区间修改清空超出范围的值,然后更新到对应位置上就行了(比如对\(v\)取\ ...

  2. BZOJ 4399: 魔法少女LJJ 线段树合并 + 对数

    Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着 ...

  3. BZOJ 4399: 魔法少女LJJ(线段树)

    传送门 解题思路 出题人真会玩..操作\(2\)线段树合并,然后每棵线段树维护元素个数和.对于\(6\)这个询问,因为乘积太大,所以要用对数.时间复杂度\(O(nlogn)\) 代码 #include ...

  4. BZOJ4399魔法少女LJJ——线段树合并+并查集

    题目描述 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着醉人的奶浆味: ...

  5. 【BZOJ4399】魔法少女LJJ 线段树合并

    [BZOJ4399]魔法少女LJJ Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的 ...

  6. bzoj4399 魔法少女LJJ 线段树合并

    只看题面绝对做不出系列.... 注意到\(c \leqslant 7\),因此不会有删边操作(那样例删边干嘛) 注意到\(2, 5\)操作十分的有趣,启示我们拿线段树合并来做 操作\(7\)很好处理 ...

  7. bzoj4399 魔法少女LJJ 线段树合并+线段树二分+并查集

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4399 题解 毒瘤题 \(9\) 种操作还有支持动态图的连通性 仔细读题 $ c<=7$. ...

  8. bzoj 4399 魔法少女LJJ

    4399: 魔法少女LJJ Time Limit: 20 Sec  Memory Limit: 162 MBhttp://www.lydsy.com/JudgeOnline/problem.php?i ...

  9. 魔法少女 LJJ——线段树

    题目 [题目描述] 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女 LJJ 已经觉得自己见过世界上的所有稀奇古怪的事情了. LJJ 感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处 ...

随机推荐

  1. .net nancy

    官网 文档 入门教程 参考

  2. npm 包下载的各种姿势

    最近在写Node程序的时候,突然对 npm install 的-save和-save-dev 这两个参数的使用比较混乱.其实博主在这之前对这两个参数的理解也是模糊的,各种查资料和实践后对它们之间的异同 ...

  3. 微信小程序制作家庭记账本之五

    第五天,开发进度没有进行,不知道该如何下手完成小程序的开发.

  4. RN与android原生开发混合后的环境报错问题

    RN与android原生开发混合后的环境报错问题 需要先安装nodejs$ yarn --version1.12.1更新当前版本yarn upgrade --latest安装 | Yarnhttps: ...

  5. Python+OpenCV图像处理(十)—— 图像二值化

    简介:图像二值化就是将图像上的像素点的灰度值设置为0或255,也就是将整个图像呈现出明显的黑白效果的过程. 一.普通图像二值化 代码如下: import cv2 as cv import numpy ...

  6. jquery 页面分页的实现

      <!DOCTYPE html> <html> <head> <title>分页</title> <link rel="s ...

  7. 怎样从外网访问内网Jupyter Notebook?

    本地安装了一个Jupyter Notebook,只能在局域网内访问,怎样从外网也能访问到本地的Jupyter Notebook呢?本文将介绍具体的实现步骤. 准备工作 安装并启动Jupyter Not ...

  8. 怎样从外网访问内网IIS?

    本地安装了一个IIS,只能在局域网内访问,怎样从外网也能访问到本地的IIS呢?本文将介绍具体的实现步骤. 准备工作 安装并启动IIS 默认安装的IIS端口是80. 实现步骤 下载并解压holer软件包 ...

  9. Logstash 安装和使用

    下载地址 https://artifacts.elastic.co/downloads/logstash/logstash-5.6.8.zip 下载后解压,测试 #将键盘内容输出到控制台 logsta ...

  10. Java连接数据库 #01# JDBC单线程适用

    官方教程(包括 javase的基础部分):JDBC Basics 重新梳理.学习一下“Java连接数据库”相关的内容. 因为最开始没有认真学多线程和JDBC,一直在自己写的多线程程序中维持下面的错误写 ...