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
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
|
--!A cross-platform build utility based on Lua
--
-- Licensed under the Apache License, Version 2.0 (the "License");
-- you may not use this file except in compliance with the License.
-- You may obtain a copy of the License at
--
-- http://www.apache.org/licenses/LICENSE-2.0
--
-- Unless required by applicable law or agreed to in writing, software
-- distributed under the License is distributed on an "AS IS" BASIS,
-- WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-- See the License for the specific language governing permissions and
-- limitations under the License.
--
-- Copyright (C) 2015-present, TBOOX Open Source Group.
--
-- @author ruki
-- @file graph.lua
--
-- load modules
local table = require("base/table")
local list = require("base/list")
local object = require("base/object")
-- define module
local graph = graph or object { _init = {"_directed"} } {true}
local edge = edge or object { _init = {"_from", "_to"} }
-- new edge, from -> to
function edge.new(from, to)
return edge {from, to}
end
function edge:from()
return self._from
end
function edge:to()
return self._to
end
function edge:other(v)
if v == self._from then
return self._to
else
return self._from
end
end
function edge:__tostring()
return string.format("<edge:%s-%s>", self:from(), self:to())
end
-- clear graph
function graph:clear()
self._vertices = {}
self._edges = {}
self._adjacent_edges = {}
self._edges_map = {}
end
-- is empty?
function graph:empty()
return #self:vertices() == 0
end
-- is directed?
function graph:is_directed()
return self._directed
end
-- get vertices
function graph:vertices()
return self._vertices
end
-- get adjacent edges of the the given vertex
function graph:adjacent_edges(v)
return self._adjacent_edges[v]
end
-- get the vertex at the given index
function graph:vertex(idx)
return self:vertices()[idx]
end
-- has the given vertex?
function graph:has_vertex(v)
return table.contains(self:vertices(), v)
end
-- remove the given vertex?
function graph:remove_vertex(v)
local contains = false
table.remove_if(self._vertices, function (_, item)
if item == v then
contains = true
return true
end
end)
if contains then
self._edges_map[v] = nil
self._adjacent_edges[v] = nil
-- remove the adjacent edge with this vertex in the other vertices
if not self:is_directed() then
for _, w in ipairs(self:vertices()) do
local edges = self:adjacent_edges(w)
if edges then
table.remove_if(edges, function (_, e)
if e:other(w) == v then
self._edges_map[w] = nil
return true
end
end)
end
end
end
end
end
-- topological sort, use Kahn's algorithm
--
-- e.g.
--
-- add_edge(a, b) -- a depend on b
-- add_edge(b, c) -- b depend on c
--
-- it will return {c, b, a}
function graph:topological_sort(opt)
opt = opt or {}
if not self:is_directed() then
return
end
-- calculate in-degree for each vertex
local in_degree = {}
for _, v in ipairs(self:vertices()) do
in_degree[v] = 0
end
-- count incoming edges for each vertex
for _, v in ipairs(self:vertices()) do
local edges = self:adjacent_edges(v)
if edges then
for _, e in ipairs(edges) do
if e:from() == v then
local w = e:to()
in_degree[w] = (in_degree[w] or 0) + 1
end
end
end
end
-- queue of vertices with no incoming edges (no dependencies)
local queue = list.new()
for _, v in ipairs(self:vertices()) do
if in_degree[v] == 0 then
queue:insert(v)
end
end
-- result list for topologically sorted vertices
local order_vertices = {}
-- process queue
while not queue:empty() do
-- remove a vertex with no incoming edges
local v = queue:remove_first()
table.insert(order_vertices, v)
-- for each outgoing edge, remove it and update in-degrees
local edges = self:adjacent_edges(v)
if edges then
for _, e in ipairs(edges) do
if e:from() == v then
local w = e:to()
in_degree[w] = in_degree[w] - 1
-- if in-degree becomes zero, add to queue
if in_degree[w] == 0 then
queue:insert(w)
end
end
end
end
end
-- if we couldn't process all vertices, there must be a cycle
local has_cycle = #order_vertices ~= #self:vertices()
return order_vertices, has_cycle
end
-- find cycle
function graph:find_cycle()
local visited = {}
local stack = {}
local cycle = {}
local function dfs(v)
visited[v] = true
stack[v] = true
table.insert(cycle, v)
local edges = self:adjacent_edges(v)
if edges then
for _, e in ipairs(edges) do
local w = e:other(v)
if not visited[w] then
if dfs(w) then
return true
elseif stack[w] then
return true
end
elseif stack[w] then
for i = #cycle, 1, -1 do
if cycle[i] == w then
cycle = table.slice(cycle, i)
return true
end
end
end
end
end
table.remove(cycle)
stack[v] = false
return false
end
for _, v in ipairs(self:vertices()) do
if not visited[v] then
if dfs(v) then
return cycle
end
end
end
end
-- get edges
function graph:edges()
return self._edges
end
-- add edge
function graph:add_edge(from, to)
local e = edge.new(from, to)
if not self:has_vertex(from) then
table.insert(self._vertices, from)
self._adjacent_edges[from] = {}
end
if not self:has_vertex(to) then
table.insert(self._vertices, to)
self._adjacent_edges[to] = {}
end
local edges_map = self._edges_map
edges_map[from] = edges_map[from] or {}
edges_map[from][to] = true
if self:is_directed() then
table.insert(self._adjacent_edges[from], e)
else
table.insert(self._adjacent_edges[from], e)
table.insert(self._adjacent_edges[to], e)
edges_map[to] = edges_map[to] or {}
edges_map[to][from] = true
end
table.insert(self._edges, e)
end
-- has the given edge?
function graph:has_edge(from, to)
local edges = self:adjacent_edges(from)
if edges then
local edges_map = self._edges_map
local from_map = edges_map[from]
if from_map and from_map[to] then
return true
else
for _, e in ipairs(edges) do
if e:to() == to then
return true
end
end
end
end
return false
end
-- clone graph
function graph:clone()
local gh = graph.new(self:is_directed())
for _, v in ipairs(self:vertices()) do
local edges = self:adjacent_edges(v)
if edges then
for _, e in ipairs(edges) do
gh:add_edge(e:from(), e:to())
end
end
end
return gh
end
-- reverse graph
function graph:reverse()
if not self:is_directed() then
return self:clone()
end
local gh = graph.new(self:is_directed())
for _, v in ipairs(self:vertices()) do
local edges = self:adjacent_edges(v)
if edges then
for _, e in ipairs(edges) do
gh:add_edge(e:to(), e:from())
end
end
end
return gh
end
-- dump graph
function graph:dump()
local vertices = self:vertices()
local edges = self:edges()
print(string.format("graph: %s, vertices: %d, edges: %d", self:is_directed() and "directed" or "not-directed", #vertices, #edges))
print("vertices: ")
for _, v in ipairs(vertices) do
print(string.format(" %s", v))
end
print("")
print("edges: ")
for _, e in ipairs(edges) do
print(string.format(" %s -> %s", e:from(), e:to()))
end
end
-- new graph
function graph.new(directed)
local gh = graph {directed}
gh:clear()
return gh
end
-- return module: graph
return graph
|