cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A179017 Odd numbers c such that c*(c^2 - 1)/4 is squarefree.

Original entry on oeis.org

3, 5, 11, 13, 21, 29, 43, 59, 61, 67, 69, 77, 83, 85, 93, 115, 123, 131, 133, 139, 141, 155, 157, 165, 173, 187, 203, 205, 211, 213, 219, 221, 227, 229, 237, 259, 267, 277, 283, 285, 291, 309, 317, 331, 347, 355, 357, 365, 371, 373, 381, 389, 403, 411, 419, 421
Offset: 1

Views

Author

Artur Jasinski, Jun 24 2010

Keywords

Comments

Original title was: "Numbers c such that (c^2-1)c is square free and gcd(c-1,c,c+1)=1", but (c^2-1)c is never squarefree for odd c, and gcd(n,n+1) is always = 1. - M. F. Hasler, Nov 03 2013
These numbers c with distribution a+b=c such that a=(c-1)/2 (see A172186) and b=(c+1)/2 (see A179019) have minimal possible values with function L(a,b,c) = log(c)/log(N(a,b,c)) = log(c)/log((c^2-1)c/4).
This function is minimal orbital in hypothesis (a,b,c).
There are no numbers or distributions which have value L less than log(c)/log((c^2-1)*c/4).
Equivalently, odd squarefree numbers c such that (c^2 - 1)/4 is also squarefree. - Jon E. Schoenfield, Feb 13 2023
The asymptotic density of this sequence is Product_{p prime} (1 - 3/p^2) = A206256 = 0.125486980905... (Tsang, 1985). - Amiram Eldar, Feb 26 2024

Crossrefs

Programs

  • Mathematica
    aa = {}; Do[If[(GCD[x, (x - 1)/2] == 1) && (GCD[x, (x + 1)/2] == 1) && (GCD[(x - 1)/2, (x + 1)/2] == 1), If[SquareFreeQ[(x^2 - 1) x/4], AppendTo[aa, x]]], {x, 2, 1000}]; aa
  • PARI
    forstep(n=3,421,2,issquarefree(n*(n^2-1)/4)&&print1(n",")) \\ M. F. Hasler, Nov 03 2013
    
  • PARI
    is(n)=n%2 && issquarefree(n) && issquarefree(n^2\4) \\ Charles R Greathouse IV, Mar 11 2014

Formula

a(n) = 2*A172186(n) + 1. - Bernard Schott, Mar 06 2023

Extensions

Edited by M. F. Hasler, Nov 03 2013