A071568 -id:A071568 - cp's OEIS frontend

A163422 Primes p such that A071568((p-1)/2) is also prime.

3, 5, 7, 11, 13, 17, 19, 31, 37, 43, 59, 61, 79, 83, 89, 97, 107, 109, 113, 139, 149, 167, 191, 233, 241, 263, 271, 293, 307, 311, 337, 359, 373, 383, 439, 443, 479, 487, 491, 523, 557, 617, 641, 647, 659, 673, 683, 701, 733, 757, 811, 829, 853, 857, 859, 877

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 27 2009

Primes p such that (p-1)^3/8+(p+1)/2 is also prime, i.e., in A095692.

			p=3 is in the sequence because (3-1)^3/8+(3+1)/2=3 is prime.
p=5 is in the sequence because (5-1)^3/8+(5+1)/2=11 is prime.

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

Cf. A162652, A163418, A163419, A163420, A163421.

[p: p in PrimesUpTo(1000) | IsPrime((p^3-3*p^2+7*p+3) div 8)]; // Vincenzo Librandi, Apr 10 2013

f[n_]:=((n-1)/2)^3+((n+1)/2); lst={}; Do[p=Prime[n]; If[PrimeQ[f[p]], AppendTo[lst, p]], {n,6!}]; lst
Select[Prime[Range[180]], PrimeQ[(#-1)^3/8+(#+1)/2]&]  (* Harvey P. Dale, Jan 05 2011 *)

Definition rewritten by R. J. Mathar, Aug 17 2009

A101220 a(n) = Sum_{k=0..n} Fibonacci(n-k)*n^k.

Original entry on oeis.org

0, 1, 3, 14, 91, 820, 9650, 140601, 2440317, 49109632, 1123595495, 28792920872, 816742025772, 25402428294801, 859492240650847, 31427791175659690, 1234928473553777403, 51893300561135516404, 2322083099525697299278

Offset: 0

Ross La Haye, Dec 14 2004

In what follows a(i,j,k) denotes a three-dimensional array, the terms a(n) are defined as a(n,n,n) in that array. - Joerg Arndt, Jan 03 2021
Previous name was: Three-dimensional array: a(i,j,k) = expansion of x*(1 + (i-j)*x)/((1-j*x)*(1-x-x^2)), read by a(n,n,n).
a(i,j,k) = the k-th value of the convolution of the Fibonacci numbers (A000045) with the powers of i = Sum_{m=0..k} a(i-1,j,m), both for i = j and i > 0; a(i,j,k) = a(i-1,j,k) + a(j,j,k-1), for i,k > 0; a(i,1,k) = Sum_{m=0..k} a(i-1,0,m), for i > 0. With F = Fibonacci and L = Lucas, then a(1,1,k) = F(k+2) - 1; a(2,1,k) = F(k+3) - 2; a(3,1,k) = L(k+2) - 3; a(4,1,k) = 4*F(k+1) + F(k) - 4; a(1,2,k) = 2^k - F(k+1); a(2,2,k) = 2^(k+1) - F(k+3); a(3,2,k) = 3(2^k - F(k+2)) + F(k); a(4,2,k) = 2^(k+2) - F(k+4) - F(k+2); a(1,3,k) = (3^k + L(k-1))/5, for k > 0; a(2,3,k) = (2 * 3^k - L(k)) /5, for k > 0; a(3,3,k) = (3^(k+1) - L(k+2))/5; a(4,3,k) = (4 * 3^k - L(k+2) - L(k+1))/5, etc..

			a(1,3,3) = 6 because a(1,3,0) = 0, a(1,3,1) = 1, a(1,3,2) = 2 and 4*2 - 2*1 - 3*0 = 6.

G. C. Greubel, Table of n, a(n) for n = 0..385
Eric Weisstein's World of Mathematics, Fibonacci Number
Eric Weisstein's World of Mathematics, Lucas Number

a(0, j, k) = A000045(k).
a(1, 2, k+1) - a(1, 2, k) = A099036(k).
a(3, 2, k+1) - a(3, 2, k) = A104004(k).
a(4, 2, k+1) - a(4, 2, k) = A027973(k).
a(1, 3, k+1) - a(1, 3, k) = A099159(k).
a(i, 0, k) = A109754(i, k).
a(i, i+1, 3) = A002522(i+1).
a(i, i+1, 4) = A071568(i+1).
a(2^i-2, 0, k+1) = A118654(i, k), for i > 0.
Sequences of the form a(n, 0, k): A000045(k+1) (n=1), A000032(k) (n=2), A000285(k-1) (n=3), A022095(k-1) (n=4), A022096(k-1) (n=5), A022097(k-1) (n=6), A022098(k-1) (n=7), A022099(k-1) (n=8), A022100(k-1) (n=9), A022101(k-1) (n=10), A022102(k-1) (n=11), A022103(k-1) (n=12), A022104(k-1) (n=13), A022105(k-1) (n=14), A022106(k-1) (n=15), A022107(k-1) (n=16), A022108(k-1) (n=17), A022109(k-1) (n=18), A022110(k-1) (n=19), A088209(k-2) (n=k-2), A007502(k) (n=k-1), A094588(k) (n=k).
Sequences of the form a(1, n, k): A000071(k+2) (n=1), A027934(k-1) (n=2), A098703(k) (n=3).
Sequences of the form a(2, n, k): A001911(k) (n=1), A008466(k+1) (n=2), A106517(k-1) (n=3).
Sequences of the form a(3, n, k): A027961(k) (n=1), A094688(k) (n=3).
Sequences of the form a(4, n, k): A053311(k-1) (n=1), A027974(k-1) (n=2).
Cf. A001622, A001891, A001924, A023537, A104161, A106517.

A101220:= func< n | (&+[n^k*Fibonacci(n-k): k in [0..n]]) >;
[A101220(n): n in [0..30]]; // G. C. Greubel, Jun 01 2025

Join[{0}, Table[Sum[Fibonacci[n-k]*n^k, {k, 0, n}], {n, 1, 20}]] (* Vaclav Kotesovec, Jan 03 2021 *)

a(n)=sum(k=0,n,fibonacci(n-k)*n^k) \\ Joerg Arndt, Jan 03 2021

def A101220(n): return sum(n^k*fibonacci(n-k) for k in range(n+1))
print([A101220(n) for n in range(31)]) # G. C. Greubel, Jun 01 2025

a(i, j, 0) = 0, a(i, j, 1) = 1, a(i, j, 2) = i+1; a(i, j, k) = ((j+1)*a(i, j, k-1)) - ((j-1)*a(i, j, k-2)) - (j*a(i, j, k-3)), for k > 2.
a(i, j, k) = Fibonacci(k) + i*a(j, j, k-1), for i, k > 0.
a(i, j, k) = (Phi^k - (-Phi)^-k + i(((j^k - Phi^k) / (j - Phi)) - ((j^k - (-Phi)^-k) / (j - (-Phi)^-1)))) / sqrt(5), where Phi denotes the golden mean/ratio (A001622).
i^k = a(i-1, i, k) + a(i-2, i, k+1).
A104161(k) = Sum_{m=0..k} a(k-m, 0, m).
a(i, j, 0) = 0, a(i, j, 1) = 1, a(i, j, 2) = i+1, a(i, j, 3) = i*(j+1) + 2; a(i, j, k) = (j+2)*a(i, j, k-1) - 2*j*a(i, j, k-2) - a(i, j, k-3) + j*a(i, j, k-4), for k > 3. a(i, j, 0) = 0, a(i, j, 1) = 1; a(i, j, k) = a(i, j, k-1) + a(i, j, k-2) + i * j^(k-2), for k > 1.
G.f.: x*(1 + (i-j)*x)/((1-j*x)*(1-x-x^2)).
a(n, n, n) = Sum_{k=0..n} Fibonacci(n-k) * n^k. - Ross La Haye, Jan 14 2006
Sum_{m=0..k} binomial(k,m)*(i-1)^m = a(i-1,i,k) + a(i-2,i,k+1), for i > 1. - Ross La Haye, May 29 2006
From Ross La Haye, Jun 03 2006: (Start)
a(3, 3, k+1) - a(3, 3, k) = A106517(k).
a(1, 1, k) = A001924(k) - A001924(k-1), for k > 0.
a(2, 1, k) = A001891(k) - A001891(k-1), for k > 0.
a(3, 1, k) = A023537(k) - A023537(k-1), for k > 0.
Sum_{j=0..i+1} a(i-j+1, 0, j) - Sum_{j=0..i} a(i-j, 0, j) = A001595(i). (End)
a(i,j,k) = a(j,j,k) + (i-j)*a(j,j,k-1), for k > 0.
a(n) ~ n^(n-1). - Vaclav Kotesovec, Jan 03 2021

New name from Joerg Arndt, Jan 03 2021

A104257 Square array T(a,n) read by antidiagonals: replace 2^i with a^i in binary representation of n, where a,n >= 2.

Original entry on oeis.org

2, 3, 3, 4, 4, 4, 5, 5, 9, 5, 6, 6, 16, 10, 6, 7, 7, 25, 17, 12, 7, 8, 8, 36, 26, 20, 13, 8, 9, 9, 49, 37, 30, 21, 27, 9, 10, 10, 64, 50, 42, 31, 64, 28, 10, 11, 11, 81, 65, 56, 43, 125, 65, 30, 11, 12, 12, 100, 82, 72, 57, 216, 126, 68, 31, 12, 13, 13, 121, 101, 90, 73, 343

Offset: 2

Ralf Stephan, Mar 05 2005

Sums of distinct powers of a. Numbers having only {0,1} in a-ary representation.

			Array begins:
  2,  3,  4,  5,  6,  7,   8,   9, ...
  3,  4,  9, 10, 12, 13,  27,  28, ...
  4,  5, 16, 17, 20, 21,  64,  65, ...
  5,  6, 25, 26, 30, 31, 125, 126, ...
  6,  7, 36, 37, 42, 43, 216, 217, ...
  7,  8, 49, 50, 56, 57, 343, 344, ...
  8,  9, 64, 65, 72, 73, 512, 513, ...
  9, 10, 81, 82, 90, 91, 729, 730, ...
  ...

John Tyler Rascoe, Antidiagonals n = 2..140, flattened

Rows include (essentially) A005836, A000695, A033042, A033043, A033044, A033045, A033046, A033047, A033048, A033049, A033050, A033051, A033052.
Columns include A000290, A002522, A002378, A000578, A001093, A034262, A071568, A011379, A098547, A027444.
Main diagonal is A104258.

T[, 0] = 0; T[2, n] := n; T[a_, 2] := a;
T[a_, n_] := T[a, n] = If[EvenQ[n], a T[a, n/2], a T[a, (n-1)/2]+1];
Table[T[a-n+2, n], {a, 2, 13}, {n, 2, a}] // Flatten (* Jean-François Alcover, Feb 09 2021 *)

T(a, n) = fromdigits(binary(n), a); \\ Michel Marcus, Aug 19 2022

def T(a, n): return n if n < 2 else (max(a, n) if min(a, n) == 2 else a*T(a, n//2) + n%2)
print([T(a-n+2, n) for a in range(2, 14) for n in range(2, a+1)]) # Michael S. Branicky, Aug 02 2022

T(a, n) = (1/(a-1))*Sum_{j>=1} floor((n+2^(j-1))/2^j) * ((a-2)*a^(j-1) + 1).
T(a, n) = (1/(a-1))*Sum_{j=1..n} ((a-2)*a^A007814(j) + 1).
G.f. of a-th row: (1/(1-x)) * Sum_{k>=0} a^k*x^2^k/(1+x^2^k).
Recurrence: T(a, 2n) = a*T(a, n), T(a, 2n+1) = a*T(a, n) + 1, T(a, 0) = 0.

A098547 a(n) = n^3 + n^2 + 1.

Original entry on oeis.org

1, 3, 13, 37, 81, 151, 253, 393, 577, 811, 1101, 1453, 1873, 2367, 2941, 3601, 4353, 5203, 6157, 7221, 8401, 9703, 11133, 12697, 14401, 16251, 18253, 20413, 22737, 25231, 27901, 30753, 33793, 37027, 40461, 44101, 47953, 52023, 56317, 60841, 65601, 70603, 75853

Offset: 0

Douglas Winston (douglas.winston(AT)srupc.com), Oct 26 2004

G. C. Greubel, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A000578, A001093, A011379, A027444, A033431, A033562, A034262, A053698, A061317, A066023, A071568.

Cf. A000217, A081423.

[(n^3+n^2+1): n in [1..60]]; // Vincenzo Librandi, Apr 06 2011

with(combinat): seq(fibonacci(3,n)+n^3, n=0..40); # Zerinvary Lajos, May 25 2008

Table[n^3+n^2+1,{n,0,100}] (* Vladimir Joseph Stephan Orlovsky, Mar 09 2011 *)

Vec(-(x^3-7*x^2+x-1)/(x-1)^4 + O(x^100)) \\ Colin Barker, Aug 29 2014

From Colin Barker, Aug 29 2014: (Start)
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
G.f.: (1 - x + 7*x^2 - x^3)/(1-x)^4. (End)
a(n) = A081423(n) + A000217(n-1). - Bruce J. Nicholson, Jan 06 2019
E.g.f.: exp(x)*(1 + 2*x + 4*x^2 + x^3). - Elmo R. Oliveira, Apr 20 2025

A099721 a(n) = n^2(2n+1).

Original entry on oeis.org

0, 3, 20, 63, 144, 275, 468, 735, 1088, 1539, 2100, 2783, 3600, 4563, 5684, 6975, 8448, 10115, 11988, 14079, 16400, 18963, 21780, 24863, 28224, 31875, 35828, 40095, 44688, 49619, 54900, 60543, 66560, 72963, 79764, 86975, 94608, 102675, 111188, 120159, 129600

Offset: 0

Douglas Winston (douglas.winston(AT)srupc.com), Nov 07 2004

For a right triangle with sides of lengths 8*n^3 + 12*n^2 + 8*n + 2, 4*n^4 + 8*n^3 + 4*n^2, and 4*n^4 + 8*n^3 + 12*n^2 + 8*n + 2, dividing the area by the perimeter gives a(n). - J. M. Bergot, Jul 30 2013
This sequence is the difference between the centered icosahedral (or cuboctahedral) numbers (A005902(n)) and the centered octagonal pyramidal numbers (A000447(n+1)). - Peter M. Chema, Jan 09 2016
a(n) is the sum of the integers in the closed interval (n-1)*n to n*(n+1). - J. M. Bergot, Apr 19 2017

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Row n=3 of A229079.
Cf. A000578, A001093, A011379, A015237, A027444, A033431, A033562, A034262, A053698, A061317, A066023, A071568, A098547.
Cf. A000290, A000447, A002378, A005408, A005902, A024196.

[n^2*(2*n+1): n in [0..50]]; // Vincenzo Librandi, May 01 2011

A099721 := proc(n) n^2*(2*n+1) ; end proc:
seq(A099721(n),n=0..10) ;

a[n_]:=2*n^3+n^2; (* Vladimir Joseph Stephan Orlovsky, Dec 21 2008 *)
LinearRecurrence[{4,-6,4,-1},{0,3,20,63},40] (* Harvey P. Dale, Aug 19 2022 *)

a(n) = ceil(sum(i=n^2-(n-1), n^2+(n-1), if(!issquare(4*i+1), (2*i+1+sqrt(4*i+1))/2, 0))); \\ Michel Marcus, Nov 14 2014, after Richard R. Forberg

G.f.: x*(3 + 8*x + x^2)/(x-1)^4.
a(n) = A024196(n) - A024196(n-1). - Philippe Deléham, May 07 2012
a(n) = ceiling(Sum_{i=n^2-(n-1)..n^2+(n-1)} s(i)), for n > 0 and integer i, where s(i) are the real solutions to x = i + sqrt(x), and the summation range excludes the integer solutions which occur where i is an oblong number (A002378). The fractional portion of the summation converges to 2/3 for large n. If s(i) is replaced with i, then the summation equals n^2*(2*n-1) = A015237. - Richard R. Forberg, Oct 15 2014
a(n) = A005902(n) - A000447(n+1). - Peter M. Chema, Jan 09 2016
From Amiram Eldar, May 17 2022: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/6 + 4*log(2) - 4.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/12 - Pi - 2*log(2) + 4. (End)
From Elmo R. Oliveira, Aug 08 2025: (Start)
E.g.f.: x*(1 + 2*x)*(3 + x)*exp(x).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) = A000290(n)*A005408(n). (End)

A292741 Number A(n,k) of partitions of n with k sorts of part 1; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 2, 2, 1, 1, 3, 5, 3, 2, 1, 4, 10, 11, 5, 2, 1, 5, 17, 31, 24, 7, 4, 1, 6, 26, 69, 95, 50, 11, 4, 1, 7, 37, 131, 278, 287, 104, 15, 7, 1, 8, 50, 223, 657, 1114, 865, 212, 22, 8, 1, 9, 65, 351, 1340, 3287, 4460, 2599, 431, 30, 12, 1, 10, 82, 521, 2459, 8042, 16439, 17844, 7804, 870, 42, 14

Offset: 0

Alois P. Heinz, Sep 22 2017

			A(1,3) = 3: 1a, 1b, 1c.
A(2,3) = 10: 2, 1a1a, 1a1b, 1a1c, 1b1a, 1b1b, 1b1c, 1c1a, 1c1b, 1c1c.
A(3,2) = 11: 3, 21a, 21b, 1a1a1a, 1a1a1b, 1a1b1a, 1a1b1b, 1b1a1a, 1b1a1b, 1b1b1a, 1b1b1b.
Square array A(n,k) begins:
  1,  1,   1,    1,     1,     1,      1,      1, ...
  0,  1,   2,    3,     4,     5,      6,      7, ...
  1,  2,   5,   10,    17,    26,     37,     50, ...
  1,  3,  11,   31,    69,   131,    223,    351, ...
  2,  5,  24,   95,   278,   657,   1340,   2459, ...
  2,  7,  50,  287,  1114,  3287,   8042,  17215, ...
  4, 11, 104,  865,  4460, 16439,  48256, 120509, ...
  4, 15, 212, 2599, 17844, 82199, 289540, 843567, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0-2 give: A002865, A000041, A090764.
Rows n=0-2 give: A000012, A001477, A002522, A071568.
Main diagonal gives A292462.
Cf. A003992, A004248, A009998, A051129, A292508, A292622, A292745.

b:= proc(n, i, k) option remember; `if`(n=0 or i<2, k^n,
      add(b(n-i*j, i-1, k), j=0..iquo(n, i)))
    end:
A:= (n, k)-> b(n$2, k):
seq(seq(A(n, d-n), n=0..d), d=0..14);

b[0, , ] = 1; b[n_, i_, k_] := b[n, i, k] = If[i < 2, k^n, Sum[b[n - i*j, i - 1, k], {j, 0, Quotient[n, i]}]];
A[n_, k_] := b[n, n, k];
Table[A[n, d - n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, May 19 2018, translated from Maple *)

G.f. of column k: 1/(1-k*x) * 1/Product_{j>=2} (1-x^j).

A(n,k) = Sum_{j=0..n} A002865(j) * k^(n-j).

A341907 T(n, k) is the result of replacing 2^e with k^e in the binary expansion of n; square array T(n, k) read by antidiagonals upwards, n, k >= 0.

Original entry on oeis.org

0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 2, 2, 1, 0, 1, 1, 3, 3, 1, 0, 0, 2, 4, 4, 4, 1, 0, 1, 2, 5, 9, 5, 5, 1, 0, 0, 3, 6, 10, 16, 6, 6, 1, 0, 1, 1, 7, 12, 17, 25, 7, 7, 1, 0, 0, 2, 8, 13, 20, 26, 36, 8, 8, 1, 0, 1, 2, 9, 27, 21, 30, 37, 49, 9, 9, 1, 0, 0, 3, 10, 28, 64, 31, 42, 50, 64, 10, 10, 1, 0

Offset: 0

Rémy Sigrist, Jun 04 2021

For any n >= 0, the n-th row, k -> T(n, k), corresponds to a polynomial in k with coefficients in {0, 1}.

For any k > 1, the k-th column, n -> T(n, k), corresponds to sums of distinct powers of k.

			Array T(n, k) begins:
  n\k|  0  1   2   3   4    5    6    7    8    9    10    11    12
  ---+-------------------------------------------------------------
    0|  0  0   0   0   0    0    0    0    0    0     0     0     0
    1|  1  1   1   1   1    1    1    1    1    1     1     1     1
    2|  0  1   2   3   4    5    6    7    8    9    10    11    12
    3|  1  2   3   4   5    6    7    8    9   10    11    12    13
    4|  0  1   4   9  16   25   36   49   64   81   100   121   144
    5|  1  2   5  10  17   26   37   50   65   82   101   122   145
    6|  0  2   6  12  20   30   42   56   72   90   110   132   156
    7|  1  3   7  13  21   31   43   57   73   91   111   133   157
    8|  0  1   8  27  64  125  216  343  512  729  1000  1331  1728
    9|  1  2   9  28  65  126  217  344  513  730  1001  1332  1729
   10|  0  2  10  30  68  130  222  350  520  738  1010  1342  1740
   11|  1  3  11  31  69  131  223  351  521  739  1011  1343  1741
   12|  0  2  12  36  80  150  252  392  576  810  1100  1452  1872

Rémy Sigrist, Table of n, a(n) for n = 0..10010
Index entries for sequences related to binary expansion of n

See A342707 for a similar sequence.

Cf. A000035, A000120, A000583, A000695, A001093, A002061, A002378, A002522, A002523, A005836, A007088, A011379, A027444, A033042, A033043, A033044, A033045, A033046, A033047, A033048, A033049, A034262, A053698, A071568, A098547, A104258.

T(n,k) = { my (v=0, e); while (n, n-=2^e=valuation(n,2); v+=k^e); v }

T(n, n) = A104258(n).
T(n, 0) = A000035(n).
T(n, 1) = A000120(n).
T(n, 2) = n.
T(n, 3) = A005836(n).
T(n, 4) = A000695(n).
T(n, 5) = A033042(n).
T(n, 6) = A033043(n).
T(n, 7) = A033044(n).
T(n, 8) = A033045(n).
T(n, 9) = A033046(n).
T(n, 10) = A007088(n).
T(n, 11) = A033047(n).
T(n, 12) = A033048(n).
T(n, 13) = A033049(n).
T(0, k) = 0.
T(1, k) = 1.
T(2, k) = k.
T(3, k) = k + 1.
T(4, k) = k^2.
T(5, k) = k^2 + 1 = A002522(k).
T(6, k) = k^2 + k = A002378(k).
T(7, k) = k^2 + k + 1 = A002061(k).
T(8, k) = k^3.
T(9, k) = k^3 + 1 = A001093(k).
T(10, k) = k^3 + k = A034262(k).
T(11, k) = k^3 + k + 1 = A071568(k).
T(12, k) = k^3 + k^2 = A011379(k).
T(13, k) = k^3 + k^2 + 1 = A098547(k).
T(14, k) = k^3 + k^2 + k = A027444(k).
T(15, k) = k^3 + k^2 + k + 1 = A053698(k).
T(16, k) = k^4 = A000583(k).
T(17, k) = k^4 + 1 = A002523(k).
T(m + n, k) = T(m, k) + T(n, k) when m AND n = 0 (where AND denotes the bitwise AND operator).

A260534 Square array read by ascending antidiagonals, T(n,k) = Sum_{j=0..k} n^j*(C(k-j,j) mod 2).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 4, 1, 3, 1, 1, 1, 5, 1, 7, 2, 1, 1, 1, 6, 1, 13, 5, 3, 1, 1, 1, 7, 1, 21, 10, 11, 1, 1, 1, 1, 8, 1, 31, 17, 31, 1, 4, 1, 1, 1, 9, 1, 43, 26, 69, 1, 23, 3, 1, 1, 1, 10, 1, 57, 37, 131, 1, 94, 21, 5, 1, 1, 1, 11

Offset: 0

Peter Luschny, Sep 20 2015

A parametrization of Stern's diatomic series (which is here T(1,k)). (For other generalizations of Dijkstra's fusc function see the Luschny link.)

			Array starts:
n\k[0, 1,  2, 3,  4,  5,   6, 7,    8,    9,    10]
[0] 1, 1,  1, 1,  1,  1,   1, 1,    1,    1,     1, ...
[1] 1, 1,  2, 1,  3,  2,   3, 1,    4,    3,     5, ... [A002487]
[2] 1, 1,  3, 1,  7,  5,  11, 1,   23,   21,    59, ... [A101624]
[3] 1, 1,  4, 1, 13, 10,  31, 1,   94,   91,   355, ...
[4] 1, 1,  5, 1, 21, 17,  69, 1,  277,  273,  1349, ... [A101625]
[5] 1, 1,  6, 1, 31, 26, 131, 1,  656,  651,  3881, ...
[6] 1, 1,  7, 1, 43, 37, 223, 1, 1339, 1333,  9295, ...
[7] 1, 1,  8, 1, 57, 50, 351, 1, 2458, 2451, 19559, ...
[8] 1, 1,  9, 1, 73, 65, 521, 1, 4169, 4161, 37385, ...
-,-,-,-,A002061,A002522,A071568,-,-,A059826,-,A002523,

Chai Wah Wu, Table of n, a(n) for n = 0..10010
Peter Luschny, Rational Trees and Binary Partitions.

Cf. A002061, A002487, A002522, A002523, A059826, A071568, A101624, A101625.

T := (n,k) -> add(modp(binomial(k-j,j),2)*n^j, j=0..k):
seq(lprint(seq(T(n,k),k=0..10)),n=0..5);

Table[If[(n - k) == 0, 1, Sum[(n - k)^j Mod[Binomial[k - j, j], 2], {j, 0, k}]], {n, 0, 10}, {k, 0, n}] (* Michael De Vlieger, Sep 21 2015 *)

def A260534_T(n,k):
    return sum(0 if ~(k-j) & j else n**j for j in range(k+1)) # Chai Wah Wu, Feb 08 2016

A300381 Numbers of the form x^3+x*y+y^3 where x,y>0.

Original entry on oeis.org

3, 11, 20, 31, 41, 63, 69, 80, 103, 131, 143, 144, 167, 209, 223, 236, 261, 275, 304, 351, 365, 371, 391, 435, 468, 503, 521, 536, 563, 601, 608, 677, 735, 739, 755, 776, 783, 829, 899, 911, 999, 1011, 1028, 1057, 1088, 1104, 1135, 1175, 1276, 1313, 1343, 1361, 1391, 1413, 1439, 1511

Offset: 1

R. J. Mathar, Mar 04 2018

nonn

Contains all terms larger than 1 from A071568 (contributed by y=1).

Cf. A003136, A266970 (prime subsequence), A024614.

A386407 a(n) = -floor(log(Integral_{x=2n^3+2n+2...oo} n^(-x^3) dx)/log(n)).

Original entry on oeis.org

10658, 238337, 2628081, 17984736, 88716544, 345948416, 1131366096, 3228667360, 8266914656, 19378404864, 42216896176, 86468159456, 168013851840, 311941643776, 556629367184, 959168691936, 1602434480416, 2604149515520, 4128340746096, 6399632545504

Offset: 2

Jason Bard, Jul 20 2025

nonn

-a(10) = -2022^3 - 8 was the solution to the final problem of the 2022 MIT Integral Bee Finals; see MIT link.

Jason Bard, Table of n, a(n) for n = 2..36
MIT, 2022 Integration Bee Finals. See Problem 5.
Prime Newtons, I solved this using Upper Incomplete Gamma, YouTube video.

Cf. A071568.

Table[-Floor[(Log[Gamma[1/3, 8 (n^3 + n + 1)^3*Log[n]]] - Log[3] - (1/3) Log[Log[n]])/Log[n]], {n, 2, 36}]

a(n) = -floor((log(Gamma(1/3, 8 * log(n) * (n^3 + n + 1)^3)) - log(3) - (1/3) * log(log(n))) / log(n)).

cp's OEIS Frontend

A163422 Primes p such that A071568((p-1)/2) is also prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

Extensions

A101220 a(n) = Sum_{k=0..n} Fibonacci(n-k)*n^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

Extensions

A104257 Square array T(a,n) read by antidiagonals: replace 2^i with a^i in binary representation of n, where a,n >= 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A098547 a(n) = n^3 + n^2 + 1.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A099721 a(n) = n^2*(2*n+1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A292741 Number A(n,k) of partitions of n with k sorts of part 1; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A341907 T(n, k) is the result of replacing 2^e with k^e in the binary expansion of n; square array T(n, k) read by antidiagonals upwards, n, k >= 0.

Views

Author

A099721 a(n) = n^2(2n+1).