[go: up one dir, main page]

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!)
A232397 a(n) = ceiling(sqrt(n^4 + n^3 + n^2 + n + 1))^2 - (n^4 + n^3 + n^2 + n + 1). 2
0, 4, 5, 0, 20, 3, 45, 8, 80, 15, 125, 24, 180, 35, 245, 48, 320, 63, 405, 80, 500, 99, 605, 120, 720, 143, 845, 168, 980, 195, 1125, 224, 1280, 255, 1445, 288, 1620, 323, 1805, 360, 2000, 399, 2205, 440, 2420, 483, 2645, 528, 2880, 575, 3125, 624, 3380, 675 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
a(n) = 0 if and only if n^4 + n^3 + n^2 + n + 1 is a perfect square.
Using formula below, we immediately prove that a(n)=0 iff n=0 or n=3. This means that all nonnegative solutions of the Diophantine equation n^4 + n^3 + n^2 + n + 1 = m^2 are n=0, m=1 and n=3, m=11.
For m >=0, if we also consider negative values of n, we obtain only one more solution: n=-1, m=1.
Indeed, if one considers sequence b(n) = ceiling(sqrt(n^4 - n^3 + n^2 - n + 1))^2 - (n^4 - n^3 + n^2 - n +1 ), then, for even n, a(n) = b(n), while for odd n>=3, a(n) = b(n-2).
LINKS
Max Alekseyev, Re: oddness of iterations of sigma(), SeqFan Mailing List.
FORMULA
a(1) = 4, for other odd n, a(n) = ((n-1)/2)^2 - 1; for even n>=0, a(n) = 5/4 * n^2.
a(n) = A068527(A053699(n)). [Straight from the description: Difference between smallest square >= (n^4 + n^3 + n^2 + n + 1) and (n^4 + n^3 + n^2 + n + 1)]. - Antti Karttunen, Nov 28 2013
a(n) = (6*n^2-2*n-3+(4*n^2+2*n+3)*(-1)^n+20*(1-(-1)^(2^abs(n-1))))/8. - Luce ETIENNE, Jan 30 2016
G.f.: 4*x+x^2*(x^5-3*x^3-5*x^2-5)/(x^2-1)^3. - Robert Israel, Feb 02 2016
MAPLE
0, 4, seq(op([5*k^2, k^2-1]), k=1..100); # Robert Israel, Feb 02 2016
MATHEMATICA
Table[Ceiling[Sqrt[n^4 + n^3 + n^2 + n + 1]]^2 - (n^4 + n^3 + n^2 + n + 1), {n, 0, 60}] (* Vincenzo Librandi, Jan 31 2016 *)
PROG
(Magma) [Ceiling(Sqrt(n^4+n^3+n^2+n+1))^2-(n^4+n^3+n^2+n+1): n in [0..60]]; // Vincenzo Librandi, Jan 31 2016
(Python)
from math import isqrt
def A232397(n): return (1+isqrt(m:=n*(n*(n*(n+1)+1)+1)))**2-m-1 # Chai Wah Wu, Jul 29 2022
CROSSREFS
Sequence in context: A200013 A178219 A322232 * A360220 A122753 A016714
KEYWORD
nonn,easy
AUTHOR
Vladimir Shevelev, Nov 23 2013
EXTENSIONS
More terms from Peter J. C. Moses
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 August 29 09:16 EDT 2024. Contains 375511 sequences. (Running on oeis4.)