A087969 Integers m such that gcd(2^m+1, m^2) is not a square.

9, 45, 50, 63, 81, 99, 117, 150, 153, 171, 207, 225, 261, 279, 315, 333, 350, 369, 387, 405, 423, 441, 450, 477, 495, 531, 549, 550, 567, 585, 603, 605, 639, 650, 657, 693, 711, 729, 747, 765, 801, 819, 850, 855, 873, 891, 909, 927, 950, 963, 981, 1014, 1017

Offset

1,1

Links

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

Example

For m = 9: gcd(513, 81) = 27 is not a square number, so 9 is a term.

Mathematica

Do[s=Sqrt[GCD[n^2, 2^n+1]]; If[ !IntegerQ[s], Print[n]], {n, 1, 1000}]

Prog

(PARI) isok(m) = !issquare(gcd(2^m+1, m^2)); \\ Michel Marcus, Aug 27 2019
(Magma) [m:m in [1..1020]|not IsSquare(Gcd(1+2^m, m^2))]; // Marius A. Burtea, Aug 27 2019

Crossrefs

Cf. A087968, A087967.

Sequence in context: A249067 A366420 A124983 * A044111 A006886 A053816

Adjacent sequences: A087966 A087967 A087968 * A087970 A087971 A087972

Keyword

nonn

Author

Labos Elemer, Sep 22 2003

Status

approved