A008805 -id:A008805 - cp's OEIS frontend

A005514 Number of n-bead bracelets (turnover necklaces) with 8 red beads and n-8 black beads.

1, 1, 5, 10, 29, 57, 126, 232, 440, 750, 1282, 2052, 3260, 4950, 7440, 10824, 15581, 21879, 30415, 41470, 56021, 74503, 98254, 127920, 165288, 211276, 268228, 337416, 421856, 523260, 645456, 790704, 963793, 1167645, 1408185

Offset: 8

N. J. A. Sloane

From Vladimir Shevelev, Apr 23 2011: (Start)
Also number of non-equivalent necklaces of 8 beads each of them painted by one of n colors.
The sequence solves the so-called Reis problem about convex k-gons in case k=8 (see our comment at A032279). (End)
From Petros Hadjicostas, Jul 14 2018: (Start)
Let (c(n): n >= 1) be a sequence of nonnegative integers and let C(x) = Sum_{n>=1} c(n)*x^n be its g.f. Let k be a positive integer. Let a_k = (a_k(n): n >= 1) be the output sequence of the DIK[k] transform of sequence (c(n): n >= 1), and let A_k(x) = Sum_{n>=1} a_k(n)*x^n be its g.f. See Bower's web link below. It can be proved that, when k is even, A_k(x) = ((1/k)*Sum_{d|k} phi(d)*C(x^d)^(k/d) + (1/2)*C(x^2)^((k/2)-1)*(C(x)^2 + C(x^2)))/2.
For this sequence, k=8, c(n) = 1 for all n >= 1, and C(x) = x/(1-x). Thus, a(n) = a_8(n) for all n >= 1. Since a_k(n) = 0  for 1 <= n <= k-1, the offset of this sequence is n = k = 8. Applying the formula for the g.f. of DIK[8] of (c(n): n >= 1) with C(x) = x/(1-x) and k=8, we get Herbert Kociemba's formula below.
Here, a(n) is defined to be the number of n-bead bracelets of two colors with 8 red beads and n-8 black beads. But it is also the number of dihedral compositions of n with 8 parts. (This statement is equivalent to Vladimir Shevelev's statement above that a(n) is the "number of non-equivalent necklaces of 8 beads each of them painted by one of n colors." By "necklaces", he means "turnover necklaces". See the second paragraph of Section 2 in his 2004 paper in the Indian Journal of Pure and Applied Mathematics.)
Two cyclic compositions of n (with k = 8 parts) belong to the same equivalence class corresponding to a dihedral composition of n if and only if one can be obtained from the other by a rotation or reversal of order. (End)

			From _Petros Hadjicostas_, Jul 14 2018: (Start)
Every n-bead bracelet of two colors such that 8 beads are red and n-8 are black can be transformed into a dihedral composition of n with 8 parts in the following way. Start with one R bead and go in one direction (say clockwise) until you reach the next R bead. Continue this process until you come back to the original R bead.
Let b_i be the number of beads from R bead i until you reach the last B bead before R bead i+1 (or R bead 1). Here, b_i = 1 iff there are no B beads between R bead i and R bead i+1 (or R bead 8 and R bead 1). Then b_1 + b_2 + ... + b_8 = n, and we get a dihedral composition of n. (Of course, b_2 + b_3 + ... + b_8 + b_1 and b_8 + b_7 + ... + b_1 belong to the same equivalence class of the dihedral composition b_1 + ... + b_8.)
For example, a(10) = 5, and we have the following bracelets with 8 R beads and 2 B beads. Next to the bracelets we list the corresponding dihedral compositions of n with k=8 parts (they must be viewed on a circle):
RRRRRRRRBB <-> 1+1+1+1+1+1+1+3
RRRRRRRBRB <-> 1+1+1+1+1+1+2+2
RRRRRRBRRB <-> 1+1+1+1+1+2+1+2
RRRRRBRRRB <-> 1+1+1+1+2+1+1+2
RRRRBRRRRB <-> 1+1+1+2+1+1+1+2
(End)

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
N. Zagaglia Salvi, Ordered partitions and colourings of cycles and necklaces, Bull. Inst. Combin. Appl., 27 (1999), 37-40.

Andrew Howroyd, Table of n, a(n) for n = 8..1000
Christian G. Bower, Transforms (2)
S. J. Cyvin, B. N. Cyvin, J. Brunvoll, I. Gutman, Chen Rong-si, S. El-Basil, and Zhang Fuji, Polygonal systems including the corannulene and coronene homologs: novel applications of Pólya's theorem, Z. Naturforsch., 52a (1997), 867-873.
Hansraj Gupta, Enumeration of incongruent cyclic k-gons, Indian J. Pure and Appl. Math., 10 (1979), no. 8, 964-999.
W. D. Hoskins and Anne Penfold Street, Twills on a given number of harnesses, J. Austral. Math. Soc. Ser. A 33 (1982), no. 1, 1-15.
W. D. Hoskins and A. P. Street, Twills on a given number of harnesses, J. Austral. Math. Soc. (Series A), 33 (1982), 1-15. (Annotated scanned copy)
Richard H. Reis, A formula for C(T) in Gupta's paper, Indian J. Pure and Appl. Math., 10 (1979), no. 8, 1000-1001.
Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc.
Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
Vladimir Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
Vladimir Shevelev, Necklaces and convex k-gons, Indian J. Pure and Appl. Math., 35 (2004), no. 5, 629-638.
Vladimir Shevelev, Spectrum of permanent's values and its extremal magnitudes in Lambda_n^3 and Lambda_n(alpha,beta,gamma), arXiv:1104.4051 [math.CO], 2011. (Cf. Section 5).
A. P. Street, Letter to N. J. A. Sloane, N.D.
Index entries for sequences related to bracelets
Index entries for linear recurrences with constant coefficients, signature (4,-4,-4,11,-8,0,8,-10,0,8,0,-10,8,0,-8,11,-4,-4,4,-1).

Column k=8 of A052307.

Cf. A008805, A032193, A032279, A032280, A032281, A032282, A119963, A292906.

k = 8; Table[(Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n + Binomial[If[OddQ[n], n - 1, n - If[OddQ[k], 2, 0]]/2, If[OddQ[k], k - 1, k]/2])/2, {n, k, 50}] (* Robert A. Russell, Sep 27 2004 *)
k=8;CoefficientList[Series[x^k*(1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])+(1+x)/(1-x^2)^Floor[(k+2)/2])/2,{x,0,50}],x] (* Herbert Kociemba, Nov 04 2016 *)

S. J. Cyvin et al. (1997) give a g.f. (See equation (18) on p. 870 of their paper. Their g.f. is the same as the one given by V. Jovovic below except for the extra x^8.) - Petros Hadjicostas, Jul 14 2018
G.f.: (x^8/16)*(1/(1 - x)^8 + 4/(1 - x^8) + 5/(1 - x^2)^4 + 2/(1 - x^4)^2 + 4/(1 - x)^2/(1 - x^2)^3) = x^8*(2*x^10 - 3*x^9 + 7*x^8 - 6*x^7 + 7*x^6 - 2*x^5 + 2*x^4 - 2*x^3 + 5*x^2 - 3*x + 1)/(1 - x)^8/(1 + x)^4/(1 + x^2)^2/(1 + x^4). - Vladeta Jovovic, Jul 17 2002
a(n) = ((n+4)/32)*s(n,0,8) + ((n-4)/32)*s(n,4,8) + (48*C(n-1,7) + (n+1)*(n-2)*(n-4)*(n-6))/768, if n is even >= 8; a(n) = (48*C(n-1,7) + (n-1)*(n-3)*(n-5)*(n-7))/768, if n odd >= 8, where s(n,k,d)=1, if n == k (mod d), and 0 otherwise. - Vladimir Shevelev, Apr 23 2011
G.f.: k=8, x^k*((1/k)*Sum_{d|k} phi(d)*(1-x^d)^(-k/d) + (1+x)/(1-x^2)^floor((k+2)/2))/2. - Herbert Kociemba, Nov 05 2016 [edited by Petros Hadjicostas, Jul 18 2018]
From Petros Hadjicostas, Jul 14 2018: (Start)
a(n) = (A032193(n) + A119963(n, 8))/2 = (A032193(n) + C(floor(n/2), 4))/2 for n >= 8.
The sequence (a(n): n >= 8) is the output sequence of Bower's "DIK[ 8 ]" (bracelet, indistinct, unlabeled, 8 parts) transform of 1, 1, 1, 1, ...
(End)

Sequence extended and description corrected by Christian G. Bower

Name edited by Petros Hadjicostas, Jul 20 2018

A278688 Triangle read by rows T(n, k) = number of non-equivalent ways to place k non-attacking ferses on an n X n board.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 3, 6, 7, 6, 2, 1, 1, 3, 17, 45, 92, 99, 76, 27, 7, 1, 6, 43, 225, 832, 2102, 3773, 4860, 4643, 3356, 1868, 795, 248, 56, 8, 1, 1, 6, 84, 709, 4500, 19987, 66201, 164423, 314224, 465230, 540247, 492206, 352300, 195717, 83247, 26083, 5754, 780, 55

Offset: 1

Heinrich Ludwig, Nov 27 2016

The triangle T(n, k) is irregularly shaped: 0 <= k <= A093005(n), which means that A093005(n) is the maximal number of non-attacking ferses that can be placed on an n X n board. First row corresponds to n = 1. First column corresponds to k = 0.
Two placements that differ by rotation or reflection are counted only once.
A fers is a fairy chess piece attacking one step ne-nw-sw-se.

			Triangle begins:
1, 1;
1, 1,  1;
1, 3,  6,  7,  6,  2,  1;
1, 3, 17, 45, 92, 99, 76, 27, 7;

Heinrich Ludwig, Table of n, a(n) for n = 1..130
Wikipedia, Fairy chess piece

Cf. A008805, A232567, A278682, A278683, A278684, A278685, A278686, (columns 2 through 8 of this sequence, respectively), A278687, A093005 (row length - 1).

A049798 a(n) = (1/2)*Sum_{k = 1..n} T(n,k), array T as in A049800.

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 2, 2, 3, 7, 2, 7, 10, 8, 8, 15, 11, 19, 16, 15, 22, 32, 19, 25, 34, 34, 33, 46, 33, 47, 47, 48, 61, 65, 45, 62, 77, 79, 68, 87, 74, 94, 97, 86, 105, 127, 98, 114, 120, 124, 129, 154, 141, 151, 142, 147, 172, 200, 151, 180

Offset: 1

Clark Kimberling

a(n) is the sum of the remainders after dividing each larger part by its corresponding smaller part for each partition of n+1 into two parts. - Wesley Ivan Hurt, Dec 20 2020

			From _Lei Zhou_, Mar 10 2014: (Start)
For n = 3, n+1 = 4, floor((n+1)/2) = 2, mod(4,2) = 0, and so a(3) = 0.
For n = 4, n+1 = 5, floor((n+1)/2) = 2, mod(5,2) = 1, and so a(4) = 1.
...
For n = 12, n+1 = 13, floor((n+1)/2) = 6, mod(13,2) = 1, mod(13,3) = 1, mod(13,4) = 1, mod(13,5) = 3, mod(13,6) = 1, and so a(12) = 1 + 1 + 1 + 3 + 1 = 7. (End)

Lei Zhou, Table of n, a(n) for n = 1..10000

Cf. A004125, A008611, A008805, A049797, A049799, A049801.

Half row sums of A049800.

List([1..60], n-> Sum([1..n], k-> (n+1) mod Int((k+1)/2))/2 ); # G. C. Greubel, Dec 09 2019

[ (&+[(n+1) mod Floor((k+1)/2): k in [1..n]])/2: n in [1..60]]; // G. C. Greubel, Dec 09 2019

seq( add( (n+1) mod floor((k+1)/2), k=1..n)/2, n=1..60); # G. C. Greubel, Dec 09 2019

Table[Sum[Mod[n+1, Floor[(k+1)/2]], {k,n}]/2, {n, 60}] (* G. C. Greubel, Dec 09 2019 *)

vector(60, n, sum(k=1,n, lift(Mod(n+1, (k+1)\2)) )/2 ) \\ G. C. Greubel, Dec 09 2019

def A049798(n): return sum((n+1)%k for k in range(2,(n+1>>1)+1)) # Chai Wah Wu, Oct 20 2023

def a(n):
    return sum([(n+1)%k for k in range(2,floor((n+3)/2))])
# Ralf Stephan, Mar 14 2014

a(n) = Sum_{k=2..floor((n+1)/2)} ((n+1) mod k). - Lei Zhou, Mar 10 2014
a(n) = A004125(n+1) - A008805(n-2), for n >= 2. - Carl Najafi, Jan 31 2013
a(n) = Sum_{i = 1..ceiling(n/2)} ((n-i+1) mod i). - Wesley Ivan Hurt, Jan 05 2017

A115514 Triangle read by rows: row n >= 1 lists first n positive terms of A004526 (integers repeated) in decreasing order.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 3, 2, 2, 1, 1, 3, 3, 2, 2, 1, 1, 4, 3, 3, 2, 2, 1, 1, 4, 4, 3, 3, 2, 2, 1, 1, 5, 4, 4, 3, 3, 2, 2, 1, 1, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1, 6, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1, 6, 6, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1, 7, 6, 6, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1, 7, 7, 6, 6, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1

Offset: 1

Roger L. Bagula, Mar 07 2006

T(n,k) = number of 2-element subsets of {1,2,...,n+2} such that the absolute difference of the elements is k+1, where 1 <= k < = n. E.g., T(7,3) = 3, the subsets are {1,5}, {2,6}, and {3,7}. - Christian Barrientos, Jun 27 2022

			Triangle begins as, for n >= 1, 1 <= k <= n,
  1;
  1, 1;
  2, 1, 1;
  2, 2, 1, 1;
  3, 2, 2, 1, 1;
  3, 3, 2, 2, 1, 1;
  4, 3, 3, 2, 2, 1, 1;
  ...

G. C. Greubel, Rows n = 1..100 of the triangle, flattened

Cf. A000012, A000073, A004526, A008805, A008967.

Cf. A002620 (row sums), A008805 (diagonal sums), A142150 (alternating row sums)

[Floor((n-k+2)/2): k in [1..n], n in [1..15]]; // G. C. Greubel, Mar 14 2024

# Assuming offset 0:
Even := n -> (1 + (-1)^n)/2: # Iverson's even.
p := n -> add(add(Even(k)*x^j, j = 0..n-k), k = 0..n):
for n from 0 to 9 do seq(coeff(p(n), x, k), k=0..n) od; # Peter Luschny, Jun 03 2021

Table[Floor[(n-k+2)/2], {n,15}, {k,n}]//Flatten (* G. C. Greubel, Mar 14 2024 *)

flatten([[(n-k+2)//2 for k in range(1,n+1)] for n in range(1,16)]) # G. C. Greubel, Mar 14 2024

Sum_{k=1..n} T(n, k) = A002620(n+1) (row sums). - Gary W. Adamson, Oct 25 2007
T(n, k) = [x^k] p(n), where p(n) are partial Gaussian polynomials (A008967) defined by p(n) = Sum_{k=0..n} Sum_{j=0..n-k} even(k)*x^j, and even(k) = 1 if k is even and otherwise 0. We assume offset 0. - Peter Luschny, Jun 03 2021
T(n, k) = floor((n+2-k)/2). - Christian Barrientos, Jun 27 2022
From G. C. Greubel, Mar 14 2024: (Start)
T(n, k) = A128623(n, k)/n.
Sum_{k=1..n} (-1)^(k-1)*T(n, k) = A142150(n+1).
Sum_{k=1..floor((n+1)/2)} T(n-k+1, k) = A008805(n-1).
Sum_{k=1..floor((n+1)/2)} (-1)^(k-1)*T(n-k+1, k) = A002265(n+3). (End)

Edited by N. J. A. Sloane, Mar 23 2008 and Dec 15 2017

A128508 Number of partitions p of n such that max(p) - min(p) = 3.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 3, 3, 7, 7, 12, 14, 20, 22, 32, 34, 45, 51, 63, 69, 87, 93, 112, 124, 144, 156, 184, 196, 225, 245, 275, 295, 335, 355, 396, 426, 468, 498, 552, 582, 637, 679, 735, 777, 847, 889, 960, 1016, 1088, 1144, 1232, 1288, 1377, 1449, 1539, 1611, 1719

Offset: 0

John W. Layman, May 07 2007

nonn

See A008805 and A049820 for the numbers of partitions p of n such that max(p)-min(p)=1 or 2, respectively.

Alois P. Heinz, Table of n, a(n) for n = 0..1000
G. E. Andrews, M. Beck and N. Robbins, Partitions with fixed differences between largest and smallest parts, arXiv:1406.3374 [math.NT], 2014

Cf. A008805, A049820, A097364, A116685, A117142, A117143.

np[n_]:=Length[Select[IntegerPartitions[n],Max[#]-Min[#]==3&]]; Array[np,60] (* Harvey P. Dale, Jul 02 2012 *)

Conjecture. a(1)=0 and, for n>1, a(n+1)=a(n)+d(n), where d(n) is defined as follows: d=0,0,0,1,0 for n=1,...,5 and, for n>5, d(n)=d(n-2)+1 if n=6k or n=6k+4, d(n)=d(n-2) if n=6k+1 or n=6k+3, d(n)=d(n-2)+2Floor[n/6] if n=6k+2 and d(n)=d(n-5) if n=6k+5.
G.f. for number of partitions p of n such that max(p)-min(p) = m is Sum_{k>0} x^(2*k+m)/Product_{i=0..m} (1-x^(k+i)). - Vladeta Jovovic, Jul 04 2007
a(n) = A097364(n,3) = A116685(n,3) = A117143(n) - A117142(n). - Alois P. Heinz, Nov 02 2012

More terms from Vladeta Jovovic, Jul 04 2007

A243716 Irregular triangle read by rows: T(n, k) = number of inequivalent (mod the dihedral group D_8 of order 8) ways to place k nonattacking knights on an n X n board.

Original entry on oeis.org

1, 1, 2, 1, 1, 3, 7, 9, 6, 2, 3, 18, 40, 66, 49, 30, 8, 3, 6, 43, 195, 609, 1244, 1767, 1710, 1148, 510, 154, 31, 6, 1, 6, 83, 618, 3375, 12329, 32524, 61731, 86748, 90059, 70128, 40770, 18053, 6089, 1643, 344, 61, 7, 1, 10, 156, 1751, 14181, 81900, 348541

Offset: 1

Heinrich Ludwig, Jun 10 2014

The triangle is irregularly shaped: 1 <= k <= A030978(n). A030978(n) is the maximal number of knights that can be placed on an n X n board.
First row corresponds to n = 1.
Counting "inequivalent ways" means: Rotations or reflections of a placement of knights on the board are considered to be the same placement.

			The triangle begins:
  1;
  1,  2,   1,   1;
  3,  7,   9,   6,    2;
  3, 18,  40,  66,   49,   30,    8,    3;
  6, 43, 195, 609, 1244, 1767, 1710, 1148, 510, 154, 31, 6, 1;
  ...

Heinrich Ludwig, Table of n, a(n) for n = 1..116

Cf. A030978, A008805 (column 1), A243717 (column 2), A243718 (column 3), A243719 (column 4), A243720 (column 5).

A232569 Triangle T(n, k) = number of non-equivalent (mod D_4) binary n X n matrices with k pairwise not adjacent 1's; k=0,...,n^2.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 0, 1, 3, 6, 6, 3, 1, 0, 0, 0, 0, 1, 3, 17, 40, 62, 45, 20, 4, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 6, 43, 210, 683, 1425, 1936, 1696, 977, 366, 101, 21, 5, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 6, 84, 681, 4015, 16149, 46472, 95838, 143657

Offset: 1

Heinrich Ludwig, Nov 29 2013

Also number of non-equivalent ways to place k non-attacking wazirs on an n X n board.
Two matrix elements are considered adjacent, if the difference of their row indices is 1 and the column indices are equal, or vice versa (von Neumann neighborhood).
Counted for this sequence are equivalence classes induced by the dihedral group D_4. If equivalent matrices are being destinguished, the corresponding numbers are A232833(n).
Row index starts from n = 1, column index k ranges from 0 to n^2.
T(n, 1) = A008805(n-1); T(n, 2) = A232567(n) for n >= 2; T(n, 3) = A232568(n) for n >= 2;
Into an n X n binary matrix there can be placed maximally A000982(n) = ceiling(n^2/2) pairwise not adjacent 1's.

			Triangle begins:
1,1;
1,1,1,0,0;
1,3,6,6,3,1,0,0,0,0;
1,3,17,40,62,45,20,4,1,0,0,0,0,0,0,0,0;
1,6,43,210,683,1425,1936,1696,977,366,101,21,5,1,0,0,0,0,0,0,0,0,0,0,0,0;
...
There are T(3, 2) = 6 non-equivalent binary 3 X 3 matrices with 2 not adjacent 1's (and no other 1's):
  [1 0 0]   [0 1 0]   [1 0 0]   [0 1 0]   [1 0 1]   [1 0 0]
  |0 0 0|   |0 0 0|   |0 1 0|   |1 0 0|   |0 0 0|   |0 0 1|
  [0 0 1]   [0 1 0]   [0 0 0]   [0 0 0]   [0 0 0]   [0 0 0]

Heinrich Ludwig, Rows n = 1..8 of irregular triangle, flattened

Cf. A232567, A232568, A239576, A008805, A000982, A201511, A232833.

A325198 Positive numbers whose maximum prime index minus minimum prime index is 2.

Original entry on oeis.org

10, 20, 21, 30, 40, 50, 55, 60, 63, 80, 90, 91, 100, 105, 120, 147, 150, 160, 180, 187, 189, 200, 240, 247, 250, 270, 275, 300, 315, 320, 360, 385, 391, 400, 441, 450, 480, 500, 525, 540, 551, 567, 600, 605, 637, 640, 713, 720, 735, 750, 800, 810, 900, 945

Offset: 1

Gus Wiseman, Apr 11 2019

nonn

Also Heinz numbers of integer partitions whose maximum minus minimum part is 2 (counted by A008805). The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			The sequence of terms together with their prime indices begins:
   10: {1,3}
   20: {1,1,3}
   21: {2,4}
   30: {1,2,3}
   40: {1,1,1,3}
   50: {1,3,3}
   55: {3,5}
   60: {1,1,2,3}
   63: {2,2,4}
   80: {1,1,1,1,3}
   90: {1,2,2,3}
   91: {4,6}
  100: {1,1,3,3}
  105: {2,3,4}
  120: {1,1,1,2,3}
  147: {2,4,4}
  150: {1,2,3,3}
  160: {1,1,1,1,1,3}
  180: {1,1,2,2,3}
  187: {5,7}

Robert Israel, Table of n, a(n) for n = 1..10000

Positions of 2's in A243055.
A061395(a(n)) - A055396((n)) = 2.
Cf. A000961, A008805, A046660, A056239, A093641, A112798, A118914, A174090, A195086, A256617, A325180, A325197.

N:= 1000: # for terms <= N
q:= 2: r:= 3:
Res:= NULL:
do
  p:= q; q:= r; r:= nextprime(r);
  if p*r > N then break fi;
  for i from 1 do
    pi:= p^i;
    if pi*r > N then break fi;
    for j from 0 do
      piqj:= pi*q^j;
      if piqj*r > N then break fi;
      Res:= Res, seq(piqj*r^k,k=1 .. floor(log[r](N/piqj)))
    od
  od
od:
sort([Res]); # Robert Israel, Apr 12 2019

Select[Range[100],PrimePi[FactorInteger[#][[-1,1]]]-PrimePi[FactorInteger[#][[1,1]]]==2&]

A330640 a(n) is the number of partitions of n with Durfee square of size <= 2.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 15, 22, 29, 40, 51, 67, 83, 105, 127, 156, 185, 222, 259, 305, 351, 407, 463, 530, 597, 676, 755, 847, 939, 1045, 1151, 1272, 1393, 1530, 1667, 1821, 1975, 2147, 2319, 2510, 2701, 2912, 3123, 3355, 3587, 3841, 4095, 4372, 4649, 4950, 5251, 5577, 5903, 6255, 6607, 6986

Offset: 0

Omar E. Pol, Dec 22 2019

This is an easy sequence since A006918 is the partial sums of A008805 (triangular numbers repeated).

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2,1,-4,1,2,-1).

Cf. A000041, A006918, A008805, A028310, A115994, A115720.

Vec((1 - x - x^2 + 2*x^3 - x^5 + x^6) / ((1 - x)^4*(1 + x)^2) + O(x^60)) \\ Colin Barker, Dec 31 2019

a(n) = A028310(n), 0 <= n <= 2.
a(n) = A028310(n) + A006918(n-3), n >= 3.
Or without A028310:
a(0) = 1, a(1) = 1, a(2) = 2.
a(n) = n + A006918(n-3), n >= 3.
From Colin Barker, Dec 31 2019: (Start)
G.f.: (1 - x - x^2 + 2*x^3 - x^5 + x^6) / ((1 - x)^4*(1 + x)^2).
a(n) = 2*a(n-1) + a(n-2) - 4*a(n-3) + a(n-4) + 2*a(n-5) - a(n-6) for n>6.
a(n) = (3 - 3*(-1)^n + (49+3*(-1)^n)*n - 6*n^2 + 2*n^3) / 48.
(End)

A343875 Array read by antidiagonals: T(n,k) is the number of n X n nonnegative integer matrices with sum of elements equal to k, up to rotations and reflections.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 3, 3, 1, 0, 1, 4, 11, 3, 1, 0, 1, 8, 31, 24, 6, 1, 0, 1, 10, 84, 113, 55, 6, 1, 0, 1, 16, 198, 528, 410, 99, 10, 1, 0, 1, 20, 440, 2003, 2710, 1091, 181, 10, 1, 0, 1, 29, 904, 6968, 15233, 10488, 2722, 288, 15, 1, 0, 1, 35, 1766, 21593, 75258, 82704, 34399, 5806, 461, 15, 1

Offset: 0

Andrew Howroyd, May 06 2021

			Array begins:
=====================================================
n\k | 0  1   2    3     4      5       6        7
----+------------------------------------------------
  0 | 1  0   0    0     0      0       0        0 ...
  1 | 1  1   1    1     1      1       1        1 ...
  2 | 1  1   3    4     8     10      16       20 ...
  3 | 1  3  11   31    84    198     440      904 ...
  4 | 1  3  24  113   528   2003    6968    21593 ...
  5 | 1  6  55  410  2710  15233   75258   331063 ...
  6 | 1  6  99 1091 10488  82704  563864  3376134 ...
  7 | 1 10 181 2722 34399 360676 3235551 25387944 ...
  ...

Andrew Howroyd, Table of n, a(n) for n = 0..1325

Rows n=0..3 are A000007, A000012, A005232, A054343.
Columns 0..1 are A000012, A008805(n-1).
Cf. A054252 (binary case), A318795, A343097, A343874.

U(n,s) = {(s(1)^(n^2) + s(1)^(n%2)*(2*s(4)^(n^2\4) + s(2)^(n^2\2)) + 2*s(1)^n*s(2)^(n*(n-1)/2) + 2*(s(1)^(n%2)*s(2)^(n\2))^n )/8}
T(n,k)={polcoef(U(n,i->1/(1-x^i) + O(x*x^k)), k)}

cp's OEIS Frontend

A005514 Number of n-bead bracelets (turnover necklaces) with 8 red beads and n-8 black beads.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A278688 Triangle read by rows T(n, k) = number of non-equivalent ways to place k non-attacking ferses on an n X n board.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A049798 a(n) = (1/2)*Sum_{k = 1..n} T(n,k), array T as in A049800.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Python

Sage

Formula

A115514 Triangle read by rows: row n >= 1 lists first n positive terms of A004526 (integers repeated) in decreasing order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

SageMath

Formula

Extensions

A128508 Number of partitions p of n such that max(p) - min(p) = 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A243716 Irregular triangle read by rows: T(n, k) = number of inequivalent (mod the dihedral group D_8 of order 8) ways to place k nonattacking knights on an n X n board.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A232569 Triangle T(n, k) = number of non-equivalent (mod D_4) binary n X n matrices with k pairwise not adjacent 1's; k=0,...,n^2.

Views

Author

Keywords

Comments

Examples

Links