Path integral formulation of centroid dynamics for systems obeying Bose–Einstein statistics