ブラーマグプタの銘言に「数学者とは

　　ｘ^2−９２ｙ^2＝１

を１年以内に解ける人のことである」とあるそうです．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【１】ブラーマグプタの恒等式

　　（ｘ1^2−Ｎｙ1^2）（ｘ2^2−Ｎｙ2^2）＝（ｘ1ｘ2＋Ｎｙ1ｙ2）^2−Ｎ（ｘ1ｙ2＋ｘ2ｙ1）^2

　そして，ブラーマグプタはこの恒等式を使って

　　ｘ^2−９２ｙ^2＝１

を解いたそうです（６２８年）．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

　たとえば，ｘ＝１０，ｙ＝１のとき

　　ｘ^2−９２ｙ^2＝８

比較的小さい値なので，これを使うことにする．

［１］ブラーマグプタの恒等式に

　　Ｎ＝９２，ｘ1＝ｘ2＝１０，ｙ1＝ｙ2＝１

を代入すると

　　８^2＝１９２^2−１９２・２０^2

　　１＝２４^2−９２・（５／２）^2

［２］ブラーマグプタの恒等式に

　　Ｎ＝９２，ｘ1＝ｘ2＝２４，ｙ1＝ｙ2＝５／２

を代入すると

　　１＝１１５１^2−９２・１２０^2

　したがって，（ｘ，ｙ）＝（１１５１，１２０）はペル方程式の整数解となる．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【１】チャクラバーラ法

チャクラバーラ法は近い方程式の解を見つけることから出発する

x=a,y=bがひとつの小さい数値ｋに対するNx^2+k=y^2の解であるとする

たとえばx=1,y=mはNx^2+(m^2-N)=y^2の解である。

このときブラーマグプタの恒等式より

N(am+b)^2+(m^2-n)k=(bm+Na)^2

x=(am+b)/k,y(bm+Na)/kはNx^2+(m^2-N)/k=y^2の解になっている。

ｍをam+bがｋで割り切れるようにとるとbm+Naもｋで割り切れることになり、

x=(am+b)/k,y(bm+Na)/kはNx^2+(m^2-N)/k=y^2の解であって、(m^2-N)/kもまた整数である。

そこで、(m^2-N)/kがｋよりも小さくなるように工夫すれば、この操作の繰り返しによりNx^2+1=y^2の整数解が得られることになる。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝