function [X Y] = init_bipartite(N)

% assumes N is even
m = N/2;

X(1:m) = 2;
X(m+1:N) = 5;
Y(1:m) = 1:m;
Y(m+1:N) = 1:m;

X = X'
Y = Y';