Saturday, June 8, 2013

On Beal's conjecture and related topics

Recently, Beal's conjecture award has raised to $1 million, and this was one of the big headlines in the math community in the past days. This is a very interesting statement in number theory inspired by, as Andy Beal himself says, on Fermat's last theorem

It says that if the positive integers $A,B,C,x,y,z$ satisfy
$$A^x+B^y=C^z$$
with $x,y,z >2$, then $A,B,C$ have a common prime factor.

Notice that this statement work on the basis that there is a solution for the equation, that is, it talks about a property that solutions of the equation must have. This then could be used to easily prove Fermat's last theorem by contradiction using the infinite decent method.

Analyzing these type of equations, it is not hard to find some nice properties arising from them. When I first saw the equation, I was doubtful about the statement, so I went to look for some examples. It was precisely this that led Beal to formulate the conjecture, having looked for numerical solutions with the aid of computer power for values of all the variables up to 1000. So my next inquire was to find a way to create solutions to  these equations.

My approach was to have $x,y,z$ fixed and try to cook up solutions $A,B,C$. The idea is to play around with the prime factors that appear in the game. Basically one has to balance the prime factors in the equation and make them to have the right power.

Grab any pair of positive integers $A,B$ and analyze $A^x+B^y$. The equation requires for this to be a $z$-power, that is, each prime factor appearing in $A^x+B^y$ must have a power multiple of $z$. The idea is to balance those powers by multiplying by powers of $p$.

Suppose that $p$ has a power of $\alpha x$ in $A^x$, $\beta y$ in $B^y$ and $\gamma$ in $A^x+B^y$. Hence, the idea is to multiply by $p^m$ to the expression, where $m$ is such that $\alpha x+m$ is multiple of $x$, $\beta y+m$ is multiple of $y$ and $\gamma+m$ is multiple of $z$. This is a classical problem of chinese remainder, which can be solve by finding solutions to the system of congruences

$$m\equiv 0\text{ mod lcm}(x,y)$$
$$m\equiv -\gamma\text{ mod }z$$.

In order to solve this, we must require that $\text{lcm}(x,y)$ and $z$ to be relative prime. Performing this, we can make $A'=p^{m/x}A$ and $B'=p^{m/y}B$ and then $A'^x+B'^y$ has the right powers for $p$. 

Following this idea for every prime factor $p$ appearing in $A^x+B^y$ leads to a solution of the equation 
$$A^x+B^y=C^z$$.

Therefore, every pair $(A,B)$ leads to a solution of the previous equation. Here then is necessary to make the assumption that $(xy,z)=1$, which makes Fermat's last theorem to fail for this methodology ($x=y=z$).

In general, something similar can be made for $(xz,y)=1$ or $(yz,x)=1$ just by rewriting the equation as 
$$C^z-A^x=B^y$$ or $$C^z-B^y=A^x$$
the only constrain is that $(C,A)$ (respectively $(C,B)$) should be such that the left hand side must be possitive.

This is totally unrelated with Beal's conjecture, as it provides a way to construct solutions, but the method agrees with its statement.   


Friday, June 7, 2013

Solución paramétrica de una elipse

Hace unos días estaba analizando una ecuación diferencial parcial de segundo grado y me topé con la necesidad de encontrar todas las soluciones de una ecuación homogénea de segundo orden

$$ax^2+bxy+cy^2=f.$$

Esta no es más que una elipse rotada, así que las soluciones $(x,y)$ pueden encontrarse de forma paramétrica por medio de las soluciones de una elipse. Un problema bastante rutinario y nada fuera de lo común. Sin embargo, estuve buscando en internet alguna solución paramétrica de una elipse, encontrando únicamente la solución para una elipse normal, no rotada. 

Desde que estaba en los primeros cursos de cálculo, recuerdo el truco de rotar los ejes para eliminar el termino cruzado $xy$ en la ecuación de una cónica, y con esto simplificar la expresión, así que me decidí a escribir la ecuación paramétrica de una elipse general

$$ax^2+bxy+cy^2+dx+ey+f=0.$$

Esta es la ecuación general de una elipse, esto es, una elipse trasladada con centro $(h,k)$ y con ejes rotados por un ángulo $\theta$. 

Sustituyendo $x\mapsto x-h$ e $y\mapsto y-k$, se tiene que la ecuación se puede expresar como

$$a(x-h)^2+b(x-h)(y-k)+c(y-k)^2$$
$$+(f-ah^2-bhk-ck^2)=0,$$

donde 
$$h=\frac{be-2cd}{4ac-b^2}$$ y $$k=\frac{bd-2ae}{4ac-b^2}.$$

Ahora, resta encontrar las soluciones de una ecuación del tipo

$$\alpha x^2+\beta xy+\gamma y^2+\delta=0.$$

Esto se logra por medio de la rotación de ejes

$$x=u \cos(\theta)-v \sin(\theta)$$ e $$y=u\sin(\theta)+v\cos(\theta)$$.

Sustituyendo esto se obtiene que para anular el factor $uv$ es necesario hacer 

$$\theta=\frac{1}{2}\arctan\left(\frac{\beta}{\alpha-\gamma}\right),$$

Luego de esto, se tiene una ecuación del tipo 

$$Au^2+Bv^2=1$$

en donde se puede suponer que $Au^2=\cos(t)$ y $Bv^2=\sin(t)$, con lo que se puede obtener entonces las soluciones paramétricas

$$u=-\frac{1}{\delta}\sqrt{\frac{2\sqrt{(\gamma-\alpha)^2+\beta^2}}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}+(\alpha-\gamma)|\gamma-\alpha|+\beta^2}}$$

$$v=-\frac{1}{\delta}\sqrt{\frac{2\sqrt{(\gamma-\alpha)^2+\beta^2}}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}-(\alpha-\gamma)|\gamma-\alpha|-\beta^2}}$$

Por lo que regresando a $(x,y)$ da

$$x=-\frac{\cos(t)}{\delta}$$
$$\times\sqrt{\frac{\sqrt{(\gamma-\alpha)^2+\beta^2}+|\gamma-\alpha|}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}+(\alpha-\gamma)|\gamma-\alpha|+\beta^2}}$$
$$+\frac{\sin(t)}{\delta}$$
$$\times\sqrt{\frac{\sqrt{(\gamma-\alpha)^2+\beta^2}-|\gamma-\alpha|}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}-(\alpha-\gamma)|\gamma-\alpha|-\beta^2}}$$

$$y=-\frac{\cos(t)}{\delta}$$
$$\times\sqrt{\frac{\sqrt{(\gamma-\alpha)^2+\beta^2}-|\gamma-\alpha|}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}+(\alpha-\gamma)|\gamma-\alpha|+\beta^2}}$$
$$-\frac{\sin(t)}{\delta}$$
$$\times\sqrt{\frac{\sqrt{(\gamma-\alpha)^2+\beta^2}+|\gamma-\alpha|}{(\gamma+\alpha)\sqrt{(\gamma-\alpha)^2+\beta^2}-(\alpha-\gamma)|\gamma-\alpha|-\beta^2}}$$.

Para finalizar, solo basta incorporar $h,k$ y hacer la correspondencia de $\alpha, \beta, \gamma, \delta$ en las ecuaciones, 

$$\alpha=a$$
$$\beta=b$$
$$\gamma=c$$
$$\delta=f-ah^2-bhk-ck^2$$

dando

$$x=-\frac{\cos(t)}{f-ah^2-bhk-ck^2}$$
$$\times\sqrt{\frac{\sqrt{(c-a)^2+b^2}+|c-a|}{(c+a)\sqrt{(c-a)^2+b^2}+(a-c)|c-a|+b^2}}$$
$$+\frac{\sin(t)}{f-ah^2-bhk-ck^2}$$
$$\times\sqrt{\frac{\sqrt{(c-a)^2+b^2}-|c-a|}{(c+a)\sqrt{(c-a)^2+b^2}-(a-c)|c-a|-b^2}}+h$$

$$y=-\frac{\cos(t)}{f-ah^2-bhk-ck^2}$$
$$\times\sqrt{\frac{\sqrt{(c-a)^2+b^2}-|c-a|}{(c+a)\sqrt{(c-a)^2+b^2}+(a-c)|c-a|+b^2}}$$
$$-\frac{\sin(t)}{f-ah^2-bhk-ck^2}$$
$$\times\sqrt{\frac{\sqrt{(c-a)^2+b^2}+|c-a|}{(c+a)\sqrt{(c-a)^2+b^2}-(a-c)|c-a|-b^2}}+k$$.