OFFSET
1,3
COMMENTS
This is "one-fifteenth of triangular numbers (integers only)". - Vladimir Joseph Stephan Orlovsky, Mar 04 2009
The sequence terms are the exponents in the expansion of Product_{n >= 1} (1 - q^n)/( (1 - q^(10*n-2))*(1 - q^(10*n-8)) ) = 1 - q - q^3 + q^7 + q^8 - q^14 - q^20 + + - - ... . - Peter Bala, Dec 26 2024
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..20000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,2,-2,0,0,-1,1).
FORMULA
a(-1 - n) = a(n). a(n) = (A057538(n) * 2 - 1) / 120.
G.f.: -x^2*(1+2*x+4*x^2+x^3+4*x^4+x^6+2*x^5) / ( (1+x)^2*(x^2+1)^2*(x-1)^3 ). - R. J. Mathar, Jun 09 2013
From Peter Bala, Dec 26 2024: (Start)
a(n) is quasi-polynomial in n
a(4*n) = n*(15*n + 1)/2; a(4*n+1) = (3*n + 1)*(5*n + 2)/2;
a(4*n+2) = (3*n + 2)*(5*n + 3)/2; a(4*n+3) = (n + 1)*(15*n + 14)/2.
For 0 <= k <= 3, a(4*n+k) = (N_k(n)^2 - 1)/120, where N_0(n) = 30*n + 1, N_1(n) = 30*n + 11, N_2(n) = 30*n + 19 and N_3(n) = 30*n + 29. (End)
MAPLE
A093722 := proc(q) local n;
for n from 0 to q do
if type(sqrt(120*n+1), integer) then print(n);
fi; od; end:
A093722(1500); # Peter Bala, Dec 26 2024
MATHEMATICA
Select[Table[(n^2-1)/120, {n, 0, 700}], IntegerQ] (* Harvey P. Dale, Nov 26 2010 *)
PROG
(PARI) {a(n) = (((n\4 * 3 + n%4) * 10 + (-1)^(n\2))^2 - 1) / 120 } /* Michael Somos, Oct 17 2006 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael Somos, Apr 13 2004
EXTENSIONS
More terms from Harvey P. Dale, Nov 26 2010
Offset corrected to 1 by Ray Chandler, Jul 29 2019
STATUS
approved