-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCloseness centrality1.cpp
153 lines (113 loc) · 3.42 KB
/
Closeness centrality1.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
// Closeness centrality1.cpp : This file contains the 'main' function. Program execution begins and ends there.
//
#include "pch.h"
#include <iostream>
#include <vector>
#include <queue>
#include <list>
#include <iomanip>
#include <fstream>
using namespace std;
vector<int> get_shortest_distance(int src, vector< list< pair<int, int> > > graph);
double closeness_centrality_node(int index, vector< list< pair<int, int> > > graph);
vector<double>closeness_centrality_graph(vector< list< pair<int, int> > > graph);
int main()
{
int n, m;
ifstream myfile;
myfile.open("input.txt");
if (myfile.is_open())
{
}
else
{
cout<<"error";
}
myfile >> n >> m;
vector< list< pair<int, int> > > graph(n);
vector<double>out1;
for (int i = 0; i < m; i++)
{
int x, y, z;
myfile >> x >> y >> z;
graph[x].push_back(make_pair(y, z));
graph[y].push_back(make_pair(x, z));
}
myfile.close();
out1 = closeness_centrality_graph(graph);
for (int j = 0; j < out1.size(); j++)
{
cout << setprecision(11) << out1[j] << '\n';
}
}
/*
a function that implement dijestra algorithm
parameters:
1-the index of the source node that we want to calculate its distance from all other nodes in the graph
2- the given graph
*/
class Comparepairs
{
public:
bool operator()(pair<int, int> p1, pair<int, int> p2) {
return p1.second > p2.second;
}
};
vector<int> get_shortest_distance(int src, vector< list< pair<int, int> > > graph)
{
list< pair<int, int> >::iterator i;
priority_queue < pair<int, int>, vector< pair<int, int> >, Comparepairs> heap;
vector<int> distance_from_nodes(graph.size(), MAX_INT);
heap.push({ src,0 });
while (!heap.empty())
{
int current_node = heap.top().first;
int current_distance = heap.top().second;
heap.pop();
if (distance_from_nodes[current_node] != MAX_INT)
continue;
distance_from_nodes[current_node] = current_distance;
for (i = graph[current_node].begin(); (i != graph[current_node].end()); i++)
{
pair<int, int>p = *i;
p.second += current_distance;
if (distance_from_nodes[p.first] != MAX_INT)
continue;
heap.push(p);
}
}
return distance_from_nodes;
}
double closeness_centrality_node(int index, vector< list< pair<int, int> > > graph)
{
vector <int> out;
out = get_shortest_distance(index, graph);
int sum = 0, n = graph.size();
double cc = 0;
int size = out.size();
for (int i = 0; i < size; i++)
{
sum += out[i];
}
cc = ((double)(n - 1)) / ((double)sum);
return cc;
}
vector<double>closeness_centrality_graph(vector< list< pair<int, int> > > graph)
{
int n = graph.size();
vector<double> out(n);
for (int i = 0; i < n; i++)
{
out[i] = closeness_centrality_node(i, graph);
}
return out;
}
// Run program: Ctrl + F5 or Debug > Start Without Debugging menu
// Debug program: F5 or Debug > Start Debugging menu
// Tips for Getting Started:
// 1. Use the Solution Explorer window to add/manage files
// 2. Use the Team Explorer window to connect to source control
// 3. Use the Output window to see build output and other messages
// 4. Use the Error List window to view errors
// 5. Go to Project > Add New Item to create new code files, or Project > Add Existing Item to add existing code files to the project
// 6. In the future, to open this project again, go to File > Open > Project and select