The equation $y_0=px_0+f(p)$ need not have a solution $p$ for all initial values $(x_0,y_0)$. In these points there is no solution of the Clairaut differential equation.
Assume now that the initial point $(x_0,y_0)$ is well-behaved, that is a $p_0$ with $y_0=x_0p_0+f(p_0)$ exists and $x_0+f'(p_0)\ne 0$. Then by the implicit function theorem for the function
$$
F(x,y,p)=px+f(p)-y
$$
at $(x,y)\approx(x_0,y_0)$ the equation $$F=0\iff y=xp+f(p)$$ is solvable for $p\approx p_0$, and the local solution $p=g(x,y)$ is as smooth as $f$. That is, as $f\in C^k$, with $k\ge 1$, giving $F\in C^k$, so is $g\in C^k$ (locally).
The solution of $y'=g(x,y)$ now in consequence is $C^{k+1}$, so at least $C^2$, as long as it stays in the neighborhood of $(x_0,y_0)$.
Note that where the linear solutions $y=cx+f(c)$ meet the singular solution of $x+f'(y')=0$, which can be parametrized as $$x(p)=-f'(p), ~~ y(p)=f(p)-f'(p)p,$$ one can switch the branches of solutions and at that point the second derivative jumps.