The line between two points $(x_0,y_0,z_0)$ and $(x_1,y_1,z_1)$ can be written as follows:
$$
\frac{X-x_0}{x_1-x_0}=\frac{Y-y_0}{y_1-y_0}=\frac{Z-z_0}{z_1-z_0}
$$
Therefore all three fractions should be equal, let's say $t$. This can be used to parametrize the line with $t$. We get:
$$
{X}=x_0+t({x_1-x_0})\\
Y=y_0+t({y_1-y_0})\\
{Z=z_0}+t({z_1-z_0})
$$
Therefore for two points $(0,0,1)$ and $(x,y,0)$ we get the following line:
$$
\frac{X}{x}=\frac{Y}{y}=\frac{Z-1}{-1}=t
$$
Therefore the paramterization of the line is $(tx,ty,1-t)$. To obtain the points on the line and on the sphere, $(tx,ty,1-t)$ should satisfy the equation of sphere. Putting this into the equation of sphere, we get:
$$
t^2(x^2+y^2)+(1-t)^2=1\implies t=0 \quad or \quad t=\frac{2}{x^2+y^2+1}
$$
which gives us two points:
$$
(0,0,1) \quad (\frac{2x}{x^2+y^2+1},\frac{2y}{x^2+y^2+1},1-\frac{2}{x^2+y^2+1})
$$