k = 15 x = round(((2+3**.5)**k-(2-3**.5)**k)/(2*3**.5)) y = round((1/6)*(3*(2-3**.5)**k-2*3**.5*(2-3**.5)**k+3*(2+3**.5)**k+2*3**.5*(2+3**.5)**k)) print(0,0) print(x,y) print(y,x)