Here is a function to maximize:
$$V(x_1, \dots , x_n) := \sum_{i=1}^n x_i a_i - b \sum_{i=1}^n d_i | x_i - c_i |$$
where $x$ is the independent variable (in an $n$-dimensional space), $a$, $d$ and $c$ are constant n-dimensional arrays and $b$ is a constant real number.
The optimization problem is subject to some linear constraints on $x_1, ..., x_n$. The objective function contains a first term, which is linear in $x$, and a second term, involving absolute values. Maximizing the first term is a standard linear programming problem, and maximizing the second term alone can be reduced to a linear programming problem introducing dummy variables (see this answer for example).
Is there a way to reduce this "mixed" optimization problem to a linear programming problem?