0,2

Two other unrelated methods give the same sequence: A047241, A086515 The lucky numbers were found by a sieve in True basic used in my previous lucky number sequences and transfered to Mathematica. It is interesting that the matrix that gives a Prime like Markov is the same as this one, so that in M<arkovian terms they are closely related.

A=3 X 3 matrix Lucky[n+m+1] B=3 X 3 matrix Lucky[n+m+2] M=B.A^(-1) a(n) = unique elements of M^n*A

a= {1, 3, 7, 9, 13, 15, 21, 25, 31, 33, 37, 43, 49, 51, 63, 67, 69, 73, 75, 79, 87, 93, 99, 105, 111, 115, 127, 129, 133, 135, 141, 151, 159, 163, 169, 171, 189, 193, 195, 201, 205, 211, 219, 223, 231, 235, 237, 241, 259, 261, 267, 273, 283, 285, 289, 297, 303, 307, 319, 321, 327, 331, 339, 349, 357, 361, 367, 385, 391, 393, 399, 409, 415, 421, 427, 429, 433, 451, 463, 475, 477, 483, 487, 489, 495, 511, 517, 519, 529, 535, 537, 541, 553, 559, 577, 579, 583, 591, 601, 613, 615 } aa = Table[a[[n + m + 1]], {n, 0, 2}, {m, 0, 2}] bb = Table[a[[n + m + 2]], {n, 0, 2}, {m, 0, 2}] M = bb.MatrixPower[aa, -1] c = Union[Flatten[Table[MatrixPower[M, n].aa, {n, 0, 200}]]]

Cf. A047241, A086515.

Sequence in context: A087550 A047241 A086515 this_sequence A132222 A111225 A032678

Adjacent sequences: A103137 A103138 A103139 this_sequence A103141 A103142 A103143

nonn,uned

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Mar 17 2005