动态树专题 WC 2006 Tube 还有范浩强的“动态树好题”
来源:互联网 发布:mac os常用软件 编辑:程序博客网 时间:2024/04/26 15:23
虽然是O(Qn)的算法可过水题,但是既然是练手的话,就不得不维护一下动态MST了。
还是倒着做,删边改为加边,利用环切性质直接找出切点,再将其切开、反向、连边就成了
但是反向是这题最蛋疼的一个操作,标记维护限于水平木有在splay中维护(反正我调不出来),直接在外维护了。做了一天半啊,中途怒删6KB的代码重打,终于。。。。。内牛满面。
裸代码题,无需多说50 1000 2000
#include<cstdio>#include<cstdlib>#include<cmath>#include<cstring>typedef int arr[100005];int n=0,m=0,task=0,e=0,top=0;arr w,px,py,next,head,poin,foll,con,maxw,othm,ww,dsf,x,y,nw,k,f,index,posw,poso,sta,ans;bool vis[1005],revese[1005];int poi[1005][1005];int son[1005][2];inline void swap(int &u,int &v){ int tmp=u; u=v; v=tmp;}inline void init(){ scanf("%d%d%d",&n,&m,&task); int i=0; memset(w,65,sizeof(w)); for (i=1;i<=m;i++) { scanf("%d%d%d",&x[i],&y[i],&w[i]); poi[x[i]][y[i]]=i; poi[y[i]][x[i]]=i; } for (i=1;i<=task;i++) { scanf("%d%d%d",&k[i],&px[i],&py[i]); if (k[i]==2){ nw[i]=w[poi[px[i]][py[i]]]; w[poi[px[i]][py[i]]]=w[0]+1;} }}inline int find(int x){ if (dsf[x]==x) return x; dsf[x]=find(dsf[x]); return dsf[x];}inline int cmp(const void *p,const void *q){ return w[*(int*)p]-w[*(int*)q];}inline void kruskal(){ int i=0; for (i=1;i<=n;i++) dsf[i]=i; for (i=1;i<=m;i++) index[i]=i; qsort(index+1,m,sizeof(index[0]),cmp); int sum=0,u=0,v=0; for (i=1;i<=m;i++) if (find(x[index[i]])!=find(y[index[i]])) {u=x[index[i]];v=y[index[i]];dsf[dsf[u]]=dsf[v];next[++e]=head[u];head[u]=e;poin[e]=v;ww[e]=w[index[i]];next[++e]=head[v];head[v]=e;poin[e]=u;ww[e]=w[index[i]];++sum;if (sum>=n-1) break; }}inline void dfs(int s){ vis[s]=1; int ne=head[s],y=0; while (ne) { if (!vis[poin[ne]]){ y=poin[ne]; dfs(y); f[y]=s; con[y]=ww[ne];} ne=next[ne]; }}inline void mark(int s){ revese[s]^=1; swap(con[s],foll[s]); swap(maxw[s],othm[s]); swap(son[s][0],son[s][1]); swap(posw[s],poso[s]);}inline void collectmark(int s){ sta[top=1]=s; while (f[s] &&(son[f[s]][0]==s || son[f[s]][1]==s)) { sta[++top]=f[s]; s=f[s]; } while (top) { if (revese[sta[top]]) revese[sta[top]]=0,mark(son[sta[top]][0]),mark(son[sta[top]][1]); top--; }}inline void update(int s){ maxw[s]=con[s],posw[s]=s; if (maxw[s]<maxw[son[s][0]]) maxw[s]=maxw[son[s][0]],posw[s]=posw[son[s][0]]; if (maxw[s]<maxw[son[s][1]]) maxw[s]=maxw[son[s][1]],posw[s]=posw[son[s][1]]; othm[s]=foll[s],poso[s]=s; if (othm[s]<othm[son[s][0]]) othm[s]=othm[son[s][0]],poso[s]=poso[son[s][0]]; if (othm[s]<othm[son[s][1]]) othm[s]=othm[son[s][1]],poso[s]=poso[son[s][1]];}inline void rotate(int s,int t){ int fs=f[s],ffs=f[fs],ch=son[s][t^1]; if (son[ffs][0]==fs) son[ffs][0]=s; else if (son[ffs][1]==fs) son[ffs][1]=s; f[s]=ffs; son[s][t^1]=fs; f[fs]=s; son[fs][t]=ch; if (ch) f[ch]=fs; update(fs);}inline void splay(int s){ collectmark(s); int fs=0,ffs=0,p=0,q=0; while (f[s] &&(son[f[s]][0]==s || son[f[s]][1]==s)) { fs=f[s],ffs=f[fs]; if (ffs &&(son[ffs][0]==fs || son[ffs][1]==fs)){ if (son[ffs][0]==fs) p=0; else p=1; if (son[fs][0]==s) q=0; else q=1; if (p==q) rotate(fs,p),rotate(s,q); else rotate(s,q),rotate(s,p);} elseif (son[fs][0]==s) rotate(s,0); else rotate(s,1); } update(s);}inline int flast(int s){ if (revese[s]) revese[s]=0,mark(son[s][0]),mark(son[s][1]); s=son[s][1]; if (revese[s]) revese[s]=0,mark(son[s][0]),mark(son[s][1]); while (son[s][0]) { s=son[s][0]; if (revese[s]) revese[s]=0,mark(son[s][0]),mark(son[s][1]); } return s;}inline int froot(int s){ if (revese[s]) revese[s]=0,mark(son[s][0]),mark(son[s][1]); while (son[s][0]) { s=son[s][0]; if (revese[s]) revese[s]=0,mark(son[s][0]),mark(son[s][1]); } return s;}inline int access(int s){ int last=0; for (;s;s=f[s]) { splay(s); son[s][1]=last; foll[s]=con[froot(last)]; last=s; update(s); } return last;}inline void cut(int s){ access(s); splay(s); con[s]=0; foll[flast(s)]=0; f[son[s][0]]=0; update(son[s][0]); son[s][0]=0; f[s]=0; update(s);}inline void evert(int s){ access(s); splay(s); mark(s);}int main(){ freopen("tube.in","r",stdin); freopen("tube.out","w",stdout); init(); kruskal(); dfs(1); int i=0,u=0,v=0,now=0,lca=0; for (i=task;i>=1;i--) { u=px[i],v=py[i]; if (k[i]==1){ access(u); lca=access(v); ans[i]=maxw[son[lca][1]]; if (u!=lca) { splay(u); ans[i]=ans[i]>maxw[u]?ans[i]:maxw[u]; }} else{ access(u); lca=access(v); ans[i]=maxw[son[lca][1]],now=posw[son[lca][1]]; if (u!=lca) { splay(u); if (ans[i]<maxw[u]){ ans[i]=maxw[u]; now=posw[u];} } if (ans[i]<=nw[i]) continue; cut(now); evert(v); f[v]=u; con[v]=nw[i];} } for (i=1;i<=task;i++) if (k[i]==1) printf("%d\n",ans[i]); return 0;}
至于范浩强的那题,其实也是比较裸的,首先一颗仙人掌只有一个环,并查集预处理找环,之后change时分类讨论是否从环上切下,是否连接时又新成了一个环,同时维护root的“父亲”,也就是范浩强所说的基点,求答案时就直接求基点的解,access(expose)时用迭代的方法维护一个sk,sb。link cut tree其实并不难实现,倒是还未学树块剖分
代码就不放了,去fhq博客就可以了
一道动态树好题 - fanhq666的日志 - 网易博客
- 动态树专题 WC 2006 Tube 还有范浩强的“动态树好题”
- you tube的架构
- [WC 2015复习](五)动态规划
- HDU动态规划专题
- 动态规划专题:beginner
- 动态规划专题:BZOJ1207
- 动态规划专题
- 动态规划专题(III)
- POJ 动态规划专题
- 动态规划专题总结
- 动态规划专题总结
- 动态规划专题总结
- 动态规划专题
- LeetCode专题----动态规划
- 【专题】基础动态规划
- 动态规划专题
- 动态规划专题总结
- 动态规划专题总结!
- 安防监控进入云时代
- JS操作GRIDVIEW多行编辑及新增行与保存,外加ACCESS数据库多行插入(事务)
- linux下vi命令大全
- 上海扬讯
- Oracle数据导入导出imp/exp命令 10g以上expdp/impdp命令
- 动态树专题 WC 2006 Tube 还有范浩强的“动态树好题”
- 界面原型设计工具使用系列(五)-CogTool
- Android开发_如何设置按钮背景透明与半透明_图片背景透明
- 在Ubuntu 11.10 编译android的改动
- Servlet
- centos6最小化安装后
- 浏览器系列:css+div如何解决文字溢出
- ptrace
- VC++中播放声音的方法