login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A226805 P_n(n+1) where P_n(x) is the polynomial of degree n-1 which satisfies P_n(i) = i^i for i = 1,...,n. 1
1, 7, 70, 877, 13316, 237799, 4885980, 113566121, 2946476764, 84417530491, 2647176188372, 90183424037293, 3316840864313484, 130985236211745959, 5528094465439087876, 248308899812296990033, 11827417687501017074876, 595470029978391175571923 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
LINKS
EXAMPLE
P_3(x) = 18 - 27*x + 10*x^2; a(3) = P_3(3+1) = 70.
MATHEMATICA
P[n_][x_] = Sum[a[i]*x^i, {i, 0, n - 1}]; ecu[n_] := Table[P[n][i] == i^i, {i, 1, n}]; PP[n_][x_] := P[n][x] /. Solve[ecu[n]][[1]]; Table[PP[i][i + 1], {i, 1, 22}]
a[n_] := InterpolatingPolynomial[Table[{i, i^i}, {i, n}], n+1]; Array[a, 20] (* Giovanni Resta, Jun 18 2013 *)
PROG
(PARI) a(n)=subst(polinterpolate(vector(n, i, i^i)), 'x, n+1) \\ Charles R Greathouse IV, Nov 19 2013
CROSSREFS
Sequence in context: A209327 A255519 A305966 * A113343 A124566 A141151
KEYWORD
nonn
AUTHOR
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:58 EDT 2024. Contains 371235 sequences. (Running on oeis4.)