Есть некая выпуклая функция. Одна переменная

целочисленная и положительная, остальные

переменных принадлежат множеству

Причем количество переменных очень велико (360). После декомпозиции удалось сократить их количество до 89. Перебор

комбинаций очень долгий. Это число даже не воспринимается компилятором. Ее требуется максимизировать по

при фиксированном

.

. По сути задачи

и

математические ожидания и дисперсии компонент нормального случайного вектора.
Помогите пожалуйста с методом решения такой задачи.