Home | History | Annotate | Line # | Download | only in tsort
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