An approximation algorithm for the maximum leaf spanning arborescence problem