login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 25 23:24 EDT 2020. Contains 334604 sequences. (Running on oeis4.)