login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A087969 Integers m such that gcd(1+2^m, m^2) is not a square. 2
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 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Table of n, a(n) for n=1..53.

EXAMPLE

n=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: A181431 A249067 A124983 * A044111 A006886 A053816

Adjacent sequences: A087966 A087967 A087968 * A087970 A087971 A087972

KEYWORD

nonn

AUTHOR

Labos Elemer, Sep 22 2003

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 6 10:18 EST 2023. Contains 360104 sequences. (Running on oeis4.)