[luogu3767]膜法
[luogu3767]膜法
luogu
神仙题
线段树分治+带权并查集
把每个操作看成点
首先这个操作的结构是一棵树
你发现每个点的对它的子树产生影响
我们可以想到用dfn序把它转成一段区间用线段树分治来做
但是还有删除操作,相当于在一个大区间里面挖掉几个小区间
可以对每个操作开一个vector记录区间搞一搞
然后带权并查集是模5意义下的,可以认为给你的操作相当于从u连向v的一条权值为1或2的边
当u,v在同一个集合时,判断是否满足条件,否则就连边
#define pb push_back
#define ls x<<1,l,mid
#define rs x<<1|1,mid+1,r
#include<bits/stdc++.h>
using namespace std;
const int _=1e5+5;
int re(){
int x=0,w=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
return x*w;
}
int n,m,ts,top,dis;
int fa[_],del[_],sz[_],dfn[_],par[_],siz[_],d[_];
bool ans[_];
struct node{int u,v;}st[_];
struct edge{int u,v,w;}e[_];
vector<int>son[_],s[_];
vector<edge>t[_<<2];
void dfs(int u){
sz[u]=1;if(u)dfn[u]=++ts;
if(del[u])s[del[u]].pb(u);
for(int i=0,j=son[u].size();i<j;i++){
int v=son[u][i];
dfs(v);sz[u]+=sz[v];
}
}
void add(int&x,int y){x=(x+y)%5;}
int find(int x){
add(dis,d[x]);
if(x==par[x])return x;
return find(par[x]);
}
void upd(int x,int l,int r,int ql,int qr,edge E){
if(ql<=l&&r<=qr){t[x].pb(E);return;}
int mid=(l+r)>>1;if(ql<=mid)upd(ls,ql,qr,E);
if(qr>mid)upd(rs,ql,qr,E);
}
void solve(int x,int l,int r,bool ok){
int pre=top;
for(int i=0,j=t[x].size();i<j;i++){
int u=t[x][i].u,v=t[x][i].v,w=t[x][i].w;
dis=0;int fu=find(u),du=dis;
dis=0;int fv=find(v),dv=dis;
if(fu==fv&&(du-dv+5)%5!=w)ok=0;
if(fu^fv){
if(siz[fu]>siz[fv]){swap(du,dv);swap(u,v);swap(fu,fv);w=-w;}
siz[fv]+=siz[fu];par[fu]=fv;
d[fu]=(w+dv-du+10)%5;st[++top]=(node){fu,fv};
}
}
if(l==r)ans[l]=ok;
else{int mid=(l+r)>>1;solve(ls,ok);solve(rs,ok);}
while(top^pre){
int u=st[top].u,v=st[top].v;top--;
siz[v]-=siz[u];par[u]=u;d[u]=0;
}
}
int main(){
n=re(),m=re();
int op,u,v;
for(int i=1;i<=m;i++){
son[fa[i]=re()].pb(i);op=re();
if(op==3)del[i]=re();
else{u=re(),v=re();e[i]=(edge){u,v,op};}
}
dfs(0);
for(int i=1;i<=m;i++){
if(del[i])continue;
int k=s[i].size(),lst=dfn[i];
for(int j=0;j<k;j++){
int u=s[i][j];
if(lst<dfn[u])upd(1,1,m,lst,dfn[u]-1,e[i]);
lst=dfn[u]+sz[u];
}
if(lst<dfn[i]+sz[i])upd(1,1,m,lst,dfn[i]+sz[i]-1,e[i]);
}
for(int i=1;i<=n;i++)par[i]=i,siz[i]=1;
solve(1,1,m,1);
for(int i=1;i<=m;i++)puts(ans[dfn[i]]?"excited":"naive");
return 0;
}
[luogu3767]膜法的更多相关文章
- 【OpenJudge3531】【背包DP】【膜法交配律】判断整除
判断整除 总时间限制: 1000ms 内存限制: 65536kB [描述] 一个给定的正整数序列,在每个数之前都插入+号或-号后计算它们的和.比如序列:1.2.4共有8种可能的序列:(+1) + (+ ...
- luogu P3767 膜法
传送门 这题如果没有删除操作,可以直接使用可持久化并查集 注意到这种可持久化的依赖关系(是这样说的把)是一棵树,然后对于一个点,自己的操作会影响自己的那棵子树,并且如果是删除操作,就会使得一个子树没有 ...
- B - 低阶入门膜法 - D-query (查询区间内有多少不同的数)
题目链接:https://cn.vjudge.net/contest/284294#problem/B 题目大意:查询区间内有多少个不相同的数. 具体思路:主席树的做法,主席树的基础做法是查询区间第k ...
- A - 低阶入门膜法 - K-th Number (主席树查询区间第k小)
题目链接:https://cn.vjudge.net/contest/284294#problem/A 题目大意:主席树查询区间第k小. 具体思路:主席树入门. AC代码: #include<i ...
- BZOJ3669 膜法森林 - LCT
Solution 非常妙的排序啊... 仔细想想好像确实能够找出最优解QUQ 先对第一关键字排序, 在$LCT$ 维护第二关键字的最大值 所在的边. 添边时如果$u, v$ 不连通 就直接加边. 如 ...
- 2018.10.16 NOIP模拟 膜法(组合数学)
传送门 原题,原题,全TM原题. 不得不说天天考原题. 其实这题我上个月做过类似的啊,加上dzyodzyodzyo之前有讲过考试直接切了. 要求的其实就是∑i=lr(ii−l+k)\sum _{i=l ...
- NOIP模拟题 膜法
题目大意 给定若干组询问求$\sum\limits_{i=l}^r \dbinom{i}{k}$. 最终输出每组询问答案的乘积. 题解 首先把$l,r$分开处理相减,只需要求$\sum\limits_ ...
- 【bzoj5055】膜法师(离散化+树状数组)
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=5055 这道题……不得不说,从标题到题面都能看出一股浓浓的膜法气息……苟…… 题意就是统计顺序 ...
- 【PCB】扫盲总结
1.PCB是什么 PCB( Printed Circuit Board),中文名称为印制电路板,又称印刷线路板,是重要的电子部件,是电子元器件的支撑体,是电子元器件电气连接的载体.由于它是采用电子印刷 ...
随机推荐
- 针对ecshop错误404页面的优化
在ecshop系统当中,比如你随意将商品详细页面的地址中的ID修改为一个不存在的商品ID,ecshop会自动跳转到首页.ecshop在这方面做得非常的差,甚至导致了很多的站不被搜索引擎收录.最模板提供 ...
- SSO单点登录系列4:cas-server登录页面自定义修改过程(jsp页面修改)
落雨 cas 单点登录 SSO单点登录系列4:cas-server登录页面自定义修改过程,全新DIY. 目标: 下面是正文: 打开cas的默认首页,映入眼帘的是满眼的中文and英文混杂体,作为一 ...
- oc 之中的 汉字字符串转化成为拼音 汉字字符串的排序
在oc 之中的字符串为汉字的时候,我们经常要进行字符串比較,可是汉字不能比較,所以就要将汉字转化成为拼音,详细步骤例如以下: //可变字符串 必须是可变字符串. NSMutableString ...
- 网络编程readn、writen和readline函数的编写
readn 在Linux中,read的声明为: ssize_t read(int fd, void *buf, size_t count); 它的返回值有以下情形: 1.大于0,代表成功读取的字节 ...
- Pro Tools安装图文教程
Pro Tools安装图文教程 Avid Pro Tools是Digidesign公司出品的一款音质最佳.音频制作强大的软件,能够在Mac或PC上为影片编曲.录制.编辑和混制高品质音乐或声音,生成 ...
- DeleteDC、ReleaseDC 、DeleteObject的使用
DeleteDC 该函数删除指定的设备上下文环境(DC). 原型: BOOL DeleteDC(HDC hdc): 参数: hdc:设备上下文环境的句柄. 返回值: 成功,返回非零值:失败,返回零.调 ...
- js中的string.format函数代码
String.prototype.format = function(args) { if (arguments.length > 0) { var result = this; if (arg ...
- int a[3];中a+1与&a+1差别 -- C
int a[3]; a 和 &a 的地址一样的. a+1 == a + 1*sizeof(int);跳跃是一个数组元素大小 &a+1 == a + 3*sizeof(int);跳跃是整 ...
- vs无法引用项目问题
vs无法引用项目问题 2017年12月13日 14:45:31 阅读数:582 开发时编译报错--项目A未被引用,展开项目的引用,发现该项目实质已经被引用了,但是该引用上有个黄色三角感叹号,遂移除该引 ...
- PHP面试题及答案解析(1)—PHP语法基础
1. strlen( )与 mb_strlen( )的作用分别是什么? strlen和mb_strlen都是用于获取字符串长度.strlen只针对单字节编码字符,也就是说它计算的是字符串的总字节数.如 ...