A003972 -id:A003972 - cp's OEIS frontend

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

1, 3, 5, 9, 7, 15, 11, 27, 25, 21, 13, 45, 17, 33, 35, 81, 19, 75, 23, 63, 55, 39, 29, 135, 49, 51, 125, 99, 31, 105, 37, 243, 65, 57, 77, 225, 41, 69, 85, 189, 43, 165, 47, 117, 175, 87, 53, 405, 121, 147, 95, 153, 59, 375, 91, 297, 115, 93, 61, 315, 67, 111, 275, 729, 119

Offset: 1

Marc LeBrun

Meyers (see Guy reference) conjectures that for all r >= 1, the least odd number not in the set {a(i): i < prime(r)} is prime(r+1). - N. J. A. Sloane, Jan 08 2021
Meyers' conjecture would be refuted if and only if for some r there were such a large gap between prime(r) and prime(r+1) that there existed a composite c for which prime(r) < c < a(c) < prime(r+1), in which case (by Bertrand's postulate) c would necessarily be a term of A246281. - Antti Karttunen, Mar 29 2021
a(n) is odd for all n and for each odd m there exists a k with a(k) = m (see A064216). a(n) > n for n > 1: bijection between the odd and all numbers. - Reinhard Zumkeller, Sep 26 2001
a(n) and n have the same number of distinct primes with (A001222) and without multiplicity (A001221). - Michel Marcus, Jun 13 2014
From Antti Karttunen, Nov 01 2019: (Start)
More generally, a(n) has the same prime signature as n, A046523(a(n)) = A046523(n). Also A246277(a(n)) = A246277(n) and A287170(a(n)) = A287170(n).
Many permutations and other sequences that employ prime factorization of n to encode either polynomials, partitions (via Heinz numbers) or multisets in general can be easily defined by using this sequence as one of their constituent functions. See the last line in the Crossrefs section for examples.
(End)

			a(12) = a(2^2 * 3) = a(prime(1)^2 * prime(2)) = prime(2)^2 * prime(3) = 3^2 * 5 = 45.
a(A002110(n)) = A002110(n + 1) / 2.

Richard K. Guy, editor, Problems From Western Number Theory Conferences, Labor Day, 1983, Problem 367 (Proposed by Leroy F. Meyers, The Ohio State U.).

Indranil Ghosh, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to Heinz numbers.

See A045965 for another version.
Row 1 of table A242378 (which gives the "k-th powers" of this sequence), row 3 of A297845 and of A306697. See also arrays A066117, A246278, A255483, A308503, A329050.
Cf. A064989 (a left inverse), A064216, A000040, A002110, A000265, A027746, A046523, A048673 (= (a(n)+1)/2), A108228 (= (a(n)-1)/2), A191002 (= a(n)*n), A252748 (= a(n)-2n), A286385 (= a(n)-sigma(n)), A283980 (= a(n)*A006519(n)), A341529 (= a(n)*sigma(n)), A326042, A049084, A001221, A001222, A122111, A225546, A260443, A245606, A244319, A246269 (= A065338(a(n))), A322361 (= gcd(n, a(n))), A305293.
Cf. A191555, A252738.
Cf. A249734, A249735 (bisections).
Cf. A246261 (a(n) is of the form 4k+1), A246263 (of the form 4k+3), A246271, A246272, A246259, A246281 (n such that a(n) < 2n), A246282 (n such that a(n) > 2n), A252742.
Cf. A275717 (a(n) > a(n-1)), A275718 (a(n) < a(n-1)).
Cf. A003972 (Möbius transform), A003973 (Inverse Möbius transform), A318321.
Cf. A300841, A305421, A322991, A250469, A269379 for analogous shift-operators in other factorization and quasi-factorization systems.
Cf. also following permutations and other sequences that can be defined with the help of this sequence: A005940, A163511, A122111, A260443, A206296, A265408, A265750, A275733, A275735, A297845, A091202 & A091203, A250245 & A250246, A302023 & A302024, A302025 & A302026.
A version for partition numbers is A003964, strict A357853.
A permutation of A005408.
Applying the same transformation again gives A357852.
Other multiplicative sequences: A064988, A357977, A357978, A357980, A357983.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000720, A076610, A296150.

a003961 1 = 1
a003961 n = product $ map (a000040 . (+ 1) . a049084) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012, Oct 09 2011
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A003961 n) (apply * (map A000040 (map 1+ (map A049084 (factor n))))))
;; Antti Karttunen, May 20 2014

a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017

a[p_?PrimeQ] := a[p] = Prime[ PrimePi[p] + 1]; a[1] = 1; a[n_] := a[n] = Times @@ (a[#1]^#2& @@@ FactorInteger[n]); Table[a[n], {n, 1, 65}] (* Jean-François Alcover, Dec 01 2011, updated Sep 20 2019 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[n == 1], {n, 65}] (* Michael De Vlieger, Mar 24 2017 *)

a(n)=local(f); if(n<1,0,f=factor(n); prod(k=1,matsize(f)[1],nextprime(1+f[k,1])^f[k,2]))

a(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Michel Marcus, May 17 2014

use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016

from sympy import factorint, prime, primepi, prod
def a(n):
    f=factorint(n)
    return 1 if n==1 else prod(prime(primepi(i) + 1)**f[i] for i in f)
[a(n) for n in range(1, 11)] # Indranil Ghosh, May 13 2017

If n = Product p(k)^e(k) then a(n) = Product p(k+1)^e(k).
Multiplicative with a(p^e) = A000040(A000720(p)+1)^e. - David W. Wilson, Aug 01 2001
a(n) = Product_{k=1..A001221(n)} A000040(A049084(A027748(n,k))+1)^A124010(n,k). - Reinhard Zumkeller, Oct 09 2011 [Corrected by Peter Munn, Nov 11 2019]
A064989(a(n)) = n and a(A064989(n)) = A000265(n). - Antti Karttunen, May 20 2014 & Nov 01 2019
A001221(a(n)) = A001221(n) and A001222(a(n)) = A001222(n). - Michel Marcus, Jun 13 2014
From Peter Munn, Oct 31 2019: (Start)
a(n) = A225546((A225546(n))^2).
a(A225546(n)) = A225546(n^2).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p))) = 2.06399637... . - Amiram Eldar, Nov 18 2022

A003973 Inverse Möbius transform of A003961; a(n) = sigma(A003961(n)), where A003961 shifts the prime factorization of n one step towards the larger primes.

Original entry on oeis.org

1, 4, 6, 13, 8, 24, 12, 40, 31, 32, 14, 78, 18, 48, 48, 121, 20, 124, 24, 104, 72, 56, 30, 240, 57, 72, 156, 156, 32, 192, 38, 364, 84, 80, 96, 403, 42, 96, 108, 320, 44, 288, 48, 182, 248, 120, 54, 726, 133, 228, 120, 234, 60, 624, 112, 480, 144, 128, 62, 624, 68

Offset: 1

Marc LeBrun

Sum of the divisors of the prime shifted n, or equally, sum of the prime shifted divisors of n. - Antti Karttunen, Aug 17 2020

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to sigma(n).

Cf. A000203, A000290 (positions of odd terms), A003961, A007814, A048673, A108228, A151800, A295664, A336840.
Permutation of A008438.
Used in the definitions of the following sequences: A326042, A336838, A336841, A336844, A336846, A336847, A336848, A336849, A336850, A336851, A336852, A336856, A336931, A336932.
Cf. also A003972.

b[1] = 1; b[p_?PrimeQ] := b[p] = Prime[ PrimePi[p] + 1]; b[n_] := b[n] = Times @@ (b[First[#]]^Last[#] &) /@ FactorInteger[n]; a[n_] := Sum[ b[d], {d, Divisors[n]}]; Table[a[n], {n, 1, 70}]  (* Jean-François Alcover, Jul 18 2013 *)

aPrime(p,e)=my(q=nextprime(p+1));(q^(e+1)-1)/(q-1)
a(n)=my(f=factor(n));prod(i=1,#f~,aPrime(f[i,1],f[i,2])) \\ Charles R Greathouse IV, Jul 18 2013

A003973(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); sigma(factorback(f)); }; \\ Antti Karttunen, Aug 06 2020

from math import prod
from sympy import factorint, nextprime
def A003973(n): return prod(((q:=nextprime(p))**(e+1)-1)//(q-1) for p,e in factorint(n).items()) # Chai Wah Wu, Jul 05 2022

Multiplicative with a(p^e) = (q^(e+1)-1)/(q-1) where q = nextPrime(p). - David W. Wilson, Sep 01 2001
From Antti Karttunen, Aug 06-12 2020: (Start)
a(n) = Sum_{d|n} A003961(d) = Sum_{d|A003961(n)} d.
a(n) = A000203(A003961(n)) = A000593(A003961(n)).
a(n) = 2*A336840(n) - A000005(n) = 2*Sum_{d|n} (A048673(d) - (1/2)).
a(n) = A008438(A108228(n)) = A008438(A048673(n)-1).
a(n) = A336838(n) * A336856(n).
a(n) is odd if and only if n is a square.
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} p^3/((p+1)*(p^2-nextprime(p))) = 3.39513795..., where nextprime is A151800. - Amiram Eldar, Dec 08 2022, May 30 2025

More terms from David W. Wilson, Aug 29 2001

Secondary name added by Antti Karttunen, Aug 06 2020

A339821 a(n) = phi(A019565(2n)), where phi is Euler totient function.

Original entry on oeis.org

1, 2, 4, 8, 6, 12, 24, 48, 10, 20, 40, 80, 60, 120, 240, 480, 12, 24, 48, 96, 72, 144, 288, 576, 120, 240, 480, 960, 720, 1440, 2880, 5760, 16, 32, 64, 128, 96, 192, 384, 768, 160, 320, 640, 1280, 960, 1920, 3840, 7680, 192, 384, 768, 1536, 1152, 2304, 4608, 9216, 1920, 3840, 7680, 15360, 11520, 23040, 46080, 92160

Offset: 0

Antti Karttunen, Dec 18 2020

nonn

Bisection of A339820.
Cf. A000010, A003961, A003972, A006093, A019565, A339822 (2-adic valuation).
Cf. also A324651.

A019565(n) = { my(m=1, p=1); while(n>0, p = nextprime(1+p); if(n%2, m *= p); n >>= 1); (m); };
A339821(n) = eulerphi(A019565(n+n));

A339821(n) = { my(m=1, p=2); while(n>0, p = nextprime(1+p); if(n%2, m *= (p-1)); n >>= 1); (m); };

If 4n = 2^e1 + 2^e2 + ... + 2^ek [e1 ... ek distinct], then a(n) = A006093(e1) * A006093(e2) * ... * A006093(ek).
a(n) = A339820(2n) = A000010(A019565(2n)) = A000010(A019565(2n+1)).
a(n) = A003972(A019565(n)) = A000010(A003961(A019565(n))).

A349387 Dirichlet convolution of A003961 with A055615 (Dirichlet inverse of n), where A003961 is fully multiplicative with a(p) = nextprime(p).

Original entry on oeis.org

1, 1, 2, 3, 2, 2, 4, 9, 10, 2, 2, 6, 4, 4, 4, 27, 2, 10, 4, 6, 8, 2, 6, 18, 14, 4, 50, 12, 2, 4, 6, 81, 4, 2, 8, 30, 4, 4, 8, 18, 2, 8, 4, 6, 20, 6, 6, 54, 44, 14, 4, 12, 6, 50, 4, 36, 8, 2, 2, 12, 6, 6, 40, 243, 8, 4, 4, 6, 12, 8, 2, 90, 6, 4, 28, 12, 8, 8, 4, 54, 250, 2, 6, 24, 4, 4, 4, 18, 8, 20, 16, 18, 12, 6

Offset: 1

Antti Karttunen, Nov 17 2021

Multiplicative because A003961 and A055615 are.
Convolving this with A000010 gives A003972, and convolving this with A000203 gives A003973.
Multiplicative with a(p^e) = nextprime(p)^e - p * nextprime(p)^(e-1), where nextprime function is A151800. - Amiram Eldar, Nov 18 2021

Cf. A000040, A001223, A003961, A055615, A151800, A349388 (Dirichlet inverse), A349389 (sum with it), A378606 (Möbius transform).

Cf. also A000010, A000203, A003972, A003973, A347236, A349381.

f[p_,e_] := (q = NextPrime[p])^e - p * q^(e-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 18 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A055615(n) = (n*moebius(n));
A349387(n) = sumdiv(n,d,A003961(n/d)*A055615(d));

a(n) = Sum_{d|n} A003961(n/d) * A055615(d).

For all n >= 1, a(A000040(n)) = A001223(n).

A339822 The exponent of the highest power of 2 dividing A339821(n).

Original entry on oeis.org

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

Offset: 0

Antti Karttunen, Dec 18 2020

nonn

Antti Karttunen, Table of n, a(n) for n = 0..65537
Index entries for sequences related to binary expansion of n

Cf. A000010, A003972, A007814, A019565, A023506, A053574, A339814, A339815, A339821.

A339822(n) = { my(s=0, p=2); while(n>0, p = nextprime(1+p); if(n%2, s += valuation((p-1),2)); n >>= 1); (s); };

If 4n = 2^e1 + 2^e2 + ... + 2^ek [e1 ... ek distinct], then a(n) = A023506(e1) + A023506(e2) + ... + A023506(ek).
a(n) = A007814(A339821(n)) = A053574(A019565(2n)).
a(n) = A007814(A003972(A019565(n))) = A007814(A000010(A019565(2n))).

A344587 Deficiency of prime-shifted n: a(n) = 2*A003961(n) - sigma(A003961(n)).

Original entry on oeis.org

1, 2, 4, 5, 6, 6, 10, 14, 19, 10, 12, 12, 16, 18, 22, 41, 18, 26, 22, 22, 38, 22, 28, 30, 41, 30, 94, 42, 30, 18, 36, 122, 46, 34, 58, 47, 40, 42, 62, 58, 42, 42, 46, 52, 102, 54, 52, 84, 109, 66, 70, 72, 58, 126, 70, 114, 86, 58, 60, 6, 66, 70, 178, 365, 94, 54, 70, 82, 110, 78, 72, 110, 78, 78, 148, 102, 118, 78

Offset: 1

Antti Karttunen, May 28 2021

First negative value occurs as a(120) = -30.

Questions: Which subsets of natural numbers generate the "cut sigmoid" graph(s) that cross the X-axis in the (lowermost) scatter plot?

Cf. A000203, A003961, A003973, A033879, A153881, A336851, A337386 (positions of terms <= 0), A346246 (Dirichlet inverse), A349387, A378216, A378231 [= a(n^2)].

Inverse Möbius transform of A337544.

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A344587(n) = { my(u=A003961(n)); (u+u - sigma(u)); };

a(n) = A033879(A003961(n)) = 2*A003961(n) - A003973(n).
a(n) = Sum_{d|n} A337544(d).
From Antti Karttunen, Nov 23 2024: (Start)
a(n) = Sum_{d|n} A003961(d)*A153881(n/d) = A003961(n) - A336851(n).
a(n) = Sum_{d|n} A033879(d)*A349387(n/d).
a(n) = Sum_{d|n} A003972(d)*A378216(n/d).
(End)

A349127 Möbius transform of A064989, where A064989 is multiplicative with a(2^e) = 1 and a(p^e) = prevprime(p)^e for odd primes p.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 4, 0, 2, 0, 6, 0, 10, 0, 2, 0, 12, 0, 16, 0, 4, 0, 18, 0, 6, 0, 4, 0, 22, 0, 28, 0, 6, 0, 8, 0, 30, 0, 10, 0, 36, 0, 40, 0, 4, 0, 42, 0, 20, 0, 12, 0, 46, 0, 12, 0, 16, 0, 52, 0, 58, 0, 8, 0, 20, 0, 60, 0, 18, 0, 66, 0, 70, 0, 6, 0, 24, 0, 72, 0, 8, 0, 78, 0, 24, 0, 22, 0, 82, 0, 40, 0, 28, 0, 32

Offset: 1

Antti Karttunen, Nov 13 2021

The multiplicative definition of this sequence ("Möbius transform of prime shift towards lesser primes") differs from otherwise similarly defined A349128 (Euler phi applied to A064989) only in that here a(2^e) = 0, while A349128(2^e) = 1.

Compare the situation with A003961 ("prime shift towards larger primes"), where A003972(n) = A000010(A003961(n)) is also the Möbius transform of A003961.

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

Agrees with A347115, A348045 and A349128 on odd numbers.
Cf. A000004, A285702 (even and odd bisection).
Cf. A000010, A008683, A064989, A151799.
Cf. also A003961, A003972, A349125, A349136.

f[p_, e_] := ((q = NextPrime[p, -1]) - 1)*q^(e - 1); a[1] = 1; a[n_] := If[EvenQ[n], 0, Times @@ f @@@ FactorInteger[n]]; Array[a, 100] (* Amiram Eldar, Nov 13 2021 *)

A349127(n) = if(!(n%2),0, my(f = factor(n), q); prod(i=1, #f~, q = precprime(f[i,1]-1); (q-1)*(q^(f[i,2]-1))));

A064989(n) = { my(f = factor(n)); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); };
A349127(n) = if(n%2, eulerphi(A064989(n)), 0);

A349127(n) = sumdiv(n,d,moebius(n/d)*A064989(d));

Multiplicative with a(2^e) = 0, and for odd primes p, a(p^e) = (q-1)*q^(e-1), where q = prevprime(p), where prevprime is A151799.
If n is odd, then a(n) = A000010(A064989(n)), and if n is even, then a(n) = 0.
a(n) = Sum_{d|n} A008683(d) * A064989(n/d).
For all n >= 1, a(2n-1) = A347115(2n-1) = A348045(2n-1) = A349128(2n-1) = A285702(n).
Sum_{k=1..n} a(k) ~ c * n^2, where c = (16/Pi^4) / Product_{p prime > 2} (1+1/p-q(p)/p^2-q(p)/p^3) = 0.1341718..., where q(p) = prevprime(p) = A151799(p). - Amiram Eldar, Dec 24 2022

A339905 Fully multiplicative with a(prime(k)) = prime(k+1) - 1.

Original entry on oeis.org

1, 2, 4, 4, 6, 8, 10, 8, 16, 12, 12, 16, 16, 20, 24, 16, 18, 32, 22, 24, 40, 24, 28, 32, 36, 32, 64, 40, 30, 48, 36, 32, 48, 36, 60, 64, 40, 44, 64, 48, 42, 80, 46, 48, 96, 56, 52, 64, 100, 72, 72, 64, 58, 128, 72, 80, 88, 60, 60, 96, 66, 72, 160, 64, 96, 96, 70, 72, 112, 120, 72, 128, 78, 80, 144, 88, 120, 128, 82

Offset: 1

Antti Karttunen, Dec 30 2020

Antti Karttunen, Table of n, a(n) for n = 1..10205
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A003958, A003961, A151800, A339903 (the odd part).

Coincides with A003972 on A005117.

A339905(n) = if(1==n,n,my(f=factor(n)); for(i=1,#f~,f[i,1] = nextprime(1+f[i,1])-1); factorback(f));

a(n) = A003958(A003961(n)).

Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p)+1)) = 0.73732173..., where nextprime is A151800. - Amiram Eldar, Dec 08 2022

A340071 a(n) = gcd(A003961(n)-1, phi(A003961(n))), where A003961 shifts the prime factorization of n one step towards larger primes.

Original entry on oeis.org

1, 2, 4, 2, 6, 2, 10, 2, 4, 4, 12, 4, 16, 4, 2, 2, 18, 2, 22, 2, 2, 2, 28, 2, 6, 2, 4, 2, 30, 8, 36, 2, 16, 4, 4, 8, 40, 4, 4, 4, 42, 4, 46, 4, 6, 2, 52, 4, 10, 2, 2, 8, 58, 2, 18, 4, 2, 4, 60, 2, 66, 2, 2, 2, 2, 2, 70, 2, 16, 10, 72, 2, 78, 2, 4, 2, 2, 2, 82, 2, 4, 4, 88, 2, 12, 4, 2, 2, 96, 4, 2, 4, 8, 2, 4, 2, 100

Offset: 1

Antti Karttunen, Dec 28 2020

nonn

Prime shifted analog of A049559.

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

Cf. A000010, A003961, A003972, A049559, A253885, A340072, A340073.

Cf. also A340081.

A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A340071(n) = { my(u=A003961(n)); gcd(u-1, eulerphi(u)); };

a(n) = A049559(A003961(n)).
a(n) = gcd(A253885(n-1), A003972(n)) = gcd(A003961(n)-1, A000010(A003961(n))).
a(n) = A003972(n) / A340072(n).
For n > 1, a(n) = (A003961(n)-1) / A340073(n).

cp's OEIS Frontend

A337549 a(n) = A003972(n) - n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Perl

Python

Formula

A003973 Inverse Möbius transform of A003961; a(n) = sigma(A003961(n)), where A003961 shifts the prime factorization of n one step towards the larger primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

Extensions

A339821 a(n) = phi(A019565(2n)), where phi is Euler totient function.

Views

Author

Keywords

Crossrefs

Programs

PARI

PARI

Formula

A349387 Dirichlet convolution of A003961 with A055615 (Dirichlet inverse of n), where A003961 is fully multiplicative with a(p) = nextprime(p).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A339822 The exponent of the highest power of 2 dividing A339821(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A344587 Deficiency of prime-shifted n: a(n) = 2*A003961(n) - sigma(A003961(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs