OFFSET
1,1
COMMENTS
See a comment on A216323 for the definition of an abc-triple, radical r and quality q which is always > 1 by this definition. There also a link is given.
LINKS
Wolfdieter Lang, Maple program for radical of a*b*(a+b) .
FORMULA
a(n) = r(1,b(n),b(n)+1) with b(n) = A216323(n), n>=1, and r(a,b,c) is the radical, also known as squarefree kernel, of a*b*c.
EXAMPLE
2*105 = 2*a(21) = r(1,4374,4375) = 1*6*35 = 210.
MAPLE
read "radabc.txt": [seq(radabc(1, A216324(n)), n=1..42)]/2;
(with the above given link with the maple text file)
MATHEMATICA
rad[n_] := Times @@ Transpose[FactorInteger[n]][[1]]; a = 1; Table[t = {}; mx = 10^n; Do[c = a + b; If[c < mx && GCD[a, b] == 1 && Log[c] > Log[rad[a*b*c]], AppendTo[t, rad[b*c]/2]], {b, a, mx - a}], {n, 5}]; t (* T. D. Noe, Sep 24 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wolfdieter Lang, Sep 24 2012
STATUS
approved