tsort.c revision 1.23 1 1.23 joerg /* $NetBSD: tsort.c,v 1.23 2011/09/06 18:34:37 joerg Exp $ */
2 1.9 jtc
3 1.1 cgd /*
4 1.9 jtc * Copyright (c) 1989, 1993, 1994
5 1.7 cgd * The Regents of the University of California. All rights reserved.
6 1.1 cgd *
7 1.1 cgd * This code is derived from software contributed to Berkeley by
8 1.1 cgd * Michael Rendell of Memorial University of Newfoundland.
9 1.1 cgd *
10 1.1 cgd * Redistribution and use in source and binary forms, with or without
11 1.1 cgd * modification, are permitted provided that the following conditions
12 1.1 cgd * are met:
13 1.1 cgd * 1. Redistributions of source code must retain the above copyright
14 1.1 cgd * notice, this list of conditions and the following disclaimer.
15 1.1 cgd * 2. Redistributions in binary form must reproduce the above copyright
16 1.1 cgd * notice, this list of conditions and the following disclaimer in the
17 1.1 cgd * documentation and/or other materials provided with the distribution.
18 1.18 agc * 3. Neither the name of the University nor the names of its contributors
19 1.1 cgd * may be used to endorse or promote products derived from this software
20 1.1 cgd * without specific prior written permission.
21 1.1 cgd *
22 1.1 cgd * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
23 1.1 cgd * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
24 1.1 cgd * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
25 1.1 cgd * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
26 1.1 cgd * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
27 1.1 cgd * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
28 1.1 cgd * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
29 1.1 cgd * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
30 1.1 cgd * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
31 1.1 cgd * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
32 1.1 cgd * SUCH DAMAGE.
33 1.1 cgd */
34 1.1 cgd
35 1.20 lukem #if HAVE_NBTOOL_CONFIG_H
36 1.20 lukem #include "nbtool_config.h"
37 1.20 lukem #endif
38 1.20 lukem
39 1.12 lukem #include <sys/cdefs.h>
40 1.20 lukem #if !defined(lint)
41 1.22 lukem __COPYRIGHT("@(#) Copyright (c) 1989, 1993, 1994\
42 1.22 lukem The Regents of the University of California. All rights reserved.");
43 1.9 jtc #if 0
44 1.10 jtc static char sccsid[] = "@(#)tsort.c 8.3 (Berkeley) 5/4/95";
45 1.9 jtc #endif
46 1.23 joerg __RCSID("$NetBSD: tsort.c,v 1.23 2011/09/06 18:34:37 joerg Exp $");
47 1.1 cgd #endif /* not lint */
48 1.17 tv
49 1.1 cgd #include <sys/types.h>
50 1.9 jtc #include <ctype.h>
51 1.9 jtc #include <db.h>
52 1.17 tv #include <err.h>
53 1.1 cgd #include <errno.h>
54 1.7 cgd #include <fcntl.h>
55 1.1 cgd #include <stdio.h>
56 1.7 cgd #include <stdlib.h>
57 1.1 cgd #include <string.h>
58 1.10 jtc #include <unistd.h>
59 1.1 cgd
60 1.1 cgd /*
61 1.7 cgd * Topological sort. Input is a list of pairs of strings separated by
62 1.1 cgd * white space (spaces, tabs, and/or newlines); strings are written to
63 1.1 cgd * standard output in sorted order, one per line.
64 1.1 cgd *
65 1.1 cgd * usage:
66 1.7 cgd * tsort [-l] [inputfile]
67 1.1 cgd * If no input file is specified, standard input is read.
68 1.1 cgd *
69 1.15 wiz * Should be compatible with AT&T tsort HOWEVER the output is not identical
70 1.1 cgd * (i.e. for most graphs there is more than one sorted order, and this tsort
71 1.1 cgd * usually generates a different one then the AT&T tsort). Also, cycle
72 1.1 cgd * reporting seems to be more accurate in this version (the AT&T tsort
73 1.1 cgd * sometimes says a node is in a cycle when it isn't).
74 1.1 cgd *
75 1.1 cgd * Michael Rendell, michael (at) stretch.cs.mun.ca - Feb 26, '90
76 1.1 cgd */
77 1.1 cgd #define HASHSIZE 53 /* doesn't need to be big */
78 1.1 cgd #define NF_MARK 0x1 /* marker for cycle detection */
79 1.1 cgd #define NF_ACYCLIC 0x2 /* this node is cycle free */
80 1.7 cgd #define NF_NODEST 0x4 /* Unreachable */
81 1.7 cgd
82 1.1 cgd typedef struct node_str NODE;
83 1.1 cgd
84 1.1 cgd struct node_str {
85 1.1 cgd NODE **n_prevp; /* pointer to previous node's n_next */
86 1.1 cgd NODE *n_next; /* next node in graph */
87 1.7 cgd NODE **n_arcs; /* array of arcs to other nodes */
88 1.1 cgd int n_narcs; /* number of arcs in n_arcs[] */
89 1.1 cgd int n_arcsize; /* size of n_arcs[] array */
90 1.1 cgd int n_refcnt; /* # of arcs pointing to this node */
91 1.1 cgd int n_flags; /* NF_* */
92 1.7 cgd char n_name[1]; /* name of this node */
93 1.1 cgd };
94 1.1 cgd
95 1.1 cgd typedef struct _buf {
96 1.1 cgd char *b_buf;
97 1.1 cgd int b_bsize;
98 1.1 cgd } BUF;
99 1.1 cgd
100 1.23 joerg static DB *db;
101 1.23 joerg static NODE *graph, **cycle_buf, **longest_cycle;
102 1.23 joerg static int debug, longest, quiet;
103 1.23 joerg
104 1.23 joerg static void add_arc(char *, char *);
105 1.23 joerg static void clear_cycle(void);
106 1.23 joerg static int find_cycle(NODE *, NODE *, int, int);
107 1.23 joerg static NODE *get_node(char *);
108 1.23 joerg static void *grow_buf(void *, int);
109 1.23 joerg static void remove_node(NODE *);
110 1.23 joerg static void tsort(void);
111 1.23 joerg __dead static void usage(void);
112 1.1 cgd
113 1.7 cgd int
114 1.23 joerg main(int argc, char *argv[])
115 1.1 cgd {
116 1.12 lukem BUF *b;
117 1.12 lukem int c, n;
118 1.1 cgd FILE *fp;
119 1.7 cgd int bsize, ch, nused;
120 1.1 cgd BUF bufs[2];
121 1.14 cgd
122 1.14 cgd setprogname(argv[0]);
123 1.1 cgd
124 1.12 lukem fp = NULL;
125 1.12 lukem while ((ch = getopt(argc, argv, "dlq")) != -1)
126 1.9 jtc switch (ch) {
127 1.7 cgd case 'd':
128 1.9 jtc debug = 1;
129 1.9 jtc break;
130 1.7 cgd case 'l':
131 1.9 jtc longest = 1;
132 1.9 jtc break;
133 1.11 mycroft case 'q':
134 1.11 mycroft quiet = 1;
135 1.11 mycroft break;
136 1.7 cgd case '?':
137 1.7 cgd default:
138 1.7 cgd usage();
139 1.7 cgd }
140 1.7 cgd argc -= optind;
141 1.7 cgd argv += optind;
142 1.7 cgd
143 1.9 jtc switch (argc) {
144 1.7 cgd case 0:
145 1.1 cgd fp = stdin;
146 1.7 cgd break;
147 1.7 cgd case 1:
148 1.7 cgd if ((fp = fopen(*argv, "r")) == NULL)
149 1.9 jtc err(1, "%s", *argv);
150 1.7 cgd break;
151 1.7 cgd default:
152 1.7 cgd usage();
153 1.1 cgd }
154 1.1 cgd
155 1.1 cgd for (b = bufs, n = 2; --n >= 0; b++)
156 1.7 cgd b->b_buf = grow_buf(NULL, b->b_bsize = 1024);
157 1.1 cgd
158 1.1 cgd /* parse input and build the graph */
159 1.1 cgd for (n = 0, c = getc(fp);;) {
160 1.1 cgd while (c != EOF && isspace(c))
161 1.1 cgd c = getc(fp);
162 1.1 cgd if (c == EOF)
163 1.1 cgd break;
164 1.1 cgd
165 1.1 cgd nused = 0;
166 1.1 cgd b = &bufs[n];
167 1.1 cgd bsize = b->b_bsize;
168 1.1 cgd do {
169 1.1 cgd b->b_buf[nused++] = c;
170 1.7 cgd if (nused == bsize)
171 1.7 cgd b->b_buf = grow_buf(b->b_buf, bsize *= 2);
172 1.1 cgd c = getc(fp);
173 1.1 cgd } while (c != EOF && !isspace(c));
174 1.1 cgd
175 1.1 cgd b->b_buf[nused] = '\0';
176 1.1 cgd b->b_bsize = bsize;
177 1.1 cgd if (n)
178 1.1 cgd add_arc(bufs[0].b_buf, bufs[1].b_buf);
179 1.1 cgd n = !n;
180 1.1 cgd }
181 1.1 cgd (void)fclose(fp);
182 1.7 cgd if (n)
183 1.8 cgd errx(1, "odd data count");
184 1.1 cgd
185 1.1 cgd /* do the sort */
186 1.1 cgd tsort();
187 1.21 christos return(0);
188 1.1 cgd }
189 1.1 cgd
190 1.1 cgd /* double the size of oldbuf and return a pointer to the new buffer. */
191 1.23 joerg static void *
192 1.23 joerg grow_buf(void *bp, int size)
193 1.1 cgd {
194 1.19 itojun void *n;
195 1.19 itojun
196 1.19 itojun if ((n = realloc(bp, (u_int)size)) == NULL)
197 1.12 lukem err(1, "realloc");
198 1.19 itojun bp = n;
199 1.7 cgd return (bp);
200 1.1 cgd }
201 1.1 cgd
202 1.1 cgd /*
203 1.1 cgd * add an arc from node s1 to node s2 in the graph. If s1 or s2 are not in
204 1.1 cgd * the graph, then add them.
205 1.1 cgd */
206 1.23 joerg static void
207 1.23 joerg add_arc(char *s1, char *s2)
208 1.1 cgd {
209 1.12 lukem NODE *n1;
210 1.1 cgd NODE *n2;
211 1.4 cgd int bsize, i;
212 1.1 cgd
213 1.7 cgd n1 = get_node(s1);
214 1.1 cgd
215 1.1 cgd if (!strcmp(s1, s2))
216 1.1 cgd return;
217 1.1 cgd
218 1.7 cgd n2 = get_node(s2);
219 1.1 cgd
220 1.1 cgd /*
221 1.3 cgd * Check if this arc is already here.
222 1.3 cgd */
223 1.3 cgd for (i = 0; i < n1->n_narcs; i++)
224 1.3 cgd if (n1->n_arcs[i] == n2)
225 1.3 cgd return;
226 1.3 cgd /*
227 1.3 cgd * Add it.
228 1.1 cgd */
229 1.1 cgd if (n1->n_narcs == n1->n_arcsize) {
230 1.1 cgd if (!n1->n_arcsize)
231 1.1 cgd n1->n_arcsize = 10;
232 1.1 cgd bsize = n1->n_arcsize * sizeof(*n1->n_arcs) * 2;
233 1.7 cgd n1->n_arcs = grow_buf(n1->n_arcs, bsize);
234 1.1 cgd n1->n_arcsize = bsize / sizeof(*n1->n_arcs);
235 1.1 cgd }
236 1.1 cgd n1->n_arcs[n1->n_narcs++] = n2;
237 1.1 cgd ++n2->n_refcnt;
238 1.1 cgd }
239 1.1 cgd
240 1.7 cgd /* Find a node in the graph (insert if not found) and return a pointer to it. */
241 1.23 joerg static NODE *
242 1.23 joerg get_node(char *name)
243 1.1 cgd {
244 1.7 cgd DBT data, key;
245 1.7 cgd NODE *n;
246 1.1 cgd
247 1.7 cgd if (db == NULL &&
248 1.7 cgd (db = dbopen(NULL, O_RDWR, 0, DB_HASH, NULL)) == NULL)
249 1.9 jtc err(1, "db: %s", name);
250 1.7 cgd
251 1.7 cgd key.data = name;
252 1.7 cgd key.size = strlen(name) + 1;
253 1.7 cgd
254 1.9 jtc switch ((*db->get)(db, &key, &data, 0)) {
255 1.7 cgd case 0:
256 1.21 christos (void)memmove(&n, data.data, sizeof(n));
257 1.7 cgd return (n);
258 1.7 cgd case 1:
259 1.7 cgd break;
260 1.7 cgd default:
261 1.7 cgd case -1:
262 1.9 jtc err(1, "db: %s", name);
263 1.7 cgd }
264 1.1 cgd
265 1.7 cgd if ((n = malloc(sizeof(NODE) + key.size)) == NULL)
266 1.12 lukem err(1, "malloc");
267 1.1 cgd
268 1.1 cgd n->n_narcs = 0;
269 1.1 cgd n->n_arcsize = 0;
270 1.7 cgd n->n_arcs = NULL;
271 1.1 cgd n->n_refcnt = 0;
272 1.1 cgd n->n_flags = 0;
273 1.21 christos (void)memmove(n->n_name, name, key.size);
274 1.1 cgd
275 1.7 cgd /* Add to linked list. */
276 1.9 jtc if ((n->n_next = graph) != NULL)
277 1.1 cgd graph->n_prevp = &n->n_next;
278 1.1 cgd n->n_prevp = &graph;
279 1.1 cgd graph = n;
280 1.1 cgd
281 1.7 cgd /* Add to hash table. */
282 1.7 cgd data.data = &n;
283 1.7 cgd data.size = sizeof(n);
284 1.7 cgd if ((*db->put)(db, &key, &data, 0))
285 1.9 jtc err(1, "db: %s", name);
286 1.7 cgd return (n);
287 1.7 cgd }
288 1.7 cgd
289 1.7 cgd
290 1.7 cgd /*
291 1.7 cgd * Clear the NODEST flag from all nodes.
292 1.7 cgd */
293 1.23 joerg static void
294 1.23 joerg clear_cycle(void)
295 1.7 cgd {
296 1.7 cgd NODE *n;
297 1.7 cgd
298 1.9 jtc for (n = graph; n != NULL; n = n->n_next)
299 1.7 cgd n->n_flags &= ~NF_NODEST;
300 1.1 cgd }
301 1.1 cgd
302 1.1 cgd /* do topological sort on graph */
303 1.23 joerg static void
304 1.23 joerg tsort(void)
305 1.1 cgd {
306 1.12 lukem NODE *n, *next;
307 1.12 lukem int cnt, i;
308 1.1 cgd
309 1.9 jtc while (graph != NULL) {
310 1.1 cgd /*
311 1.7 cgd * Keep getting rid of simple cases until there are none left,
312 1.1 cgd * if there are any nodes still in the graph, then there is
313 1.1 cgd * a cycle in it.
314 1.1 cgd */
315 1.1 cgd do {
316 1.9 jtc for (cnt = 0, n = graph; n != NULL; n = next) {
317 1.1 cgd next = n->n_next;
318 1.1 cgd if (n->n_refcnt == 0) {
319 1.1 cgd remove_node(n);
320 1.1 cgd ++cnt;
321 1.1 cgd }
322 1.1 cgd }
323 1.9 jtc } while (graph != NULL && cnt);
324 1.1 cgd
325 1.9 jtc if (graph == NULL)
326 1.1 cgd break;
327 1.1 cgd
328 1.1 cgd if (!cycle_buf) {
329 1.1 cgd /*
330 1.7 cgd * Allocate space for two cycle logs - one to be used
331 1.1 cgd * as scratch space, the other to save the longest
332 1.1 cgd * cycle.
333 1.1 cgd */
334 1.9 jtc for (cnt = 0, n = graph; n != NULL; n = n->n_next)
335 1.1 cgd ++cnt;
336 1.7 cgd cycle_buf = malloc((u_int)sizeof(NODE *) * cnt);
337 1.7 cgd longest_cycle = malloc((u_int)sizeof(NODE *) * cnt);
338 1.7 cgd if (cycle_buf == NULL || longest_cycle == NULL)
339 1.12 lukem err(1, "malloc");
340 1.1 cgd }
341 1.13 ross for (n = graph; n != NULL; n = n->n_next) {
342 1.13 ross if (!(n->n_flags & NF_ACYCLIC)) {
343 1.12 lukem if ((cnt = find_cycle(n, n, 0, 0)) != 0) {
344 1.11 mycroft if (!quiet) {
345 1.11 mycroft warnx("cycle in data");
346 1.11 mycroft for (i = 0; i < cnt; i++)
347 1.11 mycroft warnx("%s",
348 1.11 mycroft longest_cycle[i]->n_name);
349 1.11 mycroft }
350 1.1 cgd remove_node(n);
351 1.7 cgd clear_cycle();
352 1.1 cgd break;
353 1.7 cgd } else {
354 1.1 cgd /* to avoid further checks */
355 1.7 cgd n->n_flags |= NF_ACYCLIC;
356 1.7 cgd clear_cycle();
357 1.7 cgd }
358 1.13 ross }
359 1.13 ross }
360 1.9 jtc if (n == NULL)
361 1.9 jtc errx(1, "internal error -- could not find cycle");
362 1.1 cgd }
363 1.1 cgd }
364 1.1 cgd
365 1.1 cgd /* print node and remove from graph (does not actually free node) */
366 1.23 joerg static void
367 1.23 joerg remove_node(NODE *n)
368 1.1 cgd {
369 1.12 lukem NODE **np;
370 1.12 lukem int i;
371 1.1 cgd
372 1.1 cgd (void)printf("%s\n", n->n_name);
373 1.1 cgd for (np = n->n_arcs, i = n->n_narcs; --i >= 0; np++)
374 1.1 cgd --(*np)->n_refcnt;
375 1.1 cgd n->n_narcs = 0;
376 1.1 cgd *n->n_prevp = n->n_next;
377 1.1 cgd if (n->n_next)
378 1.1 cgd n->n_next->n_prevp = n->n_prevp;
379 1.1 cgd }
380 1.1 cgd
381 1.7 cgd
382 1.7 cgd /* look for the longest? cycle from node from to node to. */
383 1.23 joerg static int
384 1.23 joerg find_cycle(NODE *from, NODE *to, int longest_len, int depth)
385 1.1 cgd {
386 1.12 lukem NODE **np;
387 1.12 lukem int i, len;
388 1.1 cgd
389 1.1 cgd /*
390 1.1 cgd * avoid infinite loops and ignore portions of the graph known
391 1.1 cgd * to be acyclic
392 1.1 cgd */
393 1.7 cgd if (from->n_flags & (NF_NODEST|NF_MARK|NF_ACYCLIC))
394 1.7 cgd return (0);
395 1.7 cgd from->n_flags |= NF_MARK;
396 1.1 cgd
397 1.1 cgd for (np = from->n_arcs, i = from->n_narcs; --i >= 0; np++) {
398 1.1 cgd cycle_buf[depth] = *np;
399 1.1 cgd if (*np == to) {
400 1.1 cgd if (depth + 1 > longest_len) {
401 1.1 cgd longest_len = depth + 1;
402 1.21 christos (void)memcpy(longest_cycle, cycle_buf,
403 1.1 cgd longest_len * sizeof(NODE *));
404 1.1 cgd }
405 1.1 cgd } else {
406 1.7 cgd if ((*np)->n_flags & (NF_MARK|NF_ACYCLIC|NF_NODEST))
407 1.7 cgd continue;
408 1.1 cgd len = find_cycle(*np, to, longest_len, depth + 1);
409 1.7 cgd
410 1.7 cgd if (debug)
411 1.9 jtc (void)printf("%*s %s->%s %d\n", depth, "",
412 1.9 jtc from->n_name, to->n_name, len);
413 1.7 cgd
414 1.7 cgd if (len == 0)
415 1.9 jtc (*np)->n_flags |= NF_NODEST;
416 1.7 cgd
417 1.7 cgd if (len > longest_len)
418 1.1 cgd longest_len = len;
419 1.7 cgd
420 1.7 cgd if (len > 0 && !longest)
421 1.5 cgd break;
422 1.1 cgd }
423 1.1 cgd }
424 1.7 cgd from->n_flags &= ~NF_MARK;
425 1.7 cgd return (longest_len);
426 1.1 cgd }
427 1.1 cgd
428 1.23 joerg static void
429 1.23 joerg usage(void)
430 1.7 cgd {
431 1.11 mycroft (void)fprintf(stderr, "usage: tsort [-lq] [file]\n");
432 1.7 cgd exit(1);
433 1.1 cgd }
434