Nota. Una pregunta natural es si realmente existen polinomios en las condiciones anteriores (para ser rigurosos, podrían no existir tales polinomios y entonces no tener ningún valor $a+b$). Planteando la ecuación coeficiente a coeficiente y suponiendo que $Q_1(x)$ y $Q_2(x)$ tienen coeficiente $1$ en $x^2$, dejamos como ejercicio ver que tiene que ser \[Q_1(x)=x^2+x-b,\qquad Q_2(x)=x^2-bx-b.\] El resto $R(x)=cx+d$ es un polinomio arbitrario de grado $1$, luego tendríamos las soluciones \[P(x)=x^5-b x^4-(2b+1)x^3+(b+2)bx^2+(b+c)x+(d-b^2)\] para cualesquiera $b,c,d\in\mathbb{R}$ con $b\not\in\{-1,0\}$.
Para ello, sólo hay que ver que $K$ y $L$ tienen la misma potencia respecto de ambas circunferencias, lo que equivale a decir que $LB\cdot LF=LE\cdot LA$ y $KD\cdot KG=KC\cdot KH$. Sin embargo, la primera igualdad se deduce directamente de que los triángulos $LFA$ y $LEB$ son semejantes, mientras que la segunda viene de que $KGC$ y $KDH$ son semejantes. En ambos casos se trata de triángulos rectángulos con un ángulo (no recto) común.

Nota. La demostración vale tanto si el cuadrilátero es cóncavo o convexo, o incluso si no hay cuadrilátero sino solo cuatro puntos dados en el plano. La única salvedad es que $ABCD$ sea un paralelogramo, en cuyo caso $I=J$ (y, por tanto, la recta $IJ$ no está definida) y $K$ y $L$ no existen ya que las rectas cuya intersección los define son paralelas.
Con todo esto ya podemos calcular el área de $AECF$ puesto que esta consta de los dos triángulos equiláteros $ADM$ y $DMC$ (en azul en la figura). Ambos tienen lado $\ell=5$, luego su área es $\frac{\ell^2\sqrt{3}}{4}=\frac{25\sqrt{3}}{4}$ (esto se comprueba fácilmente usando el teorema de Pitágoras). También tenemos que contar los dos triángulos congruentes $DEC$ y $CMF$ (en amarillo en la figura), que tienen la misma altura que los equiláteros pero un tercio de su base, luego cada uno aporta un tercio del área ya calculada $\frac{25\sqrt{3}}{4}$. En resumen, tenemos que
\begin{align*}
\mathrm{Area}(AECF)&=\mathrm{Area}(ADM)+\mathrm{Area}(DMC)+\mathrm{Area}(DEC)+\mathrm{Area}(CMF)\\
&=\frac{25\sqrt{3}}{4}+\frac{25\sqrt{3}}{4}+\frac{1}{3}\cdot\frac{25\sqrt{3}}{4}+\frac{1}{3}\cdot\frac{25\sqrt{3}}{4}=\frac{50\sqrt{3}}{3}.
\end{align*}
Como los enemigos de cualquier persona de $E$ no pueden estar en $E$, el número de enemigos de cada una debe estar entre $1$ (contando a $M$) y $100 - k$. Para que existan $k$ valores distintos en ese rango, debe cumplirse que $k \leq 100 - k$, es decir, $1\leq k \leq 50$.
Lo anterior nos dice que a cada persona se le puede asignar un número entre $0$ y $50$, que es su número de enemigos, y no puede haber tres personas con el mismo número asignado (entre esas tres habría al menos un par de amigos según la regla de enemistad). Más observaciones:
Entonces, el único escollo que nos queda por salvar es que de las $100$ personas haya una con $0$ enemigos, otra con $1$ enemigo y dos con cada número de enemigos entre $2$ y $50$, pero esto también es imposible ya que en tal caso habría un número total de enemigos impar, contradiciendo que si A es enemigo de B, B también es enemigo de A.
Si $p$ divide a $a$, entonces $a \equiv 0 \pmod p$ y la ecuación se escribe $bn + c \equiv x_n^2 \pmod p$. Si $b\not\equiv 0\pmod p$, $bn+c$ toma todos los valores módulo $p$ pero $x_n^2$ no (hay $\frac{p+1}{2}\lt p$ residuos cuadráticos módulo $p$), luego este caso es imposible. Por lo tanto, $b\equiv 0\pmod p$, lluego $\Delta = b^2 - 4ac \equiv 0 \pmod p$ y hemos terminado.
Supongamos entonces que $p$ no divide a $a$. Completando el cuadrado y multiplicando por $4a$, podemos reescribir la ecuación como \[(2an + b)^2 - (b^2 - 4ac) \equiv 4ax_n^2 \pmod p.\] Notemos que $2an+b$ recorre todos los restos módulo $p$, luego escribiendo $u = 2an + b$ y definiendo $A \equiv (4a)^{-1}\not\equiv 0\pmod p$ y $B \equiv -(b^2 - 4ac)(4a)^{-1}\pmod p$, la condición que tenemos se escribe como $Au^2 + B=y_u^2 \in S$ para todo $0\leq u\lt p$, siendo $S$ el conjunto de cuadrados módulo $p$.
La aplicación $f: S \to \mathbb{Z}_p$ dada por $f(s) = As + B$ es inyectiva. Como $f(S) \subseteq S$ y $S$ es finito, se tiene que $f(S) = S$. Sumando la igualdad $Au^2+B=y_u^2$ sobre todos los residuos $u$, obtenemos en el miembro de la izquierda \[\sum_{u=0}^{p-1} (Au^2 + B) = A \sum_{u=0}^{p-1} u^2 + pB \equiv 0 \pmod p.\] Aquí se ha usado que que $\sum_{u=0}^{p-1} u^2\equiv \frac{1}{6}(p-1)p(2p-1)\equiv 0\pmod p$ ya que $p \geq 5$. Para el miembro de la derecha, observamos que, al variar $u$, el cuadrado $u^2$ toma el valor $0$ una vez y cada $s \in S \setminus \{0\}$ dos veces. Dado que $f(S)=S$, nos queda \begin{align*} \sum_{u=0}^{p-1} x_u^2 &=\sum_{u=0}^{p-1} f(u^2)= f(0) + 2 \sum_{s \in S \setminus \{0\}} f(s)\\ & = B + 2 \sum_{s \in S \setminus \{B\}} s = 2 \sum_{s \in S} s - B \equiv -B \pmod p. \end{align*} Igualando ambos resultados, llegamos a que $0 \equiv -B \pmod p$. Dado que $B = -(b^2 - 4ac)(4a)^{-1}$, concluimos que $b^2 - 4ac \equiv 0 \pmod p$, como queríamos demostrar.
Nota. Un atajo que nos puede hacer entender mejor el problema es darse cuenta de que el miembro de la izquierda $f(x)=\sqrt{a-\sqrt{a+x}}$ es una función continua y estrictamente decreciente definida en el intervalo $[0,a^2-a]$ con $f(0)>0$ y $f(a^2-a)=0$ y que el de la derecha $g(x)=x$ es también continua y estrictamente creciente con $g(0)=0$ y $\lim_{x\to\infty}g(x)=+\infty$, luego la ecuación $f(x)=g(x)$ tiene necesariamente una única solución.
Como las cevianas $AP$, $BQ$ y $CD$ son concurrentes en el punto $E$, por el Teorema de Ceva se cumple: \[ \frac{AD}{DB} \cdot \frac{BP}{PC} \cdot \frac{CQ}{QA} = 1.\] Vamos a calcular los cocientes anteriores de otra forma para obtener el resultado.
Utilizando que la razón entre los segmentos de la base de un triángulo es igual a la razón de las áreas de los triángulos que comparten un vértice común, tenemos que \[ \frac{BP}{PC} = \frac{\text{Área}(BDP)}{\text{Área}(CDP)} = \frac{\frac{1}{2} BD \cdot DP \operatorname{sen}(90^\circ - \alpha)}{\frac{1}{2} CD \cdot DP \operatorname{sen} \alpha} = \frac{BD \cos \alpha}{CD \operatorname{sen} \alpha} = \frac{BD}{CD\tan \alpha}. \] Análogamente, \[ \frac{CQ}{QA} = \frac{\text{Área}(CDQ)}{\text{Área}(ADQ)} = \frac{\frac{1}{2} CD \cdot DQ \operatorname{sen} \beta}{\frac{1}{2} AD \cdot DQ \operatorname{sen}(90^\circ - \beta)} = \frac{CD \operatorname{sen} \beta}{AD \cos \beta} = \frac{CD}{AD} \tan \beta. \] Por lo tanto, la igualdad de Ceva nos queda \[ \frac{AD}{BD} \cdot \left( \frac{BD}{CD\tan \alpha} \right) \cdot \left( \frac{CD}{AD} \tan \beta \right) = 1\ \Longleftrightarrow \tan \alpha = \tan \beta.\] Dado que $\alpha$ y $\beta$ son ángulos agudos en un triángulo acutángulo, concluimos que $\alpha = \beta$.
Nota. A la hora de formular las rectas, nos hemos basado en la conocida como ecuación continua de una recta. Esta nos dice que la recta del plano que pasa por los puntos $(a,0)$ y $(0,b)$ tiene ecuación $\frac{x}{a}+\frac{y}{b}=1$ (como puede comprobarse fácilmente).
Sin embargo, tenemos que encontrar un ejemplo en el que no haya filas ni columnas constantes para responder a la segunda parte. Para ello, inspirados por lo anterior, consideraremos $1000$ números diferentes $u_1,\ldots, u_{1000}$ cuya suma sea $1$ y $2024$ números diferentes $v_1,\ldots, v_{2024}$ cuya suma también sea $1$ y pondremos en la casilla de coordenadas $(i,j)$ el número $u_iv_j$. En este caso, la fila $i$-ésima y la columna $j$-ésima suman respectivamente \begin{align*} A_i&=u_iv_1+u_iv_2+\ldots u_iv_{2024}=u_i(v_1+\ldots+v_{2024})=u_i,\\ B_j&=u_1v_j+u_2v_j+\ldots u_{1000}v_{j}=(u_1+\ldots+v_{1000})v_j=v_j, \end{align*} luego $A_iB_j=u_iv_j$ es el elemento colocado en la casilla $(i,j)$.