added new function expand_polytomies to improve robinson_foulds calculation