Open3
ガリレオ列
初項1,公比5のガリレオ列を求める。
十進basicを使っています。メニューよりオプション---文法---microsoft basic互換を選択。
DIM g(70)
LET g(1)=1
LET g(2)=5
FOR a=3 TO 20
IF MOD(a,2)=0 THEN LET g(a)=INT (6*g(a/2)/2)+1: PRINT g(a)
IF MOD(a,2)=1 THEN LET g(a)=INT ((6*g(INT(a/2+1))-1)/2): PRINT g(a)
NEXT A
END
上の式中6は公比5+1=6です。
公比4の場合はガリレオ列は1,4,9,11,22,23,27,28,54,56となります。
ガリレオ列を生成する公式は
a(2n-1)=int(((k+1)-1/2))
a(2*n)=int((k+1)*a(n)/2)+1
(kは公比)
追記
a(2n-1)=int(((k+1)-1/2))は
a(2n-1)=int(((k+1)-1)/2)に訂正。