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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A226165 Squarefree part of A077425(n) (numbers 4*k+1, k>=0, not a square). 3
 5, 13, 17, 21, 29, 33, 37, 41, 5, 53, 57, 61, 65, 69, 73, 77, 85, 89, 93, 97, 101, 105, 109, 113, 13, 5, 129, 133, 137, 141, 145, 149, 17, 157, 161, 165, 173, 177, 181, 185, 21, 193, 197, 201, 205, 209, 213, 217, 221, 229, 233, 237, 241, 5, 249, 253, 257, 29 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS a(n) == 1 (mod 4), n >= 1. This is because 4*k+1, k>=0, not a square, can only have an even number of odd primes of the type 3 (mod 4) with odd exponents in the prime number factorization. The squarefree part of 4*k+1 has then an even number (maybe 0) of primes of the type 3 (mod 4). Examples: a(4) = 21 = 3*7, a(6) = 33 = 3*11. D(n) = A077425(n) are the 1 (mod 4) discriminants of indefinite binary quadratic forms (they are the odd numbers from A079896). sqrt(D(n)) becomes then, up to an integer factor, sqrt(a(n)), which defines a real quadratic number field Q(sqrt(a(n))) with a basis <1, omega(a(n))> for the ring of integers in this field, where omega(a(n)) = (1 + sqrt(a(n))/2. Example: sqrt(D(9)) = sqrt(45) = 3*sqrt(a(9)) = 3*sqrt(5), with omega(5) = (1 + sqrt(5))/2 (the golden section) for Q(sqrt(5)) = Q(omega(5)). LINKS Vincenzo Librandi, Table of n, a(n) for n = 1..1000 FORMULA a(n) = A007913(A077425(n)). MATHEMATICA SquareFreePart[n_] := Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 2]} & /@ FactorInteger[n]); SquareFreePart /@ (4*Range[65] + 1) // DeleteCases[#, 1] & (* Jean-François Alcover, Jun 14 2013 *) PROG (PARI) [core(n) | n <- vector(100, n, 4*n+1), !issquare(n)] \\ Charles R Greathouse IV, Mar 11 2014 CROSSREFS Cf. A077425. Sequence in context: A082700 A212287 A174361 * A166409 A077425 A039955 Adjacent sequences: A226162 A226163 A226164 * A226166 A226167 A226168 KEYWORD nonn,easy AUTHOR Wolfdieter Lang, Jun 14 2013 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.

Last modified January 27 02:41 EST 2023. Contains 359836 sequences. (Running on oeis4.)