Denote by E(a,b,c) the given expression. First, we find the maximum value of E(a,b,c) under the assumption that a+b−c, b+c−a, c+a−b are non-negative.
Since a2+b2−ab≥ab, it follows that a3+b3=(a+b)(a2+b2−ab)≥ab(a+b), so a3+b3+abc≥ab(a+b+c). Analogously, we have b3+c3+abc≥bc(a+b+c) and c3+a3+abc≥ca(a+b+c). Therefore:
E(a,b,c)=∑a3+b3+abca+b−c≤∑ab(a+b+c)a+b−c=a+b+c1⋅∑aba+b−c
Since ∑aba+b−c=abc2(ab+bc+ca)−(a2+b2+c2)≤abc2(ab+bc+ca)−(ab+bc+ca) (due to the well-known inequality a2+b2+c2≥ab+bc+ca), we infer that
E(a,b,c)≤a+b+c1⋅abcab+bc+ca=a+b+c1(a1+b1+c1)≤1.
Consider now the case when at least one of the numbers a+b−c, b+c−a, c+a−b is negative. Without loss of generality, suppose that a+b−c<0. Then c>a+b>∣a−b∣, so both numbers b+c−a and c+a−b are positive. It follows that:
∑a3+b3+abca+b−c≤b3+c3+abcb+c−a+c3+a3+abcc+a−b≤bc(a+b+c)b+c−a+ca(a+b+c)c+a−b=a+b+c1⋅abcac+bc−(a−b)2<a+b+c1⋅abcab+bc+ca=a+b+c1(a1+b1+c1)≤1.
In summary, the maximum value of E(a,b,c) is equal to 1 and it is attained if and only if a=b=c=1.