The evolution of gene regulatory networks in variable environments poses Multi-objective Optimization Problem (MOP), where the expression levels of genes must be tuned to meet the demands of each en- vironment. When formalized in the context of monotone systems, this problem falls into a sub-class of linear MOPs. Here, the constraints are partial orders and the objectives consist of either the minimization or maximization of single variables, but their number can be very large. To efficiently and exhaustively find Pareto optimal solutions, we introduce a mapping between colored Hasse diagrams and polytopes associated with an ideal point. A dynamic program based on edge contractions yields an exact closed- form description of the Pareto optimal set, in polynomial time of the number of objectives relative to the number of faces of the Pareto front. We additionally discuss the special case of series-parallel graphs with monochromatic connected components of bounded size, for which the running time and the representation of solutions can in principle be linear in the number of objectives.