[SeqFan]: Repeated iterations of INVERT starting from A019590 ?

Antti Karttunen antti.karttunen at gmail.com
Wed Jun 7 17:38:01 CEST 2006


Cheers,

concerning the repeated invocations of INVERT-transform
(see http://www.research.att.com/~njas/doc/eigen.ps or 
http://www.research.att.com/~njas/doc/eigen.pdf
for the exact definition of INVERT) to which I already
delved 2 and half years ago:

John Layman wrote Dec 8, 2003:

 >Subject: Some INVERT transforms.
 >From: John Layman <layman at calvin.math.vt.edu>
 >To: seqfan at ext.jussieu.fr
 >Date:   Mon, 8 Dec 2003 11:33:58 -0500 (EST)
 >Cc: layman at calvin.math.vt.edu (John Layman)
 >Message-Id: <20031208163403Z10491-608+64 at calvin.math.vt.edu>

 >Antti Karttunen wrote (Dec 8, 2003): ...
 >>could you compute for me the following transforms:
 >>INVERT([1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]);
...

 >John: Here are some results that I obtained using my own ISAP (Integer 
Sequence Analysis Program),
 >written in Pascal (US denotes the entered USer sequence):

 >US =1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 (25)
 >(INV)US = 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987 1597 2584 4181 
6765 10946 17711 28657 46368 75025 121393 (25)

I am now proposing that if we start from INVERT(A019590),
and then repeatedly apply INVERT(RIGHT(...)) to the result,
(here RIGHT will just insert 1 to the front of sequence)
the resulting sequences will converge towards A000108, Catalan numbers.
This purely on combinatorial grounds, as a side product of "Gatomorphology".
Note that A000108 is an eigen-sequence regarding to this
operation, i.e. it stays exactly same after INVERT(RIGHT(...)).
(Stronger conjecture is that such an iteration procedure, starting from
_almost any_ sequence, eventually converges towards A000108.)


I.e. INVERT(A019590) produces Fibonacci numbers beginning from A000045(2),
and when we apply INVERT(RIGHT(...)) to them, i.e.
INVERT([1,1,2,3,5,8,13,21,34,55,89,...])
we get Pell numbers, A000129
1, 2, 5, 12, 29, 70, 169, 408, 985, 2378, 5741,
and again, prepending one, and taking INVERT:
INVERT([1, 1, 2, 5, 12, 29, 70, 169, 408, 985, 2378, 5741, ...]
produces what?

Could somebody compute that, and continue few iterations more?
(Remember the prepended 1 and be careful with those indices,
as INVERT-transform doesn't/shouldn't actually work on any 0-indiced terms!)

This might make a nice table, with some other nice properties as well.

Also, if anybody knows a pointer to a free INVERT-procedure/algorithm,
then let me know, then I can compute it by myself.


Terveisin,

Antti






More information about the SeqFan mailing list