OFFSET
1,2
COMMENTS
This gives a sparse subsequence of the sequence A006521 of all integers n such that n | 2^n+1. An element of A006521 is said to be *primitive* if it is not divisible by any smaller element of A006521 having the same prime divisors and further it is not the lcm of any two smaller elements of A006521. See Proposition 1 of the link.
Every element of this sequence apart from 1 and 3 is divisible either by 27 or by 171. Stronger results hold. For instance, every element of this sequence apart from 1 and 3 is divisible either by 171 or 243 or 13203 or 2354697 or 10970073 or 22032887841. See the link or A136475 for more details about such results. These alternative factors enable the sequence to be generated much more quickly than by the short Maple program given below.
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..41
Toby Bailey and Chris Smyth, Primitive solutions of n|2^n+1.
EXAMPLE
MAPLE
L:=1: S:={}: for j from 3 by 6 to 10^7 do if not 2&^j+1 mod j = 0 then next end if; if not (j in S) then L := L, j end if; S := S union map( ilcm, S, j ) union {j}; S := S union map(`*`, {map2( op, 1, ifactors(j)[2] )[]}, j); end do: L;
CROSSREFS
KEYWORD
nonn
AUTHOR
Toby Bailey and Christopher J. Smyth, Jan 13 2008
EXTENSIONS
More terms from Max Alekseyev, Aug 04 2011
STATUS
approved