ER.py

def ErdosRenyi(N, m=N):
    net = np.zeros(2*m, dtype='int')

    for i in xrange(m):
        while True:
            node1, node2 = np.random.randint(0, N, 2)

            # No self loops allowed
            if node1 != node2:
                break

        net[2*i] = node1
        net[2*i+1] = node2

    return net