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

Author

Labos Elemer, Sep 22 2003

Keywords

• nonn

Examples

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

Links

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

Crossrefs

Cf. A087968, A087967.

Programs

Magma

[m:m in [1..1020]|not IsSquare(Gcd(1+2^m,m^2))]; // Marius A. Burtea, Aug 27 2019

Mathematica

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

PARI

isok(m) = !issquare(gcd(2^m+1, m^2)); \\ Michel Marcus, Aug 27 2019