nbody = int(input("nbody")) m = int((nbody/4)**(1./3.)) unitbox = 1./m r = [[0.0]*nbody for _ in range(3)] b = [[0.0]*4 for _ in range(3)] b[0][0] = 0.0 b[1][0] = 0.0 b[2][0] = 0.0 b[0][1] = unitbox/2.0 b[1][1] = unitbox/2.0 b[2][1] = 0.0 b[0][2] = unitbox/2.0 b[1][2] = 0.0 b[2][2] = unitbox/2.0 b[0][3] = 0.0 b[1][3] = unitbox/2.0 b[2][3] = unitbox/2.0 l = 1 for i in range(m): for j in range(m): for k in range(m): for n in range(0, 4): r[0][l-1] = b[0][n] + i*unitbox r[1][l-1] = b[1][n] + j*unitbox r[2][l-1] = b[2][n] + k*unitbox l += 1 for l in range(nbody): print(r[0][l], r[1][l], r[2][l])