Your formula is inconsistent. You probably meant:
$$½ \left(τ\left(0, 0, 0, t\right) + τ\left(0, 0, 0, t + \frac{x}{c + v} + \frac{x}{c - v}\right)\right) = τ\left(x, 0, 0, t + \frac{x}{c - v}\right),$$
and you can drop the primes here and just write $x'$ as $x$ since the prime is not relevant to the discussion here. You probably forgot that extra $t$ on the right-hand side.
A couple things you should be aware of (and that includes you, too, Einstein wherever you are!) The Taylor Series, properly done, is not an approximation formula, but an exact representation. So, when we say
$$f(x, y) = f(0, 0) + x f_x(0, 0) + y f_y(0, 0) + ⋯$$
what we're actually saying is that
$$f(x, y) = f(0, 0) + x f_x(0, 0) + y f_y(0, 0) + \frac{x^2}{2} A(x, y) + x y B(x, y) + \frac{y^2}{2} C(x, y),$$
for some continuous functions $A$, $B$, $C$ of $(x,y)$, provided that the second derivative of $f$ (i.e. its Hessian matrix) is absolutely continuous in some connected domain that includes both the points $(0, 0)$ and $(x, y)$. Normally, in the Physics literature, one assumes that functions are "smooth" (i.e. that they have continuous derivatives to all orders), so that the "minimum condition" technicalities can be collated together in the simpler "smoothness" condition. In that case, the functions $A$, $B$ and $C$ are also smooth in $(x, y)$. In addition, the theorem will also provide you an integral formula for these remainder functions in terms of the second order partial differentials of $f$.
Ignore the other coordinates and just write $τ$ as a function of $(x,t)$ only. The second order Taylor Theorem is
$$τ(x,t) = τ(0,0) + x τ_x(0,0) + t τ_t(0,0) + \frac{x^2}{2} A(x,t) + x t B(x,t) + \frac{t^2}{2} C(x,t)$$
for all $(x,t)$ is some connected domain that includes $(0,0)$, where $A$, $B$ and $C$ are smooth functions of $(x,t)$. We assume that the domain is for all $(x,t)$ in some neighborhood of $(0,0)$.
Substitute this into your condition (in the form which I modified it as) to write:
$$½ \left(τ(0,0) + t τ_t(0,0) + \frac{t^2}{2} C(0,t) + τ(0,0) + t_2 τ_t(0,0) + \frac{{t_2}^2}{2} C\left(0,t_2\right)\right) \\
= τ(0,0) + x τ_x(0,0) + t_1 τ_t(0,0) + \frac{x^2}{2} A(x,t_1) + x t_1 B(x,t_1) + \frac{{t_1}^2}{2} C(x,t_1),$$
where, for convenience, I've set
$$t_1 = t + \frac{x}{c - v}, \hspace 1em t_2 = t_1 + \frac{x}{c + v} = t + \frac{x}{c - v} + \frac{x}{c + v}.$$
That's not an approximation. That's an exact formula - as they are for all applications of Taylor's Theorem, when the theorem is used right. In reality, you're supposed to pay attention to those remainder terms, because the expressions involving them contain structural information about the form of the remainder (e.g. the multipliers of $A$, $B$ and $C$) that may be relevant to the physics at hand.
Simplify to get the equation:
$$\left(t - 2 t_1 + t_2\right) τ_t(0,0) + \frac{t^2}{2} C(0,t) + \frac{{t_2}^2}{2} C\left(0,t_2\right) \\
= 2 x τ_x(0,0) + x^2 A(x,t_1) + 2 x t_1 B(x,t_1) + {t_1}^2 C(x,t_1),$$
and note that
$$t - 2 t_1 + t_2 = \frac{x}{c + v} - \frac{x}{c - v} = -\frac{2vx}{c^2 - v^2}.$$
Thus
$$-\frac{2vx}{c^2 - v^2} τ_t(0,0) + \frac{t^2}{2} C(0,t) + \frac{{t_2}^2}{2} C\left(0,t_2\right) \\
= 2 x τ_x(0,0) + x^2 A(x,t_1) + 2 x t_1 B(x,t_1) + {t_1}^2 C(x,t_1).$$
This holds for all $(x,t)$ in the said neighborhood of $(0,0)$, which includes $t = 0$, provided $x$ is small enough. So, set $t = 0$, and assume $x$ is small enough that $(x,0)$ is in the neighborhood of $(0,0)$. The result is:
$$-\frac{2vx}{c^2 - v^2} τ_t(0,0) + \frac{{x_2}^2}{2} C\left(0,x_2\right) = 2 x τ_x(0,0) + x^2 A(x,x_1) + 2 x x_1 B(x,x_1) + {x_1}^2 C(x,x_1),$$
where
$$x_1 = \frac{x}{c - v}, \hspace 1em x_2 = \frac{x}{c - v} + \frac{x}{c + v}.$$
Factor out $x$:
$$-\frac{2v}{c^2 - v^2} τ_t(0,0) + x \frac{{z_2}^2}{2} C\left(0,x_2\right) = 2 τ_x(0,0) + x A(x,x_1) + 2 x_1 B(x,x_1) + x {z_1}^2 C(x,x_1),$$
where
$$z_1 = \frac{1}{c - v}, \hspace 1em z_2 = \frac{1}{c - v} + \frac{1}{c + v}.$$
Set $x = 0$, factor out a $2$ and move everything over to one side:
$$τ_x(0,0) + \frac{v}{c^2 - v^2} τ_t(0,0) = 0.$$
What Einstein didn't realize is that you don't need any calculus at all to derive the relevant identities. The reason is that Minkowski geometry is so utterly rigid that from the relation $A ⇔ B$ defined as "space-time point $A$ is at a light-like separation from $B$", and from the primitive notion of space-time point, with a suitable set of purely geometric axioms concerning the "$⇔$" relation (and nothing more), it is possible to define all geometric relations and concepts(!) (including parallelism, orthogonality, angles, lines, planes, congruence, similarity, coordinate axes, etc.) and to characterize Minkowski space exactly up to the choice of a scale unit: the "second", and an orientation of time-direction (i.e. the selection of which direction points to the future and which points to the past).
Thus, no transformation of the desired type is possible that is not linear (or else unbounded). I'm not even sure if continuity needs to be assumed. Moreover, it must not only be linear, but must produce the result corresponding to the Lorentz transforms underlying Minkowski geometry.
For the less rigid family of curved Lorentzian (i.e. 3+1 dimensional) space-time geometries, the specification of the "$⇔$" relation lays out the structure of the light-cones. From the light-cones, alone, one has a space-time metric defined, uniquely, up to a conformal transformation.
Among other things, this uniquely specifies the Weyl Tensor in geometries of 4 or more dimensions, since the Weyl Tensor is conformally invariant. The importance of that is that the wave modes of the Weyl Tensor are the modes of gravitational radiation.
So, those who advance the idea that gravitational radiation is to be quantized are caught in the trap of also requiring the Weyl Tensor to be a quantum object - subject to quantum fluctuations. That means: light cone fluctuations, the ability to do light-cone tunneling and to violate causality and potentially to even do time travel. That's a heavy millstone to hang on anyone who wants to do quantum gravity and who thinks it can be done, at all; and, as I see it, it's one small step shy of a conditional impossibility proof: that either no quantum gravity or else causality-violation.