login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A167190 a(n) = 1 + 85*n + 2232*n^2 + 15276*n^3. 3
17594, 131307, 432796, 1013717, 1965726, 3380479, 5349632, 7964841, 11317762, 15500051, 20603364, 26719357, 33939686, 42356007, 52059976, 63143249, 75697482, 89814331, 105585452, 123102501, 142457134, 163741007, 187045776 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

As mentioned in A166957, polynomials in two variables, not necessarily homogeneous, also have a property similar to that in a single variable (cf. A165806, A165808 and A165809) viz f(x+k*f(x,y), y + k*f(x,y)) is congruent to 0 (mod(f(x,y)). The quotient has two parts: a rational integer and a rational integer coefficient of sqrt(-1), when x belongs to Z(x = 5) and y is complex (sqrt(-1)). The polynomial considered is identical with that in A166957 viz x^3 + 2xy + y^2. The present sequence is only that of the rational integers and seq A167191 will consist of rational integer coefficients of sqrt(-1). Note: k belongs to N.

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

FORMULA

G.f.: x*(17594 + 60931*x + 13132*x^2 - x^3)/(1-x)^4 . - R. J. Mathar, Sep 02 2011

a(n) = +4*a(n-1) -6*a(n-2) +4*a(n-3) -a(n-4). - R. J. Mathar, Sep 02 2011

E.g.f.: (1 + 17593*x + 48060*x^2 + 15276*x^3)*exp(x) -1. - G. C. Greubel, Apr 09 2016

EXAMPLE

When x = 5 and y = i, f(x,y) = x^3 + 2xy + y^2 = 124 + 10i. The quotient of f(x + f(x,y), y + f(x,y))/(124 + 10i) is 17594 + 2664i.

MAPLE

seq(1 + 85*n + 2232*n^2 + 15276*n^3, n=1..40); # G. C. Greubel, Sep 01 2019

MATHEMATICA

CoefficientList[Series[(17594+60931*x+13132*x^2-x^3)/(x-1)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 02 2012 *)

Table[1 +85*n +2232*n^2 +15276*n^3, {n, 40}] (* G. C. Greubel, Sep 01 2019 *)

PROG

(MAGMA)  I:=[17594, 131307, 432796, 1013717]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jul 02 2012

(PARI) vector(40, n, 1 +85*n +2232*n^2 +15276*n^3) \\ G. C. Greubel, Sep 01 2019

(Sage) [1 +85*n +2232*n^2 +15276*n^3 for n in (0..40)] # G. C. Greubel, Sep 01 2019

(GAP) List([0..40], n-> 1 +85*n +2232*n^2 +15276*n^3); # G. C. Greubel, Sep 01 2019

CROSSREFS

Cf. A165806, A165808, A165809, A166957.

Sequence in context: A203963 A253403 A213107 * A203353 A273035 A036377

Adjacent sequences:  A167187 A167188 A167189 * A167191 A167192 A167193

KEYWORD

nonn,easy

AUTHOR

A.K. Devaraj, Oct 30 2009

EXTENSIONS

Extended beyond a(6) by R. J. Mathar, Nov 17 2009

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 28 00:32 EST 2020. Contains 331313 sequences. (Running on oeis4.)