OFFSET
1,4
COMMENTS
Here a proper divisor d of n is a divisor of n such that 1 < d < n.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = 0 if n is not composite, p if n is a proper power of prime p, and 1 otherwise. - Franklin T. Adams-Watters, Mar 22 2011
MAPLE
A182936 := n -> igcd(op(numtheory[divisors](n) minus {1, n}));
seq(A182936(i), i=1..79); # Peter Luschny, Mar 22 2011
MATHEMATICA
Join[{0}, Table[GCD@@Most[Rest[Divisors[n]]], {n, 2, 110}]] (* Harvey P. Dale, May 04 2018 *)
PROG
(PARI) A182936(n) = { my(divs=divisors(n)); if(#divs<3, 0, gcd(vector(numdiv(n)-2, k, divs[k+1]))); }; \\ Antti Karttunen, Sep 23 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Mar 22 2011
EXTENSIONS
More terms from Antti Karttunen, Sep 23 2017
STATUS
approved