[seqfan] Re: A new sequence

Maximilian Hasler maximilian.hasler at gmail.com
Tue Feb 5 20:41:07 CET 2013


> Wikipedia says a "characteristic polynomial" is connected to a matrix.
> Do you use the matrix to generate the sequence (and of so, how)?

Here, that one is meant:

http://en.wikipedia.org/wiki/Recurrence_relation#Linear_homogeneous_recurrence_relations_with_constant_coefficients

although both notions can be linked...

> I guess the more general question is: Do you have software that can
> automatically generate the 1000 or so divisibility sequences?

Here, a script which appends n values to
the r initial values  of a lin.rec. of order r,
already helps:

(PARI)
{linrec(c,a,n=10)=for(k=1,n,a=concat(a,-vecextract(a,Str(-#c,".."))*c~));a}

linrec([1,-40,206,-40],[1,0,1,48])

=> [1, 0, 1, 48, 1713, 58672, 1995921, 67818880, 2303940641,
78266714528, 2658767568353, 90319849348112, 3068216219484241,...]
\\ NOT IN OEIS ;
 but since highly composite, it is quite probable that it is the
product of two or more sequences in the OEIS

gp > linrec([1,-40,206,-40],[-1,0,1,6])
=>
-1, 0, 1, 6, 35, 204, 1189, 6930, 40391, 235416, 1372105, 7997214,
46611179, 271669860,...
\\ this is A001109

gp > linrec([1,-40,206,-40],[-1,0,1,34])
 =>
-1, 0, 1, 34, 1155, 39236, 1332869, 45278310, 1538129671, 52251130504,
1775000307465, 60297759323306, ...

\\this is A091761 (from the 0 on) resp. A029547 (from the 1 on)

gp > linrec([1,-40,206,-40],[-1,0,1,48])
 ms.
-1, 0, 1, 48, 1715, 58752, 1998709, 67914000, 2307174311, 78376578048,
2662499775145, 90446634986352,

\\ NOT in the OEIS, but same remark as above.

Maximilian



More information about the SeqFan mailing list