A112798 -id:A112798 - cp's OEIS frontend

A195017 If n = Product_{k >= 1} (p_k)^(c_k) where p_k is k-th prime and c_k >= 0 then a(n) = Sum_{k >= 1} c_k*((-1)^(k-1)).

0, 1, -1, 2, 1, 0, -1, 3, -2, 2, 1, 1, -1, 0, 0, 4, 1, -1, -1, 3, -2, 2, 1, 2, 2, 0, -3, 1, -1, 1, 1, 5, 0, 2, 0, 0, -1, 0, -2, 4, 1, -1, -1, 3, -1, 2, 1, 3, -2, 3, 0, 1, -1, -2, 2, 2, -2, 0, 1, 2, -1, 2, -3, 6, 0, 1, 1, 3, 0, 1, -1, 1, 1, 0, 1, 1, 0, -1, -1, 5, -4, 2, 1, 0, 2, 0, -2, 4, -1, 0, -2, 3, 0, 2, 0, 4, 1, -1, -1, 4, -1, 1, 1, 2, -1

Offset: 1

Clark Kimberling, Feb 06 2012

Let p(n,x) be the completely additive polynomial-valued function such that p(1,x) = 0 and p(prime(n),x) = x^(n-1), like is defined in A206284 (although here we are not limited to just irreducible polynomials). Then a(n) is the value of the polynomial encoded in such a manner by n, when it is evaluated at x=-1. - The original definition rewritten and clarified by Antti Karttunen, Oct 03 2018
Positions of 0 give the values of n for which the polynomial p(n,x) is divisible by x+1. For related sequences, see the Mathematica section.
Also the number of odd prime indices of n minus the number of even prime indices of n (both counted with multiplicity), where 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. - Gus Wiseman, Oct 24 2023

			The sequence can be read from a list of the polynomials:
  p(n,x)      with x = -1, gives a(n)
------------------------------------------
  p(1,x) = 0           0
  p(2,x) = 1x^0        1
  p(3,x) = x          -1
  p(4,x) = 2x^0        2
  p(5,x) = x^2         1
  p(6,x) = 1+x         0
  p(7,x) = x^3        -1
  p(8,x) = 3x^0        3
  p(9,x) = 2x         -2
  p(10,x) = x^2 + 1    2.
(The list runs through all the polynomials whose coefficients are nonnegative integers.)

Antti Karttunen, Table of n, a(n) for n = 1..65537

Cf. A206284, A277322, A284010.
For other evaluation functions of such encoded polynomials, see A001222, A048675, A056239, A090880, A248663.
Cf. A006093, A036689, A078437.
Zeros are A325698, distinct A325700.
For sum instead of count we have A366749 = A366531 - A366528.
A000009 counts partitions into odd parts, ranked by A066208.
A035363 counts partitions into even parts, ranked by A066207.
A112798 lists prime indices, reverse A296150, sum A056239.
A257991 counts odd prime indices, even A257992.
A300061 lists numbers with even sum of prime indices, odd A300063.
Cf. A019507, A106529, A239241, A239261, A369180.

b[n_] := Table[x^k, {k, 0, n}];
f[n_] := f[n] = FactorInteger[n]; z = 200;
t[n_, m_, k_] := If[PrimeQ[f[n][[m, 1]]] && f[n][[m, 1]]
== Prime[k], f[n][[m, 2]], 0];
u = Table[Apply[Plus,
    Table[Table[t[n, m, k], {k, 1, PrimePi[n]}], {m, 1,
      Length[f[n]]}]], {n, 1, z}];
p[n_, x_] := u[[n]].b[-1 + Length[u[[n]]]]
Table[p[n, x] /. x -> 0, {n, 1, z/2}]   (* A007814 *)
Table[p[2 n, x] /. x -> 0, {n, 1, z/2}] (* A001511 *)
Table[p[n, x] /. x -> 1, {n, 1, z}]     (* A001222 *)
Table[p[n, x] /. x -> 2, {n, 1, z}]     (* A048675 *)
Table[p[n, x] /. x -> 3, {n, 1, z}]     (* A090880 *)
Table[p[n, x] /. x -> -1, {n, 1, z}]    (* A195017 *)
z = 100; Sum[-(-1)^k IntegerExponent[Range[z], Prime[k]], {k, 1, PrimePi[z]}] (* Friedjof Tellkamp, Aug 05 2024 *)

A195017(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * (-1)^(1+primepi(f[i,1])))); } \\ Antti Karttunen, Oct 03 2018

Totally additive with a(p^e) = e * (-1)^(1+PrimePi(p)), where PrimePi(n) = A000720(n). - Antti Karttunen, Oct 03 2018
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} = (-1)^(primepi(p)+1)/(p-1) = Sum_{k>=1} (-1)^(k+1)/A006093(k) = A078437 +  Sum_{k>=1} (-1)^(k+1)/A036689(k) = 0.6339266524059... . - Amiram Eldar, Sep 29 2023
a(n) = A257991(n) - A257992(n). - Gus Wiseman, Oct 24 2023
a(n) = -Sum_{k=1..pi(n)} (-1)^k * valuation(n, prime(k)). - Friedjof Tellkamp, Aug 05 2024

More terms, name changed and example-section edited by Antti Karttunen, Oct 03 2018

A320628 Products of primes of nonprime index.

Original entry on oeis.org

1, 2, 4, 7, 8, 13, 14, 16, 19, 23, 26, 28, 29, 32, 37, 38, 43, 46, 47, 49, 52, 53, 56, 58, 61, 64, 71, 73, 74, 76, 79, 86, 89, 91, 92, 94, 97, 98, 101, 103, 104, 106, 107, 112, 113, 116, 122, 128, 131, 133, 137, 139, 142, 146, 148, 149, 151, 152, 158, 161, 163

Offset: 1

Gus Wiseman, Oct 18 2018

nonn

The index of a prime number n is the number m such that n is the m-th prime.

The asymptotic density of this sequence is Product_{p in A006450} (1 - 1/p) = 1/(Sum_{n>=1} 1/A076610(n)) < 1/3. - Amiram Eldar, Feb 02 2021

			The sequence of terms begins:
   1 = 1
   2 = prime(1)
   4 = prime(1)^2
   7 = prime(4)
   8 = prime(1)^3
  13 = prime(6)
  14 = prime(1)*prime(4)
  16 = prime(1)^4
  19 = prime(8)
  23 = prime(9)
  26 = prime(1)*prime(6)
  28 = prime(1)^2*prime(4)
  29 = prime(10)
  32 = prime(1)^5
  37 = prime(12)
  38 = prime(1)*prime(8)
  43 = prime(14)
  46 = prime(1)*prime(9)
  47 = prime(15)
  49 = prime(4)^2
  52 = prime(1)^2*prime(6)
  53 = prime(16)
  56 = prime(1)^3*prime(4)
  58 = prime(1)*prime(10)
  61 = prime(18)
  64 = prime(1)^6
  71 = prime(20)
  73 = prime(21)
  74 = prime(1)*prime(12)
  76 = prime(1)^2*prime(8)
  79 = prime(22)
  86 = prime(1)*prime(14)
  89 = prime(24)
  91 = prime(4)*prime(6)
  92 = prime(1)^2*prime(9)
  94 = prime(1)*prime(15)
  97 = prime(25)
  98 = prime(1)*prime(4)^2

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A018252, A056239, A290103, A302242, A302478, A320533, A320629, A320630, A320631, A320633.
Complement of A331386.
Positions of zeros in A257994.
Primes of prime index are A006450.
Primes of nonprime index are A007821.
Products of primes of prime index are A076610.
Products of primes of nonprime index are this sequence.
The number of prime prime indices is given by A257994.
The number of nonprime prime indices is given by A330944.
Cf. A000040, A000720, A001222, A112798, A295665.

Select[Range[100],And@@Not/@PrimeQ/@PrimePi/@First/@FactorInteger[#]&]

A109082 Depth of rooted tree having Matula-Goebel number n.

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, 2, 4, 3, 2, 3, 3, 2, 2, 4, 3, 5, 1, 4, 3, 3, 2, 3, 2, 3, 3, 4, 2, 3, 4, 3, 3, 4, 2, 2, 3, 3, 3, 2, 2, 4, 2, 2, 4, 4, 3, 3, 5, 2, 1, 3, 4, 3, 3, 3, 3, 4, 2, 3, 3, 3, 2, 4, 3, 5, 3, 2, 4, 4, 2, 3, 3, 4, 4, 3, 3, 3, 3, 5, 4, 3, 2, 4, 2, 4, 3

Offset: 1

Keith Briggs, Aug 17 2005

nonn

Another term for depth is height.

Starting with n, a(n) is the number of times one must take the product of prime indices (A003963) to reach 1. - Gus Wiseman, Mar 27 2019

			a(7) = 2 because the rooted tree with Matula-Goebel number 7 is the 3-edge rooted tree Y of height 2.

François Marques, Table of n, a(n) for n = 1..10000 (terms 1..5381 from Alois P. Heinz).
Emeric Deutsch, Tree statistics from Matula numbers, arXiv preprint arXiv:1111.4288 [math.CO], 2011.
F. Goebel, On a 1-1-correspondence between rooted trees and natural numbers, J. Combin. Theory, B 29 (1980), 141-143.
I. Gutman and A. Ivic, On Matula numbers, Discrete Math., 150, 1996, 131-142.
I. Gutman and Yeong-Nan Yeh, Deducing properties of trees from their Matula numbers, Publ. Inst. Math., 53 (67), 1993, 17-22.
D. W. Matula, A natural rooted tree enumeration by prime factorization, SIAM Rev. 10 (1968) 273.
Index entries for sequences related to Matula-Goebel numbers

A left inverse of A007097.
Cf. A003963, A061775, A091233.
Cf. A000081, A000720, A001222, A109129, A112798, A196050, A290822, A317713, A320325, A324927 (positions of 2), A324928 (positions of 3), A325032.
This statistic is counted by A034781, ordered A080936.
The ordered version is A358379.
For node-height instead of edge-height we have A358552.
Cf. A000108, A055277, A056239, A342507, A358576.

with(numtheory): a := proc(n) option remember; if n = 1 then 0 elif isprime(n) then 1+a(pi(n)) else max((map (p->a(p), factorset(n)))[]) end if end proc: seq(a(n), n = 1 .. 100); # Emeric Deutsch, Sep 16 2011

a [n_] := a[n] = If[n == 1, 0, If[PrimeQ[n], 1+a[PrimePi[n]], Max[Map[a, FactorInteger[n][[All, 1]]]]]]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, May 06 2014, after Emeric Deutsch *)

a(n) = my(v=factor(n)[,1],d=0); while(#v,d++; v=fold(setunion, apply(p->factor(primepi(p))[,1]~, v))); d; \\ Kevin Ryde, Sep 21 2020

from functools import lru_cache
from sympy import isprime, primepi, primefactors
@lru_cache(maxsize=None)
def A109082(n):
    if n == 1 : return 0
    if isprime(n): return 1+A109082(primepi(n))
    return max(A109082(p) for p in primefactors(n)) # Chai Wah Wu, Mar 19 2022

a(1)=0; if n is the t-th prime, then a(n) = 1 + a(t); if n is composite, n=t*s, then a(n) = max(a(t),a(s)). The Maple program is based on this.
a(A007097(n)) = n.
a(n) = A358552(n) - 1. - Gus Wiseman, Nov 27 2022

Edited by Emeric Deutsch, Sep 16 2011

A344606 Number of alternating permutations of the prime factors of n, counting multiplicity, including twins (x,x).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 0, 1, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 1, 2, 0, 1, 1, 4, 1, 0, 2, 2, 2, 2, 1, 2, 2, 0, 1, 4, 1, 1, 1, 2, 1, 0, 1, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 4, 1, 2, 1, 0, 2, 4, 1, 1, 2, 4, 1, 1, 1, 2, 1, 1, 2, 4, 1, 0, 0, 2, 1, 4, 2, 2, 2

Offset: 1

Gus Wiseman, May 28 2021

nonn

Differs from A335448 in having a(x^2) = 0 and a(270) = 0.
These are permutations of the prime factors of n, counting multiplicity, with no adjacent triples (..., x, y, z, ...) where x <= y <= z or x >= y >= z.
The version without twins (x,x) is A345164, which is identical to this sequence except when n is the square of a prime.

			The permutations for n = 2, 6, 30, 180, 210, 300, 420, 720, 840:
  2   23   253   23253   2537   25253   23275   2323252   232527
      32   325   32325   2735   25352   25273   2325232   232725
           352   32523   3275   32525   25372   2523232   252327
           523   35232   3527   35252   27253             252723
                 52323   3725   52325   27352             272325
                         5273   52523   32527             272523
                         5372           32725             325272
                         5723           35272             327252
                         7253           37252             523272
                         7352           52327             527232
                                        52723             723252
                                        57232             725232
                                        72325
                                        72523
For example, there are no alternating permutations of the prime factors of 270 because the only anti-runs are {3,2,3,5,3} and {3,5,3,2,3}, neither of which is alternating, so a(270) = 0.

The version for permutations is A001250.
The extension to anti-run permutations is A335452.
The version for compositions is A344604.
The version for patterns is A344605.
Positions of zeros are A344653 (counted by A344654).
Not including twins (x,x) gives A345164.
A008480 counts permutations of prime indices (strict: A335489, rank: A333221).
A056239 adds up prime indices,  row sums of A112798.
A071321 and A071322 are signed sums of prime factors.
A316523 is a signed sum of prime multiplicities.
A316524 and A344616 are signed sums of prime indices.
A325534 counts separable partitions (ranked by A335433).
A325535 counts inseparable partitions (ranked by A335448).
A344740 counts partitions with an alternating permutation or twin (x,x).
Cf. A000041, A000607, A000961, A001222, A003242, A026424, A028260, A049774, A103919, A181796, A343938, A344652.

Table[Length[Select[Permutations[Flatten[ConstantArray@@@FactorInteger[n]]],!MatchQ[#,{_,x_,y_,z_,_}/;x<=y<=z||x>=y>=z]&]],{n,100}]

A362611 Number of modes in the prime factorization of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 2, 2

Offset: 1

Gus Wiseman, May 05 2023

nonn

A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes of {a,a,b,b,b,c,d,d,d} are {b,d}.

a(n) depends only on the prime signature of n. - Andrew Howroyd, May 08 2023

			The factorization of 450 is 2*3*3*5*5, modes {3,5}, so a(450) = 2.
The factorization of 900 is 2*2*3*3*5*5, modes {2,3,5}, so a(900) = 3.
The factorization of 1500 is 2*2*3*5*5*5, modes {5}, so a(1500) = 1.
The factorization of 8820 is 2*2*3*3*5*7*7, modes {2,3,7}, so a(8820) = 3.

Andrew Howroyd, Table of n, a(n) for n = 1..10000

Positions of first appearances are A002110.
Positions of 1's are A356862, counted by A362608.
Positions of terms > 1 are A362605, counted by A362607.
For co-mode we have A362613, counted by A362615.
This statistic (mode-count) has triangular form A362614.
A027746 lists prime factors (with multiplicity).
A112798 lists prime indices, length A001222, sum A056239.
A359178 ranks partitions with a unique co-mode, counted by A362610.
A362606 ranks partitions with more than one co-mode, counted by A362609.
Cf. A000040, A000720, A001221, A002865, A072774, A215366, A327473, A327476, A359908.

Table[x=Last/@If[n==1,0,FactorInteger[n]];Count[x,Max@@x],{n,100}]

a(n) = if(n==1, 0, my(f=factor(n)[,2], m=vecmax(f)); #select(v->v==m, f)) \\ Andrew Howroyd, May 08 2023

from sympy import factorint
def A362611(n): return list(v:=factorint(n).values()).count(max(v,default=0)) # Chai Wah Wu, May 08 2023

For n > 1, 1 <= a(n) << log n. - Charles R Greathouse IV, May 09 2023

a(n) <= A001221(n), with equality if and only if n is a power of a squarefree number (A072774). - Amiram Eldar, Mar 02 2025

A182857 Smallest number that requires exactly n iterations to reach a fixed point under the x -> A181819(x) map.

Original entry on oeis.org

1, 3, 4, 6, 12, 60, 2520, 1286485200, 35933692027611398678865941374040400000

Offset: 0

Matthew Vandermast, Jan 05 2011

nonn

a(9) has 296 digits.
Related to Levine's sequence (A011784): A011784(n) = A001222(a(n)) = A001221(a(n+1)) = A051903(a(n+2)) = A071625(a(n+2)). Also see A182858.
Values of n where A182850(n) increases to a record.
The multiplicity of prime(k) in a(n+1) is the k-th largest prime index of a(n), which is A296150(a(n),k). - Gus Wiseman, May 13 2018

			From _Gus Wiseman_, May 13 2018: (Start)
Like A001462 the following sequence of multisets whose Heinz numbers belong to this sequence is a run-length describing sequence, as the number of k's in row n + 1 is equal to the k-th term of row n.
{2}
{1,1}
{1,2}
{1,1,2}
{1,1,2,3}
{1,1,1,2,2,3,4}
{1,1,1,1,2,2,2,3,3,4,4,5,6,7}
{1,1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,7,7,7,8,8,9,9,10,10,11,12,13,14}
(End)

Gus Wiseman, Table of n, a(n) for n = 0..9

Cf. A001462, A007755, A007916, A009287, A012257, A112798, A181819, A182850-A182858, A296150, A304455, A304464, A304465.

Prepend[Function[m,Times@@Prime/@m]/@NestList[Join@@Table[Table[i,{Reverse[#][[i]]}],{i,Length[#]}]&,{2},8],1] (* Gus Wiseman, May 13 2018 *)

For n > 0, a(n) = A181819(a(n+1)). For n > 1, a(n) = A181821(a(n-1)).

A353833 Numbers whose multiset of prime indices has all equal run-sums.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, 40, 41, 43, 47, 49, 53, 59, 61, 63, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 112, 113, 121, 125, 127, 128, 131, 137, 139, 144, 149, 151, 157, 163, 167, 169, 173, 179

Offset: 1

Gus Wiseman, May 23 2022

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 sequence of runs of a sequence consists of its maximal consecutive constant subsequences when read left-to-right. For example, the runs of (2,2,1,1,1,3,2,2) are (2,2), (1,1,1), (3), (2,2), with sums (4,3,3,4).

			The prime indices of 12 are {1,1,2}, with run-sums (2,2), so 12 is in the sequence.

Mathematics Stack Exchange, What is a sequence run? (answered 2011-12-01)

For parts instead of run-sums we have A000961, counted by A000005.
For run-lengths instead of run-sums we have A072774, counted by A047966.
These partitions are counted by A304442.
These are the positions of powers of primes in A353832.
The restriction to nonprimes is A353834.
For distinct instead of equal run-sums we have A353838, counted by A353837.
The version for compositions is A353848, counted by A353851.
A001222 counts prime factors, distinct A001221.
A005811 counts runs in binary expansion, distinct run-lengths A165413.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A300273 ranks collapsible partitions, counted by A275870.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353840-A353846 deal with iterated run-sums for partitions.
A353862 gives greatest run-sum of prime indices, least A353931.
A353866 ranks rucksack partitions, counted by A353864.
Cf. A007947, A071625, A073093, A181819, A238279, A304660, A323014, A333755, A353839.

Select[Range[100],SameQ@@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]*k]&]

A257993 Least gap in the partition having Heinz number n; index of the least prime not dividing n.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3

Offset: 1

Emeric Deutsch, May 18 2015

nonn

The "least gap" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (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.
In the Maple program the subprogram B yields the partition with Heinz number n.
Sum of least gaps of all partitions of m = A022567(m).
From Antti Karttunen, Aug 22 2016: (Start)
Index of the least prime not dividing n. (After a formula given by Heinz.)
Least k such that A002110(k) does not divide n.
One more than the number of trailing zeros in primorial base representation of n, A049345.
(End)
The least gap is also called the mex (minimal excludant) of the partition. - Gus Wiseman, Apr 20 2021

			a(18) = 3 because the partition having Heinz number 18 = 2*3*3 is [1,2,2], having least gap equal to 3.

G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
Miklós Bóna, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Alois P. Heinz, Table of n, a(n) for n = 1..20000
George E. Andrews and David Newman, Partitions and the Minimal Excludant, Annals of Combinatorics, Volume 23, May 2019, Pages 249-254.
P. J. Grabner and A. Knopfmacher, Analysis of some new partition statistics, Ramanujan J., 12, 2006, 439-454.
Brian Hopkins, James A. Sellers, and Dennis Stanton, Dyson's Crank and the Mex of Integer Partitions, arXiv:2009.10873 [math.CO], 2020.
Wikipedia, Mex (mathematics).
Index entries for sequences related to primorial base.

Cf. A215366, A002110, A022567, A049345, A053669, A055396, A064648, A276086, A276094, A276152.
Positions of 1's are A005408.
Positions of 2's are A047235.
The number of gaps is A079067.
The version for crank is A257989.
The triangle counting partitions by this statistic is A264401.
One more than A276084.
The version for greatest difference is A286469 or A286470.
A maximal instead of minimal version is A339662.
Positions of even terms are A342050.
Positions of odd terms are A342051.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339737 counts partitions by sum and greatest gap.
Cf. A001223, A001511, A005117, A026794, A029707, A072233, A079068, A098743, A124010, A279945, A325351, A325352.

with(numtheory): a := proc (n) local B, q: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: for q while member(q, B(n)) = true do  end do: q end proc: seq(a(n), n = 1 .. 150);
# second Maple program:
a:= n-> `if`(n=1, 1, (s-> min({$1..(max(s)+1)} minus s))(
        {map(x-> numtheory[pi](x[1]), ifactors(n)[2])[]})):
seq(a(n), n=1..100);  # Alois P. Heinz, May 09 2016
# faster:
A257993 := proc(n) local p, c; c := 1; p := 2;
while n mod p = 0 do p := nextprime(p); c := c + 1 od: c end:
seq(A257993(n), n=1..100); # Peter Luschny, Jun 04 2017

A053669[n_] := For[p = 2, True, p = NextPrime[p], If[CoprimeQ[p, n], Return[p]]]; a[n_] := PrimePi[A053669[n]]; Array[a, 100] (* Jean-François Alcover, Nov 28 2016 *)
Table[k = 1; While[! CoprimeQ[Prime@ k, n], k++]; k, {n, 100}] (* Michael De Vlieger, Jun 22 2017 *)

a(n) = forprime(p=2,, if (n % p, return(primepi(p)))); \\ Michel Marcus, Jun 22 2017

from sympy import nextprime, primepi
def a053669(n):
    p = 2
    while True:
        if n%p!=0: return p
        else: p=nextprime(p)
def a(n): return primepi(a053669(n)) # Indranil Ghosh, May 12 2017

(define (A257993 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) i (loop (/ (- n d) p) (+ 1 i))))))
;; Antti Karttunen, Aug 22 2016

a(n) = A000720(A053669(n)). - Alois P. Heinz, May 18 2015
From Antti Karttunen, Aug 22-30 2016: (Start)
a(n) = 1 + A276084(n).
a(n) = A055396(A276086(n)).
A276152(n) = A002110(a(n)).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} 1/A002110(k) = 1.705230... (1 + A064648). - Amiram Eldar, Jul 23 2022
a(n) << log n/log log n. - Charles R Greathouse IV, Dec 03 2022

A simpler description added to the name by Antti Karttunen, Aug 22 2016

A300273 Sorted list of Heinz numbers of collapsible integer partitions.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 36, 37, 40, 41, 43, 47, 48, 49, 53, 59, 61, 63, 64, 67, 71, 73, 79, 81, 83, 84, 89, 97, 101, 103, 107, 108, 109, 112, 113, 121, 125, 127, 128, 131, 137, 139, 144, 149, 151, 157, 163, 167, 169

Offset: 1

Gus Wiseman, Mar 01 2018

nonn

A positive integer is in this sequence iff it can be reduced to a prime number by a sequence of collapses, where a collapse is a replacement of prime(n)^k with prime(n*k) in a number's prime factorization (k > 1).

			A sequence of collapses is 84 -> 63 -> 49 -> 19 corresponding to the sequence of partitions (4211) -> (422) -> (44) -> (8). Hence 84 is in the sequence.

Cf. A000041, A001222, A002577, A005117, A018818, A056239, A094457, A112798, A145515, A215366, A275870, A289078, A289079, A291441, A296150, A299202.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
repcaps[q_]:=Union[{q},If[SquareFreeQ[q],{},Union@@repcaps/@Union[Times[q/#,Prime[Plus@@primeMS[#]]]&/@Select[Rest[Divisors[q]],!PrimeQ[#]&&PrimePowerQ[#]&]]]];
Select[Range[200],MemberQ[repcaps[#],_?PrimeQ]&]

A325781 Heinz numbers of complete integer partitions.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 30, 32, 36, 40, 42, 48, 54, 56, 60, 64, 72, 80, 84, 90, 96, 100, 108, 112, 120, 126, 128, 132, 140, 144, 150, 160, 162, 168, 176, 180, 192, 198, 200, 210, 216, 220, 224, 234, 240, 252, 256, 260, 264, 270, 280, 288, 294, 300

Offset: 1

Gus Wiseman, May 21 2019

nonn

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

The sum of prime indices of n is A056239(n). A number is in this sequence iff its divisors have sums of prime indices covering an initial interval of nonnegative integers. For example, the divisors of 60 are {1, 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, 60}, with respective sums of prime indices {0, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 7}, so 60 is in the sequence.

			The sequence of terms together with their prime indices begins:
     1: {}
     2: {1}
     4: {1,1}
     6: {1,2}
     8: {1,1,1}
    12: {1,1,2}
    16: {1,1,1,1}
    18: {1,2,2}
    20: {1,1,3}
    24: {1,1,1,2}
    30: {1,2,3}
    32: {1,1,1,1,1}
    36: {1,1,2,2}
    40: {1,1,1,3}
    42: {1,2,4}
    48: {1,1,1,1,2}
    54: {1,2,2,2}
    56: {1,1,1,4}
    60: {1,1,2,3}
    64: {1,1,1,1,1,1}

Cf. A002033, A056239, A103295, A112798, A126796, A188431, A299702, A325684, A325770, A325780, A325782.

normQ[m_]:=Or[m=={},Union[m]==Range[Max[m]]];
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Select[Range[1000],normQ[hwt/@Rest[Divisors[#]]]&]

cp's OEIS Frontend

A195017 If n = Product_{k >= 1} (p_k)^(c_k) where p_k is k-th prime and c_k >= 0 then a(n) = Sum_{k >= 1} c_k*((-1)^(k-1)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A320628 Products of primes of nonprime index.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A109082 Depth of rooted tree having Matula-Goebel number n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A344606 Number of alternating permutations of the prime factors of n, counting multiplicity, including twins (x,x).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A362611 Number of modes in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A182857 Smallest number that requires exactly n iterations to reach a fixed point under the x -> A181819(x) map.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A353833 Numbers whose multiset of prime indices has all equal run-sums.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs