maxp_greedy function

A greedy algorithm to solve the max-p-region problem