cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A093519 Numbers with no representation as the sum of two (not necessarily distinct) generalized pentagonal numbers.

Original entry on oeis.org

11, 18, 21, 25, 32, 39, 43, 46, 49, 54, 60, 65, 67, 68, 74, 76, 81, 87, 88, 90, 95, 98, 106, 109, 111, 113, 116, 120, 123, 125, 130, 136, 137, 142, 144, 153, 158, 159, 163, 164, 165, 172, 173, 175, 179, 182, 186, 193, 197, 201, 204, 205, 207, 208, 214, 219, 220
Offset: 1

Views

Author

Jon Perry, Mar 29 2004

Keywords

Links

Crossrefs

Cf. A001318 (generalized pentagonal numbers), A093518.
Cf. A204382.

Programs

  • Maple
    GP:= [0, seq(op([m*(3*m-1)/2, m*(3*m+1)/2]), m=1..50)]:
N:= GP[-1]:
V:= Array(0..N, datatype=integer[4]):
for i from 1 to nops(GP) do
for j from 1 to i do
   r:= GP[i]+GP[j];
   if r > N then break fi;
   V[r]:= V[r]+1
od od:
select(t -> V[t] = 0, [$0..N]); # Robert Israel, Feb 26 2025

Extensions

Definition clarified by Robert Israel, Feb 26 2025

A204384 G.f.: Product_{n>=1} (1 - A002203(n)*x^n + (-x^2)^n) / (1 + A002203(n)*x^n + (-x^2)^n) where A002203(n) is the companion Pell numbers.

Original entry on oeis.org

1, -4, -4, 0, 68, 56, 0, 0, 4, -5572, -4616, 0, 0, -328, 0, 0, 2663428, 2206456, -4, 0, 156808, 0, 0, 0, 0, -7420309452, -6147187208, 0, 0, -436867144, 0, 0, 4, 0, -5326856, 0, 120491016385604, 99818026262072, 0, 0, 7093848711176, -11144, 0, 0, 0, 86497488056, 0, 0, 0
Offset: 0

Views

Author

Paul D. Hanna, Jan 14 2012

Keywords

Comments

a(A022544(n)) = 0 where A022544 lists numbers that are not the sum of 2 squares.
Compare to: Product_{n>=1} (1-q^k)/(1+q^k) = 1 + 2*Sum_{n>=1} (-1)^n*q^(n^2), the Jacobi theta_4 function, which has the g.f: exp( Sum_{n>=1} -(sigma(2*k)-sigma(k)) * x^n/n ).

Examples

			G.f.: A(x) = 1 - 4*x - 4*x^2 + 68*x^4 + 56*x^5 + 4*x^8 - 5572*x^9 - 4616*x^10 +...
-log(A(x)) = 2*2*x + 4*6*x^2/2 + 8*14*x^3/3 + 8*34*x^4/4 + 12*82*x^5/5 + 16*198*x^6/6 +...+ (sigma(2*n)-sigma(n))*A002203(n)*x^n/n +...
Compare to the logarithm of Jacobi theta4 H(x) = 1 + 2*Sum_{n>=1} (-1)^n*q^(n^2):
-log(H(x)) = 2*x + 4*x^2/2 + 8*x^3/3 + 8*x^4/4 + 12*x^5/5 + 16*x^6/6 + 16*x^7/7 +...+ (sigma(2*n)-sigma(n))*x^n/n +...
The g.f. equals the products:
A(x) = (1-2*x-x^2)/(1+2*x-x^2) * (1-6*x^2+x^4)/(1+6*x^2+x^4) * (1-14*x^3-x^6)/(1+14*x^3-x^6) * (1-34*x^4+x^8)/(1+34*x^4+x^8) * (1-82*x^5-x^10)/(1+82*x^5-x^10) *...* (1 - A002203(n)*x^n + (-x^2)^n)/(1 + A002203(n)*x^n + (-x^2)^n) *...
A(x) = (1-2*x-x^2)^2 * (1-6*x^2+x^4) * (1-14*x^3-x^6)^2 * (1-34*x^4+x^8) * (1-82*x^5-x^10)^2 *(1-198*x^6+x^12) * (1-478*x^7-x^14)^2 * (1-1154*x^8+x^16) *...
Positions of zeros form A022544:
[3,6,7,11,12,14,15,19,21,22,23,24,27,28,30,31,33,35,38,39,42,43,44,...]
which are numbers that are not the sum of 2 squares.

Links

Crossrefs

Cf. A203850, A204382, A204383, A022544.

Programs

  • PARI
    /* Subroutine used in PARI programs below: */
{A002203(n)=polcoeff(2*(1-x)/(1-2*x-x^2+x*O(x^n)), n)}
  • PARI
    {a(n)=polcoeff(prod(m=1, n, 1 - A002203(m)*x^m + (-1)^m*x^(2*m) +x*O(x^n))/prod(m=1, n, 1 + A002203(m)*x^m + (-1)^m*x^(2*m) +x*O(x^n)), n)}
  • PARI
    {a(n)=polcoeff(prod(m=1, n\2+1, (1 - A002203(2*m-1)*x^(2*m-1) - x^(4*m-2))^2*(1 - A002203(2*m)*x^(2*m) + x^(4*m) +x*O(x^n))), n)}
  • PARI
    {a(n)=polcoeff(exp(sum(k=1, n,-(sigma(2*k)-sigma(k))*A002203(k)*x^k/k)+x*O(x^n)), n)}

Formula

G.f.: Product_{n>=1} (1 - A002203(2*n-1)*x^(2*n-1) - x^(4*n-2))^2 * (1 - A002203(2*n)*x^(2*n) + x^(4*n)).
G.f.: exp( Sum_{n>=1} -(sigma(2*n)-sigma(n)) * A002203(n) * x^n/n ) where A002203(n) is the companion Pell numbers.

A204383 G.f.: Product_{n>=1} (1 - A002203(n)*x^n + (-1)^n*x^(2*n))^3 where A002203(n) is the companion Pell numbers.

Original entry on oeis.org

1, -6, -9, 70, 90, 0, -1411, -1722, 0, 490, 60534, 75222, 49, -21510, 0, -6067754, -7542180, 0, 2156110, 0, 81, 1420032740, 1764323886, 0, -504516870, -8118, 0, -50196874, -783087782910, -973096740630, -121, 278263575996, 0, 0, 27685627830, 0, 1024173639305948
Offset: 0

Views

Author

Paul D. Hanna, Jan 14 2012

Keywords

Comments

a(A020757(n)) = 0 where A020757 lists numbers that are not the sum of two triangular numbers.

Examples

			G.f.: A(x) = 1 - 6*x - 9*x^2 + 70*x^3 + 90*x^4 - 1411*x^6 - 1722*x^7 +...
-log(A(x))/3 = 1*2*x + 3*6*x^2/2 + 4*14*x^3/3 + 7*34*x^4/4 + 6*82*x^5/5 + 12*198*x^6/6 +...+ sigma(n)*A002203(n)*x^n/n +...
The g.f. equals the product:
A(x) = (1-2*x-x^2)^3 * (1-6*x^2+x^4)^3 * (1-14*x^3-x^6)^3 * (1-34*x^4+x^8)^3 * (1-82*x^5-x^10)^3 * (1-198*x^6+x^12)^3 *...* (1 - A002203(n)*x^n + (-1)^n*x^(2*n))^3 *...
Positions of zeros form A020757:
[5,8,14,17,19,23,26,32,33,35,40,41,44,47,50,52,53,54,59,62,63,...].

Links

Crossrefs

Cf. A203861, A204382, A204384, A020757.

Programs

  • PARI
    /* Subroutine used in PARI programs below: */
{A002203(n)=polcoeff(2*(1-x)/(1-2*x-x^2+x*O(x^n)), n)}
  • PARI
    {a(n)=polcoeff(exp(sum(k=1, n, -3*sigma(k)*A002203(k)*x^k/k)+x*O(x^n)), n)}
  • PARI
    {a(n)=polcoeff(prod(m=1, n, 1 - A002203(m)*x^m + (-1)^m*x^(2*m) +x*O(x^n))^3, n)}

Formula

G.f.: exp( Sum_{n>=1} -3 * sigma(n) * A002203(n) * x^n/n ).
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.