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 21-30 of 65 results. Next

A324922 a(n) = unique m such that m/A003963(m) = n, where A003963 is product of prime indices.

Original entry on oeis.org

1, 2, 6, 4, 30, 12, 28, 8, 36, 60, 330, 24, 156, 56, 180, 16, 476, 72, 152, 120, 168, 660, 828, 48, 900, 312, 216, 112, 1740, 360, 10230, 32, 1980, 952, 840, 144, 888, 304, 936, 240, 6396, 336, 2408, 1320, 1080, 1656, 8460, 96, 784, 1800, 2856, 624, 848, 432
Offset: 1

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

Every positive integer has a unique factorization into factors q(i) = prime(i)/i, i > 0 given by the rows of A324924. Then a(n) is the number obtained by encoding this factorization as a standard factorization into prime numbers (A112798).

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    difac[n_]:=If[n==1,{},With[{m=Product[Prime[i]/i,{i,primeMS[n]}]},Sort[Join[primeMS[n],difac[n/m]]]]];
    Table[Times@@Prime/@difac[n],{n,30}]
  • PARI
    a(n) = my (f=factor(n)); prod (i=1, #f~, (f[i,1] * a(primepi(f[i,1])))^f[i,2]) \\ Rémy Sigrist, Jul 18 2019

Formula

a(n) = Product_t mg(t) where the product is over all (not necessarily distinct) terminal subtrees of the rooted tree with Matula-Goebel number n, and mg(t) is the Matula-Goebel number of t.
Completely multiplicative with a(prime(n)) = prime(n) * a(n). - Rémy Sigrist, Jul 18 2019

Extensions

Keyword mult added by Rémy Sigrist, Jul 18 2019

A324849 Positive integers divisible by none of their prime indices > 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 16, 17, 19, 20, 21, 22, 23, 25, 26, 27, 29, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 43, 44, 46, 47, 49, 50, 51, 52, 53, 57, 58, 59, 61, 62, 63, 64, 65, 67, 68, 69, 70, 71, 73, 74, 76, 77, 79, 80, 81, 82, 83, 85, 86, 87
Offset: 1

Views

Author

Gus Wiseman, Mar 18 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

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}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  16: {1,1,1,1}
  17: {7}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  23: {9}
  25: {3,3}
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) andmap(t -> not ((n/numtheory:-pi(t))::integer), numtheory:-factorset(n) minus {2}) end proc:
    select(filter, [$1..200]); # Robert Israel, Mar 20 2019
  • Mathematica
    Select[Range[100],!Or@@Cases[If[#==1,{},FactorInteger[#]],{p_,_}:>If[p==2,False,Divisible[#,PrimePi[p]]]]&]
  • PARI
    is(n) = my(f=factor(n)[, 1]~, idc=[]); for(k=1, #f, idc=concat(idc, [primepi(f[k])])); for(t=1, #idc, if(idc[t]==1, next); if(n%idc[t]==0, return(0))); 1 \\ Felix Fröhlich, Mar 21 2019

A324847 Numbers divisible by at least one of their prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 18 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.
If n is in the sequence, then so are all multiples of n. - Robert Israel, Mar 19 2019

Examples

			The sequence of terms together with their prime indices begins:
   2: {1}
   4: {1,1}
   6: {1,2}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  16: {1,1,1,1}
  18: {1,2,2}
  20: {1,1,3}
  22: {1,5}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
  30: {1,2,3}
  32: {1,1,1,1,1}
  34: {1,7}
  36: {1,1,2,2}
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local F;
      F:= map(numtheory:-pi, numtheory:-factorset(n));
      ormap(t -> n mod t = 0, F);
    end proc:
    select(filter, [$1..200]); # Robert Israel, Mar 19 2019
  • Mathematica
    Select[Range[100],Or@@Cases[If[#==1,{},FactorInteger[#]],{p_,_}:>Divisible[#,PrimePi[p]]]&]
  • PARI
    isok(n) = {my(f = factor(n)[,1]); for (k=1, #f, if (!(n % primepi(f[k])), return (1));); return (0);} \\ Michel Marcus, Mar 19 2019

A330953 Number of integer partitions of n whose Heinz number (product of primes of parts) is divisible by their sum of primes of parts.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 3, 4, 6, 3, 12, 10, 12, 14, 27, 38, 44, 52, 48, 77, 101, 106, 127, 206, 268, 377, 392, 496, 602, 671, 821, 1090, 1318, 1568, 1926, 2260, 2703, 3258, 3942, 4858, 5923, 6891, 8286, 9728, 11676, 13775, 16314, 19749, 23474, 27793, 32989, 38775
Offset: 1

Views

Author

Gus Wiseman, Jan 15 2020

Keywords

Comments

The Heinz number of an integer 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 a(1) = 1 through a(11) = 12 partitions: (A = 10, B = 11):
  1  2   3  4     5  6    7      8         9        A         B
     11     1111     222  3211   431       432      5311      542
                     321  22111  4211      3321     22111111  5411
                                 11111111  32211              33221
                                           321111             42221
                                           2211111            53111
                                                              322211
                                                              431111
                                                              521111
                                                              2222111
                                                              3311111
                                                              32111111
For example, the partition (3,3,2,2,1) is counted under a(11) because 5*5*3*3*2 = 450 is divisible by 5+5+3+3+2 = 18.
		

Crossrefs

The Heinz numbers of these partitions are given by A036844.
Numbers divisible by the sum of their prime indices are A324851.
Partitions whose product is divisible by their sum are A057568.
Partitions whose Heinz number is divisible by all parts are A330952.
Partitions whose Heinz number is divisible by their product are A324925.
Partitions whose Heinz number is divisible by their sum are A330950.
Partitions whose product is divisible by their sum of primes are A330954.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Divisible[Times@@Prime/@#,Plus@@Prime/@#]&]],{n,30}]

A324848 Number of prime indices of n (counted with multiplicity) that divide n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 18 2019

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

			The prime indices of 6776 are {1,1,1,4,5,5}, four of which {1,1,1,4} divide 6776, so a(6776) = 4.
		

Crossrefs

The version for distinct prime indices is A324852.
Positions of zeros are A324846.
Positions of ones are A324856.

Programs

  • Mathematica
    Table[Total[Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>k/;Divisible[n,PrimePi[p]]]],{n,100}]

A324931 Integers in the list of quotients of positive integers by their product of prime indices.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 12, 7, 5, 32, 9, 24, 14, 10, 64, 18, 48, 28, 20, 128, 36, 19, 13, 21, 15, 96, 27, 56, 40, 256, 72, 38, 26, 11, 42, 30, 192, 54, 112, 17, 80, 512, 144, 76, 52, 22, 84, 60, 384, 49, 23, 35, 53, 108, 37, 224, 25, 57, 39, 34, 160, 63, 1024
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2019

Keywords

Comments

These quotients are given by A324932(n)/A324933(n).
This is a permutation of the positive integers, with inverse A324934.

Examples

			The sequence of quotients n/A003963(n) begins: 1, 2, 3/2, 4, 5/3, 3, 7/4, 8, 9/4, 10/3, 11/5, 6, 13/6, 7/2, 5/2, 16, ...
		

Crossrefs

Programs

  • Mathematica
    Select[Table[n/Times@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>PrimePi[p]^k],{n,100}],IntegerQ]

Formula

a(n) = A324850(n)/A003963(A324850(n)).

A355735 Number of ways to choose a divisor of each prime index of n (taken in weakly increasing order) such that the result is weakly increasing.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 16 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

			The a(15) = 3 ways are: (1,1), (1,3), (2,3).
The a(18) = 3 ways are: (1,1,1), (1,1,2), (1,2,2).
The a(2) = 1 through a(19) = 4 ways:
  1  1  11  1  11  1  111  11  11  1  111  1  11  11  1111  1  111  1
     2      3  12  2       12  13  5  112  2  12  13        7  112  2
                   4       22              3  14  23           122  4
                                           6                        8
		

Crossrefs

Allowing any choice of divisors gives A355731, firsts A355732.
Choosing a multiset instead of sequence gives A355733, firsts A355734.
Positions of first appearances are A355736.
Choosing only prime divisors gives A355745, variations A355741, A355744.
The reverse version is A355749.
A000005 counts divisors.
A001414 adds up distinct prime divisors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061395 selects the maximum prime index.
A120383 lists numbers divisible by all of their prime indices.
A324850 lists numbers divisible by the product of their prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Tuples[Divisors/@primeMS[n]],LessEqual@@#&]],{n,100}]

A352487 Excedance set of A122111. Numbers k < A122111(k), where A122111 represents partition conjugation using Heinz numbers.

Original entry on oeis.org

3, 5, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93, 94
Offset: 1

Views

Author

Gus Wiseman, Mar 19 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The sequence lists all Heinz numbers of partitions whose Heinz number is less than that of their conjugate.

Examples

			The terms together with their prime indices begin:
   3: (2)
   5: (3)
   7: (4)
  10: (3,1)
  11: (5)
  13: (6)
  14: (4,1)
  15: (3,2)
  17: (7)
  19: (8)
  21: (4,2)
  22: (5,1)
  23: (9)
  25: (3,3)
  26: (6,1)
  28: (4,1,1)
For example, the partition (4,1,1) has Heinz number 28 and its conjugate (3,1,1,1) has Heinz number 40, and 28 < 40, so 28 is in the sequence, and 40 is not.
		

Crossrefs

These partitions are counted by A000701.
The weak version is A352489, counted by A046682.
The opposite version is A352490, weak A352488.
These are the positions of negative terms in A352491.
A000041 counts integer partitions, strict A000009.
A000700 counts self-conjugate partitions, ranked by A088902 (cf. A258116).
A003963 = product of prime indices, conjugate A329382.
A008292 is the triangle of Eulerian numbers (version without zeros).
A008480 counts permutations of prime indices, conjugate A321648.
A056239 adds up prime indices, row sums of A112798 and A296150.
A122111 = partition conjugation using Heinz numbers, parts A321649/A321650.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A173018 counts permutations by excedances, weak A123125.
A238744 = partition conjugate of prime signature, ranked by A238745.
A330644 counts non-self-conjugate partitions, ranked by A352486.
A352521 counts compositions by subdiagonals, rank statistic A352514.

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]}]];
    Select[Range[100],#
    				

Formula

a(n) < A122111(a(n)).

A355734 Least k such that there are exactly n multisets that can be obtained by choosing a divisor of each prime index of k.

Original entry on oeis.org

1, 3, 7, 13, 21, 35, 39, 89, 133, 105, 91, 195, 351, 285, 247, 333, 273, 481, 455, 555, 623, 801, 791, 741, 1359, 1157, 1281, 1335, 1365, 1443, 1977, 1729, 1967, 1869, 2109, 3185, 2373, 2769, 2639, 4361, 3367, 3653, 3885, 3471, 4613, 5883, 5187, 5551, 6327
Offset: 1

Views

Author

Gus Wiseman, Jul 21 2022

Keywords

Comments

This is the position of first appearance of n in A355733.
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:
    1: {}
    3: {2}
    7: {4}
   13: {6}
   21: {2,4}
   35: {3,4}
   39: {2,6}
   89: {24}
  133: {4,8}
  105: {2,3,4}
   91: {4,6}
  195: {2,3,6}
  351: {2,2,2,6}
For example, the choices for a(12) = 195 are:
  {1,1,1}  {1,2,2}  {1,3,6}
  {1,1,2}  {1,2,3}  {2,2,3}
  {1,1,3}  {1,2,6}  {2,3,3}
  {1,1,6}  {1,3,3}  {2,3,6}
		

Crossrefs

Counting all choices of divisors gives A355732, firsts of A355731.
Positions of first appearances in A355733.
Choosing weakly increasing divisors gives A355736, firsts of A355735.
A000005 counts divisors.
A001414 adds up distinct prime divisors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A120383 lists numbers divisible by all of their prime indices.
A324850 lists numbers divisible by the product of their prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    mnrm[s_]:=If[Min@@s==1,mnrm[DeleteCases[s-1,0]]+1,0];
    az=Table[Length[Union[Sort/@Tuples[Divisors/@primeMS[n]]]],{n,1000}];
    Table[Position[az,k][[1,1]],{k,mnrm[az]}]

A355535 Odd numbers of which it is not possible to choose a different prime factor of each prime index.

Original entry on oeis.org

9, 21, 25, 27, 45, 49, 57, 63, 75, 81, 99, 105, 115, 117, 121, 125, 133, 135, 147, 153, 159, 171, 175, 189, 195, 207, 225, 231, 243, 245, 261, 273, 275, 279, 285, 289, 297, 315, 325, 333, 343, 345, 351, 357, 361, 363, 369, 371, 375, 387, 393, 399, 405, 423
Offset: 1

Views

Author

Gus Wiseman, Jul 22 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

			The terms together with their prime indices begin:
    9: {2,2}
   21: {2,4}
   25: {3,3}
   27: {2,2,2}
   45: {2,2,3}
   49: {4,4}
   57: {2,8}
   63: {2,2,4}
   75: {2,3,3}
   81: {2,2,2,2}
   99: {2,2,5}
  105: {2,3,4}
For example, the prime indices of 897 are {2,6,9}, of which we can choose prime factors in two ways: (2,2,3) or (2,3,3); but neither of these has all distinct elements, so 897 is in the sequence.
		

Crossrefs

Including evens gives A355529.
The version for all divisors including evens is A355740, zeros of A355739.
Choices of a prime factor of each prime index: A355741, unordered A355744.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A120383 lists numbers divisible by all of their prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],OddQ[#]&&Select[Tuples[primeMS/@primeMS[#]],UnsameQ@@#&]=={}&]
Previous Showing 21-30 of 65 results. Next