On computing the maximum parsimony score of a phylogenetic network