If you're multiplying a polynomial P by a linear factor L = (x-a) with a being a nonreal complex number then at least one of these statements is true
1. P does not have real coefficients
2. PL does not have real coefficients
I think you try to avoid this in what you say by multiplying P by a linear term with a complex root whose conjugate is already in P, but if the conjugate is in P then P wasn't real to start with.
The term with the complex root has to "balance" with its conjugate for things to work out. E.g. (x-a)^n (x-a*)^k has real coefficients only when n=k or a is real.
That’s why you have to do a double induction (I’m not a mathematician so my terminology may be wrong) and also keep proving that those special polynomials that are a product of a complex linear term with a polynomial of real coefficients also have a higher (or equal) to 50% chance of having the largest root be real. For degree 2, I think one can prove this directly, and the rest is also an induction.
I don't think this extra claim about polynomials which are a product of a complex linear term and a polynomial with real coefficients is true.
Let a be distributed uniformly in (-1,1) and let b be distributed uniformly in the unit disk (i.e its a complex number of absolute value less than or equal to 1).
The roots of (x-a)(x-b) are (obviously) a and b. The average absolute value of b is 2/3 while the average absolute value of a is 1/2, you can also work out that the probability that |a|<|b| is 2/3, so with probability 2/3 the largest root will be complex.
This just gets worse if you let b be uniformly distributed in the unit square instead of disk, so you need some quite funky distribution on b to make this work which seems to be unjustified.
Thanks. The funky distribution of b has to be exactly the one that results in uniform (or the desired) random distribution of the real coefficients of the 3rd order real polynomial that is the resulting polynomial when you multiply your example polynomial by the linear term with conjugate b. And if this hold for the resulting 3rd order polynomial with real (random) coefficients then I think that the double induction works.
(Edit: the distribution of b has be to such that the distribution of the resulting 3rd order polynomials together with the distribution of 3rd order polynomials with 3 real roots have the desired random distribution of the coefficients.)
We know the desired random distribution of the 3rd order real polynomials (a product of uniform random distributions say), and I believe we can work out the distribution of the subset that have linear roots only and the fraction of random real polynomials that have real roots, so the distribution for the polynomials with roots of type a, b, and conj b, comes from the (weighted) difference of these two, and leads to the joint distribution of a, b (because conj b is fixed given b).
I didn’t want to claim that a full proof is not going to be messy in the details, I just think that this type of construct makes me feel not surprised that the largest root is more likely to be real.
It feels like this problem must have been solved previously, probably in some physics context (maybe in areas related to applications of random matrix theory or dynamical systems).
It is possible that this second recursion doesn’t hold at N=3 (because there is enough of the density covered by the fraction of real polynomials), but then only starts at higher degrees (it will eventually be easier because the N-2 degree will have a higher chance to have a root higher than the new complex b). So this makes is certainly more complicated to prove but still seems intuitively correct to me.
If you're multiplying a polynomial P by a linear factor L = (x-a) with a being a nonreal complex number then at least one of these statements is true
1. P does not have real coefficients
2. PL does not have real coefficients
I think you try to avoid this in what you say by multiplying P by a linear term with a complex root whose conjugate is already in P, but if the conjugate is in P then P wasn't real to start with.
The term with the complex root has to "balance" with its conjugate for things to work out. E.g. (x-a)^n (x-a*)^k has real coefficients only when n=k or a is real.