  The $k$-NN hard core lattice gas model on a square lattice, in which the
first $k$ next nearest neighbor sites of a particle are excluded from being
occupied by another particle, is the lattice version of the hard disc model in
two dimensional continuum. It has been conjectured that the lattice model, like
its continuum counterpart, will show multiple entropy-driven transitions with
increasing density if the high density phase has columnar or striped order.
Here, we determine the nature of the phase at full packing for $k$ up to
$820302$. We show that there are only eighteen values of $k$, all less than
$k=4134$, that show columnar order, while the others show solid-like sublattice
order.
