for Ri, Rj in ((i, j) for i in range(NR) for j in range(i + 1, NR)):