Combining Lagrangian relaxation and a two-set column generation model for integrated railway freight planning