Let $*$ be an associative binary law on the interval $(0,\infty)$ of the positive reals. Suppose that for every $a,b,c \in (0, \infty)$ we know $$ a * b * c = \frac{abc}{ab+bc+ca}.$$ Can we show that for all $a,b \in (0, \infty)$ we have $a * b = \frac{ab}{a+b}$?
I suspect that the answer to this question is yes, but was not able to make any meaningful progress towards a full proof in the last few days. Any ideas would be appreciated.
The following is an idea that won't fit into a comment. I'm not sure whether it is helpful, or whether it is merely the kind of thing that people are encouraged to put into the question body to indicate "what you tried."
What if we use the binary operation $*$ to define a new binary operation $\oplus$, and this time our set has an identity element (say "e")?
If for $a > 0$ and $b > 0$ we have $a \oplus b = a * b$ and if $$(a \oplus b \oplus e) = (a \oplus b)$$ then we could use a formula for $$(a \oplus b \oplus c)$$ (that is similar to our formula for $(a * b * c)$ ) and we would have some hope of obtaining a formula for $a \oplus b$, which should itself simply be equal to $(a*b)$, and then we will be done.
Now, one difficulty is that the original expression provided for $(a * b * c)$ has the product $a b$ in the denominator, and we need a generalization of that expression so that, when c is our identity element, the term $a b$ will disappear.
Consider the following:
$$(a \oplus b \oplus c) = \frac{abc}{((a b)(c - 1)) + b c + c a}$$
If the above equation could be relied upon, then we would be able to obtain in particular $$(a \oplus b \oplus 1) = \frac{ab}{b + a}$$ and in that scenario it would seem that 1 serves as an identity element.
Now, one difficulty with the above proposal is that we need to continue to have our original formula:
$$(a * b * c) = \frac{abc}{(a b + b c + c a)}$$ Furthermore, we could run into trouble if $$((a b)(c - 1)) + b c + c a) < 0$$ So it seems that we need to try a somewhat different approach.
Observe that $$\frac{abc}{(a b + b c + c a)} = \frac{1}{\frac{1}{a} + \frac{1}{b} + \frac{1}{c}}$$
So when $c$ is equal to the identity element, we could evaluate our expression by considering the limit as $c$ approaches infinity. In other words, we could define:
$$(a \oplus b \oplus e) = \lim\limits_{c \to \infty} (a * b * c)$$
Technically, that's not going to be a definition, but a consequence of some other statement that will be a definition, because $\oplus$ is simply going to be a binary operation, and thus it needs to be defined as such.
Now, the above isn't an answer, but perhaps it is the beginning of an attempt to fulfill the request for ideas.