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
动态开点线段树+启发式合并
将所有输入数据读入后离散化
添加操作 在从上往下找的过程中顺带加上信息 ，而不是找到之后在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][];
int hash[N],cnt;
double logg[N];
struct node
{
int l,r,size;
double sum;
}tr[];
int read()
{
int x=,f=; char c=getchar();
while(c<''||c>'') { if(c=='-') f=-; c=getchar(); }
while(c>=''&&c<='') { x=x*+c-''; 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>>;
if(pos<=mid) add(tr[k].l,l,mid,w);
else add(tr[k].r,mid+,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>>;
tr[x].l=merge(tr[x].l,tr[y].l,l,mid);
tr[x].r=merge(tr[x].r,tr[y].r,mid+,r);
up(x);
return x;
}
int query(int x,int l,int r,int k)
{
if(l==r) return l;
int mid=l+r>>,t=tr[tr[x].l].size;
if(t>=k) return query(tr[x].l,l,mid,k);
else return query(tr[x].r,mid+,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=;
return;
}
int mid=l+r>>;
if(opl<=mid) del(tr[x].l,l,mid,opl,opr);
if(opr>mid) del(tr[x].r,mid+,r,opl,opr);
up(x);
}
int main()
{
freopen("data","r",stdin);
freopen("std.out","w",stdout);
m=read();
for(int i=;i<=m;i++)
{
scanf("%d",&z);
data[i][]=z;
data[i][]=read();
if(z!=&&z!=) data[i][]=read();
if(z==) hash[++cnt]=data[i][],data[i][]=++id;
else if(z==||z==) hash[++cnt]=data[i][];
}
sort(hash+,hash+cnt+);
cnt=unique(hash+,hash+cnt+)-(hash+);
for(int i=;i<=m;i++)
switch(data[i][])
{ case :
{
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
add(root[data[i][]],,cnt,);
fa[data[i][]]=data[i][];
break;
}
case :
{
r1=find(data[i][]); r2=find(data[i][]);
if(r1!=r2)
{
fa[r2]=r1;
merge(root[r1],root[r2],,cnt);
}
break;
}
case :
{
s=; r1=find(data[i][]);
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
del(root[r1],,cnt,,max(,pos-));
if(s) add(root[r1],,cnt,s);
break;
}
case :
{
s=; r1=find(data[i][]);
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
del(root[r1],,cnt,min(pos+,cnt),cnt);
if(s) add(root[r1],,cnt,s);
break;
}
case :
{
r1=find(data[i][]);
printf("%d\n",hash[query(root[r1],,cnt,data[i][])]);
break;
}
case :
{
r1=find(data[i][]); r2=find(data[i][]);
if(tr[root[r1]].sum>tr[root[r2]].sum) puts("");
else puts("");
break;
}
case :
{
r1=find(data[i][]);
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][];
int hash[N],cnt;
double logg[N];
struct node
{
int l,r,size;
bool f;
double sum;
}tr[];
int read()
{
int x=,f=; char c=getchar();
while(c<''||c>'') { if(c=='-') f=-; c=getchar(); }
while(c>=''&&c<='') { x=x*+c-''; 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=;
tr[tr[k].l].f=tr[tr[k].r].f=true;
tr[k].f=;
}
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>>;
if(pos<=mid) add(tr[k].l,l,mid,w);
else add(tr[k].r,mid+,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>>;
tr[x].l=merge(tr[x].l,tr[y].l,l,mid);
tr[x].r=merge(tr[x].r,tr[y].r,mid+,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>>,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+,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=;
tr[x].f=;
return;
}
if(tr[x].f) down(x);
int mid=l+r>>;
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+,r,opl,opr);
up(x);
}
int main()
{
freopen("data","r",stdin);
freopen("my.out","w",stdout);
m=read();
for(int i=;i<=m;i++)
{
scanf("%d",&z);
data[i][]=z;
data[i][]=read();
if(z!=&&z!=) data[i][]=read();
if(z==) hash[++cnt]=data[i][],data[i][]=++id;
else if(z==||z==) hash[++cnt]=data[i][];
}
sort(hash+,hash+cnt+);
cnt=unique(hash+,hash+cnt+)-(hash+);
for(int i=;i<=m;i++)
switch(data[i][])
{
case :
{
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
add(root[data[i][]],,cnt,);
fa[data[i][]]=data[i][];
break;
}
case :
{
r1=find(data[i][]); r2=find(data[i][]);
if(r1!=r2)
{
fa[r2]=r1;
merge(root[r1],root[r2],,cnt);
}
break;
}
case :
{
s=; r1=find(data[i][]);
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
del(root[r1],,cnt,,max(,pos-));
if(s) add(root[r1],,cnt,s);
break;
}
case :
{
s=; r1=find(data[i][]);
pos=lower_bound(hash+,hash+cnt+,data[i][])-hash;
logg[pos]=log(data[i][]);
del(root[r1],,cnt,min(pos+,cnt),cnt);
if(s) add(root[r1],,cnt,s);
break;
}
case :
{
r1=find(data[i][]);
printf("%d\n",hash[query(root[r1],,cnt,data[i][])]);
break;
}
case :
{
r1=find(data[i][]); r2=find(data[i][]);
if(tr[root[r1]].sum>tr[root[r2]].sum) puts("");
else puts("");
break;
}
case :
{
r1=find(data[i][]);
printf("%d\n",tr[root[r1]].size);
break;
}
}
}
神奇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;
}
bzoj 4399 魔法少女LJJ的更多相关文章
- BZOJ 4399: 魔法少女LJJ 线段树合并 + 对数
Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着 ...
- BZOJ.4399.魔法少女LJJ(线段树合并)
BZOJ 注意\(c\leq7\)→_→ 然后就是裸的权值线段树+线段树合并了. 对于取\(\max/\min\)操作可以直接区间修改清空超出范围的值,然后更新到对应位置上就行了(比如对\(v\)取\ ...
- BZOJ 4399 魔法少女LJJ（线段树合并）
题意 https://www.lydsy.com/JudgeOnline/problem.php?id=4399 思路 码农题,需要一定代码功底.方法很暴力,先将权值离散,表示在线段树里储存的位置,每 ...
- BZOJ 4399: 魔法少女LJJ(线段树)
传送门 解题思路 出题人真会玩..操作\(2\)线段树合并,然后每棵线段树维护元素个数和.对于\(6\)这个询问,因为乘积太大,所以要用对数.时间复杂度\(O(nlogn)\) 代码 #include ...
- AC日记——魔法少女LJJ bzoj 4399
魔法少女LJJ 思路: 动态开点权值线段树+启发式合并: 来,上代码: #include <cmath> #include <cstdio> #include <cstr ...
- BZOJ4399魔法少女LJJ——线段树合并+并查集
题目描述 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着醉人的奶浆味: ...
- BZOJ4399 魔法少女LJJ【线段树合并】【并查集】
Description 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了 LJJ感叹道"这里真是个迷人的绿色世界,空气清新.淡雅 ...
- 魔法少女 LJJ——线段树
题目 [题目描述] 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女 LJJ 已经觉得自己见过世界上的所有稀奇古怪的事情了. LJJ 感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处 ...
- 【bzoj4399】魔法少女LJJ 并查集+权值线段树合并
题目描述 在森林中见过会动的树,在沙漠中见过会动的仙人掌过后,魔法少女LJJ已经觉得自己见过世界上的所有稀奇古怪的事情了LJJ感叹道“这里真是个迷人的绿色世界,空气清新.淡雅,到处散发着醉人的奶浆味: ...
随机推荐
- 【BZOJ2006】超级钢琴（主席树，优先队列）
[BZOJ2006]超级钢琴(主席树,优先队列) 题面 BZOJ 题解 既然是一段区间 首先就要变成单点 所以求一个前缀和 这个时候贪心很明显了: 枚举每一个点和可以和它组成一段的可行的点 全部丢进一 ...
- 【THUWC 2017】在美妙的数学王国中畅游
数学王国里有n座城市,每座城市有三个参数\(f\),\(a\),\(b\),一个智商为\(x\)的人经过一座城市的获益\(f(x)\)是 若\(f=1\),则\(f(x)=\sin(ax+b)\): ...
- Java interview questions(No1)
1.什么是构造和析构方法?功能是? 答: 构造方法: 每个类至少有一个构造方法,类初始化时调用的方法 1.方法名和类名相同 2.无返回值类型 格式:访问权限 类名(参数列表) {}; 1.自己定义构造 ...
- 22.jQuery(实例)
1.开关灯效果 <!DOCTYPE html> <html lang="en"> <head> <meta charset="U ...
- LeetCode之Easy篇 ——（1）Two Sum
1.Two Sum Given an array of integers, return indices of the two numbers such that they add up to a s ...
- ISAPI和CGI限制中没有ASP.NET v4.0
[服务器搭建]ISAPI和CGI限制中没有ASP.NET v4.0解决方式: 1.确保安装IIS时确实安装了ASP.NET,如果没有的话,勾上重新装一下,一般出现404.2时这么干 2.如果你是先装了 ...
- CentOS7安装最新版git教程
下载编译工具 yum -y groupinstall "Development Tools" 下载依赖包 yum -y install zlib-devel perl-ExtUti ...
- Django+xadmin打造在线教育平台（二）
三.xadmin后台管理 3.1.xadmin的安装 django2.0的安装(源码安装方式): https://github.com/sshwsfc/xadmin/tree/django2 把zip ...
- Unity3D项目程序加密-VirboxProtector加壳工具
各位Unity3D的开发者,你还为你的代码被反编译而头疼, 混淆和加密已经失效,为内存dump代码而烦恼?是否辛苦制作的游戏被盗版被抄袭而烦恼? 是否害怕算法被别人参考要把算法写成C++而费劲周折? ...
- business expressions(一）
一. 1.24/7 24/7 :twenty four hours a day, seven days a week I work 24/7.——meaing i work all the time. ...