BZOJ 4129 树上带修莫队+线段树
思路:
可以先做做BZOJ3585 是序列上的mex
考虑莫队的转移 如果当前数字出现过 线段树上把它置成1
对于询问 二分ans 线段树上查 0到ans的和 是不是ans+1
本题就是把它搞到了序列上 带了个修改…
麻烦一点 本质上是一样的
//By SiriusRen
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=55555;
int n,m,cnt=1,Block,block[N],op,xx,yy,a[N],vis[N],last[N],sum[N],tree[N*8];
int first[N],next[N*2],v[N*2],tot,s[N],top,fa[N][20],cnt1,cnt2,deep[N],Ans[N];
void add(int x,int y){v[tot]=y,next[tot]=first[x],first[x]=tot++;}
void dfs(int x){
for(int i=first[x];~i;i=next[i])
if(v[i]!=fa[x][0])
deep[v[i]]=deep[x]+1,fa[v[i]][0]=x,dfs(v[i]);
s[++top]=x;
if(top==Block){
for(int i=1;i<=top;i++)block[s[i]]=cnt;
cnt++,top=0;
}
}
struct Ask{
int time,id,l,r;Ask(){}
Ask(int T,int I,int L,int R){time=T,id=I,l=L,r=R;}
friend bool operator<(Ask a,Ask b){
if(block[a.l]==block[b.l]){
if(block[a.r]==block[b.r])return a.time<b.time;
return block[a.r]<block[b.r];
}return block[a.l]<block[b.l];
}
}ask[N];
struct Change{
int last,position,num;Change(){}
Change(int L,int P,int nn){last=L,position=P,num=nn;}
}change[N];
void insert(int l,int r,int pos,int num,int f){
if(l==r){tree[pos]=f;return;}
int mid=(l+r)>>1,lson=pos<<1,rson=pos<<1|1;
if(mid<num)insert(mid+1,r,rson,num,f);
else insert(l,mid,lson,num,f);
tree[pos]=tree[lson]+tree[rson];
}
int query(int l,int r,int pos,int L,int R){
if(l>=L&&r<=R)return tree[pos];
int mid=(l+r)>>1,lson=pos<<1,rson=pos<<1|1;
if(mid<L)return query(mid+1,r,rson,L,R);
else if(mid>=R)return query(l,mid,lson,L,R);
else return query(l,mid,lson,L,R)+query(mid+1,r,rson,L,R);
}
void reverse(int x){
vis[x]^=1;
if(a[x]<=n){
if(vis[x]){
sum[a[x]]++;
if(sum[a[x]]==1)insert(0,n,1,a[x],1);
}
else{
sum[a[x]]--;
if(!sum[a[x]])insert(0,n,1,a[x],0);
}
}
}
void change_color(int x,int y){
if(vis[x])reverse(x),a[x]=y,reverse(x);
else a[x]=y;
}
void work(int x,int y){
while(x!=y){
if(deep[x]<deep[y])swap(x,y);
reverse(x),x=fa[x][0];
}
}
int lca(int x,int y){
if(deep[x]<deep[y])swap(x,y);
for(int i=19;i>=0;i--)if(deep[x]-(1<<i)>=deep[y])x=fa[x][i];
if(x==y)return x;
for(int i=19;i>=0;i--)if(fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i];
return fa[x][0];
}
bool check(int x){
return query(0,n,1,0,x)==x+1;
}
int main(){
memset(first,-1,sizeof(first));
scanf("%d%d",&n,&m),Block=sqrt(n);
for(int i=1;i<=n;i++)scanf("%d",&a[i]),last[i]=a[i];
for(int i=1;i<n;i++)scanf("%d%d",&xx,&yy),add(xx,yy),add(yy,xx);
dfs(1);
for(int i=1;i<=top;i++)block[s[i]]=cnt;
for(int j=1;j<=19;j++)
for(int i=1;i<=n;i++)
fa[i][j]=fa[fa[i][j-1]][j-1];
while(m--){
scanf("%d%d%d",&op,&xx,&yy);
if(op==1){
if(block[xx]>block[yy])swap(xx,yy);
ask[++cnt1]=Ask(cnt2,cnt1,xx,yy);
}
else change[++cnt2]=Change(last[xx],xx,yy),last[xx]=yy;
}
sort(ask+1,ask+1+cnt1);
for(int i=1,T=0;i<=cnt1;i++){
for(;T<ask[i].time;T++){
change_color(change[T+1].position,change[T+1].num);
a[change[T+1].position]=change[T+1].num;
}
for(;T>ask[i].time;T--){
change_color(change[T].position,change[T].last);
a[change[T].position]=change[T].last;
}
if(i!=1)work(ask[i-1].l,ask[i].l),work(ask[i-1].r,ask[i].r);
else work(ask[i].l,ask[i].r);
reverse(lca(ask[i].l,ask[i].r));
int l=0,r=n,ans=0;
while(l<=r){
int mid=(l+r)>>1;
if(check(mid))ans=mid+1,l=mid+1;
else r=mid-1;
}
Ans[ask[i].id]=ans;
reverse(lca(ask[i].l,ask[i].r));
}
for(int i=1;i<=cnt1;i++)printf("%d\n",Ans[i]);
}
BZOJ 4129 树上带修莫队+线段树的更多相关文章
- BZOJ 3052 树上带修莫队
思路: 就是把带修莫队移到了树上 块的大小开到(n^2/3)/2 比较好- 这是一个卡OJ好题 //By SiriusRen #include <cmath> #include <c ...
- BZOJ 4129 Haruna’s Breakfast ( 树上带修莫队 )
题面 求树上某路径上最小的没出现过的权值,有单点修改 添加链接描述 分析 树上带修莫队板题,问题是怎么求最小的没出现过的权值. 因为只有nnn个点,所以没出现过的最小值一定在[0,n][0,n][0, ...
- BZOJ 3052/Luogu P4074 [wc2013]糖果公园 (树上带修莫队)
题面 中文题面,难得解释了 BZOJ传送门 Luogu传送门 分析 树上带修莫队板子题... 开始没给分块大小赋初值T了好一会... CODE #include <bits/stdc++.h&g ...
- bzoj4129 Haruna’s Breakfast 树上带修莫队+分块
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4129 题解 考虑没有修改的序列上的版本应该怎么做: 弱化的题目应该是这样的: 给定一个序列,每 ...
- 【BZOJ-3052】糖果公园 树上带修莫队算法
3052: [wc2013]糖果公园 Time Limit: 200 Sec Memory Limit: 512 MBSubmit: 883 Solved: 419[Submit][Status] ...
- LUOGU P4074 [WC2013]糖果公园 (树上带修莫队)
传送门 解题思路 树上带修莫队,搞了两天..终于开O2+卡常大法贴边过了...bzoj上跑了183s..其实就是把树上莫队和带修莫队结合到一起,首先求出括号序,就是进一次出一次那种的,然后如果求两个点 ...
- BZOJ3052: [wc2013]糖果公园【树上带修莫队】
Description Input Output Sample Input Sample Input Sample Output 84 131 27 84 HINT 思路 非常模板的树上带修莫队 真的 ...
- luogu4074 [WC2013]糖果公园(树上带修莫队)
link 题目大意:给一个树,树上每个点都有一种颜色,每个颜色都有一个收益 每次修改一个点上的颜色 或询问一条链上所有颜色第i次遇到颜色j可以获得w[i]*v[j]的价值,求链上价值和 题解:树上带修 ...
- 「洛谷1903」「BZOJ2120」「国家集训队」数颜色【带修莫队,树套树】
题目链接 [BZOJ传送门] [洛谷传送门] 题目大意 单点修改,区间查询有多少种数字. 解法1--树套树 可以直接暴力树套树,我比较懒,不想写. 稍微口胡一下,可以直接来一个树状数组套主席树,也就是 ...
随机推荐
- js-字符串方法
字符串 遍历字符串 方法:(类似数组) 使用for 或 for… in 结果:得到字符串中的每个字符 查找字符 ² charAt(索引值) 注: 超出索引值范围时,则返回空字符 ² ch ...
- Django[pronounced dʒ] installation on windows
1.Install python, download python windows installer from http://www.python.org/download/ and do inst ...
- 用python(2.7)自定义实现SQL的集合操作
有的时候需要在不同的数据库实例之间做集合操作,这就无法直接使用SQL语句的join,left join了.相同类型的数据库之间虽然也有类似于DBLINK和FEDERATED之类的东西,但一来这些东西不 ...
- R 连接DB2数据库,并制作词图
#写在前面的话:此教程主要是用R连接了DB2数据库,并进行文本分析,制作了词图 #教程为markdown编写 ---title: "网站留言分析"output: html_docu ...
- UVa 1585 待解决
是在遇到第一个ooxx的时候会出错,会少算一个1 #include<stdio.h> int main() { int i,k=0,sum=0; char a[100]={"oo ...
- RabbitMQ学习笔记(6)----RabbitMQ 持久化和非持久化
持久化:将交换机或队列数据保存到磁盘,服务器宕机或重启之后依然存在. 非持久化:将交换机或队列的数据保存到内存中,服务器宕机或重启之后数据将不存在. 在RabbitMQ中也提供了持久化和非持久化方式. ...
- CorelDRAW X6、X7&2018,500现金返利等你拿!
购物狂欢节要来了,你准备好了么? 不知何时起,四根神棍的日子却成了大家拼爹.拼钱包.拼手速.拼网速.拼钱包...各种火拼日子 你是从哪年关注并重视双11的,记得小编我第一次邂逅双11真的只是凑凑热闹 ...
- 打包c++项目
InstallShield Limited Edition for Visual Studio 2013 图文教程(教你如何打包.NET程序) InstallShield 2015 Limited E ...
- python笔记4----字典
1.哈希: 输入任意长度,输出固定长度. 即判断是否哈希,即判断可不可变. 2.创建字典 (1)直接创建:dic={1:'a',2:'b',3:'c'} (2)dict函数创建: list=[(1,' ...
- CSS - Span 下的width设置不可用?
解决:Span 下的width设置不可用? 内联元素-span有根据内容自动伸缩的能力,当需要对其宽度设定时,出现无效的情况. Demo:http://jsfiddle.net/JSDavi/ad62 ...