C++實現(xiàn)基于不相交集合的O(mlgn)復雜度的kruskal算法
C++實現(xiàn)基于不相交集合的O(mlgn)復雜度的kruskal算法
本文實現(xiàn)完全參考<<Introduction to Algorithms Third edition>>,
不相交集合的數(shù)據(jù)結(jié)構(gòu)
我們采用森林的方式實現(xiàn)不相交集合。這個森林是極簡化的,每個節(jié)點只有一個指向父親的指針,而且森林中的每一顆樹都是一個集合,我們?nèi)涞母?jié)點為這個集合的代表元。
int rank[505]; int father[505]; void make_set(int x) { father[x]=x; rank[x]=0; } int find_set(int x) { if (x!=father[x]) { father[x]=find_set(father[x]); } return father[x]; } void simply_union_set(int u,int v) { u=find_set(u); v=find_set(v); father[u]=v; } void perfect_union_set(int u,int v) { u=find_set(u); v=find_set(v); if (rank[u]>rank[v]) { father[v]=u; } else { father[u]=v; if(rank[u]==rank[v]) rank[v]++; } }
可以看到在find_set()函數(shù)中采用了兩趟遍歷的思想,第一趟遍歷找的根節(jié)點,第二趟遍歷將路徑上的節(jié)點全部指向根節(jié)點,完成了壓縮樹高。
在實現(xiàn)集合合并的時候,我們采用了兩種方法:一種方法是直接合并simply_union_set,另一種是采用按秩合并的思想perfect_union_set,即總是讓秩小合并到秩大的集合中,這是一種減少樹高的有效策略;
當我們采用按秩合并時時,上述每一個操作的最差時間復雜度,都約等于O(1)
詳情見<<Introduction to Algorithms Third edition>>中證明
kruskal 算法
void kruskal() { for(int i=0;i<num_v;i++)make_set(i); sort(arr_edge.begin(),arr_edge.end(),mycompare); for(int i=0;i<arr_edge.size();i++) { int fr=arr_edge[i].fr; int to=arr_edge[i].to; int w=arr_edge[i].w; if( find_set(fr)!=find_set(to)) { result+=w; perfect_union_set(fr,to); } } }
kruskal 算法是一種基于貪心策略的算法,它的時間復雜度的最大開銷就是排序算法,即O(mlgm)=O(mlgn),這里m表示邊數(shù),n表示頂點數(shù)
知識補充
乘勝追擊一下,通過一個例題再深入了解一下kruskal 算法吧
題目:http://poj.org/problem?id=2485
思路:就是最小生成樹啊
代碼
#include<iostream> #include<algorithm> #include<cstring> #include<cstdio> #include<vector> using namespace std; #define INTMAX 0x3f3f3f3f typedef pair<int,int> pii; typedef long long ll; #define x first #define y second int rank[505]; int father[505]; int find_set(int x) { if (x!=father[x]) { father[x]=find_set(father[x]); } return father[x]; } void simply_union_set(int u,int v) { u=find_set(u); v=find_set(v); father[u]=v; } void perfect_union_set(int u,int v) { u=find_set(u); v=find_set(v); if (rank[u]>rank[v]) { father[v]=u; } else { father[u]=v; if(rank[u]==rank[v]) rank[v]++; } } struct edge { int fr,to,w; }; int num_case,num_v,result; vector<edge> arr_edge; void debug() { for(int i=0;i<arr_edge.size();i++) { cout<<arr_edge[i].fr<<" to "<<arr_edge[i].to<<"="<<arr_edge[i].w<<endl; } } void init() { arr_edge.clear(); result=0; } void input() { int w; scanf("%d",&num_v); for(int i=0;i<num_v;i++) { for(int j=0;j<num_v;j++) { scanf("%d",&w); if(i<j) { edge temp; temp.fr=i; temp.to=j; temp.w=w; arr_edge.push_back(temp); } } } } bool mycompare(const edge& x,const edge &y) { return x.w<y.w; } void kruskal() { for(int i=0;i<num_v;i++)father[i]=i; sort(arr_edge.begin(),arr_edge.end(),mycompare); for(int i=0;i<arr_edge.size();i++) { int fr=arr_edge[i].fr; int to=arr_edge[i].to; int w=arr_edge[i].w; if( find_set(fr)!=find_set(to)) { result=max(result,w); simply_union_set(fr,to); } } } void solve() { init(); input(); //debug(); kruskal(); cout<<result<<endl; } int main() { scanf("%d",&num_case); while(num_case--) { solve(); } return 0; }
到此這篇關(guān)于C++實現(xiàn)基于不相交集合的O(mlgn)復雜度的kruskal算法的文章就介紹到這了,更多相關(guān)C++ kruskal算法內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
VC++實現(xiàn)文件與應用程序關(guān)聯(lián)的方法(注冊表修改)
這篇文章主要介紹了VC++實現(xiàn)文件與應用程序關(guān)聯(lián)的方法,涉及VC++針對注冊表的相關(guān)操作技巧,需要的朋友可以參考下2016-08-08