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.

Previous Showing 11-19 of 19 results.

A299960 a(n) = (4^(2*n+1) + 1) / 5.

Original entry on oeis.org

1, 13, 205, 3277, 52429, 838861, 13421773, 214748365, 3435973837, 54975581389, 879609302221, 14073748835533, 225179981368525, 3602879701896397, 57646075230342349, 922337203685477581, 14757395258967641293, 236118324143482260685, 3777893186295716170957
Offset: 0

Views

Author

M. F. Hasler, Feb 22 2018

Keywords

Comments

It is easily seen that 4^(2n+1)+1 is divisible by 5 for all n, since 4 = -1 (mod 5). For even powers this does not hold.
The aerated sequence 1, 0, 13, 0, 205, 0, 3277, ... is a linear divisibility sequence of order 4. It is the case P1 = 0, P2 = -5^2, Q = 4 of the 3-parameter family of 4th-order linear divisibility sequences found by Williams and Guy. Cf. A007583, A095372 and A100706. - Peter Bala, Aug 28 2019
Let G be a sequence satisfying G(i) = 2*G(i-1) + G(i-2) - 2*G(i-3) for arbitrary integers i and without regard to the initial values of G. Then a(n) = (G(i)*2^(4*n+2) + G(i+8*n+4))/(5*G(i+4*n+2)) as long as G(i+4*n+2) != 0. - Klaus Purath, Feb 02 2021
Ch. Gerbr asks (personal comm.) whether we can prove that 13 is the only prime in this sequence. We can prove divisibility conditions for many residue classes of the index n (cf. formulas), but have not yet found a complete covering set. - M. F. Hasler, Jan 07 2025

Examples

			For n = 0, a(0) = (4^1+1)/5 = 5/5 = 1.
For n = 1, a(1) = (4^3+1)/5 = 65/5 = 13.
		

Crossrefs

Cf. A299959 for the smallest prime factor.

Programs

Formula

a(n) = A052539(2*n+1)/5 = A015521(2*n+1) = A014985(2*n+1) = A007910(4*n+1) = A007909(4*n+1) = A207262(n+1)/5.
O.g.f.: (1 - 4*x)/(1 - 17*x + 16*x^2). - Peter Bala, Aug 28 2019
a(n) = 17*a(n-1) - 16*a(n-2). - Wesley Ivan Hurt, Oct 02 2020
From Klaus Purath, Feb 02 2021: (Start)
a(n) = (2^(4*n+2)+1)/5.
a(n) = (A061654(n) + A001025(n))/2.
a(n) = A091881(n+1) + 7*A131865(n-1) for n > 0.
(End)
E.g.f.: (exp(x) + 4*exp(16*x))/5. - Stefano Spezia, Feb 02 2021
We have d | a(n) for all n in R, for the following pairs (d, R) of divisors d and residue classes R: (13, 1 + 3Z), (5, 2 + 5Z), (29, 3 + 7Z), (397, 5 + 11Z),
(53, 6 + 13Z), (137, 8 + 17Z), (229, 9 + 19Z), (277, 11 + 23Z),
(107367629, 14 + 29Z), (5581, 15 + 31Z), (149, 18 + 27Z), (10169, 20 + 41Z),
(173, 21 + 43Z), (3761, 23 + 47Z), (15358129, 26 + 53Z), (1181, 29 + 59Z),
(733, 30 + 61Z), (269, 33 + 67Z), (569, 35 + 71Z),(293, 36 + 73Z), (317, 39 + 79Z),
(997, 41 + 83Z), (1069, 44 + 89Z), (389, 48 + 97Z), (809, 50 + 101Z),
(41201, 51 + 103Z), (857, 53 + 107Z), (5669, 54 + 109Z), (58309, 56 + 113Z),
(509, 63 + 127Z), (269665073, 65 + 131Z), (189061, 68 + 137Z), (557, 69 + 139Z),
(1789, 74 + 149Z), (653, 81 + 163Z), (9413, 90 + 181Z), (3821, 95 + 191Z),
(773, 96 + 193Z), (4729, 98 + 197Z), (797, 99 + 199Z), ... - M. F. Hasler, Jan 07 2025

A014992 a(n) = (1 - (-10)^n)/11.

Original entry on oeis.org

1, -9, 91, -909, 9091, -90909, 909091, -9090909, 90909091, -909090909, 9090909091, -90909090909, 909090909091, -9090909090909, 90909090909091, -909090909090909, 9090909090909091, -90909090909090909
Offset: 1

Views

Author

Keywords

Comments

q-integers for q = -10.

Crossrefs

Programs

  • Magma
    I:=[1, -9]; [n le 2 select I[n] else -9*Self(n-1) +10*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 22 2012
    
  • Maple
    a:=n->sum ((-10)^j, j=0..n): seq(a(n), n=0..25); # Zerinvary Lajos, Dec 16 2008
  • Mathematica
    CoefficientList[Series[1/((1 - x)*(1 + 10*x)), {x, 0, 30}], x] (* Vincenzo Librandi, Oct 22 2012 *)
  • PARI
    for(n=1, 30, print1((1-(-10)^n)/11, ", ")) \\ G. C. Greubel, May 26 2018
  • Sage
    [gaussian_binomial(n,1,-10) for n in range(1,19)] # Zerinvary Lajos, May 28 2009
    

Formula

a(n) = a(n-1) + q^(n-1) = (q^n - 1) / (q - 1).
G.f.: x/((1 - x)*(1 + 10*x)). - Vincenzo Librandi, Oct 22 2012
a(n) = -9*a(n-1) + 10*a(n-2). - Vincenzo Librandi, Oct 22 2012
a(n) = (-1)^(n+1)*A015585(n). - R. J. Mathar, Oct 26 2015
E.g.f.: (exp(x) - exp(-10*x))/11. - G. C. Greubel, May 26 2018

Extensions

Better name from Ralf Stephan, Jul 14 2013

A239473 Triangle read by rows: signed version of A059260: coefficients for expansion of partial sums of sequences a(n,x) in terms of their binomial transforms (1+a(.,x))^n ; Laguerre polynomial expansion of the truncated exponential.

Original entry on oeis.org

1, 0, 1, 1, -1, 1, 0, 2, -2, 1, 1, -2, 4, -3, 1, 0, 3, -6, 7, -4, 1, 1, -3, 9, -13, 11, -5, 1, 0, 4, -12, 22, -24, 16, -6, 1, 1, -4, 16, -34, 46, -40, 22, -7, 1, 0, 5, -20, 50, -80, 86, -62, 29, -8, 1, 1, -5, 25, -70, 130, -166, 148, -91, 37, -9, 1, 0, 6, -30, 95, -200, 296, -314, 239, -128, 46, -10, 1
Offset: 0

Views

Author

Tom Copeland, Mar 19 2014

Keywords

Comments

With T the lower triangular array above and the Laguerre polynomials L(k,x) = Sum_{j=0..k} (-1)^j binomial(k, j) x^j/j!, the following identities hold:
(A) Sum_{k=0..n} (-1)^k L(k,x) = Sum_{k=0..n} T(n,k) x^k/k!;
(B) Sum_{k=0..n} x^k/k! = Sum_{k=0..n} T(n,k) L(k,-x);
(C) Sum_{k=0..n} x^k = Sum_{k=0..n} T(n,k) (1+x)^k = (1-x^(n+1))/(1-x).
More generally, for polynomial sequences,
(D) Sum_{k=0..n} P(k,x) = Sum_{k=0..n} T(n,k) (1+P(.,x))^k,
where, e.g., for an Appell sequence, such as the Bernoulli polynomials, umbrally, (1+ Ber(.,x))^k = Ber(k,x+1).
Identity B follows from A through umbral substitution of j!L(j,-x) for x^j in A. Identity C, related to the cyclotomic polynomials for prime index, follows from B through the Laplace transform.
Integrating C gives Sum_{k=0..n} T(n,k) (2^(k+1)-1)/(k+1) = H(n+1), the harmonic numbers.
Identity A >= 0 for x >= 0 (see MathOverflow link for evaluation in terms of Hermite polynomials).
From identity C, W(m,n) = (-1)^n Sum_{k=0..n} T(n,k) (2-m)^k = number of walks of length n+1 between any two distinct vertices of the complete graph K_m for m > 2.
Equals A112468 with the first column of ones removed. - Georg Fischer, Jul 26 2023

Examples

			Triangle begins:
   1
   0    1
   1   -1    1
   0    2   -2    1
   1   -2    4   -3    1
   0    3   -6    7   -4    1
   1   -3    9  -13   11   -5    1
   0    4  -12   22  -24   16   -6    1
   1   -4   16  -34   46  -40   22   -7    1
   0    5  -20   50  -80   86  -62   29   -8    1
   1   -5   25  -70  130 -166  148  -91   37   -9    1
		

Crossrefs

For column 2: A001057, A004526, A008619, A140106.
Column 3: A002620, A087811.
Column 4: A002623, A173196.
Column 5: A001752.
Column 6: A001753.
Cf. Bottomley's cross-references in A059260.
Embedded in alternating antidiagonals of T are the reversals of arrays A071921 (A225010) and A210220.

Programs

  • Magma
    [[(&+[(-1)^(j+k)*Binomial(j,k): j in [0..n]]): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Feb 06 2018
    
  • Maple
    A239473 := proc(n,k)
        add(binomial(j,k)*(-1)^(j+k),j=k..n) ;
    end proc; # R. J. Mathar, Jul 21 2016
  • Mathematica
    Table[Sum[(-1)^(j+k)*Binomial[j,k], {j,0,n}], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 06 2018 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(sum(j=0,n, (-1)^(j+k)*binomial(j, k)), ", "))) \\ G. C. Greubel, Feb 06 2018
    
  • Sage
    Trow = lambda n: sum((x-1)^j for j in (0..n)).list()
    for n in (0..10): print(Trow(n)) # Peter Luschny, Jul 09 2019

Formula

T(n, k) = Sum_{j=0..n} (-1)^(j+k) * binomial(j, k).
E.g.f: (exp(t) - (x-1)*exp((x-1)*t))/(2-x).
O.g.f. (n-th row): (1-(x-1)^(n+1))/(2-x).
Associated operator identities:
With D=d/dx, :xD:^n=x^n*D^n, and :Dx:^n=D^n*x^n, then bin(xD,n)= binomial(xD,n)=:xD:^n/n! and L(n,-:xD:)=:Dx:^n/n!=bin(xD+n,n)=(-1)^n bin(-xD-1,n),
A-o) Sum_{k=0..n} (-1)^k L(k,-:xD:) = Sum_{k=0..n} :-Dx:^k/k!
= Sum_{k=0..n} T(n,k) :-xD:^k/k! = Sum_{k=0..n} (-1)^k T(n,k)bin(xD,k)
B-o) Sum_{k=0..n} :xD:^k/k! = Sum_{k=0..n}, T(n,k) L(k,-:xD:)
= Sum_{k=0..n} T(n,k) :Dx:^k/k! = Sum_{k=0..n}, bin(xD,k).
Associated binomial identities:
A-b) Sum_{k=0..n} (-1)^k bin(s+k,k) = Sum_{k=0..n} (-1)^k T(n,k) bin(s,k)
= Sum_{k=0..n} bin(-s-1,k) = Sum{k=0..n} T(n,k) bin(-s-1+k,k)
B-b) Sum_{k=0..n} bin(s,k) = Sum_{k=0..n} T(n,k) bin(s+k,k)
= Sum_{k=0..n} (-1)^k bin(-s-1+k,k)
= Sum_{k=0..n} (-1)^k T(n,k) bin(-s-1,k).
In particular, from B-b with s=n, Sum_{k=0..n} T(n,k) bin(n+k,k) = 2^n. From B-b with s=0, row sums are all 1.
From identity C with x=-2, the unsigned row sums are the Jacobsthal sequence, i.e., Sum_{k=0..n} T(n,k) (1+(-2))^k = (-1)^n A001045(n+1); for x=2, the Mersenne numbers A000225; for x=-3, A014983 or signed A015518; for x=3, A003462; for x=-4, A014985 or signed A015521; for x=4, A002450; for x=-5, A014986 or signed A015531; and for x=5, A003463; for x=-6, A014987 or signed A015540; and for x=6, A003464.
With -s-1 = m = 0,1,2,..., B-b gives finite differences (recursions):
Sum_{k=0..n} (-1)^k T(n,k) bin(m,k) = Sum_{k=0..n} (-1)^k bin(m+k,k) = T(n+m,m), i.e., finite differences of the columns of T generate shifted columns of T. The columns of T are signed, shifted versions of sequences listed in the cross-references. Since the finite difference is an involution, T(n,k) = Sum_{j=0..k} (-1)^j T(n+j,j) bin(k,j)}. Gauss-Newton interpolation can be applied to give a generalized T(n,s) for s noninteger.
From identity C, S(n,m) = Sum_{k=0..n} T(n,k) bin(k,m) = 1 for m < n+1 and 0 otherwise, i.e., S = T*P, where S = A000012, as a lower triangular matrix and P = Pascal = A007318, so T = S*P^(-1), where P^(-1) = A130595, the signed Pascal array (see A132440), the inverse of P, and T^(-1) = P*S^(-1) = P*A167374 = A156644.
U(n,cos(x)) = e^(-n*i*x)*Sum_{k=0..n} T(n,k)*(1+e^(2*i*x))^k = sin((n+1)x)/sin(x), where U is the Chebyschev polynomial of the second kind A053117 and i^2 = -1. - Tom Copeland, Oct 18 2014
From Tom Copeland, Dec 26 2015: (Start)
With a(n,x) = e^(nx), the partial sums are 1+e^x+...+e^(nx) = Sum_{k=0..n} T(n,k) (1+e^x)^k = [ x / (e^x-1) ] [ e^((n+1)x) -1 ] / x = [ (x / (e^x-1)) e^((n+1)x) - (x / (e^x-1)) ] / x = Sum_{k>=0} [ (Ber(k+1,n+1) - Ber(k+1,0)) / (k+1) ] * x^k/k!, where Ber(n,x) are the Bernoulli polynomials (cf. Adams p. 140). Evaluating (d/dx)^m at x=0 of these expressions gives relations among the partial sums of the m-th powers of the integers, their binomial transforms, and the Bernoulli polynomials.
With a(n,x) = (-1)^n e^(nx), the partial sums are 1-e^x+...+(-1)^n e^(nx) = Sum_{k=0..n} T(n,k) (1-e^x)^k = [ (-1)^n e^((n+1)x) + 1 ] / (e^x+1) = [ (-1)^n (2 / (e^x+1)) e^((n+1)x) + (2 / (e^x+1)) ] / 2 = (1/2) Sum_{k>=0} [ (-1)^n Eul(k,n+1) + Eul(k,0) ] * x^k/k!, where Eul(n,x) are the Euler polynomials. Evaluating (d/dx)^m at x=0 of these expressions gives relations among the partial sums of signed m-th powers of the integers; their binomial transforms, related to the Stirling numbers of the second kind and face numbers of the permutahedra; and the Euler polynomials. (End)
As in A059260, a generator in terms of bivariate polynomials with the coefficients of this entry is given by (1/(1-y))*1/(1 + (y/(1-y))*x - (1/(1-y))*x^2) = 1 + y + (x^2 - x*y + y^2) + (2*x^2*y - 2*x*y^2 + y^3) + (x^4 - 2*x^3*y + 4*x^2*y^2 - 3*x*y^3 + y^4) + ... . This is of the form -h2 * 1 / (1 + h1*x + h2*x^2), related to the bivariate generator of A049310 with h1 = y/(1-y) and h2 = -1/(1-y) = -(1+h1). - Tom Copeland, Feb 16 2016
From Tom Copeland, Sep 05 2016: (Start)
Letting P(k,x) = x in D gives Sum_{k=0..n} T(n,k)*Sum_{j=0..k} binomial(k,j) = Sum_{k=0..n} T(n,k) 2^k = n + 1.
The quantum integers [n+1]q = (q^(n+1) - q^(-n-1)) / (q - q^(-1)) = q^(-n)*(1 - q^(2*(n+1))) / (1 - q^2) = q^(-n)*Sum{k=0..n} q^(2k) = q^(-n)*Sum_{k=0..n} T(n,k)*(1 + q^2)^k. (End)
T(n, k) = [x^k] Sum_{j=0..n} (x-1)^j. - Peter Luschny, Jul 09 2019
a(n) = -n + Sum_{k=0..n} A341091(k). - Thomas Scheuerle, Jun 17 2022

Extensions

Inverse array added by Tom Copeland, Mar 26 2014
Formula re Euler polynomials corrected by Tom Copeland, Mar 08 2024

A014991 a(n) = (1 - (-9)^n)/10.

Original entry on oeis.org

1, -8, 73, -656, 5905, -53144, 478297, -4304672, 38742049, -348678440, 3138105961, -28242953648, 254186582833, -2287679245496, 20589113209465, -185302018885184, 1667718169966657, -15009463529699912
Offset: 1

Views

Author

Keywords

Comments

q-integers for q = -9.

Crossrefs

Programs

  • Magma
    I:=[1,-8]; [n le 2 select I[n] else -8*Self(n-1)+9*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 22 2012
    
  • Maple
    a:=n->sum ((-9)^j, j=0..n): seq(a(n), n=0..25); # Zerinvary Lajos, Dec 16 2008
  • Mathematica
    ((-9)^Range[30]-1)/-10 (* or *) LinearRecurrence[{-8,9},{1,-8},30] (* Harvey P. Dale, Aug 08 2011 *)
    CoefficientList[Series[1/((1 - x)*(1 + 9*x)), {x, 0, 30}], x]; (* Vincenzo Librandi, Oct 22 2012 *)
  • PARI
    for(n=1,30, print1((1-(-9)^n)/10, ", ")) \\ G. C. Greubel, May 26 2018
  • Sage
    [gaussian_binomial(n,1,-9) for n in range(1,19)] # Zerinvary Lajos, May 28 2009
    

Formula

a(n) = a(n-1) + q^(n-1) = (q^n - 1) / (q - 1).
a(0)=1, a(1)=-8, a(n) = -8*a(n-1) + 9*a(n-2). - Harvey P. Dale, Aug 08 2011
G.f.: x/((1 - x)*(1 + 9*x)). - Vincenzo Librandi, Oct 22 2012
E.g.f.: (exp(x) - exp(-9*x))/10. - G. C. Greubel, May 26 2018

Extensions

Better name from Ralf Stephan, Jul 14 2013

A014993 a(n) = (1 - (-11)^n)/12.

Original entry on oeis.org

1, -10, 111, -1220, 13421, -147630, 1623931, -17863240, 196495641, -2161452050, 23775972551, -261535698060, 2876892678661, -31645819465270, 348104014117971, -3829144155297680, 42120585708274481
Offset: 1

Views

Author

Keywords

Comments

q-integers for q = -11.

Crossrefs

Programs

  • Magma
    I:=[1, -10]; [n le 2 select I[n] else -10*Self(n-1) +11*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 22 2012
    
  • Maple
    a:=n->sum ((-11)^j, j=0..n): seq(a(n), n=0..25); # Zerinvary Lajos, Dec 16 2008
  • Mathematica
    LinearRecurrence[{-10, 11}, {1, -10}, 40] (* Vincenzo Librandi, Oct 22 2012 *)
  • PARI
    for(n=1,30, print1((1-(-11)^n)/12, ", ")) \\ G. C. Greubel, May 26 2018
  • Sage
    [gaussian_binomial(n,1,-11) for n in range(1,18)] # Zerinvary Lajos, May 28 2009
    

Formula

a(n) = a(n-1) + q^{(n-1)} = {(q^n - 1) / (q - 1)}.
G.f.: x/((1 - x)*(1 + 11*x)). - Vincenzo Librandi, Oct 22 2012
a(n) = -10*a(n-1) + 11*a(n-2). - Vincenzo Librandi, Oct 22 2012
E.g.f.: (exp(x) - exp(-11*x))/12. - G. C. Greubel, May 26 2018

Extensions

Better name from Ralf Stephan, Jul 14 2013

A015000 q-integers for q=-13.

Original entry on oeis.org

1, -12, 157, -2040, 26521, -344772, 4482037, -58266480, 757464241, -9847035132, 128011456717, -1664148937320, 21633936185161, -281241170407092, 3656135215292197, -47529757798798560, 617886851384381281
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    I:=[1,-12]; [n le 2 select I[n] else -12*Self(n-1)+13*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 22 2012
    
  • Maple
    a:=n->sum ((-13)^j, j=0..n-1): seq(a(n), n=0..20); # Zerinvary Lajos, Dec 16 2008
  • Mathematica
    QBinomial[Range[20],1,-13] (* Harvey P. Dale, May 02 2012 *)
    LinearRecurrence[{-12, 13}, {1, -12}, 30] (* Vincenzo Librandi, Oct 22 2012 *)
  • PARI
    for(n=1, 30, print1((1-(-13)^n)/14, ", ")) \\ G. C. Greubel, May 26 2018
  • Sage
    [gaussian_binomial(n,1,-13) for n in range(1,18)] # Zerinvary Lajos, May 28 2009
    

Formula

a(n) = a(n-1) + q^(n-1) = (q^n - 1) / (q - 1), with q=-13.
a(n) = Sum_{j=0..n-1} (-13)^j. - Zerinvary Lajos, Dec 16 2008
G.f.: x/((1 - x)*(1 + 13*x)). - Vincenzo Librandi, Oct 22 2012
a(n) = -12*a(n-1) + 13*a(n-2). - Vincenzo Librandi, Oct 22 2012
From G. C. Greubel, May 26 2018: (Start)
a(n) = (1 - (-13)^n)/14.
E.g.f.: (exp(x) - exp(-13*x))/14. (End)

Extensions

Edited by N. J. A. Sloane, Jan 13 2009 at the suggestion of R. J. Mathar

A103440 a(n) = Sum[d|n, d==1 (mod 3), d^2] - Sum[d|n, d==2 (mod 3), d^2].

Original entry on oeis.org

1, -3, 1, 13, -24, -3, 50, -51, 1, 72, -120, 13, 170, -150, -24, 205, -288, -3, 362, -312, 50, 360, -528, -51, 601, -510, 1, 650, -840, 72, 962, -819, -120, 864, -1200, 13, 1370, -1086, 170, 1224, -1680, -150, 1850, -1560, -24, 1584, -2208, 205, 2451, -1803, -288, 2210, -2808, -3, 2880, -2550
Offset: 1

Views

Author

Ralf Stephan, Feb 11 2005

Keywords

Examples

			G.f. = q - 3*q^2 + q^3 + 13*q^4 - 24*q^5 - 3*q^6 + 50*q^7 - 51*q^8 + q^9 + ...
		

Crossrefs

Equals A103637(n) - A103638(n). Cf. A002173.
A109041(n) = -9 * a(n) unless n=0. A014985(n) = a(2^n). -24 * A134340(n) = a(6*n+5).

Programs

  • Maple
    f:= proc(n) local D,d;
      D:= numtheory:-divisors(n/3^padic:-ordp(n,3));
      -add((-1)^(d mod 3)*d^2, d = D)
    end proc:
    map(f, [$1..100]); # Robert Israel, Aug 16 2018
  • Mathematica
    a[n_] := Sum[m=Mod[d, 3]; (Boole[m==1]-Boole[m==2]) d^2, {d, Divisors[n]}];
    Array[a, 56] (* Jean-François Alcover, Aug 16 2018 *)
    a[ n_] := SeriesCoefficient[ (1 - QPochhammer[ x]^9 / QPochhammer[ x^3]^3) / 9, {x, 0, n}]; (* Michael Somos, Sep 07 2018 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, d^2 * kronecker( -3, d)))}; /* Michael Somos, Oct 21 2007 */
    
  • PARI
    {a(n) = my(A, p, e, a0, a1, x, y, z); if( n<1, 0, A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 1, z = kronecker( -3, p) * p^2; a0 = 1; a1 = y = z + 1; for(i=2, e, x = y * a1 - z * a0; a0 = a1; a1 = x); a1)))}; /* Michael Somos, Oct 21 2007 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - eta(x + A)^9 / eta(x^3 + A)^3) / 9, n))}; /* Michael Somos, Oct 21 2007 */

Formula

G.f.: F(q) = Sum_{n>=1} A049347(n-1) * n^2 * q^n / (1 - q^n).
G.f.: F(q) = -q * G'(q) / (9*G(q)), with G(q) = Product_{n>=1} (1 - q^n)^(9*n * A049347(n-1)).
a(n) is multiplicative with a(3^e) = 1, a(p^e) = a(p) * a(p^(e-1)) - z * a(p^(e-2)) where z = Kronecker(-3, p) * p^2 and a(p) = z + 1.
a(3*n) = a(n).
G.f.: Sum_{k>0} x^k * (1 - x^k - 6*x^(2*k) - x^(3*k) + x^(4*k)) / (1 + x^k + x^(2*k))^3. - Michael Somos, Oct 21 2007
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2 - v + w + 3*v^2 - 8*w^2 + 6*v*w - 8*u*w + 6*u*v - 9*v^3 - 54*u*v*w + 72*u*w^2 - 9*u^2*w. - Michael Somos, Dec 23 2007

A230368 A strong divisibility sequence associated with the algebraic integer 1 + i.

Original entry on oeis.org

1, 1, 1, 5, 1, 1, 1, 15, 1, 1, 1, 65, 1, 1, 1, 255, 1, 1, 1, 1025, 1, 1, 1, 4095, 1, 1, 1, 16385, 1, 1, 1, 65535, 1, 1, 1, 262145, 1, 1, 1, 1048575, 1, 1, 1, 4194305, 1, 1, 1, 16777215, 1, 1, 1, 67108865, 1, 1, 1, 268435455, 1, 1, 1, 1073741825
Offset: 1

Views

Author

Peter Bala, Jan 10 2014

Keywords

Comments

Let alpha be an algebraic integer and define a sequence of integers a(n) by the condition a(n) = max { integer d : alpha^n == 1 (mod d)}. Silverman shows that a(n) is a strong divisibility sequence, that is gcd(a(n), a(m)) = a(gcd(n, m)) for all n and m in N; in particular, if n divides m then a(n) divides a(m). For the present sequence we take alpha = 1 + i. For other examples see A230369, A235450 and (conjecturally) A082630.

Crossrefs

Programs

  • Maple
    seq( gcd( 1/2*((1 - I)^n + (1 + I)^n - 2), I/2*((1 + I)^n - (1 - I )^n ) ), n = 1..80);

Formula

a(4*n) = |(-4)^n - 1| otherwise a(n) = 1.
a(4*n) = 5*A015521(n).
O.g.f.: 1/(1 - 4*x^4) - 1/(1 + x^4) + 1/(1 - x) - 1/(1 - x^4) = x*(-1 -x -x^2 -5*x^3 +3*x^4 +3*x^5 +3*x^6 +5*x^7 +4*x^8 +4*x^9 +4*x^10) / ( (1-x) *(1+x) *(2*x^2+1) *(2*x^2-1) *(x^2+1) *(x^4+1) ).
Recurrence equation: a(n) = 4*a(n-4) + a(n-8) - 4*a(n-12).

A268413 a(n) = Sum_{k = 0..n} (-1)^k*14^k.

Original entry on oeis.org

1, -13, 183, -2561, 35855, -501969, 7027567, -98385937, 1377403119, -19283643665, 269971011311, -3779594158353, 52914318216943, -740800455037201, 10371206370520815, -145196889187291409, 2032756448622079727, -28458590280709116177, 398420263929927626479
Offset: 0

Views

Author

Ilya Gutkovskiy, Feb 04 2016

Keywords

Comments

Alternating sum of powers of 14.
More generally, the ordinary generating function for the Sum_{k = 0..n} (-1)^k*m^k is 1/(1 + (m - 1)*x - m*x^2). Also, Sum_{k = 0..n} (-1)^k*m^k = ((-1)^n*m^(n + 1) + 1)/(m + 1).

Crossrefs

Cf. similar sequences of the type Sum_{k=0..n} (-1)^k*m^k: A059841 (m=1), A077925 (m=2), A014983 (m=3), A014985 (m=4), A014986 (m=5), A014987 (m=6), A014989 (m=7), A014990 (m=8), A014991 (m=9), A014992 (m=10), A014993 (m=11), A014994 (m=12), A015000 (m=13), this sequence (m=14), A239284 (m=15).

Programs

  • Magma
    I:=[1,-19]; [n le 2 select I[n] else -13*Self(n-1) +14*Self(n-2): n in [1..30]]; // G. C. Greubel, May 26 2018
  • Mathematica
    Table[((-1)^n 14^(n + 1) + 1)/15, {n, 0, 18}]
    LinearRecurrence[{-13, 14}, {1, -13}, 19]
    Table[Sum[(-1)^k*14^k, {k, 0, n}], {n, 0, 18}]
  • PARI
    x='x+O('x^30); Vec(1/(1 + 13*x - 14*x^2)) \\ G. C. Greubel, May 26 2018
    

Formula

G.f.: 1/(1 + 13*x - 14*x^2).
a(n) = ((-1)^n*14^(n + 1) + 1)/15.
a(n) = 1 - 14*a(n - 1) for n>0 and a(0)=1.
a(n) = Sum_{k = 0..n} A033999(k)*A001023(k).
Lim_{n -> infinity} a(n)/a(n + 1) = - 1/14.
Previous Showing 11-19 of 19 results.