87 lines
1.3 KiB
C
87 lines
1.3 KiB
C
|
/* $NetBSD: warshall.c,v 1.6 2011/09/10 21:29:04 christos Exp $ */
|
||
|
|
/* Id: warshall.c,v 1.7 2010/06/06 22:48:51 tom Exp */
|
||
|
|
|
||
|
|
#include "defs.h"
|
||
|
|
|
||
|
|
#include <sys/cdefs.h>
|
||
|
|
__RCSID("$NetBSD: warshall.c,v 1.6 2011/09/10 21:29:04 christos Exp $");
|
||
|
|
|
||
|
|
static void
|
||
|
|
transitive_closure(unsigned *R, int n)
|
||
|
|
{
|
||
|
|
int rowsize;
|
||
|
|
unsigned i;
|
||
|
|
unsigned *rowj;
|
||
|
|
unsigned *rp;
|
||
|
|
unsigned *rend;
|
||
|
|
unsigned *ccol;
|
||
|
|
unsigned *relend;
|
||
|
|
unsigned *cword;
|
||
|
|
unsigned *rowi;
|
||
|
|
|
||
|
|
rowsize = WORDSIZE(n);
|
||
|
|
relend = R + n * rowsize;
|
||
|
|
|
||
|
|
cword = R;
|
||
|
|
i = 0;
|
||
|
|
rowi = R;
|
||
|
|
while (rowi < relend)
|
||
|
|
{
|
||
|
|
ccol = cword;
|
||
|
|
rowj = R;
|
||
|
|
|
||
|
|
while (rowj < relend)
|
||
|
|
{
|
||
|
|
if (*ccol & (unsigned)(1 << i))
|
||
|
|
{
|
||
|
|
rp = rowi;
|
||
|
|
rend = rowj + rowsize;
|
||
|
|
while (rowj < rend)
|
||
|
|
*rowj++ |= *rp++;
|
||
|
|
}
|
||
|
|
else
|
||
|
|
{
|
||
|
|
rowj += rowsize;
|
||
|
|
}
|
||
|
|
|
||
|
|
ccol += rowsize;
|
||
|
|
}
|
||
|
|
|
||
|
|
if (++i >= BITS_PER_WORD)
|
||
|
|
{
|
||
|
|
i = 0;
|
||
|
|
cword++;
|
||
|
|
}
|
||
|
|
|
||
|
|
rowi += rowsize;
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
void
|
||
|
|
reflexive_transitive_closure(unsigned *R, int n)
|
||
|
|
{
|
||
|
|
int rowsize;
|
||
|
|
unsigned i;
|
||
|
|
unsigned *rp;
|
||
|
|
unsigned *relend;
|
||
|
|
|
||
|
|
transitive_closure(R, n);
|
||
|
|
|
||
|
|
rowsize = WORDSIZE(n);
|
||
|
|
relend = R + n * rowsize;
|
||
|
|
|
||
|
|
i = 0;
|
||
|
|
rp = R;
|
||
|
|
while (rp < relend)
|
||
|
|
{
|
||
|
|
*rp |= (unsigned)(1 << i);
|
||
|
|
if (++i >= BITS_PER_WORD)
|
||
|
|
{
|
||
|
|
i = 0;
|
||
|
|
rp++;
|
||
|
|
}
|
||
|
|
|
||
|
|
rp += rowsize;
|
||
|
|
}
|
||
|
|
}
|