A179119 -id:A179119 - cp's OEIS frontend

A136141 Decimal expansion of Sum_{p prime} 1/(p*(p-1)).

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

Offset: 0

R. J. Mathar, Mar 09 2008

Excess of prime factors with multiplicity over distinct prime factors for random (large) integers. - Charles R Greathouse IV, Sep 06 2011
Sum of reciprocals of (proper) prime powers. The sum of reciprocals of all proper powers is A072102. - Charles R Greathouse IV, Apr 24 2012
Decimal expansion of the infinite sum of the reciprocals of the prime powers which are not prime (A246547). - Robert G. Wilson v, May 13 2019
See the second 'Applications' example under the Mathematica help file for the function PrimePowerQ. - Robert G. Wilson v, May 13 2019
It easy to prove that this constant < 1 (Sum_{p prime} 1/(p*(p-1)) < Sum_{k>=2} 1/(k*(k-1)) = 1). Luthar (1969) asks for a better upper bound. The solution shows that this constant is < 3/2 - log(2) = 0.80685... . - Amiram Eldar, Feb 14 2025

			Equals 1/2 + 1/(3*2) + 1/(5*4) + 1/(7*6) + ...
= 0.7731566690497951278643674598559423956187413360831860483110060673567...

Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
Steven R. Finch, Mathematical Constants, Cambridge Univ. Press, 2003, Meissel-Mertens constants, p. 94.

Robert G. Wilson v, Table of n, a(n) for n = 0..10000 (first 1002 terms from Jason Kimberley).
Henri Cohen, High-precision computation of Hardy-Littlewood constants, preprint 1991.
Henri Cohen, High-precision computation of Hardy-Littlewood constants. [pdf copy, with permission]
Stephan Ramon Garcia and Ethan Simpson Lee, Explicit conditional bounds for the residue of a Dedekind zeta-function at s = 1, arXiv:2506.17416 [math.NT], 2025. See p. 3.
R. S. Luthar, Problem E 2192, Elementary Problems, The American Mathematical Monthly, Vol. 76, No. 8 (1969), p. 938; An Upper Bound, Solution to Problem E 2192, by O. P. Lossers, ibid., Vol. 77, No. 7 (1970), pp. 769-770.
R. J. Mathar, Series of reciprocal powers of k-almost primes, arXiv:0803.0900 [math.NT], 2008-2009. Tables 8 and 10.
Michael Ian Shamos, Property Enumerators and a Partial Sum Theorem, 2011; alternative link.
Index to constants which are prime zeta sums {1,1,0}

Cf. A152447 (over the semiprimes), A000040, A000720, A001248, A046660 (excess, see first comment), A072102, A077761, A083342, A179119, A246547.

Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), A085964 (at 4) to A085969 (at 9).

R := RealField(105);
c := &+[R|(EulerPhi(n)-MoebiusMu(n))/n*Log(ZetaFunction(R,n)):n in[2..360]];
Reverse(IntegerToSequence(Floor(c*10^103))); // Jason Kimberley, Jan 12 2017

digits = 103; sp = NSum[PrimeZetaP[n], {n, 2, Infinity}, WorkingPrecision -> digits + 10, NSumTerms -> 2*digits]; RealDigits[sp, 10, digits] // First (* Jean-François Alcover, Sep 02 2015 *)

W(x)=solve(y=log(x)/2,max(1,log(x)),y*exp(y)-x)
eps()=2. >> (32*ceil(default(realprecision)/9.63))
primezeta(s)=my(t=s*log(2),iter=W(t/eps())\t);sum(k=1,iter, moebius(k)/k*log(abs(zeta(k*s))))
a(lim,e)={ \\ choose parameters to maximize speed and precision
    my(x,y=exp(W(lim)-.5));
    x=lim^e*(e*log(y))^e*(y*log(y))^-e*incgam(-e,e*log(y));
    forprime(p=2,lim,x+=1/((p*1.)^e*(p-1)));
    x+sum(n=2,e,primezeta(n))
}; \\ Charles R Greathouse IV, Sep 07 2011

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

Equals Sum_{n>=1} 1/A036689(n).
Equals Sum_{s>=2} P(s), where P is the prime zeta function. - Charles R Greathouse IV, Sep 06 2011
Equals A083342 - A077761, that is, Sum_{n>=2} ((EulerPhi(n) - MoebiusMu(n))/n) * log(zeta(n)). - Jean-François Alcover, Sep 02 2015
Equals 2 * Sum_{k>=2} pi(k)/(k^3-k), where pi(k) = A000720(k) (Shamos, 2011, p. 8). - Amiram Eldar, Mar 12 2024

More terms from D. S. McNeil, Sep 06 2011

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

A060687 Numbers k such that there exist exactly 2 Abelian groups of order k, i.e., A000688(k) = 2.

Original entry on oeis.org

4, 9, 12, 18, 20, 25, 28, 44, 45, 49, 50, 52, 60, 63, 68, 75, 76, 84, 90, 92, 98, 99, 116, 117, 121, 124, 126, 132, 140, 147, 148, 150, 153, 156, 164, 169, 171, 172, 175, 188, 198, 204, 207, 212, 220, 228, 234, 236, 242, 244, 245, 260, 261, 268, 275, 276, 279

Offset: 1

Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 19 2001

k belongs to this sequence iff exactly one prime in its factorization into prime powers has exponent 2 and all the other primes in the factorization have exponent 1, for example 60 = 2^2 * 3 * 5.
Numbers k such that A046660(k) = 1. - Zak Seidov, Nov 14 2012
Numbers that have twice as many unitary divisors as nonunitary divisors, the largest possible ratio for nonsquarefree numbers (i.e., numbers that have nonunitary divisors). - Amiram Eldar, Nov 01 2024

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5000 from Enrique Pérez Herrero)
Eckford Cohen, Arithmetical notes. VIII. An asymptotic formula of Rényi, Proc. Amer. Math. Soc. 13 (1962), pp. 536-539.
Alfred Rényi, On the density of certain sequences of integers, Publications de l'Institut Mathématique, Vol. 8 (1955), pp. 157-162.
Index entries for sequences related to groups.

Cf. A001221, A001222, A000688, A046660, A271971, A013661, A179119.

a060687 n = a060687_list !! (n-1)
a060687_list = filter ((== 1) . a046660) [1..]
-- Reinhard Zumkeller, Nov 29 2015

Select[Range[500], PrimeOmega[#] - PrimeNu[#] == 1 &] (* Harvey P. Dale, Sep 08 2011 *)

for(n=1,279,if(bigomega(n)-omega(n)==1,print1(n,",")))

is(n)=factorback(factor(n)[,2])==2 \\ Charles R Greathouse IV, Sep 18 2015

list(lim)=my(s=lim\4,v=List(),u=vectorsmall(s,i,1),t,x); forprime(k=2,sqrtint(s), t=k^2; forstep(i=t,s,t, u[i]=0)); forprime(k=2,sqrtint(lim\1), t=k^2; for(i=1,#u, if(u[i] && gcd(k,i)==1, x=t*i; if(x>lim, break); listput(v,x)))); Set(v) \\ Charles R Greathouse IV, Aug 02 2016

from math import isqrt
from sympy import mobius, primerange
def A060687(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def g(x): return sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
    def f(x): return int(n+x+sum(sum(g(x//p**j) if j&1 else -g(x//p**j) for j in range(2,x.bit_length())) for p in primerange(isqrt(x)+1)))
    return bisection(f,n,n) # Chai Wah Wu, Feb 24 2025

from sympy import factorint
def is_A060687(n): return sum(v := factorint(n).values()) == len(v) + 1 # David Radcliffe, Jul 28 2025

k such that A001222(k)-A001221(k) = 1.

Cohen proved that a(n) = kn + O(sqrt(n) log log n), where k = A013661/A179119 = 1/A271971 = 4.981178... - Charles R Greathouse IV, Aug 02 2016

Corrected and extended by Vladeta Jovovic, Jul 05 2001

A034953 Triangular numbers (A000217) with prime indices.

Original entry on oeis.org

3, 6, 15, 28, 66, 91, 153, 190, 276, 435, 496, 703, 861, 946, 1128, 1431, 1770, 1891, 2278, 2556, 2701, 3160, 3486, 4005, 4753, 5151, 5356, 5778, 5995, 6441, 8128, 8646, 9453, 9730, 11175, 11476, 12403, 13366, 14028, 15051, 16110, 16471, 18336, 18721, 19503

Offset: 1

Patrick De Geest, Oct 15 1998

The following sequences (allowing offset of first term) all appear to have the same parity: A034953, triangular numbers with prime indices; A054269, length of period of continued fraction for sqrt(p), p prime; A082749, difference between the sum of next prime(n) natural numbers and the sum of next n primes; A006254, numbers n such that 2n-1 is prime; A067076, 2n+3 is a prime. - Jeremy Gardiner, Sep 10 2004
Given a rectangular prism with sides 1, p, p^2 for p = n-th prime (n > 1), the area of the six sides divided by the volume gives a remainder which is 4*a(n). - J. M. Bergot, Sep 12 2011
The infinite sum over the reciprocals is given by 2*A179119. - Wolfdieter Lang, Jul 10 2019

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Triangular Number.

Cf. A000217, A034954, A034955, A011756, A179119, A195678, A307868.

Cf. A054269, A067076, A082749.

a034953 n = a034953_list !! (n-1)
a034953_list = map a000217 a000040_list
-- Reinhard Zumkeller, Sep 23 2011

a:= n-> (p-> p*(p+1)/2)(ithprime(n)):
seq(a(n), n=1..65);  # Alois P. Heinz, Apr 20 2022

t[n_] := n(n + 1)/2; Table[t[Prime[n]], {n, 44}] (* Robert G. Wilson v, Aug 12 2004 *)
(#(# + 1))/2&/@Prime[Range[50]] (* Harvey P. Dale, Feb 27 2012 *)
With[{nn=200},Pick[Accumulate[Range[nn]],Table[If[PrimeQ[n],1,0],{n,nn}],1]] (* Harvey P. Dale, Mar 05 2023 *)

forprime(p=2,1e3,print1(binomial(p+1,2)", ")) \\ Charles R Greathouse IV, Jul 19 2011

apply(n->binomial(n+1,2),primes(100)) \\ Charles R Greathouse IV, Jun 04 2013

a(n) = A000217(A000040(n)). - Omar E. Pol, Jul 27 2009
a(n) = Sum_{k=1..prime(n)} k. - Wesley Ivan Hurt, Apr 27 2021
Product_{n>=1} (1 - 1/a(n)) = A307868. - Amiram Eldar, Nov 07 2022

A162642 Number of odd exponents in the canonical prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Reinhard Zumkeller, Jul 08 2009

a(n) is also known as the squarefree rank of n. - Jason Kimberley, Jul 08 2017

The number of primes that are infinitary divisors of n. - Amiram Eldar, Oct 01 2023

Jason Kimberley, Table of n, a(n) for n = 1..20000
R. B. Eggleton, J. S. Kimberley and J. A. MacDougall, Square-free rank of integers, Journal of Combinatorial Mathematics and Combinatorial Computing, Vol. 106 (2018).
Index entries for sequences computed from exponents in factorization of n.

Cf. A000290 (positions of zeros), A001221, A001620, A002035, A007913, A056169, A162641, A295316, A295659, A295662, A295664.

A162642:=func;
[A162642(n):n in[1..105]]; // Jason Kimberley, Dec 30 2015

A162642 := proc(n) add ( op(2,f) mod 2 ,f=ifactors(n)[2]) ; end proc: # R. J. Mathar, Mar 30 2011

{0}~Join~Table[Count[Last /@ FactorInteger@ n, e_ /; OddQ@ e], {n, 2, 105}] (* Michael De Vlieger, Jan 06 2016 *)

a(n) = {my(f = factor(n)); sum(k=1, #f~, f[k,2] % 2);} \\ Michel Marcus, Jan 08 2016

;; With memoization-macro definec.
(definec (A162642 n) (if (= 1 n) 0 (+ (A000035 (A067029 n)) (A162642 (A028234 n))))) ;; Antti Karttunen, Nov 28 2017

a(n) = A001221(n) - A162641(n).
a(n) = A001221(A007913(n)). - Jason Kimberley, Jan 06 2016
a(A000290(n)) = 0, n > 0. - Michel Marcus, Jan 08 2016
G.f.: Sum_{i>=1} Sum_{j>=1} (-1)^j x^(prime(i)^j)/(x^(prime(i)^j) - 1). - Robert Israel, Jan 15 2016
From Antti Karttunen, Nov 28 2017: (Start)
Additive with a(p^e) = A000035(e).
a(n) = A056169(n) + A295662(n).
A056169(n) <= a(n) <= A056169(n) + A295659(n).
a(n) <= A295664(n).
(End)
Sum_{k=1..n} a(k) = n * log(log(n)) + c * n + O(n/log(n)), where c = gamma + Sum_{p prime} (log(1-1/p) + 1/(p+1)) = A077761 - A179119 = -0.0687327134... and gamma is Euler's constant (A001620). - Amiram Eldar, Dec 25 2021

A023998 Number of block permutations on an n-set which are uniform, i.e., corresponding blocks have same size.

Original entry on oeis.org

1, 1, 3, 16, 131, 1496, 22482, 426833, 9934563, 277006192, 9085194458, 345322038293, 15024619744202, 740552967629021, 40984758230303149, 2527342803112928081, 172490568947825135203, 12952575262915522547136, 1064521056888312620947794, 95305764621957309071404877

Offset: 0

Des FitzGerald (D.FitzGerald(AT)utas.edu.au)

Number of games of simple patience with n cards. Take a shuffled deck of n cards labeled 1..n; as each card is dealt it is placed either on a higher-numbered card or starts a new pile to the right. Cards are not moved once they are placed. Suggested by reading Aldous and Diaconis. - N. J. A. Sloane, Dec 19 1999

Number of set partitions of [2n] such that within each block the numbers of odd and even elements are equal. a(2) = 3: 1234, 12|34, 14|23; a(3) = 16: 123456, 1234|56, 1236|45, 1245|36, 1256|34, 12|3456, 12|34|56, 12|36|45, 1346|25, 1456|23, 14|2356, 14|23|56, 16|2345, 16|23|45, 14|25|36, 16|25|34. - Alois P. Heinz, Jul 14 2016

			For n=3 there are 25 block permutations, of which 9 of the form ({1} maps to {1,2}; {2,3} maps to {3}), are not uniform. Hence a(3) = 25 - 9 = 16.
Alternatively, for n=3 the 6 permutations of 3 cards produce 16 games, as follows: 123 -> {1,2,3}; 132 -> {1,32}, {1,3,2}; 213 -> {21,3}, {2,1,3}; 231 -> {21,3}, {2,31}, {2,3,1}; 312 -> {31,2}, {32,1}, {3,1,2}; 321 -> {321}, {32,1}, {31,2}, {3,21}, {3,2,1}.
G.f.: A(x) = 1 + x + 3*x^2/2!^2 + 16*x^3/3!^2 + 131*x^4/4!^2 + 1496*x^5/5!^2 + ...
log(A(x)) = x + x^2/2!^2 + x^3/3!^2 + x^4/4!^2 + x^5/5!^2 + ...

Reinhard Zumkeller, Table of n, a(n) for n = 0..300
M. Aguiar and R. C. Orellana, The Hopf algebra of uniform block permutations, J. Algebr. Comb. 28 (2008) 115-138
D. Aldous and P. Diaconis, Longest increasing subsequences: from patience sorting to the Baik-Deift-Johansson theorem, Bull. Amer. Math. Soc. 36 (1999), 413-432.
H. Cheballah, S. Giraudo, and R. Maurice, Combinatorial Hopf algebra structure on packed square matrices, arXiv preprint arXiv:1306.6605 [math.CO], 2013.
Fabian Faulstich, Bernd Sturmfels, and Svala Sverrisdóttir, Algebraic Varieties in Quantum Chemistry, arXiv:2308.05258 [math.AG], 2023.
D. G. FitzGerald and Jonathan Leech, Dual symmetric inverse monoids and representation theory, J. Australian Mathematical Society (Series A), Vol. 64 (1998), pp. 345-367.
Raúl E. González-Torres, A geometric study of cores of idempotent stochastic matrices, Linear Algebra Appl. 527, 87-127 (2017).
Rosa Orellana, Franco Saliola, Anne Schilling, and Mike Zabrocki, The lattice of submonoids of the uniform block permutations containing the symmetric group, arXiv:2405.09710 [math.CO], 2024. See p. 3.
Sebastian Volz, Design and Implementation of Efficient Algorithms for Operations on Partitions of Sets, Bachelor Thesis, Saarland Univ. (Germany, 2023). See p. 45.

Cf. A023997, A002720, A061691.
Cf. A132813.
Column k=2 of A275043.
Main diagonal of A321296 and of A322670.

a023998 n = a023998_list !! n
a023998_list = 1 : f 2 [1] a132813_tabl where
   f x ys (zs:zss) = y : f (x + 1) (ys ++ [y]) zss where
                     y = sum $ zipWith (*) ys zs
-- Reinhard Zumkeller, Apr 04 2014

b:= proc(n) option remember; `if`(n=0, 1,
      add(b(n-i)*binomial(n-1, i-1)/i!, i=1..n))
    end:
a:= n-> b(n)*n!:
seq(a(n), n=0..25);  # Alois P. Heinz, May 11 2016

a[0] = 1; a[n_] := a[n] = Sum[Binomial[n, k] Binomial[n-1, k] a[k], {k, 0, n-1}];
Array[a, 25, 0] (* Jean-François Alcover, Jul 28 2016 *)
nmax = 20; CoefficientList[Series[E^(-1 + BesselI[0, 2*Sqrt[x]]), {x, 0, nmax}], x]*Range[0, nmax]!^2 (* Vaclav Kotesovec, Jun 09 2019 *)

a(n)=if(n==0,1,sum(k=0,n-1,binomial(n,k)*binomial(n-1,k)*a(k))) \\ Paul D. Hanna, Aug 15 2007

{a(n)=n!^2*polcoeff(exp(sum(m=1, n, x^m/m!^2)+x*O(x^n)), n)} /* Paul D. Hanna */

N=66; x='x+O('x^N); /* that many terms */
Vec(serlaplace(serlaplace(exp(sum(n=1, N, x^n/n!^2))))) /* show terms */
/* Joerg Arndt, Jul 12 2011 */

v=vector(N); v[1]=1;
for (n=1,N-1, v[n+1]=sum(k=0,n-1, binomial(n,k)*binomial(n-1,k)*v[k+1]) );
v /* show terms */
/* Joerg Arndt, Jul 12 2011 */

a(n) = Sum_{k=0..n-1} C(n,k)*C(n-1,k)*a(k) for n>0 with a(0)=1. - Paul D. Hanna, Aug 15 2007
G.f.: Sum_{n>=0} a(n)*x^n/n!^2 = exp( Sum_{n>=1} x^n/n!^2 ). [Paul D. Hanna, Jan 04 2011; merged from duplicate entry A179119]
Row sums of A061691.
Generating function: Let J(z) = Sum_{n>=0} z^n/n!^2. Then exp(J(z)-1) = Sum_{n>=0} a(n)*z^n/n!^2 = 1 + z + 3*z^2/2!^2 + 16*z^3/3!^2 + .... - Peter Bala, Jul 11 2011

More terms from Vladeta Jovovic, Sep 03 2002

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

A316523 Number of odd multiplicities minus number of even multiplicities in the canonical prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 05 2018

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

Cf. A000040, A000607, A071321, A100118, A112798.
Cf. A187039 (where a(n)=0). - Michel Marcus, Jul 08 2018
Cf. A077761, A162641, A162642, A179119.

f:= proc(n) local F;
  F:= map(t -> t[2],ifactors(n)[2]);
  2*nops(select(type,F,odd))-nops(F);
end proc:
map(f, [$1..100]); # Robert Israel, Aug 27 2018

Table[Total[-(-1)^If[n==1,{},FactorInteger[n][[All,2]]]],{n,100}]

a(n) = my(f=factor(n)); -sum(k=1, #f~, (-1)^(f[k,2])); \\ Michel Marcus, Jul 08 2018; corrected Jun 13 2022

If i and j are coprime, a(i*j) = a(i)+a(j). - Robert Israel, Aug 27 2018
From Amiram Eldar, Oct 05 2023: (Start)
Additive with a(p^e) = (-1)^(e+1).
a(n) = A162642(n) - A162641(n).
Sum_{k=1..n} a(k) = n * log(log(n)) + c * n + O(n/log(n)), where c = A077761 - 2*A179119 = -0.398962... . (End)

A162641 Number of even exponents in canonical prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Reinhard Zumkeller, Jul 08 2009

nonn

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

Cf. A001221, A003557, A002035, A007913, A025487, A059841, A061742 (where records occur), A072587, A162642, A179119.

Cf. A268335 (positions of zeros), A295316.

Table[Count[FactorInteger[n][[All, -1]], ?EvenQ], {n, 105}] (* _Michael De Vlieger, Jul 23 2017 *)

A162641(n) = omega(n) - omega(core(n)); \\ Antti Karttunen, Jul 23 2017

(define (A162641 n) (if (= 1 n) 0 (+ (A059841 (A067029 n)) (A162641 (A028234 n))))) ;; Antti Karttunen, Jul 23 2017

a(n) = A001221(n) - A162642(n).
a(A002035(n)) = 0.
a(A072587(n)) > 0.
Additive with a(p^e) = A059841(e). - Antti Karttunen, Jul 23 2017
From Antti Karttunen, Nov 28 2017: (Start)
a(n) = A162642(A003557(n)).
a(n) <= A056170(n).
(End)
Asymptotic mean: lim_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} 1/(p*(p+1)) = 0.3302299262... (A179119). - Amiram Eldar, Dec 25 2021

A380840 Decimal expansion of Sum_{p prime} 1/(p-1)^3.

Original entry on oeis.org

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

Offset: 1

Artur Jasinski, Mar 19 2025

nonn

			1.1475290977585800469332838..,

Cf. A085541, A086242, A085548, A136141, A152441, A154945, A179119, A369632, A324833.

PARI
```
sumeulerrat(1/(p-1)^3)
```

A036690 Product of a prime and the following number.

Original entry on oeis.org

6, 12, 30, 56, 132, 182, 306, 380, 552, 870, 992, 1406, 1722, 1892, 2256, 2862, 3540, 3782, 4556, 5112, 5402, 6320, 6972, 8010, 9506, 10302, 10712, 11556, 11990, 12882, 16256, 17292, 18906, 19460, 22350, 22952, 24806, 26732, 28056, 30102

Offset: 1

Felice Russo

The infinite sum over the reciprocals is given in A179119. - Wolfdieter Lang, Jul 10 2019

1/a(n) is the asymptotic density of numbers whose prime(n)-adic valuation is positive and even. - Amiram Eldar, Jan 23 2021

			a(3)=30 because prime(3)=5 and prime(3)+1=6, hence 5*6 = 30.

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

Cf. A036689, A034953, A065463, A179119, A306633.

[p^2+p: p in PrimesUpTo(250)]; // Vincenzo Librandi, Dec 19 2010

Table[(Prime[n] + 1) Prime[n], {n, 1, 100}] (* Artur Jasinski, Feb 06 2007 *)

a(n)=my(p=prime(n)); p*(p+1) \\ Charles R Greathouse IV, Mar 27 2020

a(n) = prime(n)*(prime(n)+1).
a(n) = A060800(n) - 1.
a(n) = 2*A034953(n). - Artur Jasinski, Feb 06 2007
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(2)/zeta(3) (A306633).
Product_{n>=1} (1 - 1/a(n)) = A065463. (End)
Sum_{n>=1} 1/a(n) = A179119. - R. J. Mathar, Mar 31 2025

cp's OEIS Frontend

A136141 Decimal expansion of Sum_{p prime} 1/(p*(p-1)).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

Extensions

A060687 Numbers k such that there exist exactly 2 Abelian groups of order k, i.e., A000688(k) = 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

PARI

Python

Python

Formula

Extensions

A034953 Triangular numbers (A000217) with prime indices.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Formula

A162642 Number of odd exponents in the canonical prime factorization of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Scheme

Formula

A023998 Number of block permutations on an n-set which are uniform, i.e., corresponding blocks have same size.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

Formula

Extensions