|
|
|
|
7, 13, 17, 31, 37, 59, 61, 65, 77, 83, 89, 97, 107, 127, 129, 131, 133, 145, 153, 165, 169, 179, 195, 197, 201, 221, 227, 235, 245, 249, 261, 269, 281, 293, 297, 303, 321, 325, 345, 369, 373, 381, 393, 399, 405, 409, 417, 421, 425, 427, 442, 443, 447, 455, 465
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MATHEMATICA
|
Block[{s, nn}, nn = 50; s = Map[FactorInteger[#][[-1, 1]] &, 2^Range[4 nn] - 1]; Select[Range[nn], s[[#]] == s[[2 #]] == s[[4 #]] &]] (* Michael De Vlieger, Dec 15 2022 *)
|
|
PROG
|
(PARI) f(n) = vecmax(factor(2^n-1)[, 1]); \\ A005420
isok(n) = my(x=f(2*n)); (f(n) == x) && (f(4*n) == x);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|