Function to maximize
$R(P) = \sum\limits_{i=1}^4 (-\frac{1}{2}p_i^2 + 1000p_i)$
Parameters
$P = (p_1, p_2, p_3, p_4)$
Constraints
$p_1, p_2, p_3, p_4 > 0$
$\sum\limits_{i=1}^4 (-\frac{1}{2}p_i + 1000) \leq 2000$
$R(P) = \sum\limits_{i=1}^4 (-\frac{1}{2}p_i^2 + 1000p_i)$
$P = (p_1, p_2, p_3, p_4)$
$p_1, p_2, p_3, p_4 > 0$
$\sum\limits_{i=1}^4 (-\frac{1}{2}p_i + 1000) \leq 2000$
This is the first time I have answered a question, so please bear with me.
Maximizing $R(P)=\sum_1^4(−1/2p^2_i+1000p_i)$ is the same as minimizing its negative. Doubling it and adding 1000000 to each summand doesn't change the minimum and we get $R(P)=\sum_1^4(p^2_i-2000p_i+1000000)$ $ =\sum_1^4(p_i-1000)^2$. Setting each $p_1=1000$ satisfies the constraints and makes $R(P)=0$ and, since it is a sum of squares this must be the minimum.
https://en.wikipedia.org/wiki/Geometric_programming I quite like it but it is basically the same. In any case, a big clue is the symmetry under permutations of the variables. – rrogers Dec 29 '20 at 21:30