pscore function

Compute the parsimony score