码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ】【3669】【NOI2014】魔法森林

时间:2015-02-24 23:20:31      阅读:216      评论:0      收藏:0      [点我收藏+]

标签:

LCT动态维护MST

  LCT动态维护MST

  我们可以枚举a,然后找从1到n的一条路径使得:这条路径上的b的最大值最小。这个路径肯定在MST上……所以枚举一遍所有的边,动态维护一个关于b值的MST即可。

 

调了半天没出解的原因:

  rotate写错了……l=c[y][1]==x 我写成了 l=c[z][1]==y sigh……

技术分享
  1 /**************************************************************
  2     Problem: 3669
  3     User: Tunix
  4     Language: C++
  5     Result: Accepted
  6     Time:4752 ms
  7     Memory:7896 kb
  8 ****************************************************************/
  9  
 10 //BZOJ 3669
 11 #include<vector>
 12 #include<cstdio>
 13 #include<cstring>
 14 #include<cstdlib>
 15 #include<iostream>
 16 #include<algorithm>
 17 #define rep(i,n) for(int i=0;i<n;++i)
 18 #define F(i,j,n) for(int i=j;i<=n;++i)
 19 #define D(i,j,n) for(int i=j;i>=n;--i)
 20 using namespace std;
 21 int getint(){
 22     int v=0,sign=1; char ch=getchar();
 23     while(ch<0||ch>9){ if (ch==-) sign=-1; ch=getchar();}
 24     while(ch>=0&&ch<=9){ v=v*10+ch-0; ch=getchar();}
 25     return v*=sign;
 26 }
 27 const int N=150010,INF=~0u>>2;
 28 /******************tamplate*********************/
 29 struct LCT{
 30     int c[N][2],fa[N],v[N],mx[N];
 31     bool rev[N];
 32     int st[N],top;
 33 #define L c[x][0]
 34 #define R c[x][1]
 35     void Push_up(int x){
 36         mx[x]=x;
 37         if (v[mx[L]]>v[mx[x]]) mx[x]=mx[L];
 38         if (v[mx[R]]>v[mx[x]]) mx[x]=mx[R];
 39     }
 40     void Push_down(int x){
 41         if (rev[x]) rev[x]=0,rev[L]^=1,rev[R]^=1,swap(L,R);
 42     }
 43     bool not_root(int x){
 44         return c[fa[x]][0]==x || c[fa[x]][1]==x;
 45     }
 46     void rotate(int x){
 47         int y=fa[x],z=fa[y],l=c[y][1]==x,r=l^1;
 48         if (not_root(y)) c[z][c[z][1]==y]=x;
 49         fa[x]=z; fa[y]=x; fa[c[x][r]]=y;
 50         c[y][l]=c[x][r]; c[x][r]=y;
 51         Push_up(y);
 52     }
 53     void preview(int x){
 54         top=0; st[++top]=x;
 55         for(;not_root(x);x=fa[x])
 56             st[++top]=fa[x];
 57         D(i,top,1) Push_down(st[i]);
 58     }
 59     void splay(int x,int y=0){
 60         for(preview(x);not_root(x);rotate(x))
 61             if (not_root(y=fa[x]))
 62                 rotate( c[y][1]==x^c[fa[y]][1]==y ? x : y);
 63         Push_up(x);
 64     }
 65     void access(int x,int y=0){
 66         for(;x;splay(x),c[x][1]=y,y=x,x=fa[x]);
 67     }
 68     void makeroot(int x){
 69         access(x); splay(x); rev[x]^=1;
 70     }
 71     void link(int x,int y){
 72         makeroot(x);fa[x]=y;
 73     }
 74     void cut(int x,int y){
 75         makeroot(x);access(y);splay(y);
 76         if (c[y][0]==x) c[y][0]=fa[x]=0;
 77     }
 78     int query(int x,int y){
 79         makeroot(x),access(y),splay(y);
 80         return mx[y];
 81     }
 82 }t;
 83 /*********************LCT***********************/
 84 struct edge{
 85     int x,y,a,b;
 86     bool operator < (const edge &e) const {
 87         return a < e.a;
 88     }
 89 }e[N];
 90 int fa[N];
 91 int find(int x){return fa[x]==x ? x : fa[x]=find(fa[x]);}
 92  
 93 int main(){
 94 #ifndef ONLINE_JUDGE
 95     freopen("3669.in","r",stdin);
 96     freopen("3669.out","w",stdout);
 97 #endif
 98     int n,m;
 99     n=getint(); m=getint();
100     F(i,1,m){
101         e[i].x=getint(); e[i].y=getint();
102         e[i].a=getint(); e[i].b=getint();
103     }
104     sort(e+1,e+m+1);
105     F(i,1,m){
106         t.v[n+i]=e[i].b;
107         t.mx[n+i]=n+i;
108     }
109     F(i,1,n) fa[i]=i;
110      
111     int ans=INF;
112     F(i,1,m){
113         int f1=find(e[i].x),f2=find(e[i].y);
114         if (f1!=f2){
115             fa[f1]=f2;
116             t.link(e[i].x,n+i); t.link(e[i].y,n+i);
117         }
118         else{
119             int tmp=t.query(e[i].x,e[i].y);
120             if (e[i].b<t.v[tmp]){//这一步即可略去自环
121                 t.cut(e[tmp-n].x,tmp); t.cut(e[tmp-n].y,tmp);
122                 t.link(e[i].x,i+n); t.link(e[i].y,i+n);
123             }
124         }
125         f1=find(1); f2=find(n);
126         if (f1==f2)
127             ans=min(ans,e[i].a+t.v[t.query(1,n)]);
128     }
129     printf("%d\n",ans==INF ? -1 : ans);
130     return 0;
131 }
View Code

 

【BZOJ】【3669】【NOI2014】魔法森林

标签:

原文地址:http://www.cnblogs.com/Tunix/p/4299031.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!