There should be infinitely many continuous solutions, one for each continuous function $g:[1,2]\to \mathbb{R}$ with $g(1)=g(2)$. After imposing appropriate boundary and differentiability conditions on $g$, we can make the function smooth.
Let $x_1=1$ and $x_{n+1}=x_n+\frac{1}{x_n}$. Then $1\le x_n\le n$ and by the divergence of the harmonic series, $x_n\to\infty$ as $n\to \infty$. Since $h:t\mapsto t+\frac{1}{t}$ is strictly increasing on $[1,\infty)$, each $x\in[1,\infty)$ belongs to exactly one $[x_{n+1},x_{n+2})$ and $x=h^n(y)$ for exactly one $y\in[1,2)$. Then we define $f(x)=g(y)$. Using the relation $f(x)=f(1/x)$, this extends to $(0,\infty)$. It is continuous since it is continuous on each $[x_n,x_{n+1}]$ and agrees at end points.