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.

A066197 Squarefree kernel of (n*prime(n))*(n+prime(n)).

Original entry on oeis.org

6, 30, 30, 154, 110, 1482, 714, 114, 138, 11310, 14322, 1554, 3198, 34314, 43710, 7314, 38114, 28914, 109478, 64610, 144102, 175538, 202354, 60342, 59170, 333502, 40170, 22470, 436218, 484770, 622046, 42706, 768570, 817598, 239890, 169422
Offset: 1

Views

Author

Reinhard Zumkeller, Dec 15 2001

Keywords

Examples

			For n=20 we have: A = n = 20, B = A000040(20) = 71, C = A + B = 20 + 71 = 91 and A*B*C = 129220 with squarefree kernel a(20) = 64610 = 2*5*7*13*71.
		

Programs

  • Haskell
    a066197 n = a007947 $ a033286 n * a014688 n
    -- Reinhard Zumkeller, Jul 24 2013
  • Mathematica
    sfk[n_] := Times @@ FactorInteger[n][[All, 1]];
    a[n_] := sfk[n Prime[n] (n+Prime[n])];
    Array[a, 40] (* Jean-François Alcover, Feb 04 2019 *)
  • PARI
    a(n)=my(p=prime(n),f=vecsort(concat(concat(p, factor(n)[,1]), factor(n+p)[,1]),,8)~); prod(i=1,#f,f[i]) \\ Charles R Greathouse IV, Jul 23 2013
    

Formula

a(n) = A007947(A033286(n) * A014688(n)).