A304360 -id:A304360 - cp's OEIS frontend

A317963 First differences of A304360.

2, 1, 3, 2, 3, 3, 1, 3, 3, 2, 1, 5, 1, 2, 3, 3, 3, 3, 1, 3, 2, 9, 1, 2, 1, 2, 1, 5, 1, 5, 1, 5, 1, 3, 3, 2, 6, 3, 1, 3, 2, 4, 2, 7, 2, 1, 3, 2, 4, 2, 1, 9, 2, 1, 2, 4, 3, 2, 3, 4, 2, 1, 2, 6, 1, 2, 3, 1, 3, 5, 4, 3, 6, 2, 6, 1, 3, 3, 2, 3, 1, 3, 3, 2, 6, 3, 4, 2, 1, 5, 1, 3, 3, 5, 1, 2, 1, 2, 3

Offset: 1

N. J. A. Sloane, Aug 26 2018

nonn

Always positive (see A304360).

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A304360.

A304360:= NULL: count:= 0:
P:= {}:
for n  from 2 while count < 201 do
  pn:= numtheory:-factorset(n);
  if pn intersect P = {} then
    count:= count+1;
    A304360:= A304360, n;
    P:= P union {ithprime(n)};
  fi
od:
A304360:= [A304360]:
A304360[2..-1]-A304360[1..-2]; # Robert Israel, Aug 26 2018

More terms from Robert Israel, Aug 26 2018

A317964 Prime numbers in the lexicographically earliest sequence of positive integers whose prime indices are not already in the sequence (A304360).

Original entry on oeis.org

2, 5, 13, 17, 23, 31, 37, 43, 47, 61, 67, 73, 79, 89, 103, 107, 109, 113, 137, 149, 151, 163, 167, 179, 181, 193, 197, 223, 227, 233, 241, 251, 257, 263, 269, 271, 277, 281, 307, 317, 347, 349, 353, 359, 379, 383, 389, 397, 419, 421, 431, 433, 449, 457, 463, 467, 487, 499, 503, 509, 521, 523, 547

Offset: 1

N. J. A. Sloane, Aug 26 2018

nonn

Also primes whose prime index is not in A304360, or is in A324696. A prime index of n is a number m such that prime(m) divides n. - Gus Wiseman, Mar 19 2019

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A000720, A001462, A007097, A060197, A079254, A112798, A276625, A277098, A290822, A304360, A306844.

Cf. A324695, A324698, A324704, A324741, A324756, A324758, A324765, A324766, A324768.

count:= 0:
P:= {}: A:= NULL:
for n from 2 while count < 100 do
  pn:= numtheory:-factorset(n);
  if pn intersect P = {} then
    P:= P union {ithprime(n)};
    if isprime(n) then A:= A, n; count:= count+1 fi;
  fi
od:
A; # Robert Israel, Aug 26 2018

aQ[n_]:=n==1||Or@@Cases[FactorInteger[n],{p_,_}:>!aQ[PrimePi[p]]];
Prime[Select[Range[100],aQ]] (* Gus Wiseman, Mar 19 2019 *)

A066328 a(n) = sum of indices of distinct prime factors of n; here, index(i-th prime) = i.

Original entry on oeis.org

0, 1, 2, 1, 3, 3, 4, 1, 2, 4, 5, 3, 6, 5, 5, 1, 7, 3, 8, 4, 6, 6, 9, 3, 3, 7, 2, 5, 10, 6, 11, 1, 7, 8, 7, 3, 12, 9, 8, 4, 13, 7, 14, 6, 5, 10, 15, 3, 4, 4, 9, 7, 16, 3, 8, 5, 10, 11, 17, 6, 18, 12, 6, 1, 9, 8, 19, 8, 11, 8, 20, 3, 21, 13, 5, 9, 9, 9, 22, 4, 2, 14, 23, 7, 10, 15, 12, 6, 24, 6, 10

Offset: 1

Leroy Quet, Jan 01 2002

nonn

Equals row sums of triangle A143542. - Gary W. Adamson, Aug 23 2008

a(n) = the sum of the distinct parts of the partition with Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} (p_j-th prime) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(75) = 5; indeed, the partition having Heinz number 75 = 3*5*5 is [2,3,3] and 2 + 3 = 5. - Emeric Deutsch, Jun 04 2015

			a(24) = 1 + 2 = 3 because 24 = 2^3 * 3 = p(1)^3 * p(2), p(k) being the k-th prime.
From _Gus Wiseman_, Mar 09 2019: (Start)
The distinct prime indices of 1..20 and their sums.
   1: () = 0
   2: (1) = 1
   3: (2) = 2
   4: (1) = 1
   5: (3) = 3
   6: (1+2) = 3
   7: (4) = 4
   8: (1) = 1
   9: (2) = 2
  10: (1+3) = 4
  11: (5) = 5
  12: (1+2) = 3
  13: (6) = 6
  14: (1+4) = 5
  15: (2+3) = 5
  16: (1) = 1
  17: (7) = 7
  18: (1+2) = 3
  19: (8) = 8
  20: (1+3) = 4
(End)

Antti Karttunen, Table of n, a(n) for n = 1..65537 (terms 1..1000 from Harry J. Smith)
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to Heinz numbers

Cf. A143542. - Gary W. Adamson, Aug 23 2008
Cf. A000720, A056239, A136565.
Cf. A001221, A046660, A112798, A114638, A116861, A304360.

with(numtheory): seq(add(pi(d), d in factorset(n)), n=1..100); # Ridouane Oudra, Aug 19 2019

PrimeFactors[n_Integer] := Flatten[ Table[ #[[1]], {1}] & /@ FactorInteger[n]]; f[n_] := (Plus @@ PrimePi[ PrimeFactors[n]]); Table[ f[n], {n, 91}] (* Robert G. Wilson v, May 04 2004 *)

{ for (n=1, 1000, f=factor(n); a=0; for (i=1, matsize(f)[1], a+=primepi(f[i, 1])); write("b066328.txt", n, " ", a) ) } \\ Harry J. Smith, Feb 10 2010

a(n)=my(f=factor(n)[,1]); sum(i=1,#f,primepi(f[i])) \\ Charles R Greathouse IV, May 11 2015

A066328(n) = vecsum(apply(primepi,(factor(n)[,1]))); \\ Antti Karttunen, Sep 06 2018

from sympy import primepi, primefactors
def A066328(n): return sum(map(primepi,primefactors(n))) # Chai Wah Wu, Mar 13 2024

G.f.: Sum_{k>=1} k*x^prime(k)/(1-x^prime(k)). - Vladeta Jovovic, Aug 11 2004
Additive with a(p^e) = PrimePi(p), where PrimePi(n) = A000720(n).
a(n) = A056239(A007947(n)). - Antti Karttunen, Sep 06 2018
a(n) = Sum_{p|n} A000720(p), where p is a prime. - Ridouane Oudra, Aug 19 2019

A324695 Lexicographically earliest sequence of positive integers whose prime indices are not already in the sequence.

Original entry on oeis.org

1, 3, 7, 9, 11, 13, 19, 21, 27, 29, 33, 37, 39, 43, 47, 49, 53, 57, 59, 61, 63, 71, 77, 79, 81, 83, 87, 89, 91, 97, 99, 101, 107, 111, 113, 117, 121, 127, 129, 131, 133, 139, 141, 143, 147, 149, 151, 159, 163, 169, 171, 173, 177, 179, 181, 183, 189, 193, 197

Offset: 1

Gus Wiseman, Mar 10 2019

nonn

A self-describing sequence, similar to A304360.

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 sequence of terms together with their prime indices begins:
   1: {}
   3: {2}
   7: {4}
   9: {2,2}
  11: {5}
  13: {6}
  19: {8}
  21: {2,4}
  27: {2,2,2}
  29: {10}
  33: {2,5}
  37: {12}
  39: {2,6}
  43: {14}
  47: {15}
  49: {4,4}
  53: {16}
  57: {2,8}
  59: {17}
  61: {18}
  63: {2,2,4}

Complement of A324694. Prime indices are A304360.
Cf. A000002, A000720, A001222, A001462, A007097, A055396, A061395, A079000, A079254, A109298, A112798, A276625, A277098.
Cf. A324696, A324697, A324698, A324699, A324700, A324701, A324702, A324703, A324704, A324705.

aQ[n_]:=And@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>!aQ[PrimePi[p]]];
Select[Range[100],aQ]

A306844 Number of anti-transitive rooted trees with n nodes.

Original entry on oeis.org

1, 1, 2, 3, 7, 14, 36, 83, 212, 532, 1379, 3577, 9444, 25019, 66943, 179994, 487031, 1323706, 3614622, 9907911

Offset: 1

Gus Wiseman, Mar 13 2019

A rooted tree is anti-transitive if the subbranches are disjoint from the branches, i.e., no branch of a branch is a branch.

			The a(1) = 1 through a(6) = 14 anti-transitive rooted trees:
  o  (o)  (oo)   (ooo)    (oooo)     (ooooo)
          ((o))  ((oo))   ((ooo))    ((oooo))
                 (((o)))  (((oo)))   (((ooo)))
                          ((o)(o))   ((o)(oo))
                          ((o(o)))   ((o(oo)))
                          (o((o)))   ((oo(o)))
                          ((((o))))  (o((oo)))
                                     (oo((o)))
                                     ((((oo))))
                                     (((o)(o)))
                                     (((o(o))))
                                     ((o((o))))
                                     (o(((o))))
                                     (((((o)))))

Gus Wiseman, The a(7) = 36 anti-transitive rooted trees.
Gus Wiseman, The a(10) = 532 anti-transitive rooted trees.

Cf. A276625, A279861, A279861, A290689, A290760, A304360.

Cf. A324694, A324695, A324738, A324741, A324743, A324751, A324754, A324756, A324758, A324759, A324764.

rtall[n_]:=Union[Sort/@Join@@(Tuples[rtall/@#]&/@IntegerPartitions[n-1])];
Table[Length[Select[rtall[n],Intersection[Union@@#,#]=={}&]],{n,10}]

a(16)-a(20) from Jinyuan Wang, Jun 20 2020

A109298 Primal codes of finite idempotent functions on positive integers.

Original entry on oeis.org

1, 2, 9, 18, 125, 250, 1125, 2250, 2401, 4802, 21609, 43218, 161051, 300125, 322102, 600250, 1449459, 2701125, 2898918, 4826809, 5402250, 9653618, 20131375, 40262750, 43441281, 86882562, 181182375, 362364750, 386683451, 410338673, 603351125, 773366902, 820677346

Offset: 1

Jon Awbrey, Jul 06 2005

nonn

Finite idempotent functions are identity maps on finite subsets, counting the empty function as the idempotent on the empty set.
From Gus Wiseman, Mar 09 2019: (Start)
Also numbers whose ordered prime signature is equal to the distinct prime indices in increasing order. A prime index of n is a number m such that prime(m) divides n. The ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base. The case where the prime indices are taken in decreasing order is A324571.
Also numbers divisible by prime(k) exactly k times for each prime index k. These are a kind of self-describing numbers (cf. A001462, A304679).
Also Heinz numbers of integer partitions where the multiplicity of m is m for all m in the support (counted by A033461). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of distinct elements of A062457. For example, 43218 = prime(1)^1 * prime(2)^2 * prime(4)^4.
(End)

			Writing (prime(i))^j as i:j, we have the following table of examples:
Primal Codes of Finite Idempotent Functions on Positive Integers
` ` ` 1 = { }
` ` ` 2 = 1:1
` ` ` 9 = ` ` 2:2
` ` `18 = 1:1 2:2
` ` 125 = ` ` ` ` 3:3
` ` 250 = 1:1 ` ` 3:3
` `1125 = ` ` 2:2 3:3
` `2250 = 1:1 2:2 3:3
` `2401 = ` ` ` ` ` ` 4:4
` `4802 = 1:1 ` ` ` ` 4:4
` 21609 = ` ` 2:2 ` ` 4:4
` 43218 = 1:1 2:2 ` ` 4:4
`161051 = ` ` ` ` ` ` ` ` 5:5
`300125 = ` ` ` ` 3:3 4:4
`322102 = 1:1 ` ` ` ` ` ` 5:5
`600250 = 1:1 ` ` 3:3 4:4
From _Gus Wiseman_, Mar 09 2019: (Start)
The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2) has prime signature {1,2} and the distinct prime indices are also {1,2}.
       1: {}
       2: {1}
       9: {2,2}
      18: {1,2,2}
     125: {3,3,3}
     250: {1,3,3,3}
    1125: {2,2,3,3,3}
    2250: {1,2,2,3,3,3}
    2401: {4,4,4,4}
    4802: {1,4,4,4,4}
   21609: {2,2,4,4,4,4}
   43218: {1,2,2,4,4,4,4}
  161051: {5,5,5,5,5}
  300125: {3,3,3,4,4,4,4}
  322102: {1,5,5,5,5,5}
  600250: {1,3,3,3,4,4,4,4}
(End)

David A. Corneth, Table of n, a(n) for n = 1..10000
J. Awbrey, Riffs and Rotes

Cf. A076954, A106177, A108352, A108371, A109297, A109301.
Cf. A001156, A033461, A056239, A062457, A112798, A118914, A124010 (ordered prime signature), A181819, A276078, A304679.
Cf. A324524, A324525, A324570, A324571, A324572, A324587, A324588.
Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Select[Range[10000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]==k]&]

is(n) = my(f = factor(n)); for(i = 1, #f~, if(prime(f[i, 2]) != f[i, 1], return(0))); 1 \\ David A. Corneth, Mar 09 2019

Sum_{n>=1} 1/a(n) = Product_{n>=1} (1 + 1/prime(n)^n) = 1.6807104966... - Amiram Eldar, Jan 03 2021

Offset set to 1, missing terms inserted and more terms added by Alois P. Heinz, Mar 08 2019

A324694 Lexicographically earliest sequence of positive integers divisible by prime(m) for some m not already in the sequence.

Original entry on oeis.org

2, 4, 5, 6, 8, 10, 12, 14, 15, 16, 17, 18, 20, 22, 23, 24, 25, 26, 28, 30, 31, 32, 34, 35, 36, 38, 40, 41, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 76, 78, 80, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95

Offset: 1

Gus Wiseman, Mar 10 2019

nonn

A self-describing sequence, similar to A304360.

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 sequence of terms together with their prime indices begins:
   2: {1}
   4: {1,1}
   5: {3}
   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}
  17: {7}
  18: {1,2,2}
  20: {1,1,3}
  22: {1,5}
  23: {9}
  24: {1,1,1,2}
  25: {3,3}
  26: {1,6}
  28: {1,1,4}
  30: {1,2,3}

Complement of A324695.
Cf. A000002, A000720, A001222, A001462, A007097, A055396, A061395, A079000, A079254, A109298, A112798, A276625, A277098, A304360.
Cf. A324696, A324697, A324698, A324699, A324700, A324701, A324702, A324703, A324704, A324705.

aQ[n_]:=!And@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>aQ[PrimePi[p]]];
Select[Range[100],aQ]

A324758 Heinz numbers of integer partitions containing no prime indices of the parts.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 16, 17, 19, 20, 21, 22, 23, 25, 27, 29, 31, 32, 33, 34, 35, 37, 40, 41, 43, 44, 46, 47, 49, 50, 51, 53, 57, 59, 61, 62, 63, 64, 65, 67, 68, 71, 73, 77, 79, 80, 81, 82, 83, 85, 87, 88, 89, 91, 92, 93, 94, 95, 97, 99, 100, 101

Offset: 1

Gus Wiseman, Mar 17 2019

nonn

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 Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
These could be described as anti-transitive numbers (cf. A290822), as they are numbers x such that if prime(y) divides x and prime(z) divides y, then prime(z) does not divide x.
Also numbers n such that A003963(n) is coprime to n.

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

The subset version is A324741, with maximal case A324743. The strict integer partition version is A324751. The integer partition version is A324756. An infinite version is A324695.

Cf. A000720, A001221, A001462, A007097, A056239, A112798, A276625, A289509, A290822, A304360, A306844, A324742, A324753, A324764.

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

A324756 Number of integer partitions of n containing no prime indices of the parts.

Original entry on oeis.org

1, 1, 2, 2, 4, 3, 7, 7, 9, 11, 16, 16, 24, 25, 34, 39, 50, 54, 70, 79, 96, 111, 135, 152, 186, 208, 249, 285, 335, 377, 448, 506, 588, 664, 777, 873, 1010, 1139, 1309, 1471, 1697, 1890, 2175, 2435, 2772, 3106, 3532, 3941, 4478, 4995, 5643, 6297, 7107, 7897

Offset: 0

Gus Wiseman, Mar 17 2019

nonn

These could be described as anti-transitive integer partitions.

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 a(1) = 1 through a(8) = 9 integer partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (311)    (33)      (43)       (44)
                    (31)    (11111)  (42)      (52)       (71)
                    (1111)           (51)      (331)      (422)
                                     (222)     (511)      (2222)
                                     (3111)    (31111)    (3311)
                                     (111111)  (1111111)  (5111)
                                                          (311111)
                                                          (11111111)

The subset version is A324741, with maximal case A324743. The strict case is A324751. The Heinz number version is A324758. An infinite version is A324695.

Cf. A000720, A000837, A001462, A051424, A112798, A276625, A304360, A306844, A324764, A324742, A324753.

Table[Length[Select[IntegerPartitions[n],Intersection[#,PrimePi/@First/@Join@@FactorInteger/@#]=={}&]],{n,0,30}]

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

Gus Wiseman, Mar 18 2019

nonn

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.

			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}

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Complement of A324847.
Cf. A000720, A001222, A003963, A056239, A112798, A120383, A289509, A290822, A304360, A306844.
Cf. A324695, A324741, A324743, A324756, A324758, A324765, A324848, A324849, A324850, A324852, A324853.

q:= n-> ormap(i-> irem(n, numtheory[pi](i[1]))=0, ifactors(n)[2]):
remove(q, [$1..200])[];  # Alois P. Heinz, Mar 19 2019

Select[Range[100],!Or@@Cases[If[#==1,{},FactorInteger[#]],{p_,_}:>Divisible[#,PrimePi[p]]]&]

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

cp's OEIS Frontend

A317963 First differences of A304360.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Extensions

A317964 Prime numbers in the lexicographically earliest sequence of positive integers whose prime indices are not already in the sequence (A304360).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

A066328 a(n) = sum of indices of distinct prime factors of n; here, index(i-th prime) = i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

A324695 Lexicographically earliest sequence of positive integers whose prime indices are not already in the sequence.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A306844 Number of anti-transitive rooted trees with n nodes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Extensions

A109298 Primal codes of finite idempotent functions on positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A324694 Lexicographically earliest sequence of positive integers divisible by prime(m) for some m not already in the sequence.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A324758 Heinz numbers of integer partitions containing no prime indices of the parts.