De Kruskal algoritmo e disjoint-set estrutura de dados:Eu preciso fazer as seguintes duas linhas de código?
-
12-11-2019 - |
Pergunta
Eu tenho implementado de Kruskal algoritmo em C++ usando o disjoint-set estrutura de dados de acordo com a Wikipedia, como este:
#include <stdio.h>
#include <algorithm>
#define MAX_EDGES 10000000
#define MAX_VERTICES 200001
using namespace std;
int num_edges,num_vertices;
int total_cost=0;
struct edge{
int v1,v2;
int cost;
};
struct comp{
bool operator()(const edge& e1,const edge& e2){
return e1.cost<e2.cost;
}
};
edge edges[MAX_EDGES];
int parent[MAX_VERTICES];
int rank[MAX_VERTICES];
int findset(int x){
if(x!=parent[x]){
parent[x]=findset(parent[x]);
}
return parent[x];
}
void merge(int x,int y){
int px=findset(x),py=findset(y);
if(rank[px]>rank[py]){
parent[py]=px;
}else{
parent[px]=py;
}
if(rank[px]==rank[py]){
++rank[py];
}
}
int main(){
FILE* in=fopen("input","r");
FILE* out=fopen("output","w");
fscanf(in,"%d %d\n",&num_vertices,&num_edges);
for(int i=1;i<=num_vertices;++i){
parent[i]=i;
rank[i]=0;
}
for(int i=0;i<num_edges;++i){
fscanf(in,"%d %d %d\n",&edges[i].v1,&edges[i].v2,&edges[i].cost);
}
sort(edges,edges+num_edges,comp());
for(int i=0;i<num_edges;++i){
int s1=findset(edges[i].v1),s2=findset(edges[i].v2);
if(s1!=s2){
merge(s1,s2);
total_cost+=edges[i].cost;
}
}
fprintf(out,"%d\n",total_cost);
}
A minha pergunta é:Eu preciso que estas duas linhas de código?Se assim for, qual é a sua importância?
int px=findset(x),py=findset(y);
na mala, em vez deint px=parent[x],py=parent[y];
parent[x]=findset(parent[x]);
no findset em vez dereturn findset(parent[x]);
Solução
1) findset(x)
retorna a canonical representativa do conjunto que x é (a raiz da sua ascendência de árvore).Você precisa deste para ser capaz de comparar se dois elementos estão no mesmo conjunto ou não (eles têm o mesmo representante), parent[x]
retorna apenas o pai de x na árvore, que pode não ser a raiz.
1a) Você esqueceu de teste para px e py sendo idênticos em merge
.
2) É uma otimização para que as futuras chamadas para findset
vai correr mais rápido.Se parent[x]
usado para apontar para seu pai, que apontavam para a raiz de sua árvore, após esta chamada parent[x]
será que aponte diretamente para a raiz.
Outras dicas
- Você precisar isto porque
x.parent
não é necessariamente o representante da classe para a qualx
pertence, portanto, o algoritmo não seria correto sem ele. - Sem a atribuição, o algoritmo teria de ser de qualidade inferior.Este é o caminho de compressão a otimização, também discutido na Wikipédia.