According to wikipedia (but this is well known), A klee-Minty cube in D-dimension is defined by
x1≤5
4x1+x2≤25
8x1+4x2+x3≤125
$16 x_1 + 8 x_2 + 4 x_3 + x_4 \leq 625\$
32x1+16x2+8x3+4x4+x5≤3125
⋮⋮
2Dx1+2D−1x2+2D−2x3+…+4xD−1+xD≤5D
x1≥0
x2≥0
x3≥0
x4≥0
x5≥0
This poloyhedron has 2D vertices. How can I find them. I have tried the following code
A= matrix(QQ, 10,5,[0,0,0,0,-1, 0,0,0,-1,-4, 0,0,-1,-4,-8, 0,-1,-4,-8,-16, -1,-4,-8,-16,-32,
1,0,0,0,0, 0,1,0,0,0, 0,0,1,0,0, 0,0,0,1,0, 0,0,0,0,1])
show(LatexExpr(r"\text{A = }"),A)
b=vector(QQ,[5,254,125,625,3125,0,0,0,0,0])
show(LatexExpr(r"\text{b = }"),b)
c=vector(QQ,[2^4,2^3,2^2,2^1,2^0])
show(LatexExpr(r"\text{c = }"),c)
AA=matrix(list(list([b])+list(transpose(A))))
show(transpose(AA))
pol = Polyhedron(ieqs = AA)
pol.Hrepresentation()
But I am far away from the account and some of the vertices are in the negative part of the hyperplan (which is impossible). Need some help to decipher my mistake. Thanks