bzoj 4399 魔法少女LJJ
4399: 魔法少女LJJ
Time Limit: 20 Sec Memory Limit: 162 MB
http://www.lydsy.com/JudgeOnline/problem.php?id=4399
Description
在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了
LJJ感叹道“这里真是个迷人的绿色世界,空气清新、淡雅,到处散发着醉人的奶浆味;小猴在枝头悠来荡去,好不自在;各式各样的鲜花争相开放,各种树枝的枝头挂满沉甸甸的野果;鸟儿的歌声婉转动听,小河里飘着落下的花瓣真是人间仙境”
SHY觉得LJJ还是太naive,一天,SHY带着自己心爱的图找到LJJ,对LJJ说:“既然你已经见识过动态树,动态仙人掌了,那么今天就来见识一下动态图吧”
LJJ:“要支持什么操作?”
SHY:“
1.新建一个节点,权值为x。
2.连接两个节点。
3.将一个节点a所属于的联通快内权值小于x的所有节点权值变成x。
4.将一个节点a所属于的联通快内权值大于x的所有节点权值变成x。
5.询问一个节点a所属于的联通块内的第k小的权值是多少。
6.询问一个节点a所属联通快内所有节点权值之积与另一个节点b所属联通快内所有节点权值之积的大小。
7.询问a所在联通快内节点的数量
8.若两个节点a,b直接相连,将这条边断开。
9.若节点a存在,将这个点删去。
”
LJJ:“我可以离线吗?”
SHY:“可以,每次操作是不加密的,”
LJJ:“我可以暴力吗?”
SHY:“自重”
LJJ很郁闷,你能帮帮他吗
Input
第一行有一个正整数m,表示操作个数。
接下来m行,每行先给出1个正整数c。
若c=1,之后一个正整数x,表示新建一个权值为x的节点,并且节点编号为n+1(当前有n个节点)。
若c=2,之后两个正整数a,b,表示在a,b之间连接一条边。
若c=3,之后两个正整数a,x,表示a联通快内原本权值小于x的节点全部变成x。
若c=4,之后两个正整数a,x,表示a联通快内原本权值大于x的节点全部变成x。
若c=5,之后两个正整数a,k,表示询问a所属于的联通块内的第k小的权值是多少。
若c=6,之后两个正整数a,b,表示询问a所属联通快内所有节点权值之积与b所属联通快内所有节点权值之积的大小,
若a所属联通快内所有节点权值之积大于b所属联通快内所有节点权值之积,输出1,否则为0。
若c=7,之后一个正整数a,表示询问a所在联通块大小
若c=8,之后两个正整数a,b,表示断开a,b所连接的边。
若c=9,之后一个正整数a,表示断开a点的所有连边
具体输出格式见样例
Output
Sample Input
1 2
1 3
1 4
1 5
1 6
2 1 2
2 2 3
2 3 4
2 4 5
9 1
3 2 5
5 3 4
Sample Output
HINT
对100%的数据 0<=m<=400000,c<=7,所有出现的数均<=1000000000,所有出现的点保证存在
【HINT】请认真阅读题面
Source
By liuchenrui
动态开点线段树+启发式合并
将所有输入数据读入后离散化
添加操作 在从上往下找的过程中顺带加上信息 ,而不是找到之后在up date会快很多
合并操作 注意合并前先判断是否在同一联通块内
修改操作 暴力一个点一个点删除后单点加,
删除的时候注意当删除区间没有数的时候及时return,否则TLE
比较乘积大小操作 利用log(a*b)=log(a)+log(b),取对数后变为比较和
再就是这题同线段树注意注意空间
#include<cmath> #include<cstdio> #include<algorithm> #define N 400001 using namespace std; int n,m,id,tot,s; int z,r1,r2,pos; int fa[N],root[N]; int data[N][3]; int hash[N],cnt; double logg[N]; struct node {int l,r,size;double sum; }tr[7000001]; int read() {int x=0,f=1; char c=getchar();while(c<'0'||c>'9') { if(c=='-') f=-1; c=getchar(); }while(c>='0'&&c<='9') { x=x*10+c-'0'; c=getchar(); }return x*f; } int find(int i) { return fa[i]==i ? i :fa[i]=find(fa[i]);} void add(int &k,int l,int r,int w) {if(!k) k=++tot;tr[k].size+=w; tr[k].sum+=w*logg[pos];if(l==r) return; int mid=l+r>>1;if(pos<=mid) add(tr[k].l,l,mid,w);else add(tr[k].r,mid+1,r,w); } void up(int x) {tr[x].size=tr[tr[x].l].size+tr[tr[x].r].size;tr[x].sum=tr[tr[x].l].sum+tr[tr[x].r].sum; } int merge(int x,int y,int l,int r) {if(!x) return y;if(!y) return x;if(l==r) {tr[x].size+=tr[y].size;tr[x].sum+=tr[y].sum;return x;}int mid=l+r>>1;tr[x].l=merge(tr[x].l,tr[y].l,l,mid);tr[x].r=merge(tr[x].r,tr[y].r,mid+1,r);up(x);return x; } int query(int x,int l,int r,int k) {if(l==r) return l;int mid=l+r>>1,t=tr[tr[x].l].size;if(t>=k) return query(tr[x].l,l,mid,k);else return query(tr[x].r,mid+1,r,k-t); } void del(int x,int l,int r,int opl,int opr) {if(!tr[x].size) return;if(l==r){s+=tr[x].size;tr[x].size=tr[x].sum=0;return;}int mid=l+r>>1;if(opl<=mid) del(tr[x].l,l,mid,opl,opr);if(opr>mid) del(tr[x].r,mid+1,r,opl,opr);up(x); } int main() {freopen("data","r",stdin);freopen("std.out","w",stdout);m=read();for(int i=1;i<=m;i++){scanf("%d",&z);data[i][0]=z;data[i][1]=read();if(z!=7&&z!=1) data[i][2]=read();if(z==1) hash[++cnt]=data[i][1],data[i][2]=++id;else if(z==3||z==4) hash[++cnt]=data[i][2];}sort(hash+1,hash+cnt+1);cnt=unique(hash+1,hash+cnt+1)-(hash+1);for(int i=1;i<=m;i++)switch(data[i][0]){case 1:{pos=lower_bound(hash+1,hash+cnt+1,data[i][1])-hash;logg[pos]=log(data[i][1]);add(root[data[i][2]],1,cnt,1); fa[data[i][2]]=data[i][2];break;}case 2:{r1=find(data[i][1]); r2=find(data[i][2]);if(r1!=r2){fa[r2]=r1;merge(root[r1],root[r2],1,cnt);}break;}case 3:{s=0; r1=find(data[i][1]);pos=lower_bound(hash+1,hash+cnt+1,data[i][2])-hash;logg[pos]=log(data[i][2]);del(root[r1],1,cnt,1,max(1,pos-1));if(s) add(root[r1],1,cnt,s);break;}case 4:{s=0; r1=find(data[i][1]);pos=lower_bound(hash+1,hash+cnt+1,data[i][2])-hash;logg[pos]=log(data[i][2]);del(root[r1],1,cnt,min(pos+1,cnt),cnt);if(s) add(root[r1],1,cnt,s);break; }case 5:{r1=find(data[i][1]);printf("%d\n",hash[query(root[r1],1,cnt,data[i][2])]);break;}case 6:{r1=find(data[i][1]); r2=find(data[i][2]);if(tr[root[r1]].sum>tr[root[r2]].sum) puts("1");else puts("0");break;}case 7:{r1=find(data[i][1]);printf("%d\n",tr[root[r1]].size);break;}} }
神奇TLE1:然后删除操作改成了一段区间一起删除,就TLE了,原因未知
难道是因为动态开点的线段树标记下传会有影响?
#include<cmath> #include<cstdio> #include<algorithm> #define N 400001 using namespace std; int n,m,id,tot,s; int z,r1,r2,pos; int fa[N],root[N]; int data[N][3]; int hash[N],cnt; double logg[N]; struct node {int l,r,size;bool f;double sum; }tr[6400001]; int read() {int x=0,f=1; char c=getchar();while(c<'0'||c>'9') { if(c=='-') f=-1; c=getchar(); }while(c>='0'&&c<='9') { x=x*10+c-'0'; c=getchar(); }return x*f; } int find(int i) { return fa[i]==i ? i :fa[i]=find(fa[i]);} void down(int k) {tr[tr[k].l].size=tr[tr[k].r].size=0;tr[tr[k].l].f=tr[tr[k].r].f=true;tr[k].f=0; } void add(int &k,int l,int r,int w) {if(!k) k=++tot;tr[k].size+=w; tr[k].sum+=w*logg[pos];if(l==r) return; if(tr[k].f) down(k);int mid=l+r>>1;if(pos<=mid) add(tr[k].l,l,mid,w);else add(tr[k].r,mid+1,r,w); } void up(int x) {tr[x].size=tr[tr[x].l].size+tr[tr[x].r].size;tr[x].sum=tr[tr[x].l].sum+tr[tr[x].r].sum; } int merge(int x,int y,int l,int r) {if(!x) return y;if(!y) return x;if(l==r) {tr[x].size+=tr[y].size;tr[x].sum+=tr[y].sum;return x;}if(tr[x].f) down(x);if(tr[y].f) down(y);int mid=l+r>>1;tr[x].l=merge(tr[x].l,tr[y].l,l,mid);tr[x].r=merge(tr[x].r,tr[y].r,mid+1,r);up(x);return x; } int query(int x,int l,int r,int k) {if(l==r) return l;if(tr[x].f) down(x);int mid=l+r>>1,t=tr[tr[x].l].size;if(t>=k&&tr[tr[x].l].size) return query(tr[x].l,l,mid,k);else if(tr[tr[x].r].size) return query(tr[x].r,mid+1,r,k-t); } void del(int x,int l,int r,int opl,int opr) {if(l>=opl&&r<=opr){s+=tr[x].size;tr[x].size=tr[x].sum=0;tr[x].f=1;return;}if(tr[x].f) down(x);int mid=l+r>>1;if(opl<=mid&&tr[tr[x].l].size) del(tr[x].l,l,mid,opl,opr);if(opr>mid&&tr[tr[x].r].size) del(tr[x].r,mid+1,r,opl,opr);up(x); } int main() {freopen("data","r",stdin);freopen("my.out","w",stdout);m=read();for(int i=1;i<=m;i++){scanf("%d",&z);data[i][0]=z;data[i][1]=read();if(z!=7&&z!=1) data[i][2]=read();if(z==1) hash[++cnt]=data[i][1],data[i][2]=++id;else if(z==3||z==4) hash[++cnt]=data[i][2];}sort(hash+1,hash+cnt+1);cnt=unique(hash+1,hash+cnt+1)-(hash+1);for(int i=1;i<=m;i++)switch(data[i][0]){case 1:{pos=lower_bound(hash+1,hash+cnt+1,data[i][1])-hash;logg[pos]=log(data[i][1]);add(root[data[i][2]],1,cnt,1); fa[data[i][2]]=data[i][2];break;}case 2:{r1=find(data[i][1]); r2=find(data[i][2]);if(r1!=r2){fa[r2]=r1;merge(root[r1],root[r2],1,cnt);}break;}case 3:{s=0; r1=find(data[i][1]);pos=lower_bound(hash+1,hash+cnt+1,data[i][2])-hash;logg[pos]=log(data[i][2]);del(root[r1],1,cnt,1,max(1,pos-1));if(s) add(root[r1],1,cnt,s);break;}case 4:{s=0; r1=find(data[i][1]);pos=lower_bound(hash+1,hash+cnt+1,data[i][2])-hash;logg[pos]=log(data[i][2]);del(root[r1],1,cnt,min(pos+1,cnt),cnt);if(s) add(root[r1],1,cnt,s);break; }case 5:{r1=find(data[i][1]);printf("%d\n",hash[query(root[r1],1,cnt,data[i][2])]);break;}case 6:{r1=find(data[i][1]); r2=find(data[i][2]);if(tr[root[r1]].sum>tr[root[r2]].sum) puts("1");else puts("0");break;}case 7:{r1=find(data[i][1]);printf("%d\n",tr[root[r1]].size);break;}} }
View Code
神奇TLE2:这样合并:难道多一层就T?
int merge(int x,int y) {if(!x) return y;if(!y) return x;tr[x].l=merge(tr[x].l,tr[y].r);tr[x].r=merge(tr[x].r,tr[y].r);up(x);return x; }
转载于:https://www.cnblogs.com/TheRoadToTheGold/p/6820749.html
bzoj 4399 魔法少女LJJ相关推荐
- BZOJ 4399: 魔法少女LJJ
4399: 魔法少女LJJ Time Limit: 20 Sec Memory Limit: 162 MB Submit: 287 Solved: 73 Description 在森林中见过会动的树, ...
- BZOJ 4399: 魔法少女LJJ 线段树合并 + 对数
Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ感叹道"这里真是个迷人的绿色世界,空气清新.淡雅 ...
- 4399: 魔法少女LJJ
注意他这个题的题面傻逼.实际上只有7种操作.可以发现是道裸的线段树合并.只是那个乘积会爆,要取log. #include<cstdio> #include<iostream> ...
- bzoj4399 魔法少女LJJ
这题不愧为4399. 4399: 魔法少女LJJ Time Limit: 20 Sec Memory Limit: 162 MB Submit: 622 Solved: 145 [Submit][ ...
- 【BZOJ 4399】魔法少女LJJ
[题目] 传送门 题目描述: 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女 LJJ 已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ 感叹道"""这里 ...
- 【bzoj4399】魔法少女LJJ 并查集+权值线段树合并
题目描述 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ感叹道"这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着醉 ...
- 魔法少女 LJJ——线段树
题目 [题目描述] 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女 LJJ 已经觉得自己见过世界上的所有稀奇古怪的事情了. LJJ 感叹道"这里真是个迷人的绿色世界,空气清新. ...
- BZOJ4399魔法少女LJJ——线段树合并+并查集
题目描述 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ感叹道"这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着醉 ...
- BZOJ4399 魔法少女LJJ【线段树合并】【并查集】
Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ感叹道"这里真是个迷人的绿色世界,空气清新.淡雅 ...
最新文章
- 数学老师必备工具,你的最爱!
- 成功解决from scipy.linalg import _fblas ImportError: DLL load failed: 找不到指定的模块。
- 如何对 string 进行Base64编码,解码?
- eclipse使用技巧---使用正则表达式查找替换
- 恒大汽车:仍在就出售新能源汽车生活项目等资产进行磋商
- SAP License:CO常见问题
- numpy tricks(二)—— 删除多维数组的行或列
- git源代码控制客户端基本配置
- Python由浅入深编程(概念)实战
- 简单C语言小程序:求根公式求一元二次方程式的根!
- java红包正态分布,微信红包算法探讨
- Android之利用NanoHttpd搭建服务器
- qlistview 自定义控件_QT中QListView中放置自定义控件并添加滚动条
- android toast防重_安卓Toast自定义及防止重复显示
- 考研英语 常见不规则动词过去式/过去分词
- 阿里云服务器迁移内容
- Bias和Variance
- Excel PivotTable 使用心得手顺分享(五)
- ffmpeg C++推流
- 【新知实验室】腾讯云音视频应用
热门文章
- 怎么把不能编辑的qlv转换成mp4
- anemometer mysql 500_【Mysql】Anemometer-可视化慢sql查看
- 15Echarts:复杂数据展示
- 麒麟系统arm64或aarch64架构下docker打包jdk1.8镜像
- 30个漂亮的网页设计案例(组图)
- 世界崩塌前要写下游戏心得(持续更新)
- html 可脱机浏览,如何脱机浏览Web页面
- 微信小程序生命周期(onLoad,onHide,onShow,onReady,onUnload)
- 设置webhook_数据采集教程_智能模式_如何设置Webhook功能_后羿采集器
- 挖财推出Android6.0版,大幅度减少跳转页面