■ブラーマグプタ・バースカラ・チャクラバーラ法(その118)

  67x^2+1=y^2の解は、x=5967,y=48842ともともと小さくない

D・1^2+(m^2-D)=m^2

Du^2+b=v^2

を合成すれば

D(mu+v)^2+b(m^2-D)=(Du+mv)^2が得られる。

このとき、mu+vがbの倍数ならば

D{(mu+v)/b}^2+(m^2-D)/b={(Du+mv)/b}^2が得られる。

u1=(mu+v)/bが整数ならば、b1=(m^2-D)/bもv1=(Du+mv)/b=+/-(Du1^2+b1)^1/2も整数

===================================

  67x^2+1=y^2

[1]任意に平方数、たとえば1^2=1をとり67・1^2+bがは平方数となるようなbを目の子で見つける。67・1^2-3=8^2

[2]b=-3とすれば67・1^2-3=64=8^2,u=1,v=8

[3]um+v=bn,この例ではm+8=-3n,m^2がなるべくDに近くなるようにmをとる。m=7とする。(m=8はm+8が3の倍数にならないのでNG)

[4]b1=(m^2-D)/b=(49-67)/(-3)=6とする,u1=|(um+v)/b|=|(1・7+8)/(-3))|=|-5|=5,v1=|(Du+mv)/b|=|123/(-3)|=41,(Du1^2+b1)^1/2=(67・25+6)^1/2=41

[5]67・5^2+6=41^2に対して、この手順を繰り返す。u=5,v=41,b=6

[6]5m+41=6n,|m^2-67|が小さいm→m=5(m=8は5m+41が6の倍数にならないのでNG)

[7]b1=(m^2-D)/b=(25-67)/(6)=-7,u1=|(um+v)/b|=|(5・5+41)/(6))|=11,v1=|(Du+mv)/b|=|(67・5+5・41)/6|=90,=((Du1^2+b1)^1/2=(67・121-7)^1/2=90^2

[8]67・11^2-7=90^2,u=11,v=90,b=-7

[9]11m+90=-7n,|m^2-67|が小さいm→m=9(m=8は11m+90が7の倍数にならないのでNG)

[10]b1=(m^2-D)/b=(81-67)/(-7)=-2,u1=|(um+v)/b|=|(11・9+90)/(-7))|=27,v1=|(Du+mv)/b|=|(67・11+9・90)/(-7)|=221,(Du1^2+b1)^1/2=(67・121-2)^1/2=221

[11]67・27^2-2=221^2

Du0^2+c0=v0^2

Du1^2+c1=v1^2のとき

D(u0v1+u1v0)^2+c0c1=(Du0u1+v0v1)^2が成り立つ

67(27・221+27・221)+4=(67・27・27+221・221)^2

67(11936)^2+4=(97684)^2

67(5967)^2+1=(48842)^2・・・x=5967,y=48842

===================================