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 11-20 of 74 results. Next

A355745 Number of ways to choose a prime factor of each prime index of n (with multiplicity, in weakly increasing order) such that the result is also weakly increasing.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 18 2022

Keywords

Comments

First differs from A355741 and A355744 at n = 35.
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 1469 are {6,30}, and there are five valid choices: (2,2), (2,3), (2,5), (3,3), (3,5), so a(1469) = 5.
		

Crossrefs

Allowing all divisors gives A355735, firsts A355736, reverse A355749.
Not requiring an increasing sequence gives A355741.
Choosing a multiset instead of sequence gives A355744.
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.
A355731 chooses of a divisor of each prime index, firsts A355732.
A355733 chooses a multiset of divisors, firsts A355734.

Programs

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

A355733 Number of multisets that can be obtained by choosing a divisor of each prime index of n.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 1, 3, 2, 2, 2, 4, 3, 4, 1, 2, 3, 4, 2, 5, 2, 3, 2, 3, 4, 4, 3, 4, 4, 2, 1, 4, 2, 6, 3, 6, 4, 7, 2, 2, 5, 4, 2, 6, 3, 4, 2, 6, 3, 4, 4, 5, 4, 4, 3, 7, 4, 2, 4, 6, 2, 7, 1, 7, 4, 2, 2, 6, 6, 6, 3, 4, 6, 6, 4, 6, 7, 4, 2, 5, 2, 2, 5, 4, 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) = 4 multisets are: {1,1}, {1,2}, {1,3}, {2,3}.
The a(18) = 3 multisets are: {1,1,1}, {1,1,2}, {1,2,2}.
		

Crossrefs

Counting all choices of divisors gives A355731, firsts A355732.
Positions of first appearances are A355734.
Choosing weakly increasing divisors gives A355735, firsts A355736.
Choosing only prime divisors gives A355744.
The version choosing a divisor of each number from 1 to n is A355747.
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.
A340852 lists numbers that can be factored into divisors of bigomega.

Programs

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

A324923 Number of distinct factors in the factorization of n into factors q(i) = prime(i)/i, i > 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

Also the number of distinct proper terminal subtrees of the rooted tree with Matula-Goebel number n. See illustrations in A061773.

Examples

			The factorization 22 = q(1)^2 q(2) q(3) q(5) has four distinct factors, so a(22) = 4.
		

Crossrefs

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[Length[Union[difac[n]]],{n,100}]
  • PARI
    A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]);
    A324923(n) = { my(lista = List([]), gpf, i); while(n > 1, gpf=A006530(n); i = primepi(gpf); n /= gpf; n *= i; listput(lista,i)); #Set(lista); }; \\ Antti Karttunen, Oct 23 2023

Formula

a(n) = A317713(n) - 1.
a(n) = A196050(n) - A366386(n). - Antti Karttunen, Oct 23 2023

Extensions

Data section extended up to a(108) by Antti Karttunen, Oct 23 2023

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 17 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(2) = 1 through a(18) = 4 choices:
  1  1  11  1  11  1  111  11  11  1  111  1  11  11  1111  1  111
               12          12  13     112     12  13           112
                           21                 14  21           121
                                                  23           122
		

Crossrefs

Dominated by A355731, firsts A355732, primes A355741, prime-powers A355742.
For weakly increasing instead of coprime we have A355735, primes A355745.
Positions of first appearances are A355738.
For strict instead of coprime we have A355739, zeros A355740.
A000005 counts divisors.
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.
A289508 gives GCD of prime indices.
A289509 ranks relatively prime partitions, odd A302697, squarefree A302796.
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]],GCD@@#==1&]],{n,100}]

A324851 Numbers > 1 divisible by the sum of their prime indices.

Original entry on oeis.org

2, 4, 6, 12, 15, 16, 20, 30, 35, 36, 42, 48, 56, 88, 99, 112, 120, 126, 130, 135, 143, 144, 160, 162, 180, 192, 210, 216, 220, 221, 228, 231, 242, 250, 256, 270, 275, 280, 288, 297, 300, 308, 322, 330, 338, 360, 396, 400, 408, 429, 435, 440, 455, 468, 480, 493
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. The sum of prime indices of n is A056239(n). For example, the prime indices of 99 are {2,2,5}, with sum 9, a divisor of 99, so 99 is in the sequence.
For any k>=2, let d be a divisor of k such that d > A056239(k). Then 2^(d-A056239(k))*k is in the sequence. Similarly if k is in the sequence with d = A056239(k), then 2^d*k is in the sequence. - Robert Israel, Mar 19 2019

Examples

			The sequence of terms together with their prime indices begins:
    2: {1}
    4: {1,1}
    6: {1,2}
   12: {1,1,2}
   15: {2,3}
   16: {1,1,1,1}
   20: {1,1,3}
   30: {1,2,3}
   35: {3,4}
   36: {1,1,2,2}
   42: {1,2,4}
   48: {1,1,1,1,2}
   56: {1,1,1,4}
   88: {1,1,1,5}
   99: {2,2,5}
  112: {1,1,1,1,4}
  120: {1,1,1,2,3}
  126: {1,2,2,4}
  130: {1,3,6}
  135: {2,2,2,3}
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local t; n mod add(numtheory:-pi(t[1])*t[2],t=ifactors(n)[2]) = 0 end proc:
    select(filter, [$1..1000]); # Robert Israel, Mar 19 2019
  • Mathematica
    Select[Range[2,100],Divisible[#,Plus@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]*k]]&]
  • PARI
    isok(n) = {my(f = factor(n)); (n!=1) && !(n % sum(k=1, #f~, primepi(f[k,1])*f[k,2]));} \\ Michel Marcus, Mar 19 2019

A324924 Irregular triangle read by rows giving the factorization of n into factors q(i) = prime(i)/i, i > 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

Row n is the multiset of Matula-Goebel numbers of all proper terminal subtrees of the rooted tree with Matula-Goebel number n. For example, the rooted tree with Matula-Goebel number 1362 is (o(o)((oo)(oo))), with proper terminal subtrees {o,o,o,o,o,o,(o),(oo),(oo),((oo)(oo))}, which have Matula-Goebel numbers {1,1,1,1,1,1,2,4,4,49}, which is row 1362, as required.

Examples

			Triangle begins:
  {}
  1
  1  2
  1  1
  1  2  3
  1  1  2
  1  1  4
  1  1  1
  1  1  2  2
  1  1  2  3
  1  2  3  5
  1  1  1  2
  1  1  2  6
  1  1  1  4
  1  1  2  2  3
  1  1  1  1
  1  1  4  7
  1  1  1  2  2
  1  1  1  8
  1  1  1  2  3
  1  1  1  2  4
  1  1  2  3  5
  1  1  2  2  9
For example, row 65 is {1,1,1,2,2,3,6} because 65 = q(1)^3 q(2)^2 q(3) q(6).
		

Crossrefs

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[difac[n],{n,30}]

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

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 3, 7, 7, 11, 11, 22, 15, 30, 42, 77, 42, 101, 56, 176, 176, 231, 135, 490, 490, 490, 792, 1002, 490, 1575, 627, 3010, 2436, 2436, 3718, 5604, 1958, 4565, 6842, 12310, 3718, 14883, 4565, 21637, 26015, 17977, 8349, 53174, 44583, 63261
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(10) = 11 partitions:
  1  11  21  211   32   321    43    5111      522      631
             1111  311  2211   421   32111     3222     3331
                        21111  4111  41111     4221     4321
                                     221111    22221    5311
                                     311111    32211    32221
                                     2111111   222111   33211
                                     11111111  2211111  43111
                                                        322111
                                                        331111
                                                        3211111
                                                        31111111
For example, the Heinz number of (3,2) is 15, which is divisible by 5, so (3,2) is counted under a(5).
		

Crossrefs

The Heinz numbers of these partitions are given by 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.

Programs

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

A324846 Positive integers divisible by none of their prime indices.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 13, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 47, 49, 51, 53, 57, 59, 61, 63, 65, 67, 69, 71, 73, 77, 79, 81, 83, 85, 87, 89, 91, 93, 95, 97, 99, 101, 103, 107, 109, 111, 113, 115, 117, 121, 123, 125, 127, 129, 131, 133, 137
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. For example, the prime indices of 5673 are {2,11,18}, none of which divides 5673, so 5673 belongs to the sequence.

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  11: {5}
  13: {6}
  17: {7}
  19: {8}
  21: {2,4}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  29: {10}
  31: {11}
  33: {2,5}
  35: {3,4}
  37: {12}
  39: {2,6}
		

Crossrefs

Programs

  • Maple
    q:= n-> ormap(i-> irem(n, numtheory[pi](i[1]))=0, ifactors(n)[2]):
    remove(q, [$1..200])[];  # Alois P. Heinz, 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 (0));); return (1);} \\ Michel Marcus, Mar 19 2019

A355742 Number of ways to choose a sequence of prime-power divisors, one of each prime index of n. Product of bigomega over the prime indices of n, with multiplicity.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 20 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 prime indices of 49 are {4,4}, and the a(49) = 4 choices are: (2,2), (2,4), (4,2), (4,4).
The prime indices of 777 are {2,4,12}, and the a(777) = 6 choices are: (2,2,2), (2,2,3), (2,2,4), (2,4,2), (2,4,3), (2,4,4).
		

Crossrefs

The unordered version is A001970, row-sums of A061260.
Positions of 1's are A076610, just primes A355743.
Positions of 0's are A299174.
Allowing all divisors (not just primes) gives A355731, firsts A355732.
Choosing only prime factors (not prime-powers) gives A355741.
Counting multisets of primes gives A355744.
The case of weakly increasing primes A355745, all divisors A355735.
A000688 counts factorizations into prime powers.
A001414 adds up distinct prime factors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.

Programs

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

Formula

Totally multiplicative with a(prime(k)) = A001222(k).

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
Previous Showing 11-20 of 74 results. Next