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 91-100 of 340 results. Next

A317089 Numbers whose prime factors span an initial interval of prime numbers and whose prime multiplicities span an initial interval of positive integers.

Original entry on oeis.org

2, 6, 12, 18, 30, 60, 90, 150, 180, 210, 300, 360, 420, 450, 540, 600, 630, 1050, 1260, 1350, 1470, 1500, 2100, 2250, 2310, 2520, 2940, 3150, 3780, 4200, 4410, 4620, 5880, 6300, 6930, 7350, 8820, 9450, 10500, 11550, 12600, 13230, 13860, 14700, 15750, 16170
Offset: 1

Views

Author

Gus Wiseman, Jul 21 2018

Keywords

Examples

			The sequence of rows of A296150 indexed by the terms of this sequence begins: (1), (21), (211), (221), (321), (3211), (3221), (3321), (32211), (4321), (33211), (322111), (43211).
		

Crossrefs

Programs

  • Mathematica
    normalQ[m_]:=Union[m]==Range[Max[m]];
    Select[Range[10000],And[normalQ[PrimePi/@FactorInteger[#][[All,1]]],normalQ[FactorInteger[#][[All,2]]]]&]
  • PARI
    ok(n)={my(f=factor(n), p=f[,1], e=vecsort(f[,2],,8)); n > 1 && #p==primepi(p[#p]) && #e==e[#e]} \\ Andrew Howroyd, Aug 26 2018

A321648 Number of permutations of the conjugate of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 4, 3, 1, 1, 2, 1, 3, 6, 5, 1, 2, 1, 6, 1, 4, 1, 6, 1, 1, 10, 7, 4, 2, 1, 8, 15, 3, 1, 12, 1, 5, 3, 9, 1, 2, 1, 3, 21, 6, 1, 2, 10, 4, 28, 10, 1, 6, 1, 11, 6, 1, 20, 20, 1, 7, 36, 12, 1, 2, 1, 12, 3, 8, 5, 30, 1, 3, 1, 13
Offset: 1

Views

Author

Gus Wiseman, Nov 15 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The a(42) = 12 permutations: (3211), (3121), (3112), (2311), (2131), (2113), (1321), (1312), (1231), (1213), (1132), (1123).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Table[Length[Permutations[conj[primeMS[n]]]],{n,50}]
  • PARI
    A008480(n) = {my(sig=factor(n)[, 2]); vecsum(sig)!/factorback(apply(k->k!, sig))}; \\ From A008480
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A122111(n) = if(1==n,n,prime(bigomega(n))*A122111(A064989(n)));
    A321648(n) = A008480(A122111(n)); \\ Antti Karttunen, Dec 23 2018

Formula

a(n) = A008480(A122111(n)).

A334442 Irregular triangle whose reversed rows are all integer partitions sorted first by sum, then by length, and finally reverse-lexicographically.

Original entry on oeis.org

1, 2, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 3, 2, 2, 1, 1, 2, 1, 1, 1, 1, 5, 1, 4, 2, 3, 1, 1, 3, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 5, 2, 4, 3, 3, 1, 1, 4, 1, 2, 3, 2, 2, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 7, 1, 6, 2, 5, 3, 4, 1, 1, 5
Offset: 0

Views

Author

Gus Wiseman, May 07 2020

Keywords

Comments

First differs from A036036 for reversed partitions of 9. Namely, this sequence has (2,2,5) before (1,4,4), while A036036 has (1,4,4) before (2,2,5).

Examples

			The sequence of all partitions begins:
  ()         (2,3)        (1,1,1,1,2)    (1,1,1,2,2)
  (1)        (1,1,3)      (1,1,1,1,1,1)  (1,1,1,1,1,2)
  (2)        (1,2,2)      (7)            (1,1,1,1,1,1,1)
  (1,1)      (1,1,1,2)    (1,6)          (8)
  (3)        (1,1,1,1,1)  (2,5)          (1,7)
  (1,2)      (6)          (3,4)          (2,6)
  (1,1,1)    (1,5)        (1,1,5)        (3,5)
  (4)        (2,4)        (1,2,4)        (4,4)
  (1,3)      (3,3)        (1,3,3)        (1,1,6)
  (2,2)      (1,1,4)      (2,2,3)        (1,2,5)
  (1,1,2)    (1,2,3)      (1,1,1,4)      (1,3,4)
  (1,1,1,1)  (2,2,2)      (1,1,2,3)      (2,2,4)
  (5)        (1,1,1,3)    (1,2,2,2)      (2,3,3)
  (1,4)      (1,1,2,2)    (1,1,1,1,3)    (1,1,1,5)
This sequence can also be interpreted as the following triangle:
                  0
                 (1)
               (2)(11)
             (3)(12)(111)
        (4)(13)(22)(112)(1111)
  (5)(14)(23)(113)(122)(1112)(11111)
Taking Heinz numbers (A334438) gives:
   1
   2
   3   4
   5   6   8
   7  10   9  12  16
  11  14  15  20  18  24  32
  13  22  21  25  28  30  27  40  36  48  64
  17  26  33  35  44  42  50  45  56  60  54  80  72  96 128
		

Crossrefs

Row lengths are A036043.
The version for reversed partitions is A334301.
The version for colex instead of revlex is A334302.
Taking Heinz numbers gives A334438.
The version with rows reversed is A334439.
Ignoring length gives A335122.
Lexicographically ordered reversed partitions are A026791.
Reversed partitions in Abramowitz-Stegun (sum/length/lex) order are A036036.
Partitions in increasing-length colex order (sum/length/colex) are A036037.
Reverse-lexicographically ordered partitions are A080577.
Lexicographically ordered partitions are A193073.
Partitions in colexicographic order (sum/colex) are A211992.
Sorting partitions by Heinz number gives A296150.

Programs

  • Mathematica
    revlensort[f_,c_]:=If[Length[f]!=Length[c],Length[f]
    				
  • PARI
    A334442_row(n)=vecsort(partitions(n),p->concat(#p,-Vecrev(p))) \\ Rows of triangle defined in EXAMPLE (all partitions of n). Wrap into [Vec(p)|p<-...] to avoid "Vecsmall". - M. F. Hasler, May 14 2020

A352872 Numbers whose weakly increasing prime indices y have a fixed point y(i) = i.

Original entry on oeis.org

2, 4, 6, 8, 9, 10, 12, 14, 16, 18, 20, 22, 24, 26, 27, 28, 30, 32, 34, 36, 38, 40, 42, 44, 45, 46, 48, 50, 52, 54, 56, 58, 60, 62, 63, 64, 66, 68, 70, 72, 74, 75, 76, 78, 80, 81, 82, 84, 86, 88, 90, 92, 94, 96, 98, 99, 100, 102, 104, 106, 108, 110, 112, 114
Offset: 1

Views

Author

Gus Wiseman, Apr 06 2022

Keywords

Comments

First differs from A118672 in having 75.
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

			The terms together with their prime indices begin:
      2: {1}           28: {1,1,4}         56: {1,1,1,4}
      4: {1,1}         30: {1,2,3}         58: {1,10}
      6: {1,2}         32: {1,1,1,1,1}     60: {1,1,2,3}
      8: {1,1,1}       34: {1,7}           62: {1,11}
      9: {2,2}         36: {1,1,2,2}       63: {2,2,4}
     10: {1,3}         38: {1,8}           64: {1,1,1,1,1,1}
     12: {1,1,2}       40: {1,1,1,3}       66: {1,2,5}
     14: {1,4}         42: {1,2,4}         68: {1,1,7}
     16: {1,1,1,1}     44: {1,1,5}         70: {1,3,4}
     18: {1,2,2}       45: {2,2,3}         72: {1,1,1,2,2}
     20: {1,1,3}       46: {1,9}           74: {1,12}
     22: {1,5}         48: {1,1,1,1,2}     75: {2,3,3}
     24: {1,1,1,2}     50: {1,3,3}         76: {1,1,8}
     26: {1,6}         52: {1,1,6}         78: {1,2,6}
     27: {2,2,2}       54: {1,2,2,2}       80: {1,1,1,1,3}
For example, the multiset {2,3,3} with Heinz number 75 has a fixed point at position 3, so 75 is in the sequence.
		

Crossrefs

* = unproved
These partitions are counted by A238395, strict A096765.
These are the nonzero positions in A352822.
*The complement reverse version is A352826, counted by A064428.
*The reverse version is A352827, counted by A001522 (strict A352829).
The complement is A352830, counted by A238394 (strict A025147).
A000700 counts self-conjugate partitions, ranked by A088902.
A001222 counts prime indices, distinct A001221.
A008290 counts permutations by fixed points, nonfixed A098825.
A056239 adds up prime indices, row sums of A112798 and A296150.
A114088 counts partitions by excedances.
A115720 and A115994 count partitions by their Durfee square.
A122111 represents partition conjugation using Heinz numbers.
A124010 gives prime signature, sorted A118914, conjugate rank A238745.
A238349 counts compositions by fixed points, complement A352523.
A238352 counts reversed partitions by fixed points.
A352833 counts partitions by fixed points.

Programs

  • Mathematica
    pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]];
    Select[Range[100],pq[Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]>0&]

A353867 Heinz numbers of integer partitions where every partial run (consecutive constant subsequence) has a different sum, and these sums include every integer from 0 to the greatest part.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 20, 30, 32, 56, 64, 90, 128, 140, 176, 210, 256, 416, 512, 616, 990, 1024, 1088, 1540, 2048, 2288, 2310, 2432, 2970, 4096, 4950, 5888, 7072, 7700, 8008, 8192, 11550, 12870, 14848, 16384, 20020, 20672, 30030, 31744, 32768, 38896, 50490, 55936
Offset: 1

Views

Author

Gus Wiseman, Jun 07 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.
Related concepts:
- A partition whose submultiset sums cover an initial interval is said to be complete (A126796, ranked by A325781).
- In a knapsack partition (A108917, ranked by A299702), every submultiset has a different sum.
- A complete partition that is also knapsack is said to be perfect (A002033, ranked by A325780).
- A partition whose partial runs have all different sums is said to be rucksack (A353864, ranked by A353866, complement A354583).

Examples

			The terms together with their prime indices begin:
    1: {}
    2: {1}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   16: {1,1,1,1}
   20: {1,1,3}
   30: {1,2,3}
   32: {1,1,1,1,1}
   56: {1,1,1,4}
   64: {1,1,1,1,1,1}
   90: {1,2,2,3}
  128: {1,1,1,1,1,1,1}
  140: {1,1,3,4}
  176: {1,1,1,1,5}
  210: {1,2,3,4}
  256: {1,1,1,1,1,1,1,1}
		

Crossrefs

Knapsack partitions are counted by A108917, ranked by A299702.
Complete partitions are counted by A126796, ranked by A325781.
These partitions are counted by A353865.
This is a special case of A353866, counted by A353864, complement A354583.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A073093 counts prime-power divisors.
A124010 gives prime signature, sorted A118914.
A300273 ranks collapsible partitions, counted by A275870.
A353832 represents the operation of taking run-sums of a partition.
A353833 ranks partitions with all equal run-sums, nonprime A353834.
A353836 counts partitions by number of distinct run-sums.
A353852 ranks compositions with all distinct run-sums, counted by A353850.
A353863 counts partitions whose weak run-sums cover an initial interval.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    norqQ[m_]:=Sort[m]==Range[0,Max[m]];
    msubs[s_]:=Join@@@Tuples[Table[Take[t,i],{t,Split[s]},{i,0,Length[t]}]];
    Select[Range[1000],norqQ[Total/@Select[msubs[primeMS[#]],SameQ@@#&]]&]

A372436 Numbers whose binary indices and prime indices have the same maximum.

Original entry on oeis.org

3, 5, 14, 22, 39, 52, 68, 85, 102, 119, 133, 152, 171, 190, 209, 228, 247, 276, 299, 322, 345, 368, 391, 414, 437, 460, 483, 506, 522, 551, 580, 609, 638, 667, 696, 725, 754, 783, 812, 841, 870, 928, 957, 986, 1015, 1054, 1085, 1116, 1178, 1209, 1240, 1302
Offset: 1

Views

Author

Gus Wiseman, May 04 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
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.
Note that a number's binary and prime indices cannot have the same minimum; see A372437.

Examples

			The binary indices of 345 are {1,4,5,7,9}, and the prime indices are {2,3,9}. Both have maximum 9, so 345 is in the sequence.
The terms together with their prime indices begin:
     3: {2}
     5: {3}
    14: {1,4}
    22: {1,5}
    39: {2,6}
    52: {1,1,6}
    68: {1,1,7}
    85: {3,7}
   102: {1,2,7}
   119: {4,7}
   133: {4,8}
   152: {1,1,1,8}
   171: {2,2,8}
The terms together with their binary expansions and binary indices begin:
     3:           11 ~ {1,2}
     5:          101 ~ {1,3}
    14:         1110 ~ {2,3,4}
    22:        10110 ~ {2,3,5}
    39:       100111 ~ {1,2,3,6}
    52:       110100 ~ {3,5,6}
    68:      1000100 ~ {3,7}
    85:      1010101 ~ {1,3,5,7}
   102:      1100110 ~ {2,3,6,7}
   119:      1110111 ~ {1,2,3,5,6,7}
   133:     10000101 ~ {1,3,8}
   152:     10011000 ~ {4,5,8}
   171:     10101011 ~ {1,2,4,6,8}
		

Crossrefs

For length instead of maximum we have A071814.
For sum instead of maximum we have A372427.
Positions of zeros in A372442, for minimum instead of maximum A372437.
A003963 gives product of prime indices.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.

Programs

  • Mathematica
    bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Max[prix[#]]==Max[bix[#]]&]

Formula

A070939(a(n)) = A061395(a(n)).

A375128 Irregular triangle read by rows where row n lists the minima of maximal strictly increasing runs in the weakly increasing prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 04 2024

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 minima of strictly increasing runs in a sequence are obtained by splitting it into maximal strictly increasing subsequences and taking the first term of each.

Examples

			The prime indices of 540 are {1,1,2,2,2,3}, with strictly increasing runs ({1},{1,2},{2},{2,3}), with minima (1,1,2,2), which is row 540.
Triangle begins:
   1:
   2:  1
   3:  2
   4:  1  1
   5:  3
   6:  1
   7:  4
   8:  1  1  1
   9:  2  2
  10:  1
  11:  5
  12:  1  1
  13:  6
  14:  1
  15:  2
  16:  1  1  1  1
		

Crossrefs

Row-minima are A055396.
Row-sums are A374706.
Row-lengths are A375136.
For leaders of constant runs we have A304038, row-sums A066328.
For compositions we have A374683, row-sums of A374684 (length A124768).
A112798 lists prime indices:
- length A001222, distinct A001221
- leader A055396
- sum A056239
- reverse A296150

Programs

  • Mathematica
    Table[If[n==1,{},First/@Split[Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]],Less]],{n,100}]

A316314 Number of distinct nonempty-subset-averages of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 4, 1, 3, 3, 1, 1, 4, 1, 4, 3, 3, 1, 5, 1, 3, 1, 4, 1, 5, 1, 1, 3, 3, 3, 5, 1, 3, 3, 5, 1, 7, 1, 4, 4, 3, 1, 6, 1, 4, 3, 4, 1, 5, 3, 5, 3, 3, 1, 8, 1, 3, 4, 1, 3, 7, 1, 4, 3, 7, 1, 7, 1, 3, 4, 4, 3, 7, 1, 6, 1, 3, 1, 8, 3, 3, 3, 5, 1, 7, 3, 4, 3, 3, 3, 7, 1, 4, 4, 5, 1, 7, 1, 5, 5
Offset: 1

Views

Author

Gus Wiseman, Jun 29 2018

Keywords

Comments

A rational number q is a nonempty-subset-average of an integer partition y if there exists a nonempty submultiset of y with average q.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The a(42) = 7 subset-averages of (4,2,1) are 1, 3/2, 2, 7/3, 5/2, 3, 4.
The a(72) = 7 subset-averages of (2,2,1,1,1) are 1, 5/4, 4/3, 7/5, 3/2, 5/3, 2.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Union[Mean/@Rest[Subsets[primeMS[n]]]]],{n,100}]
  • PARI
    up_to = 65537;
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
    v056239 = vector(up_to,n,A056239(n));
    A316314(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((d>1)&&!mapisdefined(m,s = v056239[d]/bigomega(d)), mapput(m,s,s); k++)); (k); }; \\ Antti Karttunen, Sep 23 2018

Formula

a(n) = A316398(n) - 1.

Extensions

More terms from Antti Karttunen, Sep 23 2018

A372428 Sum of binary indices of n minus sum of prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 2, 3, 2, 4, 5, 1, -1, 2, 0, 3, 3, 4, 2, 4, 4, 4, 6, 6, 3, 8, 4, 1, 0, 0, 2, 3, -2, 2, 4, 4, -2, 5, -1, 6, 7, 5, 1, 5, 4, 6, 5, 6, -1, 9, 9, 8, 6, 6, 1, 11, 1, 8, 13, 1, -1, 1, -9, 1, 0, 4, -7, 4, -9, 0, 6, 4, 6, 7, -5, 5, 5, 0, -8
Offset: 1

Views

Author

Gus Wiseman, May 02 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
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

			The binary indices of 65 are {1,7}, and the prime indices are {3,6}, so a(65) = 8 - 9 = -1.
		

Crossrefs

Positions of zeros are A372427.
For minimum instead of sum we have A372437.
For length instead of sum we have A372441, zeros A071814.
For maximum instead of sum we have A372442, zeros A372436.
Positions of odd terms are A372586, even A372587.
A003963 gives product of prime indices.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
A326031 gives weight of the set-system with BII-number n.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Table[Total[bix[n]]-Total[prix[n]],{n,100}]
  • Python
    from itertools import count, islice
    from sympy import sieve, factorint
    def a_gen():
        for n in count(1):
            b = sum((i+1) for i, x in enumerate(bin(n)[2:][::-1]) if x =='1')
            p = sum(sieve.search(i)[0] for i in factorint(n, multiple=True))
            yield(b-p)
    A372428_list = list(islice(a_gen(), 83)) # John Tyler Rascoe, May 04 2024
    
  • Python
    from sympy import primepi, factorint
    def A372428(n): return int(sum(i for i, j in enumerate(bin(n)[:1:-1],1) if j=='1')-sum(primepi(p)*e for p, e in factorint(n).items())) # Chai Wah Wu, Oct 18 2024

Formula

a(n) = A029931(n) - A056239(n).

A316431 Least common multiple divided by greatest common divisor of the integer partition with Heinz number n > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 4, 6, 1, 1, 2, 1, 3, 2, 5, 1, 2, 1, 6, 1, 4, 1, 6, 1, 1, 10, 7, 12, 2, 1, 8, 3, 3, 1, 4, 1, 5, 6, 9, 1, 2, 1, 3, 14, 6, 1, 2, 15, 4, 4, 10, 1, 6, 1, 11, 2, 1, 2, 10, 1, 7, 18, 12, 1, 2, 1, 12, 6, 8, 20, 6, 1, 3, 1, 13, 1, 4, 21, 14, 5, 5, 1, 6, 6, 9, 22, 15, 24, 2, 1, 4, 10, 3, 1, 14, 1, 6, 12
Offset: 2

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			63 is the Heinz number of (4,2,2), which has LCM 4 and GCD 2, so a(63) = 4/2 = 2.
91 is the Heinz number of (6,4), which has LCM 12 and GCD 2, so a(91) = 12/2 = 6.
		

Crossrefs

Programs

  • Mathematica
    Table[With[{pms=Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]},LCM@@pms/GCD@@pms],{n,2,100}]
  • PARI
    A316431(n) = if(1==n,1,my(pis = apply(p -> primepi(p), factor(n)[, 1]~)); lcm(pis)/gcd(pis)); \\ Antti Karttunen, Sep 06 2018

Formula

a(n) = A290103(n)/A289508(n).
a(n) = a(A005117(n)). - David A. Corneth, Sep 06 2018

Extensions

More terms from Antti Karttunen, Sep 06 2018
Previous Showing 91-100 of 340 results. Next