MCM2016 This is our outstanding winner award paper for Problem C in 2016 COMAP's Mathematical Contest in Modeling