code: def _adjust_probability_vec_best(population, fitnesses, probability_vec,
                                   adjust_rate):
      
      best_solution = max(zip(fitnesses, population))[1]
  
      # Shift probabilities towards best solution
      return _adjust(probability_vec, best_solution, adjust_rate)
docs: Shift probabilities towards the best solution.