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 18 results. Next

A007294 Number of partitions of n into nonzero triangular numbers.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 4, 4, 4, 6, 7, 7, 10, 11, 11, 15, 17, 17, 22, 24, 25, 32, 35, 36, 44, 48, 50, 60, 66, 68, 81, 89, 92, 107, 117, 121, 141, 153, 159, 181, 197, 205, 233, 252, 262, 295, 320, 332, 372, 401, 417, 465, 501, 520, 575, 619, 645, 710, 763
Offset: 0

Views

Author

Keywords

Comments

Also number of decreasing integer sequences l(1) >= l(2) >= l(3) >= .. 0 such that sum('i*l(i)','i'=1..infinity)=n.
a(n) is also the number of partitions of n such that #{parts equal to i} >= #{parts equal to j} if i <= j.
Also the number of partitions of n (necessarily into distinct parts) where the part sizes are monotonically decreasing (including the last part, which is the difference between the last part and a "part" of size 0). These partitions are the conjugates of the partitions with number of parts of size i increasing. - Franklin T. Adams-Watters, Apr 08 2008
Also partitions with condition as in A179255, and additionally, if more than one part, first difference >= first part: for example, a(10)=7 as there are 7 such partitions of 10: 1+2+3+4 = 1+2+7 = 1+3+6 = 1+9 = 2+8 = 3+7 = 10. - Joerg Arndt, Mar 22 2011
Number of members of A181818 with a bigomega value of n (cf. A001222). - Matthew Vandermast, May 19 2012

Examples

			6 = 3+3 = 3+1+1+1 = 1+1+1+1+1+1 so a(6) = 4.
a(7)=4: Four sequences as above are (7,0,..), (5,1,0,..), (3,2,0,..),(2,1,1,0,..). They correspond to the partitions 1^7, 2 1^5, 2^2 1^3, 3 2 1^2 of seven or in the main description to the partitions 1^7, 3 1^4, 3^2 1, 6 1.
From _Gus Wiseman_, May 03 2019: (Start)
The a(1) = 1 through a(9) = 6 partitions using nonzero triangular numbers are the following. The Heinz numbers of these partitions are given by A325363.
  1   11   3     31     311     6        61        611        63
           111   1111   11111   33       331       3311       333
                                3111     31111     311111     6111
                                111111   1111111   11111111   33111
                                                              3111111
                                                              111111111
The a(1) = 1 through a(10) = 7 partitions with weakly decreasing multiplicities are the following. Equivalent to Matthew Vandermast's comment, the Heinz numbers of these partitions are given by A025487 (products of primorial numbers).
  1  11  21   211   2111   321     3211     32111     32211      4321
         111  1111  11111  2211    22111    221111    222111     322111
                           21111   211111   2111111   321111     2221111
                           111111  1111111  11111111  2211111    3211111
                                                      21111111   22111111
                                                      111111111  211111111
                                                                 1111111111
The a(1) = 1 through a(11) = 7 partitions with weakly increasing differences (where the last part is taken to be zero) are the following. The Heinz numbers of these partitions are given by A325362 (A = 10, B = 11).
  (1)  (2)  (3)   (4)   (5)   (6)    (7)    (8)    (9)    (A)     (B)
            (21)  (31)  (41)  (42)   (52)   (62)   (63)   (73)    (83)
                              (51)   (61)   (71)   (72)   (82)    (92)
                              (321)  (421)  (521)  (81)   (91)    (A1)
                                                   (531)  (631)   (731)
                                                   (621)  (721)   (821)
                                                          (4321)  (5321)
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A102462.
Row sums of array A176723 and triangle A176724. - Wolfdieter Lang, Jul 19 2010
Cf. A179255 (condition only on differences), A179269 (parts strictly increasing instead of nondecreasing). - Joerg Arndt, Mar 22 2011
Row sums of A319797.

Programs

  • Haskell
    a007294 = p $ tail a000217_list where
       p _      0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Jun 28 2013
    
  • Maple
    b:= proc(n,i) option remember;
          if n<0 then 0
        elif n=0 then 1
        elif i=0 then 0
        else b(n, i-1) +b(n-i*(i+1)/2, i)
          fi
        end:
    a:= n-> b(n, floor(sqrt(2*n))):
    seq(a(n), n=0..100);  # Alois P. Heinz, Mar 22 2011
    isNondecrP :=proc(L) slp := DIFF(DIFF(L)) ; min(op(%)) >= 0 ; end proc:
    A007294 := proc(n) local a, p; a := 0 ; if n = 0 then return 1 ; end if; for p in combinat[partition](n) do if nops(p) = nops(convert(p, set)) then if isNondecrP(p) then if nops(p) =1 then a := a+1 ; elif op(2, p) >= 2*op(1, p) then a := a+1; end if; end if; end if; end do; a ; end proc:
    seq(A007294(n), n=0..30) ; # R. J. Mathar, Jan 07 2011
  • Mathematica
    CoefficientList[ Series[ 1/Product[1 - x^(i(i + 1)/2), {i, 1, 50}], {x, 0, 70}], x]
    (* also *)
    t = Table[n (n + 1)/2, {n, 1, 200}] ; p[n_] := IntegerPartitions[n, All, t]; Table[p[n], {n, 0, 12}] (*shows partitions*)
    a[n_] := Length@p@n; a /@Range[0, 80]
    (* Clark Kimberling, Mar 09 2014 *)
    b[n_, i_] := b[n, i] = Which[n < 0, 0, n == 0, 1, i == 0, 0, True, b[n, i-1]+b[n-i*(i+1)/2, i]]; a[n_] := b[n, Floor[Sqrt[2*n]]]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Apr 09 2014, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],OrderedQ[Differences[Append[#,0]]]&]],{n,0,30}] (* Gus Wiseman, May 03 2019 *)
    nmax = 58; t = Table[PolygonalNumber[n], {n, nmax}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[t, x]], {n, 0, nmax}] (* Robert Price, Aug 02 2020 *)
  • PARI
    N=66; Vec(1/prod(k=1,N,1-x^(k*(k+1)\2))+O(x^N)) \\ Joerg Arndt, Apr 14 2013
    
  • Python
    from functools import lru_cache
    from sympy import divisors
    from sympy.ntheory.primetest import is_square
    @lru_cache(maxsize=None)
    def A007294(n):
        @lru_cache(maxsize=None)
        def a(n): return is_square((n<<3)+1)
        @lru_cache(maxsize=None)
        def c(n): return sum(d for d in divisors(n,generator=True) if a(d))
        return (c(n)+sum(c(k)*A007294(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024
  • Sage
    def A007294(n):
        has_nondecreasing_diffs = lambda x: min(differences(x, 2)) >= 0
        special = lambda x: (x[1]-x[0]) >= x[0]
        allowed = lambda x: (len(x) < 2 or special(x)) and (len(x) < 3 or has_nondecreasing_diffs(x))
        return len([1 for x in Partitions(n, max_slope=-1) if allowed(x[::-1])]) # D. S. McNeil, Jan 06 2011
    

Formula

G.f.: 1/Product_{k>=2} (1-z^binomial(k, 2)).
For n>0: a(n) = b(n, 1) where b(n, k) = if n>k*(k+1)/2 then b(n-k*(k+1)/2, k) + b(n, k+1) else (if n=k*(k+1)/2 then 1 else 0). - Reinhard Zumkeller, Aug 26 2003
For n>0, a(n) is Euler Transform of [1,0,1,0,0,1,0,0,0,1,0,0,0,0,1,...], i.e A010054, n>0. - Benedict W. J. Irwin, Jul 29 2016
a(n) ~ exp(3*Pi^(1/3) * Zeta(3/2)^(2/3) * n^(1/3) / 2) * Zeta(3/2) / (2^(7/2) * sqrt(3) * Pi * n^(3/2)) [Brigham 1950 (exponential part), Almkvist 2006]. - Vaclav Kotesovec, Dec 31 2016
G.f.: Sum_{i>=0} x^(i*(i+1)/2) / Product_{j=1..i} (1 - x^(j*(j+1)/2)). - Ilya Gutkovskiy, May 07 2017

Extensions

Additional comments from Roland Bacher, Jun 17 2001

A358194 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with partial sums summing to k, where k ranges from n to n(n+1)/2.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 2, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1
Offset: 0

Views

Author

Gus Wiseman, Dec 31 2022

Keywords

Comments

The partial sums of a sequence (a, b, c, ...) are (a, a+b, a+b+c, ...).

Examples

			Triangle begins:
  1
  1
  1 1
  1 0 1 1
  1 0 1 1 0 1 1
  1 0 0 1 1 0 1 1 0 1 1
  1 0 0 1 1 1 1 0 1 1 0 1 1 0 1 1
  1 0 0 0 1 1 1 0 1 1 1 1 1 0 1 1 0 1 1 0 1 1
  1 0 0 0 1 1 1 1 0 1 1 1 2 1 0 1 1 1 1 1 0 1 1 0 1 1 0 1 1
For example, the T(15,59) = 5 partitions are: (8,2,2,2,1), (7,3,3,1,1), (6,5,2,1,1), (4,3,2,2,2,2), (3,3,3,3,2,1).
		

Crossrefs

Row sums are A000041.
The version for compositions is A053632.
Row lengths are A152947.
The version for reversed partitions is A264034.
A048793 = partial sums of reversed standard compositions, sum A029931.
A358134 = partial sums of standard compositions, sum A359042.
A358136 = partial sums of prime indices, sum A318283.
A359361 = partial sums of reversed prime indices, sum A304818.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Total[Accumulate[#]]==k&]],{n,0,8},{k,n,n*(n+1)/2}]

A358137 Heinz number of the partial sums of the prime indices of n.

Original entry on oeis.org

1, 2, 3, 6, 5, 10, 7, 30, 21, 14, 11, 42, 13, 22, 33, 210, 17, 110, 19, 66, 39, 26, 23, 330, 65, 34, 273, 78, 29, 130, 31, 2310, 51, 38, 85, 546, 37, 46, 57, 390, 41, 170, 43, 102, 357, 58, 47, 2730, 133, 238, 69, 114, 53, 1870, 95, 510, 87, 62, 59, 714, 61
Offset: 1

Views

Author

Gus Wiseman, Oct 31 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The terms together with their prime indices begin:
      1: {}
      2: {1}
      3: {2}
      6: {1,2}
      5: {3}
     10: {1,3}
      7: {4}
     30: {1,2,3}
     21: {2,4}
     14: {1,4}
     11: {5}
     42: {1,2,4}
     13: {6}
     22: {1,5}
     33: {2,5}
    210: {1,2,3,4}
     17: {7}
    110: {1,3,5}
		

Crossrefs

The sorted version is A325362.
The prime indices are rows of A358136, partial sums of rows of A112798.
A000040 lists the primes.
A000041 counts partitions, strict A000009.
A003963 multiplies prime indices.
A056239 adds up prime indices.

Programs

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

Formula

A001222(a(n)) = A001222(n).

A358136 Irregular triangle read by rows whose n-th row lists the partial sums of the prime indices of n (row n of A112798).

Original entry on oeis.org

1, 2, 1, 2, 3, 1, 3, 4, 1, 2, 3, 2, 4, 1, 4, 5, 1, 2, 4, 6, 1, 5, 2, 5, 1, 2, 3, 4, 7, 1, 3, 5, 8, 1, 2, 5, 2, 6, 1, 6, 9, 1, 2, 3, 5, 3, 6, 1, 7, 2, 4, 6, 1, 2, 6, 10, 1, 3, 6, 11, 1, 2, 3, 4, 5, 2, 7, 1, 8, 3, 7, 1, 2, 4, 6, 12, 1, 9, 2, 8, 1, 2, 3, 6, 13
Offset: 2

Views

Author

Gus Wiseman, Oct 31 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.

Examples

			Triangle begins:
   2: 1
   3: 2
   4: 1 2
   5: 3
   6: 1 3
   7: 4
   8: 1 2 3
   9: 2 4
  10: 1 4
  11: 5
  12: 1 2 4
  13: 6
  14: 1 5
  15: 2 5
  16: 1 2 3 4
  17: 7
  18: 1 3 5
		

Crossrefs

Row-lengths are A001222.
First element in each row is A055396.
Last element in each row is A056239.
Rows are the partial sums of rows of A112798.
Row-sums are A318283.
Sorted Heinz numbers of the rows are A325362.
The version for standard compositions is A358134.
Rows are ranked by A358137.
A000041 counts partitions, strict A000009.
A003963 multiplies prime indices.
A056239 adds up prime indices.

Programs

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

A325367 Heinz numbers of integer partitions with distinct differences between successive parts (with the last part taken to be zero).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 10, 11, 13, 14, 15, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 82, 83, 85, 86, 87, 89, 91, 92, 93, 94, 95, 97
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 A325324.

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   4: {1,1}
   5: {3}
   7: {4}
   9: {2,2}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  15: {2,3}
  17: {7}
  19: {8}
  20: {1,1,3}
  22: {1,5}
  23: {9}
  25: {3,3}
  26: {1,6}
  28: {1,1,4}
		

Crossrefs

Programs

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

A359361 Irregular triangle read by rows whose n-th row lists the partial sums of the integer partition with Heinz number n.

Original entry on oeis.org

1, 2, 1, 2, 3, 2, 3, 4, 1, 2, 3, 2, 4, 3, 4, 5, 2, 3, 4, 6, 4, 5, 3, 5, 1, 2, 3, 4, 7, 2, 4, 5, 8, 3, 4, 5, 4, 6, 5, 6, 9, 2, 3, 4, 5, 3, 6, 6, 7, 2, 4, 6, 4, 5, 6, 10, 3, 5, 6, 11, 1, 2, 3, 4, 5, 5, 7, 7, 8, 4, 7, 2, 4, 5, 6, 12, 8, 9, 6, 8, 3, 4, 5, 6, 13
Offset: 2

Views

Author

Gus Wiseman, Dec 30 2022

Keywords

Comments

The partial sums of a sequence (a, b, c, ...) are (a, a+b, a+b+c, ...).
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The partition with Heinz number n is the reversed n-th row of A112798.

Examples

			Triangle begins:
   2: 1
   3: 2
   4: 1 2
   5: 3
   6: 2 3
   7: 4
   8: 1 2 3
   9: 2 4
  10: 3 4
  11: 5
  12: 2 3 4
  13: 6
  14: 4 5
  15: 3 5
  16: 1 2 3 4
For example, the integer partition with Heinz number 90 is (3,2,2,1), so row n = 90 is (3,5,7,8).
		

Crossrefs

Row-lengths are A001222.
The version for standard compositions is A048793, non-reversed A358134.
Last element in each row is A056239.
First element in each row is A061395
Rows are the partial sums of rows of A296150.
Row-sums are A304818.
A reverse version is A358136, row sums A318283, Heinz numbers A358137.
The sorted Heinz numbers of rows are A359397.
A000041 counts partitions, strict A000009.
A112798 lists prime indices, product A003963.
A355536 lists differences of prime indices.

Programs

  • Maple
    T:= n-> ListTools[PartialSums](sort([seq(numtheory
           [pi](i[1])$i[2], i=ifactors(n)[2])], `>`))[]:
    seq(T(n), n=2..50);  # Alois P. Heinz, Jan 01 2023
  • Mathematica
    Table[Accumulate[Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]],{n,2,30}]

A325394 Heinz numbers of integer partitions whose augmented differences are weakly increasing.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 15, 16, 17, 19, 23, 25, 27, 29, 31, 32, 35, 37, 41, 43, 47, 49, 53, 55, 59, 61, 64, 67, 71, 73, 75, 77, 79, 81, 83, 89, 91, 97, 101, 103, 105, 107, 109, 113, 119, 121, 125, 127, 128, 131, 137, 139, 143, 149, 151, 157, 163, 167
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 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 A325356.

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}
    9: {2,2}
   11: {5}
   13: {6}
   15: {2,3}
   16: {1,1,1,1}
   17: {7}
   19: {8}
   23: {9}
   25: {3,3}
   27: {2,2,2}
   29: {10}
   31: {11}
   32: {1,1,1,1,1}
		

Crossrefs

Programs

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

A325388 Heinz numbers of strict integer partitions with distinct differences (with the last part taken to be 0).

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 11, 13, 14, 15, 17, 19, 22, 23, 26, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 67, 69, 71, 73, 74, 77, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 115, 118, 119, 122
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 (6,3,1) (with the last part taken to be 0) are (-3,-2,-1).
The enumeration of these partitions by sum is given by A320348.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    3: {2}
    5: {3}
    7: {4}
   10: {1,3}
   11: {5}
   13: {6}
   14: {1,4}
   15: {2,3}
   17: {7}
   19: {8}
   22: {1,5}
   23: {9}
   26: {1,6}
   29: {10}
   31: {11}
   33: {2,5}
   34: {1,7}
   35: {3,4}
		

Crossrefs

Programs

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

A261079 Sum of index differences between prime factors of n, summed over all unordered pairs of primes present (with multiplicity) in the prime factorization of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 2, 0, 3, 1, 0, 0, 2, 0, 4, 2, 4, 0, 3, 0, 5, 0, 6, 0, 4, 0, 0, 3, 6, 1, 4, 0, 7, 4, 6, 0, 6, 0, 8, 2, 8, 0, 4, 0, 4, 5, 10, 0, 3, 2, 9, 6, 9, 0, 7, 0, 10, 4, 0, 3, 8, 0, 12, 7, 6, 0, 6, 0, 11, 2, 14, 1, 10, 0, 8, 0, 12, 0, 10, 4, 13, 8, 12, 0, 6, 2, 16, 9, 14, 5, 5, 0, 6, 6, 8, 0, 12, 0, 15, 4, 15, 0, 6, 0, 8, 10, 12, 0, 14, 6, 18, 8, 16, 3, 10
Offset: 1

Views

Author

Antti Karttunen, Sep 23 2015

Keywords

Examples

			For n = 1 the prime factorization is empty, thus there is nothing to sum, so a(1) = 0.
For n = 6 = 2*3 = prime(1) * prime(2), a(6) = 1 because the (absolute value of) difference between prime indices of 2 and 3 is 1.
For n = 10 = 2*5 = prime(1) * prime(3), a(10) = 2 because the difference between prime indices of 2 and 5 is 2.
For n = 12 = 2*2*3 = prime(1) * prime(1) * prime(2), a(12) = 2 because the difference between prime indices of 2 and 3 is 1, and the pair (2,3) occurs twice as one can pick either one of the two 2's present in the prime factorization to be a pair of a single 3. Note that the index difference between 2 and 2 is 0, thus the pair (2,2) of prime divisors does not contribute to the sum.
For n = 36 = 2*2*3*3, a(36) = 4 because the index difference between 2 and 3 is 1, and the prime factor pair (2,3) occurs 2^2 = four times in total. As the index difference is zero between 2 and 2 as well as between 3 and 3, the pairs (2,2) and (3,3) do not contribute to the sum.
		

Crossrefs

Cf. A000720.
Cf. A000961 (positions of zeros), A006094 (positions of ones).
Cf. also A260737.
A055396 gives minimum prime index, maximum A061395.
A112798 list prime indices, length A001222, sum A056239.
A304818 adds up partial sums of reversed prime indices, row sums of A359361.
A318283 adds up partial sums of prime indices, row sums of A358136.

Programs

  • Mathematica
    Table[Function[p, Total@ Map[Function[b, Times @@ {First@ Differences@ PrimePi@ b, Count[Subsets[p, {2}], c_ /; SameQ[c, b]]}], Subsets[Union@ p, {2}]]][Flatten@ Replace[FactorInteger@ n, {p_, e_} :> ConstantArray[p, e], 2]], {n, 120}] (* Michael De Vlieger, Mar 08 2017 *)

Formula

a(n) = A304818(n) - A318283(n). - Gus Wiseman, Jan 09 2023
a(n) = 2*A304818(n) - A359362(n). - Gus Wiseman, Jan 09 2023

A325390 Heinz number of the negated differences plus one of the integer partition with Heinz number n (with the last part taken to be 0).

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 11, 12, 10, 15, 13, 18, 17, 21, 15, 24, 19, 18, 23, 30, 25, 33, 29, 36, 14, 39, 20, 42, 31, 27, 37, 48, 35, 51, 21, 36, 41, 57, 55, 60, 43, 45, 47, 66, 30, 69, 53, 72, 22, 30, 65, 78, 59, 36, 35, 84, 85, 87, 61, 54, 67, 93, 50, 96, 49, 63, 71
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The Heinz number of a positive integer sequence (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 (6,3,1) (with the last part taken to be 0) are (-3,-2,-1).

Examples

			The Heinz number of (6,3,1) is 130, and its negated differences plus one are (4,3,2), which has Heinz number 105, so a(130) = 105.
		

Crossrefs

Number of appearances of n is A325392(n).
Positions of squarefree numbers are A325367.

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Times@@Prime/@(1-Differences[Append[primeptn[n],0]]),{n,100}]
Showing 1-10 of 18 results. Next