|
|
A066197
|
|
Squarefree kernel of (n*prime(n))*(n+prime(n)).
|
|
1
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
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.
|
|
MATHEMATICA
|
sfk[n_] := Times @@ FactorInteger[n][[All, 1]];
a[n_] := sfk[n Prime[n] (n+Prime[n])];
|
|
PROG
|
(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
(Haskell)
a066197 n = a007947 $ a033286 n * a014688 n
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|