A257993 -id:A257993 - cp's OEIS frontend

A181819 Prime shadow of n: a(1) = 1; for n>1, if n = Product prime(i)^e(i), then a(n) = Product prime(e(i)).

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

Offset: 1

Matthew Vandermast, Dec 07 2010

a(n) depends only on prime signature of n (cf. A025487). a(m) = a(n) iff m and n have the same prime signature, i.e., iff A046523(m) = A046523(n).

Because A046523 (the smallest representative of prime signature of n) and this sequence are functions of each other as A046523(n) = A181821(a(n)) and a(n) = a(A046523(n)), it implies that for all i, j: a(i) = a(j) <=> A046523(i) = A046523(j) <=> A101296(i) = A101296(j), i.e., that equivalence-class-wise this is equal to A101296, and furthermore, applying any function f on this sequence gives us a sequence b(n) = f(a(n)) whose equivalence class partitioning is equal to or coarser than that of A101296, i.e., b is then a sequence that depends only on the prime signature of n (the multiset of exponents of its prime factors), although not necessarily in a very intuitive way. - Antti Karttunen, Apr 28 2022

			20 = 2^2*5 has the exponents (2,1) in its prime factorization. Accordingly, a(20) = prime(2)*prime(1) = A000040(2)*A000040(1) = 3*2 = 6.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from exponents in factorization of n

Column 1 of A353510 (see also A325239 and A325277).
Cf. A000040, A000265, A001511, A001222, A003963, A005361, A007814, A008578, A028234, A046523, A056239, A064553, A064989, A067029, A101296 (restricted growth sequence transform), A108951, A122111, A124010, A124859, A156552, A181820, A181821, A182850, A182855, A182857 (also A323014), A115621, A101296, A238690, A238745, A238747, A238748, A246029, A304465, A304647, A305732, A305733, A320118, A323022, A325501, A325502, A325507, A325508, A325755 (A353566), A325756, A328830 [= a(a(n))], A328835, A351564 (characteristic function of A130091), A351944, A351946, A353379.
Left inverse of A304660.

a181819 = product . map a000040 . a124010_row
-- Reinhard Zumkeller, Mar 26 2012

A181819 := proc(n)
    local a;
    a := 1;
    for pf in ifactors(n)[2] do
        a := a*ithprime(pf[2]) ;
    end do:
    a ;
end proc:
seq(A181819(n),n=1..80) ; # R. J. Mathar, Jan 09 2019
# second Maple program:
a:= n-> mul(ithprime(i[2]), i=ifactors(n)[2]):
seq(a(n), n=1..105);  # Alois P. Heinz, Nov 26 2024

{1}~Join~Table[Times @@ Prime@ Map[Last, FactorInteger@ n], {n, 2, 120}] (* Michael De Vlieger, Feb 07 2016 *)

a(n) = {my(f=factor(n)); prod(k=1, #f~, prime(f[k,2]));} \\ Michel Marcus, Nov 16 2015

;; With memoization-macro definec.
(definec (A181819 n) (cond ((= 1 n) 1) (else (* (A000040 (A067029 n)) (A181819 (A028234 n))))))
;; Antti Karttunen, Feb 05 2016

;; With memoization-macro definec.
(definec (A181819 n) (cond ((= 1 n) 1) ((even? n) (* (A000040 (A007814 n)) (A181819 (A000265 n)))) (else (A181819 (A064989 n)))))
;; Antti Karttunen, Feb 07 2016

From Antti Karttunen, Feb 07 2016: (Start)
a(1) = 1; for n > 1, a(n) = A000040(A067029(n)) * a(A028234(n)).
a(1) = 1; for n > 1, a(n) = A008578(A001511(n)) * a(A064989(n)).
Other identities. For all n >= 1:
a(A124859(n)) = A122111(a(n)) = A238745(n). - from Matthew Vandermast's formulas for the latter sequence.
(End)
a(n) = A246029(A156552(n)). - Antti Karttunen, Oct 15 2016
From Antti Karttunen, Apr 28 & Apr 30 2022: (Start)
A181821(a(n)) = A046523(n) and a(A046523(n)) = a(n). [See comments]
a(n) = A329900(A124859(n)) = A319626(A124859(n)).
a(n) = A246029(A156552(n)).
a(a(n)) = A328830(n).
a(A304660(n)) = n.
a(A108951(n)) = A122111(n).
a(A185633(n)) = A322312(n).
a(A025487(n)) = A181820(n).
a(A276076(n)) = A275735(n) and a(A276086(n)) = A328835(n).
As the sequence converts prime exponents to prime indices, it effects the following mappings:
A001221(a(n)) = A071625(n). [Number of distinct indices --> Number of distinct exponents]
A001222(a(n)) = A001221(n). [Number of indices (i.e., the number of prime factors with multiplicity) --> Number of exponents (i.e., the number of distinct prime factors)]
A056239(a(n)) = A001222(n). [Sum of indices --> Sum of exponents]
A066328(a(n)) = A136565(n). [Sum of distinct indices --> Sum of distinct exponents]
A003963(a(n)) = A005361(n). [Product of indices --> Product of exponents]
A290103(a(n)) = A072411(n). [LCM of indices --> LCM of exponents]
A156061(a(n)) = A290107(n). [Product of distinct indices --> Product of distinct exponents]
A257993(a(n)) = A134193(n). [Index of the least prime not dividing n --> The least number not among the exponents]
A055396(a(n)) = A051904(n). [Index of the least prime dividing n --> Minimal exponent]
A061395(a(n)) = A051903(n). [Index of the greatest prime dividing n --> Maximal exponent]
A008966(a(n)) = A351564(n). [All indices are distinct (i.e., n is squarefree) --> All exponents are distinct]
A007814(a(n)) = A056169(n). [Number of occurrences of index 1 (i.e., the 2-adic valuation of n) --> Number of occurrences of exponent 1]
A056169(a(n)) = A136567(n). [Number of unitary prime divisors --> Number of exponents occurring only once]
A064989(a(n)) = a(A003557(n)) = A295879(n). [Indices decremented after <--> Exponents decremented before]
Other mappings:
A007947(a(n)) = a(A328400(n)) = A329601(n).
A181821(A007947(a(n))) = A328400(n).
A064553(a(n)) = A000005(n) and A000005(a(n)) = A182860(n).
A051903(a(n)) = A351946(n).
A003557(a(n)) = A351944(n).
A258851(a(n)) = A353379(n).
A008480(a(n)) = A309004(n).
a(A325501(n)) = A325507(n) and a(A325502(n)) = A038754(n+1).
a(n!) = A325508(n).
(End)

Name "Prime shadow" (coined by Gus Wiseman in A325755) prefixed to the definition by Antti Karttunen, Apr 27 2022

A053669 Smallest prime not dividing n.

Original entry on oeis.org

2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 7, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 7, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 7, 2, 3, 2, 3, 2, 5, 2, 3, 2, 3, 2, 5, 2, 3, 2

Offset: 1

Henry Bottomley, Feb 15 2000

Smallest prime coprime to n.
Smallest k >= 2 coprime to n.
a(#(p-1)) = a(A034386(p-1)) = p is the first appearance of prime p in sequence.
a(A005408(n)) = 2; for n > 2: a(n) = A112484(n,1). - Reinhard Zumkeller, Sep 23 2011
Average value is 2.920050977316134... = A249270. - Charles R Greathouse IV, Nov 02 2013
Differs from A236454, "smallest number not dividing n^2", for the first time at n=210, where a(210)=11 while A236454(210)=8. A235921 lists all n for which a(n) differs from A236454. - Antti Karttunen, Jan 26 2014
For k >= 0, a(A002110(k)) is the first occurrence of p = prime(k+1). Thereafter p occurs whenever A007947(n) = A002110(k). Thus every prime appears in this sequence infinitely many times. - David James Sycamore, Dec 04 2024

			a(60) = 7, since all primes smaller than 7 divide 60 but 7 does not.
a(90) = a(120) = a(150) = a(180) = 7 because 90,120,150,180 all have same squarefree kernel = 30 = A002110(3), and 7 is the smallest prime which does not divide 30. - _David James Sycamore_, Dec 04 2024

T. D. Noe, Table of n, a(n) for n = 1..10000
Dylan Fridman, Juli Garbulsky, Bruno Glecer, James Grime and Massi Tron Florentin, A Prime-Representing Constant, The American Mathematical Monthly, Vol. 126, No. 1 (2019), pp. 70-73; ResearchGate link.
James Grime and Brady Haran, 2.920050977316, Numberphile video (2020)
Igor Rivin, Geodesics with one self-intersection, and other stories, Advances in Mathematics, Vol. 231, No. 5 (2012), pp. 2391-2412; arXiv preprint, arXiv:0901.2543 [math.GT], 2009-2011.

One more than A071222(n-1).
Cf. A000040, A020639, A053670, A053671, A053672, A053673, A053674, A055874, A079578, A087560, A096014, A235921, A236454, A249270, A257993 (the indices of these primes), A276086, A324895, A353526, A353528, A353529, A358754, A358755, A370125.
Cf. A002110, A007947.

a053669 n = head $ dropWhile ((== 0) . (mod n)) a000040_list
-- Reinhard Zumkeller, Nov 11 2012

f:= proc(n) local p;
p:= 2;
while n mod p = 0 do p:= nextprime(p) od:
p
end proc:
map(f, [$1..100]); # Robert Israel, May 18 2016

Table[k := 1; While[Not[GCD[n, Prime[k]] == 1], k++ ]; Prime[k], {n, 1, 60}] (* Stefan Steinerberger, Apr 01 2006 *)
With[{prs=Prime[Range[10]]},Flatten[Table[Select[prs,!Divisible[ n,#]&,1],{n,110}]]] (* Harvey P. Dale, May 03 2012 *)

a(n)=forprime(p=2,,if(n%p,return(p))) \\ Charles R Greathouse IV, Nov 20 2012

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

# using standard library functions only
import math
def a(n):
    k = 2
    while math.gcd(n,k) > 1: k += 1
    return k # Ely Golden, Nov 26 2020

(define (A053669 n) (let loop ((i 1)) (cond ((zero? (modulo n (A000040 i))) (loop (+ i 1))) (else (A000040 i))))) ;; Antti Karttunen, Jan 26 2014

a(n) = A071222(n-1)+1. [Because the right hand side computes the smallest k >= 2 such that gcd(n,k) = gcd(n-1,k-1) which is equal to the smallest k >= 2 coprime to n] - Antti Karttunen, Jan 26 2014
a(n) = 1 + Sum_{k=1..n}(floor((n^k)/k!)-floor(((n^k)-1)/k!)) = 2 + Sum_{k=1..n} A001223(k)*( floor(n/A002110(k))-floor((n-1)/A002110(k)) ). - Anthony Browne, May 11 2016
a(n!) = A151800(n). - Anthony Browne, May 11 2016
a(2k+1) = 2. - Bernard Schott, Jun 03 2019
Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = A249270. - Amiram Eldar, Oct 29 2020
a(n) = A000040(A257993(n)) = A020639(A276086(n)) = A276086(n) / A324895(n). - Antti Karttunen, Apr 24 2022
a(n) << log n. For every e > 0, there is some N such that for all n > N, a(n) < (1 + e)*log n. - Charles R Greathouse IV, Dec 03 2022
A007947(n) = A002110(k) ==> a(n) = prime(k+1). - David James Sycamore, Dec 04 2024

More terms from Andrew Gacek (andrew(AT)dgi.net), Feb 21 2000 and James Sellers, Feb 22 2000

Entry revised by David W. Wilson, Nov 25 2006

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(n)).

Original entry on oeis.org

0, 1, 2, 3, 6, 7, 8, 9, 30, 31, 32, 33, 36, 37, 38, 39, 210, 211, 212, 213, 216, 217, 218, 219, 240, 241, 242, 243, 246, 247, 248, 249, 2310, 2311, 2312, 2313, 2316, 2317, 2318, 2319, 2340, 2341, 2342, 2343, 2346, 2347, 2348, 2349, 2520, 2521, 2522, 2523, 2526, 2527, 2528, 2529, 2550, 2551, 2552, 2553, 2556, 2557, 2558, 2559, 30030, 30031

Offset: 0

Antti Karttunen, Aug 24 2016

Numbers that are sums of distinct primorial numbers, A002110.

Numbers with no digits larger than one in primorial base, A049345.

Antti Karttunen, Table of n, a(n) for n = 0..8191
Index entries for sequences related to primorial base

Cf. A000040, A001511, A002110, A007088, A007814, A019565, A049345, A257993, A276084, A276085, A276154, A351073, A328461, A328473, A328474, A328571, A328831, A328836, A341518, A344591.
Complement of A177711.
Subsequences: A328233, A328832, A328462 (odd bisection).
Conjectured subsequences: A328110, A380527.
Fixed points of A328841, positions of zeros in A328828, A328842, and A329032, positions of ones in A328581, A328582, and A381032.
Positions of terms < 2 in A328114.
Indices where A327860 and A329029 coincide.
Cf. also table A328464 (and its rows).
Cf. also A059590, A283985, A290249, A342921.

nn = 65; b = MixedRadix[Reverse@ Prime@ Range[IntegerLength[nn, 2] - 1]]; Table[FromDigits[IntegerDigits[n, 2], b], {n, 0, 65}] (* Version 10.2, or *)
Table[Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ # - 1]], Reverse@ #}] &@ IntegerDigits[n, 2], {n, 0, 65}] (* Michael De Vlieger, Aug 26 2016 *)

A276156(n) = { my(s=0, p=1, r=1); while(n, if(n%2, s += r); n>>=1; p = nextprime(1+p); r *= p); (s); }; \\ Antti Karttunen, Feb 03 2022

from sympy import prime, primorial, primepi, factorint
from operator import mul
def a002110(n): return 1 if n<1 else primorial(n)
def a276085(n):
    f=factorint(n)
    return sum([f[i]*a002110(primepi(i) - 1) for i in f])
def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) # after Chai Wah Wu
def a(n): return 0 if n==0 else a276085(a019565(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 23 2017

a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1+A276154(a(n)).
Other identities. For all n >= 0:
a(n) = A276085(A019565(n)).
A049345(a(n)) = A007088(n).
A257993(a(n)) = A001511(n).
A276084(a(n)) = A007814(n).
A051903(a(n)) = A351073(n).

A264401 Triangle read by rows: T(n,k) is the number of partitions of n having least gap k.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 1, 1, 2, 2, 1, 2, 3, 2, 4, 4, 2, 1, 4, 6, 4, 1, 7, 8, 5, 2, 8, 11, 8, 3, 12, 15, 10, 4, 1, 14, 20, 15, 6, 1, 21, 26, 19, 9, 2, 24, 35, 27, 12, 3, 34, 45, 34, 17, 5, 41, 58, 47, 23, 6, 1, 55, 75, 59, 31, 10, 1, 66, 96, 79, 41, 13, 2

Offset: 0

Emeric Deutsch, Nov 21 2015

The "least gap" or "mex" 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.
Sum of entries in row n is A000041(n).
T(n,1) = A002865(n).
Sum_{k>=1} k*T(n,k) = A022567(n).

			Row n=5 is 2,3,2; indeed, the least gaps of [5], [4,1], [3,2], [3,1,1], [2,2,1], [2,1,1,1], and [1,1,1,1,1] are 1, 2, 1, 2, 3, 3, and 2, respectively (i.e., two 1s, three 2s, and two 3s).
Triangle begins:
   1
   0   1
   1   1
   1   1   1
   2   2   1
   2   3   2
   4   4   2   1
   4   6   4   1
   7   8   5   2
   8  11   8   3
  12  15  10   4   1
  14  20  15   6   1
  21  26  19   9   2

Alois P. Heinz, Rows n = 0..1000, flattened
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)

Row sums are A000041.
Row lengths are A002024.
Column k = 1 is A002865.
Column k = 2 is A027336.
The strict case is A343348.
A000009 counts strict partitions.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A257993 gives the least gap of the partition with Heinz number n.
A339564 counts factorizations with a selected factor.
A342050 ranks partitions with even least gap.
A342051 ranks partitions with odd least gap.
Cf. A003242, A022567, A048004, A083710, A098743, A130689, A338470, A343341.

g := (sum(t^j*x^((1/2)*j*(j-1))*(1-x^j), j = 1 .. 80))/(product(1-x^i, i = 1 .. 80)): gser := simplify(series(g, x = 0, 23)): for n from 0 to 30 do P[n] := sort(coeff(gser, x, n)) end do: for n from 0 to 25 do seq(coeff(P[n], t, j), j = 1 .. degree(P[n])) end do; # yields sequence in triangular form
# second Maple program:
b:= proc(n, i) option remember; `if`(n=0, `if`(i=0, [1, 0],
      [0, x]), `if`(i<1, 0, (p-> [0, p[2] +p[1]*x^i])(
      b(n, i-1)) +add(b(n-i*j, i-1), j=1..n/i)))
    end:
T:= n->(p->seq(coeff(p, x, i), i=1..degree(p)))(b(n, n+1)[2]):
seq(T(n), n=0..20);  # Alois P. Heinz, Nov 29 2015

Needs["Combinatorica`"]; {1, 0}~Join~Flatten[Table[Count[Map[If[# == {}, 0, First@ #] &@ Complement[Range@ n, #] &, Combinatorica`Partitions@ n], n_ /; n == k], {n, 17}, {k, n}] /. 0 -> Nothing] (* Michael De Vlieger, Nov 21 2015 *)
mingap[q_]:=Min@@Complement[Range[If[q=={},0,Max[q]]+1],q];Table[Length[Select[IntegerPartitions[n],mingap[#]==k&]],{n,0,15},{k,Round[Sqrt[2*(n+1)]]}] (* Gus Wiseman, Apr 19 2021 *)
b[n_, i_] := b[n, i] = If[n == 0, If[i == 0, {1, 0}, {0, x}], If[i<1, {0, 0}, {0, #[[2]] + #[[1]]*x^i}&[b[n, i-1]] + Sum[b[n-i*j, i - 1], {j, 1, n/i}]]];
T[n_] := CoefficientList[b[n, n + 1], x][[2]] // Rest;
T /@ Range[0, 20] // Flatten (* Jean-François Alcover, May 21 2021, after Alois P. Heinz *)

G.f.: G(t,x) = Sum_{j>=1} (t^j*x^{j(j-1)/2}*(1-x^j))/Product_{i>=1}(1-x^i).

A276084 a(n) = Number of trailing zeros in primorial base representation of n (A049345); largest k such that A002110(k) divides n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Aug 22 2016

Terms begin from a(1)=0 because for zero the count is ambiguous.
From Amiram Eldar, Mar 10 2021: (Start)
The asymptotic density of the occurrences of k is (prime(k+1)-1)/A002110(k+1).
The asymptotic mean of this sequence is Sum_{k>=1} 1/A002110(k) = 0.705230... (A064648). (End)

			For n=24, which is "400" in primorial base (as 24 = 4*(3*2*1) + 0*(2*1) + 0*1, see A049345), there are two trailing zeros, thus a(24) = 2.

Antti Karttunen, Table of n, a(n) for n = 1..2310
Index entries for sequences related to primorial base

Cf. A000040, A001221, A002110, A049345, A053589, A064648, A340346.
One less than A257993.
Differs from the related A230403 for the first time at n=24.

Table[If[# == 0, 0, j = #; While[! Divisible[n, Times @@ Prime@ Range@ j], j--]; j] &@ If[OddQ@ n, 0, k = 1; While[Times @@ Prime@ Range[k + 1] <= n, k++]; k], {n, 120}] (* or *)
nn = 120; b = MixedRadix[Reverse@ Prime@ Range@ PrimePi[nn + 1]]; Table[Length@ TakeWhile[Reverse@ IntegerDigits[n, b], # == 0 &], {n, nn}] (* Version 10.2, or *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Length@ TakeWhile[Reverse@ f@ n, # == 0 &], {n, 120}] (* Michael De Vlieger, Aug 30 2016 *)

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)) - 1 # Indranil Ghosh, May 12 2017

(define (A276084 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) (- i 1) (loop (/ (- n d) p) (+ 1 i))))))

a(n) = A257993(n)-1.
Other identities. For all n >= 1:
A053589(n) = A002110(a(n)).
a(n) = A001221(A053589(n)) = A001221(A340346(n)). - Peter Munn, Jan 14 2021

A342050 Numbers k which have an odd number of trailing zeros in their primorial base representation A049345(k).

Original entry on oeis.org

2, 4, 8, 10, 14, 16, 20, 22, 26, 28, 30, 32, 34, 38, 40, 44, 46, 50, 52, 56, 58, 60, 62, 64, 68, 70, 74, 76, 80, 82, 86, 88, 90, 92, 94, 98, 100, 104, 106, 110, 112, 116, 118, 120, 122, 124, 128, 130, 134, 136, 140, 142, 146, 148, 150, 152, 154, 158, 160, 164, 166, 170, 172, 176, 178, 180, 182, 184, 188, 190, 194, 196, 200, 202, 206, 208, 212

Offset: 1

Amiram Eldar, Feb 26 2021

Numbers k such that A276084(k) is odd.
All the terms are even since odd numbers have 0 trailing zeros, and 0 is not odd.
The number of terms not exceeding A002110(m) for m>=1 is A002110(m) * Sum_{k=1..m}(-1)^k/A002110(k) = 1, 2, 11, 76, 837, 10880, 184961, ...
The asymptotic density of this sequence is Sum_{k>=1} (-1)^(k+1)/A002110(k) = 0.362306... (A132120).
Also Heinz numbers of partitions with even least gap. The least gap (mex or minimal excludant) of a partition is the least positive integer that is not a part. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions. - Gus Wiseman, Apr 23 2021
Numbers k such that A000720(A053669(k)) is even. Differences from the related A353531 seem to be terms that are multiples of 210, but not all of them, for example primorial 30030 (= 143*210) is in neither sequence. Consider also A038698. - Antti Karttunen, Apr 25 2022

			2 is a term since A049345(2) = 10 has 1 trailing zero.
4 is a term since A049345(2) = 20 has 1 trailing zero.
30 is a term since A049345(2) = 1000 has 3 trailing zeros.
From _Gus Wiseman_, Apr 23 2021: (Start)
The sequence of terms together with their prime indices begins:
      2: {1}             46: {1,9}             90: {1,2,2,3}
      4: {1,1}           50: {1,3,3}           92: {1,1,9}
      8: {1,1,1}         52: {1,1,6}           94: {1,15}
     10: {1,3}           56: {1,1,1,4}         98: {1,4,4}
     14: {1,4}           58: {1,10}           100: {1,1,3,3}
     16: {1,1,1,1}       60: {1,1,2,3}        104: {1,1,1,6}
     20: {1,1,3}         62: {1,11}           106: {1,16}
     22: {1,5}           64: {1,1,1,1,1,1}    110: {1,3,5}
     26: {1,6}           68: {1,1,7}          112: {1,1,1,1,4}
     28: {1,1,4}         70: {1,3,4}          116: {1,1,10}
     30: {1,2,3}         74: {1,12}           118: {1,17}
     32: {1,1,1,1,1}     76: {1,1,8}          120: {1,1,1,2,3}
     34: {1,7}           80: {1,1,1,1,3}      122: {1,18}
     38: {1,8}           82: {1,13}           124: {1,1,11}
     40: {1,1,1,3}       86: {1,14}           128: {1,1,1,1,1,1,1}
     44: {1,1,5}         88: {1,1,1,5}        130: {1,3,6}
(End)

Amiram Eldar, Table of n, a(n) for n = 1..10000
George E. Andrews and David Newman, Partitions and the Minimal Excludant, Annals of Combinatorics, Volume 23, May 2019, Pages 249-254.
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 computed from indices in prime factorization
Index entries for sequences related to primorial base

Cf. A002110, A049345, A053669, A132120, A276084.
Complement of A342051.
A099800 is subsequence.
Analogous sequences: A001950 (Zeckendorf representation), A036554 (binary), A145204 (ternary), A217319 (base 4), A232745 (factorial base).
The version for reversed binary expansion is A079523.
Positions of even terms in A257993.
A000070 counts partitions with a selected part.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A079067 counts gaps in prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339662 gives greatest gap in prime indices.
Cf. A000720, A001223, A005408, A026794, A029707, A038698, A047235, A079068, A121539, A286469, A286470, A325351, A353525 (characteristic function).
Differs from A353531 for the first time at n=77, where a(77) = 212, as this sequence misses A353531(77) = 210.

seq[max_] := Module[{bases = Prime@Range[max, 1, -1], nmax}, nmax = Times @@ bases - 1; Select[Range[nmax], OddQ @ LengthWhile[Reverse @ IntegerDigits[#, MixedRadix[bases]], #1 == 0 &] &]]; seq[4]
Select[Range[100],EvenQ[Min@@Complement[Range[PrimeNu[#]+1],PrimePi/@First/@FactorInteger[#]]]&] (* Gus Wiseman, Apr 23 2021 *)

A353525(n) = { for(i=1,oo,if(n%prime(i),return((i+1)%2))); }
isA342050(n) = A353525(n);
k=0; n=0; while(k<77, n++; if(isA342050(n), k++; print1(n,", "))); \\ Antti Karttunen, Apr 25 2022

More terms added (to differentiate from A353531) by Antti Karttunen, Apr 25 2022

A342193 Numbers with no prime index dividing all the other prime indices.

Original entry on oeis.org

1, 15, 33, 35, 45, 51, 55, 69, 75, 77, 85, 91, 93, 95, 99, 105, 119, 123, 135, 141, 143, 145, 153, 155, 161, 165, 175, 177, 187, 195, 201, 203, 205, 207, 209, 215, 217, 219, 221, 225, 231, 245, 247, 249, 253, 255, 265, 275, 279, 285, 287, 291, 295, 297, 299

Offset: 1

Gus Wiseman, Apr 11 2021

nonn

Alternative name: 1 and numbers with smallest prime index not dividing all the other prime indices.
First differs from A339562 in having 45.
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.
Also 1 and Heinz numbers of integer partitions with smallest part not dividing all the others (counted by A338470). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

			The sequence of terms together with their prime indices begins:
      1: {}         105: {2,3,4}      201: {2,19}
     15: {2,3}      119: {4,7}        203: {4,10}
     33: {2,5}      123: {2,13}       205: {3,13}
     35: {3,4}      135: {2,2,2,3}    207: {2,2,9}
     45: {2,2,3}    141: {2,15}       209: {5,8}
     51: {2,7}      143: {5,6}        215: {3,14}
     55: {3,5}      145: {3,10}       217: {4,11}
     69: {2,9}      153: {2,2,7}      219: {2,21}
     75: {2,3,3}    155: {3,11}       221: {6,7}
     77: {4,5}      161: {4,9}        225: {2,2,3,3}
     85: {3,7}      165: {2,3,5}      231: {2,4,5}
     91: {4,6}      175: {3,3,4}      245: {3,4,4}
     93: {2,11}     177: {2,17}       247: {6,8}
     95: {3,8}      187: {5,7}        249: {2,23}
     99: {2,2,5}    195: {2,3,6}      253: {5,9}

The complement is counted by A083710 (strict: A097986).
The complement with no 1's is A083711 (strict: A098965).
These partitions are counted by A338470 (strict: A341450).
The squarefree case is A339562, with squarefree complement A339563.
The case with maximum prime index not divisible by all others is A343338.
The case with maximum prime index divisible by all others is A343339.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001221 counts distinct prime factors.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices, row sums of A112798.
A299702 lists Heinz numbers of knapsack partitions.
A339564 counts factorizations with a selected factor.
Cf. A066637, A072774, A098743, A253249, A264401, A257993, A342050, A342051, A343344.

Select[Range[100],#==1||With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(p/Min@@p)]&]

A343337 Numbers with no prime index divisible by all the other prime indices.

Original entry on oeis.org

1, 15, 30, 33, 35, 45, 51, 55, 60, 66, 69, 70, 75, 77, 85, 90, 91, 93, 95, 99, 102, 105, 110, 119, 120, 123, 132, 135, 138, 140, 141, 143, 145, 150, 153, 154, 155, 161, 165, 170, 175, 177, 180, 182, 186, 187, 190, 198, 201, 203, 204, 205, 207, 209, 210, 215

Offset: 1

Gus Wiseman, Apr 13 2021

nonn

Alternative name: 1 and numbers whose greatest prime index is not divisible by all the other prime indices.
First differs from A318992 in lacking 195.
First differs from A343343 in lacking 195.
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.
Also Heinz numbers of partitions with greatest part not divisible by all the others (counted by A343341). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

			The sequence of terms together with their prime indices begins:
      1: {}            90: {1,2,2,3}      141: {2,15}
     15: {2,3}         91: {4,6}          143: {5,6}
     30: {1,2,3}       93: {2,11}         145: {3,10}
     33: {2,5}         95: {3,8}          150: {1,2,3,3}
     35: {3,4}         99: {2,2,5}        153: {2,2,7}
     45: {2,2,3}      102: {1,2,7}        154: {1,4,5}
     51: {2,7}        105: {2,3,4}        155: {3,11}
     55: {3,5}        110: {1,3,5}        161: {4,9}
     60: {1,1,2,3}    119: {4,7}          165: {2,3,5}
     66: {1,2,5}      120: {1,1,1,2,3}    170: {1,3,7}
     69: {2,9}        123: {2,13}         175: {3,3,4}
     70: {1,3,4}      132: {1,1,2,5}      177: {2,17}
     75: {2,3,3}      135: {2,2,2,3}      180: {1,1,2,2,3}
     77: {4,5}        138: {1,2,9}        182: {1,4,6}
     85: {3,7}        140: {1,1,3,4}      186: {1,2,11}
For example, 195 has prime indices {2,3,6}, and 6 is divisible by both 2 and 3, so 195 does not belong to the sequence.

The complement is counted by A130689.
The dual version is A342193.
The case with smallest prime index not dividing all the others is A343338.
The case with smallest prime index dividing by all the others is A343340.
These are the Heinz numbers of the partitions counted by A343341.
Including the dual version gives A343343.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A067824 counts strict chains of divisors starting with n.
A253249 counts strict chains of divisors.
A339564 counts factorizations with a selected factor.
Cf. A083710, A257993, A338470, A341450, A343342, A343345, A343346, A343377, A343379, A343381, A343382.

Select[Range[1000],#==1||With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(Max@@p/p)]&]

A339563 Squarefree numbers > 1 whose smallest prime index divides all the other prime indices.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 30, 31, 34, 37, 38, 39, 41, 42, 43, 46, 47, 53, 57, 58, 59, 61, 62, 65, 66, 67, 70, 71, 73, 74, 78, 79, 82, 83, 86, 87, 89, 94, 97, 101, 102, 103, 106, 107, 109, 110, 111, 113, 114, 115, 118, 122, 127

Offset: 1

Gus Wiseman, Apr 10 2021

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.

Also Heinz numbers of strict integer partitions whose smallest part divides all the others (counted by A097986). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

			The sequence of terms together with their prime indices begins:
      2: {1}       29: {10}        59: {17}
      3: {2}       30: {1,2,3}     61: {18}
      5: {3}       31: {11}        62: {1,11}
      6: {1,2}     34: {1,7}       65: {3,6}
      7: {4}       37: {12}        66: {1,2,5}
     10: {1,3}     38: {1,8}       67: {19}
     11: {5}       39: {2,6}       70: {1,3,4}
     13: {6}       41: {13}        71: {20}
     14: {1,4}     42: {1,2,4}     73: {21}
     17: {7}       43: {14}        74: {1,12}
     19: {8}       46: {1,9}       78: {1,2,6}
     21: {2,4}     47: {15}        79: {22}
     22: {1,5}     53: {16}        82: {1,13}
     23: {9}       57: {2,8}       83: {23}
     26: {1,6}     58: {1,10}      86: {1,14}

These partitions are counted by A097986 (non-strict: A083710).
The case with no 1's is counted by A098965 (non-strict: A083711).
The squarefree complement is A339562, ranked by A341450.
The complement of the not necessarily squarefree version is A342193.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001055 counts factorizations.
A001221 counts distinct prime factors.
A005117 lists squarefree numbers.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices, row sums of A112798.
A338470 counts partitions with no dividing part.
Cf. A130714, A253249, A257993, A299702, A339564, A343340, A343378.

Select[Range[2,100],SquareFreeQ[#]&&With[{p=PrimePi/@First/@FactorInteger[#]},And@@IntegerQ/@(p/Min@@p)]&]

A276088 The least significant nonzero digit in primorial base representation of n: a(n) = A276094(n) / A002110(A276084(n)) (with a(0) = 0).

Original entry on oeis.org

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

Offset: 0

Antti Karttunen, Aug 22 2016

For any n >= 1, start from k = n and repeatedly try to divide as many successive primes as possible out of k, iterating as k/2 -> k, k/3 -> k, k/5 -> k, until a nonzero remainder is encountered, which is then the value of a(n). (See the last example).
Note that the sequence has been defined so that it will eventually include also "digits" (actually: value holders) > 9 that occur as the least significant nonzero digits in primorial base representation. Thus any eventual decimal corruption of A049345 will not affect these values.
The sums of the first 10^k terms (starting from n=1), for k = 1, 2, ..., are 12, 138, 1441, 14565, 145950, 1459992, 14600211, 146002438, 1460025336, 14600254674, ... . Apparently, the asymptotic mean of this sequence is limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1.460025... . - Amiram Eldar, Sep 10 2022
The asymptotic density of the occurrences of k = 1, 2, ..., is d(1) = c = A064648 for k = 1, and d(k) = c - Sum_{i = 1..PrimePi(k)} 1/primorial(i) for k >= 2, where primorial(i) = A002110(i). The asymptotic mean of this sequence is Sum_{k>=1} k * d(k) = c + Sum_{i>=1} (prime(i) * gap(i) + t(gap(i)-1)) * (c - Sum_{j = 1..i} 1/primorial(i)) = 1.460025488658067356046281458556..., where t(i) = A000217(i) and gap(i) = A001223(i). - Amiram Eldar, Feb 20 2025

			   n   A049345  the rightmost nonzero = a(n)
---------------------------------------------------------
   0       0             0
   1       1             1
   2      10             1
   3      11             1
   4      20             2
   5      21             1
   6     100             1
   7     101             1
   8     110             1
   9     111             1
  10     120             2
  11     121             1
  12     200             2
  13     201             1
  14     210             1
  15     211             1
  16     220             2
.
For n=48 according to the iteration interpretation, we obtain first 48/2 = 24, and the remainder is zero, so we continue: 24/3 = 8 and here the remainder is zero as well, so we try next 8/5, but this gives the nonzero remainder 3, thus a(48)=3.
For n=2100, which could be written "A0000" in primorial base (where A stands for digit "ten", as 2100 = 10*A002110(4)), the least significant nonzero value holder (also the most significant) is thus 10 and a(2100) = 10. (The first point where this sequence attains a value larger than 9).

Antti Karttunen, Table of n, a(n) for n = 0..2310
Index entries for sequences related to primorial base.

Cf. A000040, A002110, A049345, A067029, A276084, A276086, A276094, A328569.

Cf. A000217, A000720, A001223, A064648.

nn = 120; b = MixedRadix[Reverse@ Prime@ Range@ PrimePi[nn + 1]]; Table[Last[IntegerDigits[n, b] /. 0 -> Nothing, 0], {n, 0, nn}] (* Version 11, or *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; {0}~Join~Table[Last@ DeleteCases[f@ n, d_ /; d == 0], {n, 120}] (* Michael De Vlieger, Aug 30 2016 *)

A276088(n) = { my(e=0, p=2); while(n && !(e=(n%p)), n = n/p; p = nextprime(1+p)); (e); }; \\ Antti Karttunen, Oct 29 2019

from sympy import nextprime, primepi, primorial
def a053669(n):
    p = 2
    while True:
        if n%p!=0: return p
        else: p=nextprime(p)
def a257993(n): return primepi(a053669(n))
def a002110(n): return 1 if n<1 else primorial(n)
def a276094(n): return 0 if n==0 else n%a002110(a257993(n))
def a(n): return 0 if n==0 else a276094(n)//a002110(a257993(n) - 1)
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 22 2017

(define (A276088 n) (if (zero? n) n (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) d (loop (/ (- n d) p) (+ 1 i)))))))

(define (A276088 n) (if (zero? n) n (/ (A276094 n) (A002110 (A276084 n)))))

a(0) = 0, and for n >= 1, a(n) = A276094(n) / A002110(A276084(n)).
From Antti Karttunen, Oct 29 2019: (Start)
a(n) = A067029(A276086(n)).
a(A276086(n)) = A328569(n).
(End).

cp's OEIS Frontend

A181819 Prime shadow of n: a(1) = 1; for n>1, if n = Product prime(i)^e(i), then a(n) = Product prime(e(i)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Scheme

Scheme

Formula

Extensions

A053669 Smallest prime not dividing n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Python

Scheme

Formula

Extensions

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A264401 Triangle read by rows: T(n,k) is the number of partitions of n having least gap k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A276084 a(n) = Number of trailing zeros in primorial base representation of n (A049345); largest k such that A002110(k) divides n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Scheme

Formula

A342050 Numbers k which have an odd number of trailing zeros in their primorial base representation A049345(k).

Views

Author

Keywords

Comments

Examples

Links