A000040 -id:A000040 - cp's OEIS frontend

A070826 One half of product of first n primes A000040.

1, 3, 15, 105, 1155, 15015, 255255, 4849845, 111546435, 3234846615, 100280245065, 3710369067405, 152125131763605, 6541380665835015, 307444891294245705, 16294579238595022365, 961380175077106319535, 58644190679703485491635, 3929160775540133527939545, 278970415063349480483707695

Offset: 1

Wolfdieter Lang, May 10 2002

Also, with offset 0, product of first n odd primes. - N. J. A. Sloane, Feb 26 2017
Identical to A002110(n)/2, n>=1.
a(n+1) is the least odd number with exactly n distinct prime divisors. - Labos Elemer, Mar 24 2003
Also, odd numbers n for which sigma(n)*phi(n)/n^2 reaches a new record low, monotonically decreasing to the lower bound 8/Pi^2. - M. F. Hasler, Jul 08 2025

Vincenzo Librandi, Table of n, a(n) for n = 1..200

Cf. A003266 (for Fibonacci), A070825 (for Lucas), A003046 (for Catalan).
Cf. also A002110, A024451, A060389, A091852, A276086, A203008 [= A003415(a(1+n))].
Range of A196529.

a:=n->mul(ithprime(j), j=2..n):seq(a(n), n=1..17); # Zerinvary Lajos, Aug 24 2008

Rest[ FoldList[ Times, 1, Prime[ Range[ 18]] ]]/2 (* Robert G. Wilson v, Feb 17 2004 *)
FoldList[Times, 1, Prime[Range[2, 18]]] (* Zak Seidov, Jan 26 2009 *)

a(n) = prod(k=2, n, prime(k)) \\ Michel Marcus, Mar 25 2017, simplified by M. F. Hasler, Jul 09 2025

from sympy import primorial
def A070826(n): return primorial(n)>>1 # Chai Wah Wu, Jul 21 2022

a(n) = A002110(n)/2.
From Antti Karttunen, Feb 06 2024: (Start)
a(1) = 1, and for n > 1, a(n) = A276086(A060389(n-1)).
a(n) = A024451(n) - 2*A203008(n-1).
(End)
a(n) = A000040(n)*a(n-1) for n > 1, a(1) = 1. - M. F. Hasler, Jul 09 2025

Formula corrected by Gary Detlefs, Dec 07 2011

A243055 Difference between the indices of the smallest and the largest prime dividing n: If n = p_i * ... * p_k, where p_i <= ... <= p_k, where p_h = A000040(h), then a(n) = (k-i), a(1) = 0 by convention.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 3, 1, 0, 0, 1, 0, 2, 2, 4, 0, 1, 0, 5, 0, 3, 0, 2, 0, 0, 3, 6, 1, 1, 0, 7, 4, 2, 0, 3, 0, 4, 1, 8, 0, 1, 0, 2, 5, 5, 0, 1, 2, 3, 6, 9, 0, 2, 0, 10, 2, 0, 3, 4, 0, 6, 7, 3, 0, 1, 0, 11, 1, 7, 1, 5, 0, 2, 0, 12, 0, 3, 4, 13, 8, 4, 0, 2

Offset: 1

Antti Karttunen, May 31 2014

nonn

For n>=1, A100484(n+1) gives the position where n occurs for the first time (setting also the records for the sequence).

a(n) = the difference between the largest and the smallest parts of the partition having Heinz number n. 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. Example: a(57) = 6; indeed, the partition having Heinz number 57 = 3*19 is [2, 8]. - Emeric Deutsch, Jun 04 2015

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from indices in prime factorization

Differs from A242411 for the first time at n=30.
A000961 gives the positions of zeros.
Cf. A243056, A241917, A241919, A049084, A027748, A055396, A061395, A100484, A215366.

with(numtheory):
a:= n-> `if`(n=1, 0, (f-> pi(max(f[]))-pi(min(f[])))(factorset(n))):
seq(a(n), n=1..100);  # Alois P. Heinz, Jun 04 2015

a[1]=0; a[n_] := Function[{f}, PrimePi[Max[f]] - PrimePi[Min[f]]][FactorInteger[n][[All, 1]]]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Jul 29 2015, after Alois P. Heinz *)

from sympy import primepi, primefactors
def A243055(n): return primepi(max(p:=primefactors(n),default=0))-primepi(min(p,default=0)) # Chai Wah Wu, Oct 10 2023

(define (A243055 n) (- (A061395 n) (A055396 n)))

If n = p_i * ... * p_k, where p_i <= ... <= p_k are not necessarily distinct primes (sorted into nondescending order) in the prime factorization of n, where p_i = A000040(i), then a(n) = (k-i).

a(n) = A061395(n) - A055396(n).

A246282 Numbers k for which A003961(k) > 2k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2n, where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 24, 27, 28, 30, 32, 35, 36, 39, 40, 42, 44, 45, 48, 49, 50, 52, 54, 56, 57, 60, 63, 64, 66, 68, 69, 70, 72, 75, 76, 78, 80, 81, 84, 88, 90, 91, 92, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 116, 117, 120, 124, 125, 126, 128, 130, 132, 135, 136, 138, 140, 144

Offset: 1

Antti Karttunen, Aug 24 2014

Numbers n such that A003961(n) > 2*n.
Numbers n such that A048673(n) > n.
The sequence grows as:
        a(10) = 18
       a(100) = 192
      a(1000) = 1830
     a(10000) = 18636
    a(100000) = 187350
   a(1000000) = 1865226
  a(10000000) = 18654333
and the powers of 10 occur at:
        a(5) = 10
       a(53) = 100
      a(536) = 1000
     a(5423) = 10000
    a(53290) = 100000
   a(535797) = 1000000
  a(5361886) = 10000000
suggesting that the ratio a(n)/n is converging to an constant and an arbitrary natural number is slightly more likely to be in this sequence than in the complement A246281. See also comments at A246351 and compare to quite a different ratio present in the "inverse" case A246362.
From Antti Karttunen, Aug 27 2020: (Start)
Any perfect number, including all odd perfect numbers (if such numbers exist), must occur in this sequence. See A286385 and A326042 for the reason why.
Like abundancy index (ratio A000203(n)/n), also ratio A003961(n)/n is multiplicative and always > 1 for all n > 1. Thus if the number has a proper divisor that is in this sequence, then the number itself also is. See A337372 for terms included here, but with no proper divisor in this sequence. (End)
For k >= 2, if m * A130789(k) is a term then m * A130789(k-1) is a term. - Peter Munn, Sep 01 2025
Could be called "primeshift-abundant numbers", in analogy with A005101. - Antti Karttunen, Sep 01 2025

			3 = p_2 (3 is the second prime, A000040(2)) is not a member, because p_3 = 5 (5 is the next prime after 3, A000040(3)) and 5/3 < 2.
4 = 2*2 = p_1 * p_1 is a member, as p_2 * p_2 = 3*3 = 9, and 9/4 > 2.
33 = 3*11 = p_2 * p_5 is not a member, as p_3 * p_6 = 5*13 = 65, and 65/33 < 2.
35 = 5*7 = p_3 * p_4 is a member, as p_4 * p_5 = 7*11 = 77, and 77/35 > 2.

Complement: A246281.
Setwise difference of A246352 and A048674.
Cf. A000040, A003961, A048673, A130789, A246362, A252742 (characteristic function), A286385, A326042, A337345.
Positions of positive terms in A252748 and in A337345.
Union of A337372 (primitive terms), A341610 (non-primitive terms).
Subsequences: A000396, (A005101 U A341614), (A023196 U A341615), A326134, A337373, A337374, A337378, A337381, A337384, A337386, A337543, A341611.
Cf. also A275717, A275718.

Select[Range[144], 2 # < Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1] &] (* Michael De Vlieger, Feb 22 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
isA246282(n) = (A003961(n) > (n+n));
n = 0; i = 0; while(i < 10000, n++; if(isA246282(n), i++; write("b246282.txt", i, " ", n)));

;; With Antti Karttunen's IntSeq-library.
(define A246282 (MATCHING-POS 1 1 (lambda (n) (> (A003961 n) (* 2 n)))))

;; With Antti Karttunen's IntSeq-library.
(define A246282 (MATCHING-POS 1 1 (lambda (n) (> (A048673 n) n))))

A new shorter version of name prepended by Antti Karttunen, Aug 27 2020

A227413 a(1)=1, a(2n)=nthprime(a(n)), a(2n+1)=nthcomposite(a(n)), where nthprime = A000040, nthcomposite = A002808.

Original entry on oeis.org

1, 2, 4, 3, 6, 7, 9, 5, 8, 13, 12, 17, 14, 23, 16, 11, 10, 19, 15, 41, 22, 37, 21, 59, 27, 43, 24, 83, 35, 53, 26, 31, 20, 29, 18, 67, 30, 47, 25, 179, 58, 79, 34, 157, 54, 73, 33, 277, 82, 103, 40, 191, 62, 89, 36, 431, 114, 149, 51, 241, 75, 101, 39, 127, 46

Offset: 1

Antti Karttunen, Jul 10 2013

Inverse permutation of A135141.

Shares with A073846 the property that the other bisection consists of just primes and the other bisection of just nonprimes.

Antti Karttunen and Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
Index entries for sequences that are permutations of the natural numbers

Similarly constructed permutations: A227402, A227404, A227410, A227412. Cf. also A073846, A209636.

import Data.List (transpose)
a227413 n = a227413_list !! (n-1)
a227413_list = 1 : concat (transpose [map a000040 a227413_list,
                                      map a002808 a227413_list])
-- Reinhard Zumkeller, Jan 29 2014

a(1)=1, a(2n) = A000040(a(n)), a(2n+1) = A002808(a(n)).

A007097(n) = a(A000079(n)).

A034602 Wolstenholme quotient W_p = (binomial(2p-1,p) - 1)/p^3 for prime p=A000040(n).

Original entry on oeis.org

1, 5, 265, 2367, 237493, 2576561, 338350897, 616410400171, 7811559753873, 17236200860123055, 3081677433937346539, 41741941495866750557, 7829195555633964779233, 21066131970056662377432067, 59296957594629000880904587621, 844326030443651782154010715715

Offset: 3

N. J. A. Sloane

nonn

Equivalently, (binomial(2p,p)-2)/(2*p^3) where p runs through the primes >=5.

The values of this sequence's terms are replicated by conjectured general formula, given in A223886 (and also added to the formula section here) for k=2, j=1 and n>=3. - Alexander R. Povolotsky, Apr 18 2013

			Binomial(10,5)-2 = 250; 5^3=125 hence a(5)=1.

R. K. Guy, Unsolved Problems in Number Theory, Sect. B31.

Robert Israel, Table of n, a(n) for n = 3..263
R. R. Aidagulov and M. A. Alekseyev. On p-adic approximation of sums of binomial coefficients. Journal of Mathematical Sciences 233:5 (2018), 626-634. doi:10.1007/s10958-018-3948-0 arXiv:1602.02632
R. J. McIntosh, On the converse of Wolstenholme's theorem, Acta Arithmetica 71:4 (1995), 381-389.
Romeo Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.
Jonathan Sondow, Extending Babbage's (non-)primality tests, in Combinatorial and Additive Number Theory II, Springer Proc. in Math. & Stat., Vol. 220, 269-277, CANT 2015 and 2016, New York, 2017; arXiv:1812.07650 [math.NT], 2018.

Cf. A177783 (alternative definition of Wolstenholme quotient), A072984, A092101, A092103, A092193, A128673, A217772, A223886, A263882.

Cf. A268512, A268589, A268590.

[(Binomial(2*p-1,p)-1) div p^3: p in PrimesInInterval(4,100)]; // Vincenzo Librandi, Nov 23 2015

f:= proc(n) local p;
p:= ithprime(n);
(binomial(2*p-1,p)-1)/p^3
end proc:
map(f, [$3..30]); # Robert Israel, Dec 19 2018

Table[(Binomial[2 Prime[n] - 1, Prime[n] - 1] - 1)/Prime[n]^3, {n, 3, 20}] (* Vincenzo Librandi, Nov 23 2015 *)

a(n) = (A088218(p)-1)/p^3 = (A001700(p-1)-1)/p^3 = (A000984(p)-2)/(2*p^3), where p=A000040(n).
a(n) = A087754(n)/2.
a(n) = (binomial(j*k*prime(n), j*prime(n)) - binomial(k*j, j)) / (k*prime(n)^3) for k=2, j=1, and n>=3. - Alexander R. Povolotsky, Apr 18 2013
a(n) = A263882(n)/prime(n) for n > 2. - Jonathan Sondow, Nov 23 2015
a(n) = numerator(tanh(Sum_{k=1..p-1} artanh(k/p)))/p^3, where p = prime(n) for n >= 3. - Thomas Ordowski, Apr 17 2025

Edited by Max Alekseyev, May 14 2010

More terms from Vincenzo Librandi, Nov 23 2015

A154945 Decimal expansion of Sum_{p} 1/(p^2-1), summed over the primes p = A000040.

Original entry on oeis.org

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

Offset: 0

R. J. Mathar, Jan 17 2009

By geometric series expansion, the same as the sum over the prime zeta function at even arguments, P(2i), i=1,2,....

(Pi^2/6)*density of A190641, the numbers divisible by exactly one prime with exponent greater than 1. - Charles R Greathouse IV, Aug 02 2016

			0.551693297656999184439731023971343578131500377778628252230...

Jacques Grah, Comportement moyen du cardinal de certains ensembles de facteurs premiers, Monatsh. Math., Vol. 118 (1994), pp. 91-109, Corollary 6.1.
Carl Pomerance and Andrzej Schinzel, Multiplicative Properties of Sets of Residues, Moscow Journal of Combinatorics and Number Theory, Vol. 1, Iss. 1 (2011), pp. 52-66. See p. 61.
Index to constants which are prime zeta sums {0,1,1}

Cf. A000040, A000961, A056798, A084920, A085548, A085964, A085966, A085968, A152447, A154932, A190641.

digits = 105; m0 = 2 digits; Clear[rd]; rd[m_] := rd[m] = RealDigits[delta1 = Sum[PrimeZetaP[2n], {n, 1, m}] , 10, digits][[1]]; rd[m0]; rd[m = 2m0];
While[rd[m] != rd[m-m0], Print[m]; m = m+m0]; Print[N[delta1, digits]]; rd[m] (* Jean-François Alcover, Sep 11 2015, updated Mar 16 2019 *)

eps()=2.>>bitprecision(1.)
primezeta(s)=my(t=s*log(2)); sum(k=1, lambertw(t/eps())\t, moebius(k)/k*log(abs(zeta(k*s))))
sumpos(n=1,primezeta(2*n)) \\ Charles R Greathouse IV, Aug 02 2016

sumeulerrat(1/(p^2-1)) \\ Amiram Eldar, Mar 18 2021

Equals Sum_{k>=1} 1/A084920(k) = Sum_{i>=1} P(2i) = A085548+A085964+A085966+A085968+... = A152447+A085548-A154932.
Equals Sum_{k>=2} 1/A000961(k)^2 = Sum_{k>=2} 1/A056798(k). - Amiram Eldar, Sep 21 2020
Equals (A136141 + A179119)/2. - Artur Jasinski, Mar 31 2025

More digits from Jean-François Alcover, Sep 11 2015

A276075 a(1) = 0, a(n) = (e1i1! + e2i2! + ... + eziz!) for n = prime(i1)^e1 prime(i2)^e2 * ... * prime(iz)^ez, where prime(k) is the k-th prime, A000040(k).

Original entry on oeis.org

0, 1, 2, 2, 6, 3, 24, 3, 4, 7, 120, 4, 720, 25, 8, 4, 5040, 5, 40320, 8, 26, 121, 362880, 5, 12, 721, 6, 26, 3628800, 9, 39916800, 5, 122, 5041, 30, 6, 479001600, 40321, 722, 9, 6227020800, 27, 87178291200, 122, 10, 362881, 1307674368000, 6, 48, 13, 5042, 722, 20922789888000, 7, 126, 27, 40322, 3628801, 355687428096000, 10, 6402373705728000, 39916801, 28, 6, 726, 123

Offset: 1

Antti Karttunen, Aug 18 2016

nonn

Additive with a(p^e) = e * (PrimePi(p)!), where PrimePi(n) = A000720(n).

a(3181) has 1001 decimal digits. - Michael De Vlieger, Dec 24 2017

Michael De Vlieger, Table of n, a(n) for n = 1..3180 (First 120 terms from Antti Karttunen).

Cf. A000040, A000142, A000720, A002110, A007489, A019565, A028234, A033312, A055396, A059590, A067029, A076954, A206296, A260443, A276080, A276081.
Left inverse of A276076.
Cf. also A048675.

Array[If[# == 1, 0, Total[FactorInteger[#] /. {p_, e_} /; p > 1 :> e PrimePi[p]!]] &, 66] (* Michael De Vlieger, Dec 24 2017 *)

from sympy import factorint, factorial as f, primepi
def a(n):
    F=factorint(n)
    return 0 if n==1 else sum(F[i]*f(primepi(i)) for i in F)
print([a(n) for n in range(1, 121)]) # Indranil Ghosh, Jun 21 2017

a(1) = 0; for n > 1, a(n) = a(A028234(n)) + (A067029(n) * A000142(A055396(n))).
Other identities.
For all n >= 0:
a(A276076(n)) = n.
a(A002110(n)) = A007489(n).
a(A019565(n)) = A059590(n).
a(A206296(n)) = A276080(n).
a(A260443(n)) = A276081(n).
For all n >= 1:
a(A000040(n)) = n! = A000142(n).
a(A076954(n-1)) = A033312(n).

A048674 Fixed points of A048673 and A064216: Numbers n such that if n = product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) = (2*n)-1, where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

1, 2, 3, 25, 26, 33, 93, 1034, 970225, 8550146, 325422273, 414690595, 1864797542, 2438037206

Offset: 1

Antti Karttunen, Jul 14 1999

Equally: after 1, numbers n such that, if the prime factorization of 2n-1 = Product_{k >= 1} (p_k)^(c_k) then Product_{k >= 1} (p_{k-1})^(c_k) = n.
Factorization of the initial terms: 1, 2, 3, 5^2, 2*13, 3*11, 3*31, 2*11*47, 5^2*197^2, 2*11*47*8269, 3*11*797*12373, 5*11^2*433*1583, 2*23*59*101*6803, 2*11*53*1201*1741.
The only 3-cycle of permutation A048673 in range 1 .. 402653184 is (2821 3460 5639).
For 2-cycles, take setwise difference of A245449 and this sequence.
Numbers k for which A336853(k) = k-1. - Antti Karttunen, Nov 26 2021

			25 is present, as 2*25 - 1 = 49 = p_4^2, and p_3^2 = 5*5 = 25.
26 is present, as 2*26 - 1 = 51 = 3*17 = p_2 * p_8, and p_1 * p_7 = 2*13 = 26.
Alternatively, as 26 = 2*13 = p_1 * p_7, and ((p_2 * p_8)+1)/2 = ((3*17)+1)/2 = 26 also, thus 26 is present.

Fixed points of permutation pair A048673/A064216.
Positions of zeros in A349573.
Subsequence of the following sequences: A245449, A269860, A319630, A349622, A378980 (see also A379216).
This sequence is also obtained as a setwise difference of the following pairs of sequences: A246281 \ A246351, A246352 \ A246282, A246361 \ A246371, A246372 \ A246362.
Cf. A000040, A003961 (A045965), A064989, A285701, A336853.
Cf. also A348514 (fixed points of map A108228, similar to A048673).

A048673 := n -> (A003961(n)+1)/2;
A048674list := proc(upto_n) local b,i; b := [ ]; for i from 1 to upto_n do if(A048673(i) = i) then b := [ op(b), i ]; fi; od: RETURN(b); end;

Join[{1}, Reap[For[n = 1, n < 10^7, n++, ff = FactorInteger[n]; If[Times @@ Power @@@ (NextPrime[ff[[All, 1]]]^ff[[All, 2]]) == 2 n - 1, Print[n]; Sow[n]]]][[2, 1]]] (* Jean-François Alcover, Mar 04 2016 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA048674(n) = ((n+n)==(1+A003961(n))); \\ Antti Karttunen, Nov 26 2021

Entry revised and the names in Maple-code cleaned by Antti Karttunen, Aug 25 2014

Terms a(11) - a(14) added by Antti Karttunen, Sep 11-13 2014

A057335 a(0) = 1, and for n > 0, a(n) = A000040(A000120(n)) * a(floor(n/2)); essentially sequence A055932 generated using A000120, hence sorted by number of factors.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 18, 30, 16, 24, 36, 60, 54, 90, 150, 210, 32, 48, 72, 120, 108, 180, 300, 420, 162, 270, 450, 630, 750, 1050, 1470, 2310, 64, 96, 144, 240, 216, 360, 600, 840, 324, 540, 900, 1260, 1500, 2100, 2940, 4620, 486, 810, 1350, 1890, 2250, 3150, 4410

Offset: 0

Alford Arnold, Aug 27 2000

Note that for n>0 the prime divisors of a(n) are consecutive primes starting with 2. All of the least prime signatures (A025487) are included; with the other values forming A056808.
Using the formula, terms of b(n)= a(n)/A057334(n) are: 1, 1, 2, 2, 4, 4, 6, 6, 8, ..., indeed a(n) repeated. - Michel Marcus, Feb 09 2014
a(n) is the unique normal number whose unsorted prime signature is the k-th composition in standard order (graded reverse-lexicographic). This composition (row k of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. A number is normal if its prime indices cover an initial interval of positive integers. Unsorted prime signature is the sequence of exponents in a number's prime factorization. - Gus Wiseman, Apr 19 2020

			From _Gus Wiseman_, Apr 19 2020: (Start)
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}
     18: {1,2,2}
     30: {1,2,3}
     16: {1,1,1,1}
     24: {1,1,1,2}
     36: {1,1,2,2}
     60: {1,1,2,3}
     54: {1,2,2,2}
     90: {1,2,2,3}
    150: {1,2,3,3}
    210: {1,2,3,4}
     32: {1,1,1,1,1}
     48: {1,1,1,1,2}
For example, the 27th composition in standard order is (1,2,1,1), and the normal number with prime signature (1,2,1,1) is 630 = 2*3*3*5*7, so a(27) = 630.
(End)

Michael De Vlieger, Table of n, a(n) for n = 0..10000

Cf. A000120, A057334, A055932 and A056808.
Cf. A324939.
Unsorted prime signature is A124010.
Numbers whose prime signature is aperiodic are A329139.
The reversed version is A334031.
A partial inverse is A334032.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Sum is A070939.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Aperiodic compositions are A328594.
- Normal compositions are A333217.
- Permutations are A333218.
- Heinz number is A333219.
Cf. A005867, A029931, A048793, A052409, A056239, A066099, A112798, A124767, A228351, A233249, A333220, A345974.
Related to A019565 via A122111 and to A000079 via A336321.

Table[Times @@ Map[If[# == 0, 1, Prime@ #] &, Accumulate@ IntegerDigits[n, 2]], {n, 0, 54}] (* Michael De Vlieger, May 23 2017 *)

mg(n) = if (n==0, 1, prime(hammingweight(n))); \\ A057334
lista(nn) = {my(v = vector(nn)); v[1] = 1; for (i=2, nn, v[i] = mg(i-1)*v[(i+1)\2];); v;} \\ Michel Marcus, Feb 09 2014

A057335(n) = if(0==n,1,prime(hammingweight(n))*A057335(n\2)); \\ Antti Karttunen, Jul 20 2020

a(n) = A057334(n) * a (repeated).
A334032(a(n)) = n; a(A334032(n)) = A071364(n). - Gus Wiseman, Apr 19 2020
a(n) = A122111(A019565(n)); A019565(n) = A122111(a(n)). - Peter Munn, Jul 18 2020
a(n) = A336321(2^n). - Peter Munn, Mar 04 2022
Sum_{n>=0} 1/a(n) = Sum_{n>=0} 1/A005867(n) = 2.648101... (A345974). - Amiram Eldar, Jun 26 2025

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Mar 29 2003

New primary name from Antti Karttunen, Jul 20 2020

A328768 The first primorial based variant of arithmetic derivative: a(prime(i)) = A002110(i-1), where prime(i) = A000040(i), a(uv) = a(u)v + u*a(v), with a(0) = a(1) = 0.

Original entry on oeis.org

0, 0, 1, 2, 4, 6, 7, 30, 12, 12, 17, 210, 20, 2310, 67, 28, 32, 30030, 33, 510510, 44, 104, 431, 9699690, 52, 60, 4633, 54, 148, 223092870, 71, 6469693230, 80, 652, 60077, 192, 84, 200560490130, 1021039, 6956, 108, 7420738134810, 229, 304250263527210, 884, 114, 19399403, 13082761331670030, 128, 420, 145, 90124, 9292, 614889782588491410, 135, 1116, 324

Offset: 0

Antti Karttunen, Oct 28 2019

nonn

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

Cf. A002110, A108951, A276085, A276150, A328771, A373982, A374211.
Cf. A042965 (indices of even terms), A016825 (of odd terms), A152822 (antiparity of terms), A373992 (indices of multiples of 3), A374212 (2-adic valuation), A374213 (3-adic valuation), A374123 [a(n) mod 360].
Cf. A374031 [gcd(a(n), A276085(n))], A374116 [gcd(a(n), A328845(n))].
For variants of the same formula, see A003415, A258851, A328769, A328845, A328846, A371192.

A002110(n) = prod(i=1,n,prime(i));
A328768(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]*A002110(primepi(f[i,1])-1)/f[i, 1]));

A328768(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]*A002110(primepi(f[i, 1]))/(f[i, 1]^2)));

a(n) = n * Sum e_j * A276085(p_j)/p_j for n = Product p_j^e_j, where for primes p, A276085(p) = A002110(A000720(p)-1).
a(n) = n * Sum e_j * (p_j)#/(p_j^2) for n = Product p_j^e_j with (p_j)# = A034386(p_j).
For all n >= 0, A276150(a(n)) = A328771(n).

cp's OEIS Frontend

A070826 One half of product of first n primes A000040.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A243055 Difference between the indices of the smallest and the largest prime dividing n: If n = p_i * ... * p_k, where p_i <= ... <= p_k, where p_h = A000040(h), then a(n) = (k-i), a(1) = 0 by convention.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Python

Scheme

Formula

A246282 Numbers k for which A003961(k) > 2*k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2*n, where p_k indicates the k-th prime, A000040(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Scheme

Extensions

A227413 a(1)=1, a(2n)=nthprime(a(n)), a(2n+1)=nthcomposite(a(n)), where nthprime = A000040, nthcomposite = A002808.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Formula

A034602 Wolstenholme quotient W_p = (binomial(2p-1,p) - 1)/p^3 for prime p=A000040(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

Extensions

A154945 Decimal expansion of Sum_{p} 1/(p^2-1), summed over the primes p = A000040.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A246282 Numbers k for which A003961(k) > 2k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2n, where p_k indicates the k-th prime, A000040(k).

A276075 a(1) = 0, a(n) = (e1i1! + e2i2! + ... + eziz!) for n = prime(i1)^e1 prime(i2)^e2 * ... * prime(iz)^ez, where prime(k) is the k-th prime, A000040(k).

A328768 The first primorial based variant of arithmetic derivative: a(prime(i)) = A002110(i-1), where prime(i) = A000040(i), a(uv) = a(u)v + u*a(v), with a(0) = a(1) = 0.