OFFSET
1,3
FORMULA
a(1) = 1, a(2) = 1, a(n) = (1/n)*((n+1)^{n-2} + sum_{k=1}^{n-1} (n+1)^{gcd(n,k)-1}) for n > 2
PROG
(Sage) [1, 1] + [((n+1)^(n-2) + sum([(n+1)^(gcd(n, k) -1) for k in [1..n-1]]))/n for n in [3..20]]
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Nikos Apostolakis, Mar 10 2009
EXTENSIONS
Corrected the formula and Sage code - Nikos Apostolakis, Feb 27 2011.
STATUS
approved