The general formula for reflection $P'(x_2,y_2)$ of a point $P(x_1,y_1)$ along a line $ax+by+c=0\ $ is
$$
{x_2-x_1\over a}={y_2-y_1\over b}=-2\left({ax_1+by_1+c\over a^2+b^2}\right)
$$
A general point on the line $px+qy+r=0$ is $\left({t\over p},{r-t\over q}\right)$
Which means the reflection would be $$\left({t\over p}-2a\left({\frac{at}{p}+b(\frac{r-t}{q})+c\over a^2+b^2}\right),{r-t\over q}-2b\left({\frac{at}{p}+b(\frac{r-t}{q})+c\over a^2+b^2}\right)\right)$$
The locus does give a line but it is a massive pain to solve
Its easier to see the reflection property by using the angle between lines
Assuming the parabola to be $y^2=4ax$ and a light ray $y=b$ the point of incidence being $\left({b^2\over4a},b\right)$
The normal at that point is $y=mx-2am-am^3$ where $m=-{y\over2a}=-{b\over2a}$
The slope of the line joining the point of incidence to the focus is $m={4ab\over b^2-4a^2}$
Using routine co-ord geo you can see that the angle bw the normal is and ray is half the angle bw the incident ray and reflected ray