Oper: new function transitive_reduction
Rand: fixed type of gnp, with an extra argument of type unit