A176258 Numbers of the form 4k+1 with greatest prime divisor of the form 4m-1.

9, 21, 33, 49, 57, 69, 77, 81, 93, 105, 121, 129, 133, 141, 161, 165, 177, 189, 201, 209, 213, 217, 237, 245, 249, 253, 285, 297, 301, 309, 321, 329, 341, 345, 361, 381, 385, 393, 413, 417, 437, 441, 453, 465, 469, 473, 489, 497, 501, 513, 517, 525, 529, 537, 553, 573

Offset

1,1

Comments

All terms of A107978 are in the sequence.

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Mathematica

Select[4 Range@ 150 + 1, Mod[#, 4] == 3 &[FactorInteger[#][[-1, 1]]] &] (* Michael De Vlieger, Feb 07 2016 *)

Prog

(PARI) isok(n) = (n != 1) && ((n % 4) == 1) && ((vecmax(factor(n)[, 1]) % 4) == 3); \\ Michel Marcus, Feb 07 2016

Crossrefs

Cf. A176257, A176256, A176255, A107978, A002148, A002145, A016813, A004767.

Sequence in context: A176256 A017629 A216240 * A107978 A226355 A253052

Adjacent sequences: A176255 A176256 A176257 * A176259 A176260 A176261

Keyword

nonn

Author

Vladimir Shevelev, Apr 13 2010

Extensions

Corrected and extended by Michel Marcus, Feb 07 2016

Status

approved