OFFSET
1,3
COMMENTS
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
EXAMPLE
a(100) = #{1,3,9,11,33,73,99} = 7.
MAPLE
f:= n -> nops(select(`<=`, numtheory:-divisors(n^4-1), n)):
1, seq(f(n), n=2..100); # Robert Israel, Sep 21 2014
MATHEMATICA
Table[Count[Divisors[n^4-1], _?(#<=n&)], {n, 90}] (* Harvey P. Dale, Aug 13 2014 *)
PROG
(PARI) a(n) = if (n==1, 1, sumdiv(n^4-1, d, d <= n)); \\ Michel Marcus, Sep 21 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 09 2007
STATUS
approved