#include "bitset.h"
#include "bipartite.h"
+#include "xmalloc.h"
struct _bipartite_t {
int n_left, n_right;
- bitset_t **adj;
+ bitset_t *adj[1];
};
bipartite_t *bipartite_new(int n_left, int n_right)
{
- int i;
- bipartite_t *gr = malloc(sizeof(*gr));
+ int i, size;
+ bipartite_t *gr;
+
+ size = n_left > 0 ? n_left - 1 : 0;
+ gr = xmalloc(sizeof(*gr) + size * sizeof(void *));
memset(gr, 0, sizeof(*gr));
gr->n_left = n_left;
gr->n_right = n_right;
- gr->adj = malloc(n_left * sizeof(void *));
for(i = 0; i < n_left; ++i)
gr->adj[i] = bitset_malloc(n_right);
int i;
for(i = 0; i < gr->n_left; ++i)
bitset_free(gr->adj[i]);
-
- free(gr->adj);
free(gr);
}
{
int left, right;
int done_something = 0;
- bitset_t *tmp = bitset_malloc(gr->n_left);
+ bitset_t *tmp = bitset_alloca(gr->n_right);
for(left = 0; left < gr->n_left; ++left) {
bitset_t *left_adj = gr->adj[left];
assert(right != -1);
- /* We have to find another left node which has the old right
- * one as a neighbor. */
+ /*
+ We have to find another left node which has the old right one as a neighbor.
+ This node must not be part of a matching
+ */
for(i = 0; i < gr->n_left; ++i)
- if(i != left && bitset_is_set(gr->adj[i], old_right))
+ if(i != left && bitset_is_set(gr->adj[i], old_right) && !bitset_is_set(matched_left, i))
break;
/* If no such node can be found, exit. */
}
}
- bitset_free(tmp);
-
return done_something;
}
void bipartite_matching(const bipartite_t *gr, int *matching)
{
- bitset_t *matched_left = bitset_malloc(gr->n_left);
- bitset_t *matched_right = bitset_malloc(gr->n_right);
+ bitset_t *matched_left = bitset_alloca(gr->n_left);
+ bitset_t *matched_right = bitset_alloca(gr->n_right);
memset(matching, -1, gr->n_left * sizeof(int));
while(apply_alternating_path(gr, matching, matched_left, matched_right));
+}
- bitset_free(matched_left);
- bitset_free(matched_right);
+void bipartite_dump_f(FILE *f, const bipartite_t *gr)
+{
+ int i;
+
+ for(i = 0; i < gr->n_left; ++i) {
+ bitset_fprint(f, gr->adj[i]);
+ fprintf(f, "\n");
+ }
+}
+
+void bipartite_dump(const char *name, const bipartite_t *gr) {
+ FILE *f = fopen(name, "w");
+
+ if (f) {
+ bipartite_dump_f(f, gr);
+ fclose(f);
+ }
}