Для оптимизации сначала binary_search можно сначала сделать, и его результат использовать как верхнюю границу для lower_bound и как нижнюю — для upper_bound