先谈一下并查集:它是一种动态维护不重叠集合的做法,并且支持查询和合并的数据结构,是超级有用的;
1、find函数:查询某一个元素属于哪个集合;
2.union函数:将两个集合合并为一个大集合;
使用一个树形结构储存每一个集合,树上的节点表示一个元素,树根是代表集合,father数组表示当前i的父亲节点,初始化数组为本身,每次递归寻找根节点;
这样,合并的时候我们可以找到两棵树的树根,并且father[root1]=root2;就完成了操作;
并查集有压缩路径和按秩合并;
1.压缩路径:在每次执行find是将访问的节点直接指向树根;均摊复杂度O(logN);
2.按秩合并:秩表示的是树的深度(未压缩路径时),秩也可以表示集合的大小,每次合并是将秩小的元素合并到秩较大的元素上,只增加查询小的结构的代价,只使用按秩合并,复杂度:O(logN);
但同时使用可以将复杂度降到常数级别(不会写);
int find(int x) { if(x==father[x]) return x; else return father[x]=find(father[x]); } void Union(int x,int y) { int p=find(x),q=find(y); if(p!=q) father[p]=q; }
当然find函数可以使用符号运算符;
int find(int x){return x==father[x]?x:father[x]=find(father[x]);}
这样你就可以写洛谷的并查集模板了;
模板;
#include<bits/stdc++.h> using namespace std; #define N 500001 template<typename T>inline void read(T &x) { x=0;T f=1,ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();} x*=f; } int n,m,f[N]; inline int find(int x) { if(x!=f[x]) f[x]=find(f[x]); return f[x]; } inline void Union(int x,int y) { int fx=find(x); int fy=find(y); if(fx!=fy) f[fx]=fy; } int main() { read(n);read(m); for(int i=1;i<=n;i++) f[i]=i; for(int i=1;i<=m;i++) { int x,y,z; read(x);read(y);read(z); if(x==1) { int p=find(y),q=find(z); f[q]=p; } else { if(find(y)==find(z)) cout<<"Y"<<endl; else cout<<"N"<<endl; } } return 0; }
今天想谈三道题需要将并查集的思想拓展:
学长的题解:http://www.cnblogs.com/Cydiater/p/5810261.htm 二分+2-SAT 先预处理出所有的v,然后离散化一下,在那个的基础上二分,对于每次二分出的值约束边权超过所二分出的边权的两点。
还有二分图,嗯,不会嗯,那我们来看看并查集怎么处理这道题;
第一种思路:
我们可以建立两个集合,那么初始化也要注意是2*n了,一个集合表示这两个人不在一个监狱里,另一个集合表示这两个人在一个监狱里面,那么我们就可以将所有的怒气值从大到小进行一个排序,然后进行枚举,查找他们所在的集合,然后判断是不是出现了冲突(在同一个集合里面),如果出现了冲突,那么这个值就是怒气值的最大值,直接输出就好,如果没有出现冲突,那么将这两个放在一个并查集(表示不在同一个监狱的并查集)里面,表示这两个人是在不同的监狱里面。然后继续循环。
其实这道题就是运用到了并查集的容斥原理,和一点小小的贪心,我们将尽量怒气值比较大的放在不同的监狱里面,然后放次大的,如果出现冲突,也就是说在前面的安排之下,没有办法将当前的两个人分开,那么这两个人只能在一个监狱里面,这也就是最小的怒气值。
看代码:
#include<bits/stdc++.h> using namespace std; #define N 500001 template<typename T>inline void read(T &x) { x=0;T f=1,ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();} x*=f; } int n,m,d[N],f[N]; struct pink { int x,y,v; }a[N<<1]; bool mycmp(pink x,pink y) { return x.v>y.v; } int find(int x) { return f[x]==x?x:f[x]=find(f[x]); } int main() { int flag=0; read(n);read(m); for(int i=1;i<=m;i++) { read(a[i].x);read(a[i].y);read(a[i].v); } for(int i=1;i<=2*n;i++) f[i]=i; sort(a+1,a+m+1,mycmp); for(int i=1;i<=m;i++) { int x,y; x=a[i].x,y=a[i].y; int xx=find(a[i].x); int yy=find(a[i].y); if(xx==yy) { cout<<a[i].v<<endl; return 0; } f[yy]=find(a[i].x+n); f[xx]=find(a[i].y+n); } puts("0"); return 0; }
第二种思路:
可能发生冲突,所以尽量属于不同的监狱,就是属于不同的集合,怒气值越大,我们可以尽量让他们让他们属于不同的集合,如果在从大到小的安排之下不得已需要让两个人在一个集合,所以直接输出当前答案;敌人的敌人是可以放到同一个集合里的;
两种思路其实是一样的,我们用的是并查集维护一种关系;
那么:
#include<bits/stdc++.h> using namespace std; #define N 500001 template<typename T>inline void read(T &x) { x=0;T f=1,ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();} x*=f; } int n,m,d[N],f[N]; struct pink { int x,y,v; }a[N<<1]; bool mycmp(pink x,pink y) { return x.v>y.v; } int find(int x) { return f[x]==x?x:f[x]=find(f[x]); } int main() { int flag=0; read(n);read(m); for(int i=1;i<=m;i++) { read(a[i].x);read(a[i].y);read(a[i].v); } for(int i=1;i<=n;i++) f[i]=i; sort(a+1,a+m+1,mycmp); for(int i=1;i<=m;i++) { int x,y; x=a[i].x,y=a[i].y; if(find(x)==find(y)) { printf("%d ",a[i].v); flag=1; break; } else { if(!d[x]) d[x]=y; else { int p=find(d[x]); f[p]=find(y); } if(!d[y]) d[y]=x; else { int q=find(d[y]); f[q]=find(x); } } } if(!flag) printf("%d ",0); return 0; }
这是一个带权的并查集啦;
我们用一个d数组表示到根节点距离,size数组表示集合大小;
find和Union稍作修改;
inline int find(int x) { if(x==father[x]) return x; int root=find(father[x]);//递归寻找集合代表; d[x]+=d[father[x]];//求出边权和;return father[x]=root;//压缩路径; } void Union(int x,int y) { int p=find(x),q=find(y); father[p]=q; d[p]=size[q]; size[q]+=size[p]; }
当接受到c命令是查询,如果根节点相同,那么就是同一列,不同则不是;d数组分别记录x,y到跟节点,然后他们之间的abs(d[x]-d[y])-1;就是距离;
如果是M命令,合并即可,我们还需要用一个size数组记录集合大小;
#include<algorithm> #include<bitset> #include<cctype> #include<cerrno> #include<clocale> #include<cmath> #include<complex> #include<cstdio> #include<cstdlib> #include<cstring> #include<ctime> #include<deque> #include<exception> #include<fstream> #include<functional> #include<limits> #include<list> #include<map> #include<iomanip> #include<ios> #include<iosfwd> #include<iostream> #include<istream> #include<ostream> #include<queue> #include<set> #include<sstream> #include<stack> #include<stdexcept> #include<streambuf> #include<string> #include<utility> #include<vector> #include<cwchar> #include<cwctype> using namespace std; #define N 1000001 template<typename T>inline void read(T &x) { x=0;T f=1,ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();} x*=f; } int n,x,y,father[N],d[N],size[N]; inline int find(int x) { if(x==father[x]) return x; int root=find(father[x]); d[x]+=d[father[x]]; father[x]=root; return father[x]=root; } void Union(int x,int y) { int p=find(x),q=find(y); father[p]=q; d[p]=size[q]; size[q]+=size[p]; } int main() { for(int i=1;i<=N;i++) father[i]=i,d[i]=0,size[i]=1; read(n); for(int i=1;i<=n;i++) { char c; cin>>c; read(x);read(y); if(c=='M') Union(x,y); else if(c=='C') { int p=find(x),q=find(y); if(p!=q) { cout<<"-1"<<endl; } else { printf("%d ",abs(d[x]-d[y])-1); } } } return 0; }
这种传递关系只用一次并查集显然维护不了,所以我们可以考虑使用并查集的拓展域;
把每个动物x拆成三个点,同类域xself,捕食域xeat,天敌域xenemy;
如果x吃y,那么换句话讲,x的捕食域是y的同类,x的同类是y的天敌,x的天敌是y的捕食域,因为为环形关系;
如果x和y是同类,x的天敌和y的天敌是一类,x的捕食域和同类和y分别都一样;
合并就完事了;
判断谎言(flam)
与x和y是同类矛盾:
1.xeat和yself是一个集合;
2.xself和yeat是一个集合;
与x吃y矛盾:
1.xself和yself在一个集合;
2.xself和yeat在一个集合;
#include<bits/stdc++.h> using namespace std; #define N 1000001 template<typename T>inline void read(T &x) { x=0;T f=1,ch=getchar(); while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();} x*=f; } int n,m,xself,yself,xenemy,yenemy,xeat,yeat,father[N],x,y,z,flam; int find(int x) { if(x==father[x]) return x; else return father[x]=find(father[x]); } void Union(int x,int y) { int p=find(x),q=find(y); if(p!=q) father[p]=q; } int main() { read(n);read(m); for(int i=1;i<=3*n;i++) father[i]=i; for(int i=1;i<=m;i++) { read(z);read(x);read(y); xself=(x-1)*3+1,xeat=(x-1)*3+2,xenemy=x*3; yself=(y-1)*3+1,yeat=(y-1)*3+2,yenemy=y*3; if(x>n||y>n) { flam++; continue; } if(z==1) { if(find(xself)==find(yeat)||find(xeat)==find(yself)) { flam++; continue; } Union(xself,yself); Union(xenemy,yenemy); Union(xeat,yeat); } if(z==2) { if(find(xself)==find(yself)||find(yeat)==find(xself)) { flam++; continue; } Union(xeat,yself); Union(xself,yenemy); Union(xenemy,yeat); } } cout<<flam<<endl; return 0; }