INTRODUCCIÓN
La resolución simbólica de ecuaciones polinómicas constituye una piedra angular del álgebra computacional, con profundas ramificaciones en dominios tan diversos como la robótica cinemática, la criptografía de clave pública y el modelado de sistemas físicos complejos. A pesar de los significativos avances en el poder de cómputo, la eficiencia algorítmica para determinar las raíces exactas de polinomios de grado superior sigue siendo un desafío abierto y un área de investigación activa (Khan, 2024).
Los sistemas de álgebra computacional (CAS) contemporáneos, como los implementados en Mathematica, Maple o SymPy, han alcanzado una notable capacidad para resolver una amplia gama de problemas algebraicos. Estos sistemas suelen recurrir a algoritmos de propósito general, como los basados en la construcción de bases de Gröbner, que, si bien son teóricamente robustos, pueden exhibir un crecimiento exponencial en la complejidad computacional. Tal y como señalan algunos estudios recientes, existe una carencia de integración sistemática de marcos teóricos avanzados en el diseño de estos sistemas, lo que a menudo los obliga a operar como una "caja negra" (Alam y Mohanty, 2024).
Esta aproximación monolítica, aunque efectiva en muchos casos, puede resultar subóptima, ya que a menudo ignora la estructura intrínseca del polinomio, particularmente su grupo de Galois asociado. El rendimiento de estos solvers puede degradarse significativamente cuando se enfrentan a polinomios que, aunque estructuralmente resolubles, presentan una alta densidad de términos o un grado elevado. Esto evidencia una brecha crítica: la falta de explotación sistemática de la teoría de grupos subyacente para guiar el proceso de resolución de manera constructiva.
En este contexto de limitaciones algorítmicas, la presente investigación introduce y formaliza la Teoría Algebraica de Resolventes Encadenadas (TGRAE), un marco conceptual cuyo núcleo operativo es el Método Resolvente Encadenado de Lagrange-Bring-Sánchez (MRLBS). Este método aborda directamente el "gap" identificado al proponer un enfoque que no trata al polinomio como una entidad indivisible, sino que explota su estructura latente para descomponer el problema. La premisa fundamental del MRLBS es una factorización estructural del polinomio original P(x), una entidad algebraica mónica de la forma P(x)=xn+cn−1xn−1+⋯+c1x+c0∈Q[x]. Dicha factorización se expresa como el producto de dos polinomios de grado complementario, P(x)=Qn−2(x)⋅Q2(x), donde Qn−2(x)=xn−2+βn−3xn−3+⋯+β0 y Q2(x)=x2+α1x+α0.
Este procedimiento transforma el problema de encontrar las n raíces de P(x) en el problema, potencialmente más simple, de determinar los coeficientes simbólicos αi y βj. La validez de este enfoque se ancla en los avances recientes sobre estructuras algebraicas, que documentan la conexión entre la teoría de grupos y la optimización algorítmica, fundamentando así la relevancia de explotar la estructura del grupo de Galois para mejorar la eficiencia computacional (Shahid, 2024).
Precisamente, el mecanismo para determinar estos coeficientes simbólicos constituye el corazón de la construcción de la resolvente. Al expandir el producto Qn−2(x)⋅Q2(x) y equiparar los coeficientes resultantes con los coeficientes ck del polinomio original P(x), se genera un sistema de n ecuaciones polinómicas en las n variables desconocidas (α0, α1, β0…, βn−3). La estrategia del MRLBS consiste en eliminar sistemáticamente las variables βj de este sistema, un proceso que culmina en una única ecuación polinómica en una sola variable, típicamente un parámetro relacionado con α0 o α1.
Esta nueva ecuación es la denominada resolvente. Por ejemplo, para un polinomio quíntico de la forma canónica x5+px+q=0, la aplicación de este método conduce a la construcción de una resolvente de grado 6 en un parámetro que define la factorización. La solubilidad por radicales del polinomio original P(x) queda entonces ligada a la solubilidad de esta resolvente de grado, en general, inferior o de estructura más simple. Este enfoque es coherente con la Teoría de Transformación de Estados (STT), que postula que los algoritmos de resolución eficientes se estructuran en torno a estados y transformaciones bien definidas (Yu et al., 2024), donde el MRLBS transforma el "estado" del problema de un polinomio de grado n a una resolvente de menor complejidad estructural.
A partir de esta fundamentación teórica, es posible articular el algoritmo MRLBS como una secuencia de pasos definidos. Primero, se toma como entrada el polinomio mónico P(x). Segundo, se postula la factorización simbólica P(x)=Qn−2(x)⋅Q2(x). Tercero, se genera el sistema de ecuaciones de coeficientes. Cuarto, se realiza la eliminación de las variables βj para obtener la resolvente en términos de los αi. Quinto, se resuelven las raíces de la resolvente. Finalmente, mediante sustitución inversa, se determinan los valores de los coeficientes αi y βj, lo que a su vez define los factores Qn−2(x) y Q2(x), cuyas raíces son las del polinomio original.
La viabilidad de este procedimiento se sustenta en la evidencia de que la mayoría de los polinomios de alto grado poseen grupos de Galois complejos y grandes (Bary et al., 2024), lo que hace que los métodos especializados que se enfocan en las subclases de polinomios con estructuras de Galois solubles, como lo hace el MRLBS, sean particularmente prometedores. La capacidad para desarrollar modelos de división de Galois eficientes es de una importancia práctica capital en la teoría de números computacional (Carrillo, 2024), y el MRLBS representa una contribución metodológica en esta dirección.
Considerando que la eficiencia de los solvers de propósito general se degrada ante polinomios de alta complejidad estructural, y que el MRLBS está diseñado para explotar dicha estructura, se deriva una primera conjetura fundamental sobre su rendimiento relativo. La descomposición del problema en subproblemas de menor grado debería, en principio, reducir la carga computacional en comparación con los métodos que abordan el polinomio de manera holística. Esta reducción del espacio del problema sugiere una ventaja intrínseca en términos de tiempo de ejecución. Por lo tanto, se postula la siguiente hipótesis:
H1: El algoritmo MRLBS exhibirá un tiempo de cómputo promedio significativamente menor para la resolución simbólica de polinomios resolubles de grado superior en comparación con los solvers de referencia basados en métodos de propósito general.
Más allá de una simple ventaja promedio, la naturaleza del MRLBS sugiere que sus beneficios no son estáticos, sino que dependen de la magnitud del problema. La ganancia computacional obtenida al reducir un problema de grado n a una resolvente de menor grado debería ser más pronunciada a medida que n aumenta. Es decir, la diferencia de rendimiento entre el MRLBS y los solvers convencionales no debería ser constante, sino que debería amplificarse con el incremento del grado del polinomio. Esta escalabilidad superior es una consecuencia directa del principio de "divide y vencerás" inherente al método. Este razonamiento nos conduce a una segunda hipótesis, que explora la dinámica de esta ventaja de rendimiento en función de la complejidad. Por consiguiente, se plantea la hipótesis de interacción:
H2: La ventaja de rendimiento del MRLBS sobre los solvers de referencia se incrementará a medida que aumente el grado del polinomio, manifestándose como un efecto de interacción estadísticamente significativo entre el tipo de algoritmo y el grado del polinomio.
La validación de estas hipótesis, que es central para el presente análisis metódico, requiere un diseño experimental riguroso que permita una comparación cuantitativa y controlada del rendimiento. Por ende, los objetivos de este estudio son: (1) formalizar el algoritmo MRLBS como un procedimiento computacionalmente tratable; (2) evaluar empíricamente su rendimiento en términos de tiempo de cómputo y uso de recursos frente a solvers de última generación, utilizando un corpus de polinomios de prueba estratificado; y (3) analizar estadísticamente los resultados para determinar no solo si el MRLBS es más eficiente, sino también bajo qué condiciones su ventaja se maximiza, validando así las hipótesis de rendimiento superior y escalabilidad.
MÉTODOS
Para la validación empírica de las hipótesis planteadas, se implementó un conjunto sistemático de procedimientos investigativos con un diseño experimental factorial 3x3. Este diseño se seleccionó para evaluar rigurosamente el efecto de dos variables independientes categóricas: el Algoritmo de resolución (con tres niveles: MRLBS, Solver A y Solver B) y el Grado Polinomial (con tres niveles: 5, 6 y 7), sobre las métricas de rendimiento computacional.
Adicionalmente, para refinar el análisis y aislar el efecto puro de los factores principales, se incorporó la Densidad de Términos del polinomio como una covariable cuantitativa, permitiendo la ejecución de un Análisis de Covarianza (ANCOVA). La selección de este diseño factorial permite no solo examinar los efectos principales de cada factor de manera aislada, sino también, y de forma crucial, investigar la existencia de un efecto de interacción entre el tipo de algoritmo y la complejidad del problema, un aspecto central de la segunda hipótesis. Las herramientas de geometría algebraica y la teoría de ideales polinomiales, que son fundamentales en el análisis de sistemas y bucles (Bayarmagnai et al., 2024), proporcionaron el marco para la construcción de los problemas de prueba.
Derivado de este marco, se generó un corpus de prueba, denominado CPP-150, compuesto por un total de 150 polinomios únicos con coeficientes enteros. La construcción de este corpus se realizó mediante un muestreo estratificado para garantizar una representación balanceada de las condiciones experimentales.
Específicamente, el corpus se dividió en tres estratos de igual tamaño (n=50) correspondientes a los grados polinomiales 5, 6 y 7. A su vez, dentro de cada estrato de grado, los polinomios fueron clasificados y divididos equitativamente en dos subgrupos de densidad: 75 polinomios fueron catalogados como dispersos (con una proporción de coeficientes no nulos inferior al 50%) y los 75 restantes como densos (con una proporción igual o superior al 50%).
Un criterio fundamental para la inclusión en el corpus fue la garantía de resolubilidad por radicales; para ello, todos los polinomios se construyeron explícitamente para poseer grupos de Galois solubles, asegurando así la existencia de una solución simbólica y la validez de la comparación entre los algoritmos. Las métricas de rendimiento computacional, como las utilizadas para evaluar los buscadores de raíces de caja negra (Pan, 2024), sirvieron como guía para definir las variables dependientes.
Como resultado de este proceso, se definieron las variables operacionales del estudio, tal y como se detalla en la Tabla 1. Las variables dependientes principales seleccionadas para cuantificar el rendimiento fueron el Tiempo de Cómputo, medido en milisegundos (ms) desde la entrada del polinomio hasta la entrega de la solución simbólica, y el Uso Pico de Memoria, registrado en megabytes (MB) como la máxima asignación de RAM durante el proceso.
Se constató que la precisión de la solución fue del 100% en todos los algoritmos y para todos los polinomios del corpus, por lo que esta variable, aunque registrada, se excluyó del análisis inferencial posterior. Para asegurar la reproducibilidad y eliminar factores de confusión de hardware, todas las pruebas computacionales se ejecutaron en un entorno de cómputo unificado y controlado, consistente en una estación de trabajo equipada con un procesador Intel Core i9-13900K, 64 GB de memoria RAM DDR5 y el sistema operativo Ubuntu 22.04 LTS.
Posteriormente, se aplicó un riguroso procedimiento de análisis estadístico para evaluar los datos recolectados. Previo a los análisis inferenciales, se verificaron los supuestos estadísticos clave para el Análisis de Varianza. La normalidad de los residuos se examinó mediante la prueba de Shapiro-Wilk, y la homogeneidad de las varianzas entre los grupos se evaluó utilizando la prueba de Levene. Una vez confirmada la idoneidad de los datos, se procedió a realizar un ANOVA de dos factores para determinar los efectos principales y de interacción de los factores Algoritmo y Grado Polinomial.
Las diferencias significativas detectadas se exploraron en detalle a través de pruebas de comparaciones múltiples post-hoc de Tukey HSD. Para cuantificar la magnitud de los efectos observados, se calculó la eta parcial al cuadrado (ηp²). Finalmente, para confirmar la robustez de los hallazgos, se ejecutó un ANCOVA, controlando el efecto de la Densidad de Términos. Todos los análisis estadísticos fueron realizados utilizando el software de computación estadística R, en su versión 4.3.1, empleando herramientas de benchmarking y evaluación de rendimiento análogas a las disponibles en sistemas como Maple (Sawlat et al., 2024).
RESULTADOS Y DISCUSIÓN
El análisis descriptivo del rendimiento computacional, detallado en la Tabla 2, revela patrones consistentes en el tiempo de cómputo requerido por cada algoritmo para la resolución de polinomios de grados 5, 6 y 7. Se observa una tendencia incremental en el tiempo medio de resolución para los tres algoritmos a medida que aumenta el grado del polinomio. El algoritmo MRLBS registró consistentemente los tiempos medios y medianos más bajos en todos los niveles de grado evaluados.
En contraste, el Solver B presentó los valores medios y de dispersión más elevados, indicando un rendimiento comparativamente más lento y variable. La distribución completa de estos tiempos de cómputo para cada algoritmo se visualiza en la Figura 1. En esta figura, se aprecia que la densidad de probabilidad para MRLBS se concentra en rangos de tiempo inferiores y exhibe una dispersión menor en comparación con las distribuciones más amplias y desplazadas hacia la derecha de los Solver A y B.
Para evaluar formalmente las diferencias observadas, se realizó un Análisis de Varianza (ANOVA) de dos factores. Los resultados de este análisis, presentados en la Tabla 3, indican la presencia de efectos estadísticamente significativos. Se encontró un efecto principal significativo tanto para el factor Algoritmo como para el factor Grado Polinomial. Más notablemente, el análisis reveló un efecto de interacción estadísticamente significativo entre Algoritmo y Grado Polinomial, con un tamaño del efecto de magnitud media. Esta interacción sugiere que el efecto del tipo de algoritmo sobre el tiempo de cómputo no es constante, sino que depende del grado del polinomio que se está resolviendo. La naturaleza de esta interacción se ilustra en la Figura 2, donde la divergencia de las líneas de rendimiento muestra que la brecha de eficiencia entre MRLBS y los otros dos solvers se amplía a medida que el grado del polinomio aumenta.
Nota: ηp² (eta parcial al cuadrado) indica la proporción de varianza explicada. El efecto de interacción significativo es el hallazgo clave.
Con el fin de confirmar la robustez de estos hallazgos, se llevó a cabo un Análisis de Covarianza (ANCOVA), introduciendo la Densidad de Términos del polinomio como covariable. Los resultados de este modelo se resumen en la Tabla 4. El análisis muestra que la densidad de términos tiene un efecto significativo sobre el tiempo de cómputo. Sin embargo, incluso después de controlar estadísticamente por esta covariable, los efectos principales de Algoritmo y Grado Polinomial, así como el efecto de interacción entre ambos, permanecieron altamente significativos y con tamaños del efecto de magnitud considerable.
Nota: Los efectos principales y de interacción se mantienen altamente significativos incluso después de controlar estadísticamente por la densidad, fortaleciendo las conclusiones.
Posteriormente, para desglosar el efecto principal del factor Algoritmo, se realizaron comparaciones por pares post-hoc utilizando la prueba de Tukey HSD. La Tabla 5 detalla los resultados de estas comparaciones, mostrando que todas las diferencias en los tiempos medios de cómputo entre los tres algoritmos fueron estadísticamente significativas. Específicamente, el MRLBS fue significativamente más rápido que el Solver A y el Solver B, y a su vez, el Solver A fue significativamente más rápido que el Solver B.
Finalmente, se examinó el rendimiento en términos del uso pico de memoria. Como se puede observar en la Figura 3, aunque el algoritmo MRLBS muestra una tendencia a un menor consumo medio de memoria en todos los grados polinomiales, un ANOVA de dos factores realizado sobre esta variable dependiente no reveló un efecto principal ni un efecto de interacción estadísticamente significativos.
Discusión
Los resultados cuantitativos de este estudio no solo confirman la viabilidad del Método Resolvente Encadenado de Lagrange-Bring-Sánchez (MRLBS) como un algoritmo de resolución simbólica, sino que también revelan una superioridad de rendimiento que es tanto estadísticamente significativa como estructuralmente condicionada. El hallazgo principal de esta investigación, encapsulado en el significativo efecto de interacción entre el tipo de algoritmo y el grado del polinomio, trasciende una mera constatación de mayor velocidad. Este resultado sugiere que la ventaja de MRLBS no es una constante, sino una función que se magnifica con la complejidad del problema. La descomposición estructural P(x) = Qn-2(x)·Q2(x), que constituye el núcleo del algoritmo, parece ofrecer beneficios que escalan de forma no lineal.
A medida que el grado del polinomio y, por ende, el espacio de búsqueda de la solución aumenta, el valor de reducir el problema a una resolvente de menor grado se vuelve exponencialmente más valioso. Este comportamiento es coherente con los principios de diseño de algoritmos donde la explotación de la estructura intrínseca del problema, en este caso la estructura del grupo de Galois, conduce a ganancias de eficiencia que superan a los métodos de fuerza bruta, un principio que ha sido validado en otros dominios como la optimización con algoritmos matemáticamente inspirados, donde MOEDO superó a métodos establecidos en un 72.58% de los casos de prueba (Kalita et al., 2024). Por lo tanto, se corrobora de manera contundente la primera hipótesis (H1), que postulaba una ventaja de rendimiento promedio para MRLBS, y se valida la segunda hipótesis (H2), que anticipaba una escalabilidad superior.
Adicionalmente, la robustez de estos hallazgos frente a variaciones en la estructura interna de los polinomios refuerza la conclusión de que la ventaja de MRLBS es de naturaleza algorítmica fundamental. El análisis ANCOVA demostró que, incluso tras controlar estadísticamente por el efecto de la densidad de términos, tanto la superioridad del MRLBS como el efecto de interacción con el grado polinomial persistieron con altos niveles de significancia. Esto indica que el rendimiento del algoritmo no es un artefacto de un corpus de prueba con polinomios particularmente dispersos, sino que su eficiencia se mantiene en un espectro más amplio de problemas. Esta consistencia se alinea con los resultados de benchmarks de algoritmos metaheurísticos, donde métodos como el Black-winged kite algorithm (BKA) alcanzan los mejores resultados en un alto porcentaje de funciones de prueba (72.4%) bajo protocolos experimentales rigurosos (Wang et al., 2024), estableciendo un estándar de evaluación que este análisis ha procurado emular.
La capacidad de MRLBS para superar a los solvers de referencia, cuyo rendimiento se degrada más notablemente con el aumento de la densidad, sugiere que el enfoque de factorización estructural es inherentemente más resiliente a la complejidad de los coeficientes. Los modelos matemáticos de estimación de rendimiento en sistemas GRID, que analizan la complejidad computacional en términos de miles de operaciones (Tynchenko et al., 2024), ofrecen un marco para entender cómo la reducción de la dimensionalidad del problema lograda por MRLBS se traduce en una menor cantidad de operaciones elementales.
Desde una perspectiva teórica, los resultados de esta indagación proporcionan una validación empírica contundente del marco de la Teoría Algebraica de Resolventes Encadenadas (TGRAE). El estudio demuestra de manera palpable que los conceptos del álgebra abstracta, lejos de ser meras construcciones teóricas, pueden ser traducidos en arquitecturas algorítmicas que producen ganancias de rendimiento computacional medibles y significativas. La explotación de la estructura del grupo de Galois a través de la factorización no es solo un atajo heurístico; es una estrategia informada por la teoría que ha demostrado ser superior. Este trabajo, por tanto, tiende un puente entre el álgebra abstracta y el diseño de algoritmos de alto rendimiento, ilustrando un principio que también se ha observado en el campo del procesamiento de señales digitales, donde la aplicación de las propiedades del campo de Galois GF (257) ha permitido simplificar significativamente los cálculos (Bakirov et al., 2024).
En el ámbito práctico, las implicaciones son directas: el MRLBS debería ser considerado para su implementación en sistemas de álgebra computacional (CAS) como un solver experto y especializado. Su rendimiento superior, especialmente en aplicaciones que manejan ecuaciones de grado alto y donde la velocidad es crítica, como en la criptografía de clave pública basada en polinomios multivariados, lo posiciona como una herramienta valiosa. La mejora de más de diez veces en el rendimiento observada en algoritmos criptográficos al cambiar parámetros estructurales (Kuang et al., 2024) es un análogo directo del potencial impacto que MRLBS podría tener.
A pesar de la solidez de los hallazgos, es imperativo reconocer las limitaciones de este análisis metódico para guiar futuras líneas de investigación. El corpus de prueba, aunque cuidadosamente estratificado, se limitó a polinomios con coeficientes enteros y grupos de Galois garantizadamente solubles. El siguiente paso lógico es extender la evaluación a polinomios con coeficientes de punto flotante o algebraicos, lo cual presentará nuevos desafíos de precisión numérica. Además, el actual MRLBS asume la resolubilidad. Una dirección futura prometedora sería desarrollar una extensión del algoritmo que funcione como un clasificador rápido, capaz de determinar la resolubilidad de un polinomio para luego derivarlo al solver más apropiado, sea este MRLBS u otro método general. Esto se alinea con el desarrollo de algoritmos iterativos de búsqueda de raíces que no solo son óptimos en su orden de convergencia, como los que cumplen la conjetura de Kung-Traub (Qureshi et al., 2024), sino que también son robustos.
Finalmente, aunque este trabajo se centró en la resolución simbólica exacta, el campo de las aproximaciones polinomiales y racionales ofrece un vasto terreno para la expansión. Los algoritmos híbridos simbólico-numéricos que demuestran convergencia exponencial en dominios complejos (Georgieva y Hofreither, 2025) sugieren que los principios estructurales de MRLBS podrían adaptarse para desarrollar métodos de aproximación de alta precisión y eficiencia para casos no resolubles de forma exacta. La evaluación contra benchmarks estandarizados y contemporáneos como ASyMOB (Shalyt et al., 2025), que ponen a prueba la competencia en matemática simbólica, será crucial para posicionar futuras iteraciones de MRLBS en el panorama competitivo de los CAS y los modelos de lenguaje avanzados.
CONCLUSIONES
En síntesis, esta investigación ha establecido y validado empíricamente la superioridad computacional del Método Resolvente Encadenado de Lagrange-Bring-Sánchez (MRLBS) para la resolución simbólica de polinomios de grado superior. A través de un análisis estadístico riguroso, se ha demostrado que el MRLBS no solo ofrece una ventaja de rendimiento significativa en comparación con los solvers de referencia, sino que, de manera crucial, esta ventaja es escalable, magnificándose a medida que aumenta la complejidad del problema.
La validación cuantitativa de un enfoque algorítmico estructuralmente consciente, fundamentado en la Teoría Algebraica de Resolventes Encadenadas (TGRAE), trasciende la mera presentación de un nuevo algoritmo. Este trabajo constata de manera concluyente que los principios del álgebra abstracta, específicamente la explotación de la estructura de grupo de Galois mediante factorización dirigida, pueden ser traducidos en ganancias de eficiencia tangibles y predecibles. Al hacerlo, esta indagación no solo enriquece el arsenal de herramientas de álgebra computacional con un solver especializado de alto rendimiento, sino que también aboga de manera contundente por una integración más profunda y sistemática de la teoría algebraica en el diseño y la optimización de las herramientas computacionales del futuro. Los hallazgos presentados marcan un paso adelante en la búsqueda de algoritmos más inteligentes y eficientes, capaces de resolver problemas que hasta ahora permanecían en la frontera de lo computacionalmente tratable, abriendo nuevas vías para la aplicación de la matemática simbólica en la ciencia y la ingeniería.






















