Not a reference. Just checking my comment.
I wanted to be sure I had the statement correct. We are asking about roots and factoring for $a x^2 + b x + c.$ Once we find that $b^2 - 4ac$ is a perfect square $N^2,$ we get roots $\frac{-b \pm N}{2a},$ or $\frac{-P}{a}, \; \frac{-Q}{a},$ with $PQ = a c$ and $P+Q = b.$ This does use the fact that
$ N \equiv b \pmod 2.$ We begin with the rational factorization $$ \frac{1}{a} ( a x + P) (a x + Q) $$ with all integers. Next, take
$$g = \gcd(a,P), \; \; a = g \alpha, \; \; P = g \pi$$
so that $\gcd(\alpha, \pi) = 1.$
Well, from $PQ = ac,$ we get $P|ac,$ then $g \pi | g \alpha c,$ then $\pi | \alpha c.$ However, from $\gcd(\alpha, \pi) = 1$ we find
$\pi | c,$ and we may let $$ c = \gamma \pi. $$
Comparing $PQ = ac,$ we find $Q = \alpha \gamma,$ so our integral factorization is
$$ a x^2 + b x + c = (\alpha x + \pi) (g x + \gamma). $$
If, in addition, we have $\gcd(a,b,c) = 1,$ it follows that $\gcd(g,\gamma) = 1.$ It also follows that $\gcd(a,Q) = \alpha.$
In this case, it is correct to say that we could have calculated the roots with the quadratic formula, written those fractions in lowest terms, and used those fractions to write the factorization directly.