Efficiently simulating the work distribution of multiple identical bosons with boson sampling