|
编程算 1+2^(1/2),(用变通的3x+1方法),发散好快.
对前面第十五段的补充。(本帖中“奇元”“偶元”(或者“奇的”“偶的”)的含义在“奇偶数概念的推广”这篇帖子中讨论过,仅供参考)。
(一).在2^(1/2)生成的整环内,以2^(1/2)作除数(即u=2^(1/2)),从( 13)+( 6)*g算起,
(注:g=2^(1/2))
(方法是乘3加1,再用2^(1/2)除成奇的,……)
有如下一个循环圈:
( 13)+( 6)*g m( 1 )= 3 ( 9)+( 10)*g m( 2 )= 3 ( 15)+( 7)*g m( 3 )= 1 ( 21)+( 23)*g m( 4 )= 1 ( 69)+( 32)*g m( 5 )=8 ( 13)+( 6)*g
此圈有5个奇元。
以下是计算过程:
3*(13+6*g)+1=(g^3)*(9+10*g)
3*(9+10*g)+1=(g^3)*(15+7*g)
3*(15+7*g)+1=g*(21+23*g)
3*(21+23*g)+1=g*(69+32*g)
3*(69+32*g)+1=(g^8)*(13+6*g)
方程组如下(g = 2^(1/2)).
3*x(1)+1=(g^m(1))*x(2)
3*x(2)+1=(g^m(2))*x(3)
3*x(3)+1=(g^m(3))*x(4)
3*x(4)+1=(g^m(4))*x(5)
3*x(5)+1=(g^m(5))*x(1)
这里x(1)解公式为
a = 3^(n-1) + 3^(n-2)*(g^(m(1))) +3^(n-3)*(g^(m(1)+m(2))) +……+3*g^(m(1)+m(2)+m(3)+……+m(n-2) ) + g^(m(1)+m(2)+……+m(n-1))
b = ( g^(m(1)+m(2)+……+m(n-1)+m(n) ) - 3^n )
x(1) =a / b
令m(1)=3,m(2)=3,m(3)=1,m(4)=1,m(5)=8,n=5
(5个奇元,故n=5)
代入解公式就算出x(1)=13+6*g
后面的解类似。
相对于前面十五段,找了个循环圈,但不知道1+2^(1/2)用推广3x+1规则(乘3加1,再用2^(1/2)除成奇的,……),会不会算到这个循环圈内。
(二). 用方程组的解公式找循环圈,对形如s+t*(2^(1/2))的代数整数,(s和t为普通整数),如果它是某个循环圈里的奇代数整数,则s和t同号,即s*t>0,已找到的几个循环圈如(-3-2*2^(1/2))和(-3-2^(1/2),-3-4*2^(1/2)),包括上面的有5个奇元的循环圈,均是同号。(证明很简单略)
而对于一个形如s+t*(2^(1/2))的代数整数,(s和t为普通整数),如果s和t异号, (如1-2^(1/2)),用推广3x+1规则(乘3加1,再用2^(1/2)除成奇的,……),算下去s和t总是异号,例如用1-2^(1/2)算下去,前6个奇元如下(g=1-2^(1/2))
( 1)+(-1)*g m( 1 )= 1 (-3)+( 2)*g m( 2 )= 3 ( 3)+(-2)*g m( 3 )= 2 ( 5)+(-3)*g m( 4 )= 1 (-9)+( 8)*g m( 5 )= 2 (-13)+( 12)*g
再算下去 s和t 全部是一正一负的,而前面已指出循环圈内奇元的s和t是同号的,故1-2^(1/2)这样算下去是不进入循环圈的。
(也许x^2-x-2=0的两个根很有意思! )
[ 本帖最后由 laodiao8014 于 2009-8-25 14:04 编辑 ] |
|