Java - Algoritmo factorización polinomios

/**

* Factor Común

* 

* @param p

*            Polinomio a sacar factor común

*

* @return int Número de factores x comunes a los monomios de p (grado de x)

*/

public static int factorcomun(float[]p){

	int factoresx=0;

	while(p[0]==0){

		for(int i=1;i<p.length;i++){

			p[i-1]=p[i];

		}

		factoresx++;

	}

	return factoresx;

}

/**

* Raíces de un polinomio

* 

* @param p

*            Polinomio a obtener raíces

*

* @return float[] Vector de reales cuyas componentes son las raíces del polinomio p

*/

public static float[] raices(float[]p){

	//GRADO POLINOMIO P

	int grado=0;

	for(int i=0;i<p.length;i++){

		if (p[i]!=0){grado=i;}

	}

	//CALCULO DIVISORES

	float divisores[]=new float[2*Math.abs((int)p[0])];

	float divisor=1;

	int i=0;

	while(divisor<Math.abs(p[0])){

		if(p[0]%divisor==0){

			divisores[i]=divisor;

			divisores[i+1]=-divisor;

		}

		divisor++;

		i=i+2;

	}

	//CALCULO RAICES

	float raices[]=new float[grado];

	int nraices=0;

	int k=0;

	for(int j=0;j<divisores.length;j++){

	if(divisores[j]!=0){

		if(Polinomio.valor(p,divisores[j])==0){

			raices[k]=divisores[j];

			nraices++;

			k++;

			}

		}

	}

	return raices;

}

/**

* Factorización de un polinomio

* @see Utiliza el método Ruffini

* @see Muestra por pantalla el proceso seguido para la factorización por ruffini

* Factorización: 

* Polinomio a factorizar = 2,0X^5+X^4-8,0X^3-X^2+6,0X

* 1. Factor Común

* X(2,0X^4+X^3-8,0X^2-X+6,0)

* 2. Ruffini

* X(2,0X^4+X^3-8,0X^2-X+6,0)

* 		|	2.0		1.0		-8.0	-1.0	6.0	

* 1	|			2.0		3.0		-5.0	-6.0	

* ___________________________________________________________

* 		|	2.0		3.0		-5.0	-6.0	0.0	

* 

* 		|	2.0		3.0		-5.0	-6.0	

* -1	|			-2.0	-1.0	6.0	

* ___________________________________________________________

* 		|	2.0		1.0		-6.0	0.0	

* 

* 		|	2.0		1.0		-6.0	

* -2	|			-4.0	6.0	

* ___________________________________________________________

* 		|	2.0		-3.0	0.0	

* 	

* 	X (X-1,0)(X+1,0)(X+2,0)(2,0X-3,0)

* @param p polinomio a factorizar

* @return float [][] array de polinomios que representan la factorizacion de p

*/

public static float[][] factorizacion(float[] p)

{	int gradomaximo=0;

	float raices[]=raices(p);

	int numeroraices=0;

	for(int i=0;i<raices.length;i++){

		if(raices[i]!=0){

			numeroraices++;

		}

	}

	float factorizacion[][]=new float[3*numeroraices][];

	int Ruffini=0;

	int fila=0;

	while(fila<numeroraices){

		//1ªFila

		factorizacion[fila]=new float[gradomaximo+1-Ruffini];

		for(int i=0;i<factorizacion[fila].length;i++){

			factorizacion[fila][i]=p[i];

		}

		fila=fila+2;

		//3ªFila

		factorizacion[fila]=new float[gradomaximo-Ruffini];

		for(int i=factorizacion[fila].length;i>0;i--){

			if(i==factorizacion[fila].length){

				factorizacion[fila][i]=factorizacion[fila-2][i+1];

			}

			else{factorizacion[fila][i]=factorizacion[fila-2][i+1]+raices[Ruffini]*factorizacion[fila][i+1];

			}

		}

		fila--;

		//2ªFila

		factorizacion[fila]=new float[gradomaximo-Ruffini];

		for(int i=factorizacion[fila].length;i>0;i--){

			if(i==0){

				factorizacion[fila][i]=0-factorizacion[fila-1][i];

			}

			else factorizacion[fila][i]=factorizacion[fila+1][i-1]-factorizacion[fila-1][i];

		}

		p=factorizacion[fila];

		fila=fila+2;

	}

	return factorizacion;

}