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-10 of 26 results. Next

A049988 Number of nondecreasing arithmetic progressions of positive integers with sum n.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 7, 5, 7, 9, 9, 7, 14, 8, 11, 16, 13, 10, 20, 11, 17, 21, 16, 13, 27, 17, 18, 26, 22, 16, 35, 17, 23, 31, 23, 25, 41, 20, 25, 36, 33, 22, 46, 23, 31, 48, 30, 25, 52, 29, 38, 47, 36, 28, 57, 37, 41, 52, 37, 31, 71, 32, 39, 62, 44, 43, 69, 35, 45, 62, 57, 37, 79, 38
Offset: 0

Views

Author

Keywords

Comments

From Gus Wiseman, May 03 2019: (Start)
a(n) is the number of integer partitions of n with equal differences. The Heinz numbers of these partitions are given by A325328. For example, the a(1) = 1 through a(9) = 9 partitions are:
1 2 3 4 5 6 7 8 9
11 21 22 32 33 43 44 54
111 31 41 42 52 53 63
1111 11111 51 61 62 72
222 1111111 71 81
321 2222 333
111111 11111111 432
531
111111111
(End)
From Petros Hadjicostas, Sep 29 2019: (Start)
We show how Leroy Quet's g.f. Sum_{n >= 0} a(n)*x^n = 1/(1-x) + Sum_{k >= 2} x^k/(1-x^(k*(k-1)/2))/(1-x^k) in the Formula section below can be derived from Graeme McRae's g.f. for A049982 (see one of the links below).
Let b(n) = A049982(n) for n >= 1. Then Graeme McRae proved that Sum_{n >= 1} b(n)*x^n = Sum_{k >= 2} x^t(k)/(x^t(k) - x^t(k-1) - x^k + 1) = Sum_{k >= 2} x^t(k)/((1 - x^k) * (1 - x^t(k-1))), where t(k) = A000217(k) = k*(k+1)/2.
Since a(n) - b(n) = A000005(n) for n >= 1, to finish the proof, we only need to show that K(x) := 1 + Sum_{n >= 1} a(n)*x^n - Sum_{n >= 1} b(n)*x^n is the g.f. of A000005 (= number of divisors). But it is easy to show that K(x) = 1 + Sum_{k >= 1} x^k/(1 - x^k) = 1 + Sum_{n >= 1} A000005(n)*x^n (Lambert series for the number of divisors function). (End)

Crossrefs

Programs

  • Mathematica
    a[n_]:=If[n==0,1,Block[{i,c=Floor[(n-1)/2]+DivisorSigma[0,n]},Do[i=1;While[i*kGus Wiseman, May 07 2019 *)
    Table[Length[Select[IntegerPartitions[n],SameQ@@Differences[#]&]],{n,0,30}] (* Gus Wiseman, May 03 2019 *)
  • PARI
    seq(n)={Vec(1/(1-x) + sum(k=2, n, x^k/(1 - x^(k*(k-1)/2))/(1-x^k) + O(x*x^n)))} \\ Andrew Howroyd, Sep 28 2019

Formula

G.f.: 1/(1-x) + Sum_{k>=2} x^k/(1-x^(k*(k-1)/2))/(1-x^k). - Leroy Quet, Apr 08 2010. [Edited by Gus Wiseman, May 03 2019]
a(n) = A049982(n) + A000005(n) = A049980(n) + A000005(n) - 1 for n >= 1. - Petros Hadjicostas, Sep 28 2019

Extensions

Edited by Max Alekseyev, May 03 2010
a(0) = 1 prepended by Gus Wiseman, May 03 2019

A175342 Number of arithmetic progressions (where the difference between adjacent terms is either positive, 0, or negative) of positive integers that sum to n.

Original entry on oeis.org

1, 2, 4, 5, 6, 10, 8, 10, 15, 14, 12, 22, 14, 18, 28, 21, 18, 34, 20, 28, 38, 28, 24, 46, 31, 32, 48, 38, 30, 62, 32, 40, 58, 42, 46, 73, 38, 46, 68, 58, 42, 84, 44, 56, 90, 56, 48, 94, 55, 70, 90, 66, 54, 106, 70, 74, 100, 70, 60, 130, 62, 74, 118, 81, 82, 130, 68, 84, 120
Offset: 1

Views

Author

Leroy Quet, Apr 17 2010

Keywords

Examples

			From _Gus Wiseman_, May 15 2019: (Start)
The a(1) = 1 through a(8) = 10 compositions with equal differences:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (12)   (13)    (14)     (15)      (16)       (17)
             (21)   (22)    (23)     (24)      (25)       (26)
             (111)  (31)    (32)     (33)      (34)       (35)
                    (1111)  (41)     (42)      (43)       (44)
                            (11111)  (51)      (52)       (53)
                                     (123)     (61)       (62)
                                     (222)     (1111111)  (71)
                                     (321)                (2222)
                                     (111111)             (11111111)
(End)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],SameQ@@Differences[#]&]],{n,0,15}] (* returns a(0) = 1, Gus Wiseman, May 15 2019*)

Formula

a(n) = 2*A049988(n) - A000005(n).
G.f.: x/(1-x) + Sum_{k>=2} x^k * (1 + x^(k(k-1)/2)) / (1 - x^(k(k-1)/2)) / (1 -x^k).

Extensions

Edited and extended by Max Alekseyev, May 03 2010

A325368 Heinz numbers of integer partitions with distinct differences between successive parts.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The enumeration of these partitions by sum is given by A325325.

Examples

			Most small numbers are in the sequence, but the sequence of non-terms together with their prime indices begins:
    8: {1,1,1}
   16: {1,1,1,1}
   24: {1,1,1,2}
   27: {2,2,2}
   30: {1,2,3}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   40: {1,1,1,3}
   48: {1,1,1,1,2}
   54: {1,2,2,2}
   56: {1,1,1,4}
   60: {1,1,2,3}
   64: {1,1,1,1,1,1}
   72: {1,1,1,2,2}
   80: {1,1,1,1,3}
   81: {2,2,2,2}
   88: {1,1,1,5}
   90: {1,2,2,3}
   96: {1,1,1,1,1,2}
  100: {1,1,3,3}
		

Crossrefs

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],UnsameQ@@Differences[primeptn[#]]&]

A355536 Irregular triangle read by rows where row n lists the differences between adjacent prime indices of n; if n is prime, row n is empty.

Original entry on oeis.org

0, 1, 0, 0, 0, 2, 0, 1, 3, 1, 0, 0, 0, 1, 0, 0, 2, 2, 4, 0, 0, 1, 0, 5, 0, 0, 0, 3, 1, 1, 0, 0, 0, 0, 3, 6, 1, 0, 1, 0, 7, 4, 0, 0, 2, 1, 2, 0, 4, 0, 1, 8, 0, 0, 0, 1, 0, 2, 0, 5, 0, 5, 1, 0, 0, 2, 0, 0, 3, 6, 9, 0, 1, 1, 10, 0, 2, 0, 0, 0, 0, 0, 3, 1, 3, 0, 6
Offset: 2

Views

Author

Gus Wiseman, Jul 12 2022

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The version where zero is prepended to the prime indices is A287352.
One could argue that row n = 1 is empty, but adding it changes only the offset, not the data.

Examples

			Triangle begins (showing n, prime indices, differences*):
   2:    (1)       .
   3:    (2)       .
   4:   (1,1)      0
   5:    (3)       .
   6:   (1,2)      1
   7:    (4)       .
   8:  (1,1,1)    0 0
   9:   (2,2)      0
  10:   (1,3)      2
  11:    (5)       .
  12:  (1,1,2)    0 1
  13:    (6)       .
  14:   (1,4)      3
  15:   (2,3)      1
  16: (1,1,1,1)  0 0 0
		

Crossrefs

Row-lengths are A001222 minus one.
The prime indices are A112798, sum A056239.
Row-sums are A243055.
Constant rows have indices A325328.
The Heinz numbers of the rows plus one are A325352.
Strict rows have indices A325368.
Row minima are A355524.
Row maxima are A286470, also A355526.
An adjusted version is A358169, reverse A355534.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Differences[primeMS[n]],{n,2,100}]

A325352 Heinz number of the differences plus one of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 5, 1, 6, 1, 7, 3, 8, 1, 6, 1, 10, 5, 11, 1, 12, 2, 13, 4, 14, 1, 9, 1, 16, 7, 17, 3, 12, 1, 19, 11, 20, 1, 15, 1, 22, 6, 23, 1, 24, 2, 10, 13, 26, 1, 12, 5, 28, 17, 29, 1, 18, 1, 31, 10, 32, 7, 21, 1, 34, 19, 15, 1, 24, 1, 37, 6, 38
Offset: 1

Views

Author

Gus Wiseman, Apr 23 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The only fixed point is 1 because otherwise the sequence decreases omega (A001222) by one.

Examples

			The partition (3,2,2,1) with Heinz number 90 has differences plus one (2,1,2) with Heinz number 18, so a(90) = 18.
		

Crossrefs

Positions of m's are A008578 (m = 1), A001248 (m = 2), A006094 (m = 3), A030078 (m = 4), A090076 (m = 5).

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    db[n_]:=Times@@Prime/@(1+Differences[primeMS[n]]);
    Table[db[n],{n,100}]

A325349 Number of integer partitions of n whose augmented differences are distinct.

Original entry on oeis.org

1, 1, 1, 2, 3, 2, 4, 5, 7, 7, 12, 10, 13, 15, 21, 21, 31, 34, 38, 45, 55, 60, 71, 80, 84, 103, 119, 134, 152, 186, 192, 228, 263, 292, 321, 377, 399, 454, 514, 565, 618, 709, 752, 840, 958, 1050, 1140, 1297, 1402, 1568, 1755, 1901, 2080, 2343, 2524, 2758, 3074
Offset: 0

Views

Author

Gus Wiseman, Apr 23 2019

Keywords

Comments

The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
The Heinz numbers of these partitions are given by A325366.

Examples

			The a(1) = 1 through a(11) = 10 partitions (A = 10, B = 11):
  (1)  (2)  (3)   (4)   (5)   (6)   (7)    (8)    (9)    (A)    (B)
            (21)  (22)  (41)  (33)  (43)   (44)   (54)   (55)   (65)
                  (31)        (42)  (52)   (62)   (63)   (64)   (83)
                              (51)  (61)   (71)   (72)   (73)   (92)
                                    (421)  (422)  (81)   (82)   (A1)
                                           (431)  (522)  (91)   (443)
                                           (521)  (621)  (433)  (641)
                                                         (442)  (722)
                                                         (541)  (731)
                                                         (622)  (821)
                                                         (631)
                                                         (721)
For example, (4,4,3) has augmented differences (1,2,3), which are distinct, so (4,4,3) is counted under a(11).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@Differences[Append[#,1]]&]],{n,0,30}]

A014405 Number of arithmetic progressions of 3 or more positive integers, strictly increasing with sum n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 2, 1, 0, 3, 0, 1, 5, 1, 0, 6, 0, 2, 7, 2, 0, 8, 2, 2, 9, 3, 0, 13, 0, 2, 11, 3, 4, 15, 0, 3, 13, 6, 0, 18, 0, 4, 20, 4, 0, 19, 2, 8, 18, 5, 0, 23, 6, 6, 20, 5, 0, 30, 0, 5, 25, 6, 7, 29, 0, 6, 24, 15, 0, 32, 0, 6, 34, 7, 4, 34, 0, 14, 31, 7, 0, 39, 9, 7, 31, 9, 0, 49, 5, 9, 33, 8, 10, 42, 0, 12
Offset: 1

Views

Author

Keywords

Examples

			E.g., 15 = 1+2+3+4+5 = 1+5+9 = 2+5+8 = 3+5+7 = 4+5+6.
		

Crossrefs

Programs

  • PARI
    a(n)= t=0; st=0; forstep(s=(n-3)\3,1,-1, st++; for(c=1,st, m=3; w=m*(s+c); while(wRick L. Shepherd, Aug 30 2006

Formula

G.f.: Sum_{k >= 3} x^t(k)/(x^t(k) - x^t(k-1) - x^k + 1) = Sum_{k >= 3} x^t(k)/((1 - x^k) * (1 - x^t(k-1))), where t(k) = k*(k+1)/2 = A000217(k) is the k-th triangular number [Graeme McRae]. - Petros Hadjicostas, Sep 29 2019
a(n) = A049992(n) - A023645(n). - Antti Karttunen, Feb 20 2023

A325360 Heinz numbers of integer partitions whose differences are weakly increasing.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 71, 73
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (x, y, z) are (y - x, z - y). We adhere to this standard for integer partitions also even though they are always weakly decreasing. For example, the differences of (6,3,1) are (-3,-2).
The enumeration of these partitions by sum is given by A240026.

Examples

			Most small numbers are in the sequence. However, the sequence of non-terms together with their prime indices begins:
   18: {1,2,2}
   36: {1,1,2,2}
   50: {1,3,3}
   54: {1,2,2,2}
   70: {1,3,4}
   72: {1,1,1,2,2}
   75: {2,3,3}
   90: {1,2,2,3}
   98: {1,4,4}
  100: {1,1,3,3}
		

Crossrefs

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],OrderedQ[Differences[primeptn[#]]]&]

A307824 Heinz numbers of integer partitions whose augmented differences are all equal.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 11, 13, 15, 16, 17, 19, 23, 29, 31, 32, 37, 41, 43, 47, 53, 55, 59, 61, 64, 67, 71, 73, 79, 83, 89, 97, 101, 103, 105, 107, 109, 113, 119, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227
Offset: 1

Views

Author

Gus Wiseman, May 03 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The augmented differences aug(y) of an integer partition y of length k are given by aug(y)i = y_i - y{i + 1} + 1 if i < k and aug(y)_k = y_k. For example, aug(6,5,5,3,3,3) = (2,1,3,1,1,3).
The enumeration of these partitions by sum is given by A129654.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    3: {2}
    4: {1,1}
    5: {3}
    7: {4}
    8: {1,1,1}
   11: {5}
   13: {6}
   15: {2,3}
   16: {1,1,1,1}
   17: {7}
   19: {8}
   23: {9}
   29: {10}
   31: {11}
   32: {1,1,1,1,1}
   37: {12}
   41: {13}
   43: {14}
		

Crossrefs

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    aug[y_]:=Table[If[i
    				

A295235 Numbers k such that the positions of the ones in the binary representation of k are in arithmetic progression.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 17, 18, 20, 21, 24, 28, 30, 31, 32, 33, 34, 36, 40, 42, 48, 56, 60, 62, 63, 64, 65, 66, 68, 72, 73, 80, 84, 85, 96, 112, 120, 124, 126, 127, 128, 129, 130, 132, 136, 144, 146, 160, 168, 170, 192, 224, 240, 248
Offset: 1

Views

Author

Rémy Sigrist, Nov 18 2017

Keywords

Comments

Also numbers k of the form Sum_{b=0..h-1} 2^(i+j*b) for some h >= 0, i >= 0, j > 0 (in fact, h = A000120(k), and if k > 0, i = A007814(k)).
There is a simple bijection between the finite sets of nonnegative integers in arithmetic progression and the terms of this sequence: s -> Sum_{i in s} 2^i; the term 0 corresponds to the empty set.
For any n > 0, A054519(n) gives the numbers of terms with n+1 digits in binary representation.
For any n >= 0, n is in the sequence iff 2*n is in the sequence.
For any n > 0, A000695(a(n)) is in the sequence.
The first prime numbers in the sequence are: 2, 3, 5, 7, 17, 31, 73, 127, 257, 8191, 65537, 131071, 262657, 524287, ...
This sequence contains the following sequences: A000051, A000079, A000225, A000668, A002450, A019434, A023001, A048645.
For any k > 0, 2^k - 2, 2^k - 1, 2^k, 2^k + 1 and 2^k + 2 are in the sequence (e.g., 14, 15, 16, 17, and 18).
Every odd term is a binary palindrome (and thus belongs to A006995).
Odd terms are A064896. - Robert Israel, Nov 20 2017

Examples

			The binary representation of the number 42 is "101010" and has ones evenly spaced, hence 42 appears in the sequence.
The first terms, alongside their binary representations, are:
   n  a(n)  a(n) in binary
  --  ----  --------------
   1    0           0
   2    1           1
   3    2          10
   4    3          11
   5    4         100
   6    5         101
   7    6         110
   8    7         111
   9    8        1000
  10    9        1001
  11   10        1010
  12   12        1100
  13   14        1110
  14   15        1111
  15   16       10000
  16   17       10001
  17   18       10010
  18   20       10100
  19   21       10101
  20   24       11000
		

Crossrefs

Cf. A029931, A048793 (binary indices triangle), A070939, A291166, A325328 (prime indices rather than binary indices), A326669, A326675.

Programs

  • Maple
    f:= proc(d) local i,j,k;
      op(sort([seq(seq(add(2^(d-j*k),k=0..m),m=1..d/j),j=1..d),2^(d+1)]))
    end proc:
    0,1,seq(f(d),d=0..10); # Robert Israel, Nov 20 2017
  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[100],SameQ@@Differences[bpe[#]]&] (* Gus Wiseman, Jul 22 2019 *)
  • PARI
    is(n) = my(h=hammingweight(n)); if(h<3, return(1), my(i=valuation(n,2),w=#binary(n)); if((w-i-1)%(h-1)==0, my(j=(w-i-1)/(h-1)); return(sum(k=0,h-1,2^(i+j*k))==n), return(0)))
Showing 1-10 of 26 results. Next