Optimizing Gurobi model (PuLP) for a symmetric MIP with binary variables

Question

I\'m trying to get into integer linear programming. For my first project, I want to replicate the results from the solution to the "8 soldiers lining up for the morning