A190272 -id:A190272 - cp's OEIS frontend

A054377 Primary pseudoperfect numbers: numbers n > 1 such that 1/n + sum 1/p = 1, where the sum is over the primes p | n.

2, 6, 42, 1806, 47058, 2214502422, 52495396602, 8490421583559688410706771261086

Offset: 1

Primary pseudoperfect numbers are the solutions of the "differential equation" n' = n-1, where n' is the arithmetic derivative of n. - Paolo P. Lava, Nov 16 2009
Same as n > 1 such that 1 + sum n/p = n (and the only known numbers n > 1 satisfying the weaker condition that 1 + sum n/p is divisible by n). Hence a(n) is squarefree, and is pseudoperfect if n > 1. Remarkably, a(n) has exactly n (distinct) prime factors for n < 9. - Jonathan Sondow, Apr 21 2013
From the Wikipedia article: it is unknown whether there are infinitely many primary pseudoperfect numbers, or whether there are any odd primary pseudoperfect numbers. - Daniel Forgues, May 27 2013
Since the arithmetic derivative of a prime p is p' = 1, 2 is obviously the only prime in the sequence. - Daniel Forgues, May 29 2013
Just as 1 is not a prime number, 1 is also not a primary pseudoperfect number, according to the original definition by Butske, Jaje, and Mayernik, as well as Wikipedia and MathWorld. - Jonathan Sondow, Dec 01 2013
Is it always true that if a primary pseudoperfect number N > 2 is adjacent to a prime N-1 or N+1, then in fact N lies between twin primes N-1, N+1? See A235139. - Jonathan Sondow, Jan 05 2014
Also, integers n > 1 such that A069359(n) = n - 1. - Jonathan Sondow, Apr 16 2014

			From _Daniel Forgues_, May 24 2013: (Start)
With a(1) = 2, we have 1/2 + 1/2 = (1 + 1)/2 = 1;
with a(2) = 6 = 2 * 3, we have
  1/2 + 1/3 + 1/6 = (3 + 2 + 1)/6 = (1*3 + 3)/(2*3) = (1 + 1)/2 = 1;
with a(3) = 42 = 6 * 7, we have
  1/2 + 1/3 + 1/7 + 1/42 = (21 + 14 + 6 + 1)/42 =
  (3*7 + 2*7 + 7)/(6*7) = (3 + 2 + 1)/6 = 1;
with a(4) = 1806 = 42 * 43, we have
  1/2 + 1/3 + 1/7 + 1/43 + 1/1806 = (903 + 602 + 258 + 42 + 1)/1806 =
  (21*43 + 14*43 + 6*43 + 43)/(42*43) = (21 + 14 + 6 + 1)/42 = 1;
with a(5) = 47058 (not oblong number), we have
  1/2 + 1/3 + 1/11 + 1/23 + 1/31 + 1/47058 =
  (23529 + 15686 + 4278 + 2046 + 1518 + 1)/47058 = 1.
For n = 1 to 8, a(n) has n prime factors:
  a(1) = 2
  a(2) = 2 * 3
  a(3) = 2 * 3 *  7
  a(4) = 2 * 3 *  7 * 43
  a(5) = 2 * 3 * 11 * 23 *  31
  a(6) = 2 * 3 * 11 * 23 *  31 * 47059
  a(7) = 2 * 3 * 11 * 17 * 101 *   149 *       3109
  a(8) = 2 * 3 * 11 * 23 *  31 * 47059 * 2217342227 * 1729101023519
If a(n)+1 is prime, then a(n)*[a(n)+1] is also primary pseudoperfect. We have the chains: a(1) -> a(2) -> a(3) -> a(4); a(5) -> a(6). (End)
A primary pseudoperfect number (greater than 2) is oblong if and only if it is not the initial member of a chain. - _Daniel Forgues_, May 29 2013
If a(n)-1 is prime, then a(n)*(a(n)-1) is a Giuga number (A007850). This occurs for a(2), a(3), and a(5). See A235139 and the link "The p-adic order . . .", Theorem 8 and Example 1. - _Jonathan Sondow_, Jan 06 2014

M. A. Alekseyev, J. M. Grau, A. M. Oller-Marcen. Computing solutions to the congruence 1^n + 2^n + ... + n^n == p (mod n). Discrete Applied Mathematics, 2018. doi:10.1016/j.dam.2018.05.022 arXiv:1602.02407 [math.NT]
W. Butske, L. M. Jaje, and D. R. Mayernik, On the Equation Sum_{p|N} 1/p + 1/N = 1, Pseudoperfect numbers and partially weighted graphs, Math. Comput., 69 (1999), 407-420. [Title corrected by _Jonathan Sondow_, Apr 11 2012]
J. M. Grau, A. M. Oller-Marcen, and J. Sondow, On the congruence 1^m + 2^m + ... + m^m == n (mod m) with n|m, arXiv:1309.7941 [math.NT], 2013.
J. M. Grau, A. M. Oller-Marcén, and D. Sadornil, On µ-Sondow Numbers, arXiv:2111.14211 [math.NT], 2021.
John Machacek, Egyptian Fractions and Prime Power Divisors, arXiv:1706.01008 [math.NT], 2017.
J. Sondow and K. MacMillan, Reducing the Erdos-Moser equation 1^n + 2^n + . . . + k^n = (k+1)^n modulo k and k^2, Integers 11 (2011), #A34.
J. Sondow and K. MacMillan, Primary pseudoperfect numbers, arithmetic progressions, and the Erdos-Moser equation, Amer. Math. Monthly, 124 (2017) 232-240; arXiv:math/1812.06566 [math.NT], 2018.
J. Sondow and E. Tsukerman, The p-adic order of power sums, the Erdos-Moser equation, and Bernoulli numbers, arXiv:1401.0322 [math.NT], 2014; see section 4.
Eric Weisstein's World of Mathematics, Primary pseudoperfect number.
Wikipedia, Primary pseudoperfect number.
OEIS Wiki, Primary pseudoperfect numbers.

Cf. A005835, A007850, A069359, A168036, A190272, A191975, A203618, A216825, A216826, A230311, A235137, A235138, A235139, A236433.

pQ[n_] := (f = FactorInteger[n]; 1/n + Sum[1/f[[i]][[1]], {i, Length[f]}] == 1)
Select[Range[2, 10^6], pQ[#] &] (* Robert Price, Mar 14 2020 *)

isok(n) = if (n > 1, my(f=factor(n)[,1]); 1/n + sum(k=1, #f, 1/f[k]) == 1); \\ Michel Marcus, Oct 05 2017

from sympy import primefactors
A054377 = [n for n in range(2,10**5) if sum([n/p for p in primefactors(n)]) +1 == n] # Chai Wah Wu, Aug 20 2014

A031971(a(n)) (mod a(n)) = A233045(n). - Jonathan Sondow, Dec 11 2013
A069359(a(n)) = a(n) - 1. - Jonathan Sondow, Apr 16 2014
a(n) == 36*(n-2) + 6 (mod 288) for n = 2,3,..,8. - Kieren MacMillan and Jonathan Sondow, Sep 20 2017

A190275 Semiprimes of the form p*(p^2 - p + 1).

Original entry on oeis.org

6, 21, 301, 2041, 296341, 486877, 2666437, 3420301, 4304341, 7152001, 38159521, 42387097, 54296677, 95235601, 158048281, 229971241, 265434901, 383712781, 454166017, 775307917, 972261181, 1063290841, 1304557801, 1392422041, 1730882401, 1863895261, 2631883561, 2879450461, 3714274297, 3845297341, 4070454361, 4256780041, 4849695001, 5328809461, 5722533337, 5838483601, 7218898681, 7841065621

Offset: 1

Giorgio Balzarotti, May 07 2011

nonn

This sequence is infinite, assuming Schinzel's Hypothesis H.
Related to Rassias Conjecture ("for any odd prime p there are primes q < r such that p*q = q + r + 1") setting p = q. Generalization can be achieved by removing semiprimality condition and accepting p^e, e >= 2.
These are semiprimes m = p*q such that 1/p + 1/q - 1/m = p/q. Cf. A326690. - Amiram Eldar and Thomas Ordowski, Jul 22 2019

			a(1) = 6 = 2*3 = 2*(2^2-2+1).
a(2) = 21 = 3*7 = 3*(3^2-3+1).
a(3) = 301 = 7*43 = 7*(7^2-7+1).

Giovanni Resta, Table of n, a(n) for n = 1..10000
For Rassias conjecture: Preda Mihăilescu, Review of Problem Solving and Selected Topics in Number Theory, Newsletter of the European Mathematical Society, March 2011, p. 46.

Cf. A065508 (primes p such that p^2-p+1 is prime).

Cf. A001358 (semiprime), A003415 (arithmetic derivative), A164643, A190272 (n'=a-1), A190273 (n'=a+1), A190274 (n'=p^2-1).

seq(`if`(isprime((ithprime(i)^2-ithprime(i)+1))=true,(ithprime(i)^2-ithprime(i)+1)*ithprime(i),NULL),i=1..300);

p = Select[Prime@ Range@ 500, PrimeQ[#^2 - # + 1] &]; p (p^2 - p + 1) (* Giovanni Resta, Jul 22 2019 *)

forprime(p=2,1e4,if(isprime(k=p^2-p+1),print1(p*k", "))) \\ Charles R Greathouse IV, May 08 2011

A190273 Numbers n such that n' = m+1, with n and m semiprimes and gcd(m,n)>1, where n' is the arithmetic derivative of n.

Original entry on oeis.org

6, 10, 21, 26, 39, 55, 57, 74, 93, 111, 122, 146, 155, 201, 203, 253, 301, 305, 314, 327, 381, 386, 417, 471, 497, 543, 554, 597, 626, 633, 689, 737, 755, 791, 794, 842, 889, 905, 914, 921, 1011, 1027, 1055, 1081, 1082, 1137, 1226, 1227, 1322, 1346, 1379, 1461, 1466, 1477, 1497, 1514, 1623, 1655, 1703, 1711, 1713, 1731, 1751, 1754, 1893, 1967, 1994

Offset: 1

Giorgio Balzarotti, May 07 2011

nonn

The sequence is related to the Rassias Conjecture ("for any prime p there are two primes p1 and p2 such that p*p1=p1+p2+1, p>2, p2>p1", see A190272-A190275), because n = p1*p2, m=p1*p -> p1*p = p1+p2-1. The sequence includes the cases with p=p1 (or p2). Generalization can be achieved by removing semiprimarity condition or accepting gcd(n,m)=1. The differential equation in its general form n'=m+1 includes Giuga Numbers, i.e., n'=b*n+1, or n'=n+1 (A007850).

These are semiprimes n = p*q such that 1/p + 1/q - 1/n = P/Q, where P <> Q are primes. Cf. A326690. - Amiram Eldar and Thomas Ordowski, Jul 25 2019

			n=6, 6'=5, m=5+1=6, gcd(6,6)=6 -> a(1)=6

For Rassias conjecture: Preda Mihăilescu, Review of Problem Solving and Selected Topics in Number Theory, Newsletter of the European Mathematical Society, March 2011, p. 46.

Cf. A001358 (semiprimes), A003415 (arithmetic derivative), A007850 (Giuga numbers), A190272 (n'=m-1), A190273, A190274, A190275.

der:=n->n*add(op(2,p)/op(1,p),p=ifactors(n)[2]);
seq(`if`(bigomega(i)=2 and bigomega(der(i)-1)=2 and gcd(i,der(i)-1)>1,i,NULL),i=1..2000);

A190274 Numbers n such that n' = p^2-1, with n = semiprime = pq, n' is the arithmetic derivative of n. Also: semiprimes of the form p(p^2-p-1).

Original entry on oeis.org

15, 95, 287, 1199, 4607, 23519, 28799, 101567, 223199, 296207, 352799, 903167, 1019999, 2032127, 2230799, 2666159, 3285599, 5896799, 7606367, 13939199, 19392479, 28839887, 36154799, 46139039, 54295919, 62412767, 68250239, 73384079, 74440799, 90316799, 95234687, 109672319, 115263647, 118129199, 214562399, 223279487, 234503807, 236792879, 262963199, 270420767, 309829727, 355897439, 422999999, 486823247, 589884959, 628687487

Offset: 1

Giorgio Balzarotti, May 07 2011

nonn

The sequence shows similarity with the Rassias Conjecture ("for any prime p there are two primes p1 and p2 such that p*p1=p1+p2+1, p>2, p2>p1") with p1=p we have p*p=p+p2-1 (see A190272). Generalization can be achieved by removing semiprimarity condition and accepting p^e, e>=2.

			n=15, 15'=8, a=8+1=9=3^2 -> a(1)=15

European Mathematical Society, Newsletter (see book reviews), March 2011, page 46.

Cf. A001358 (semiprime), A003415 (arithmetic derivative), A190273 (n'=a-1), A190273 (n'=a+1).

seq(`if`(isprime((ithprime(i)^2-ithprime(i)-1))=true,(ithprime(i)^2-ithprime(i)-1)*ithprime(i),NULL),i=1..300);

cp's OEIS Frontend

A054377 Primary pseudoperfect numbers: numbers n > 1 such that 1/n + sum 1/p = 1, where the sum is over the primes p | n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A190275 Semiprimes of the form p*(p^2 - p + 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A190273 Numbers n such that n' = m+1, with n and m semiprimes and gcd(m,n)>1, where n' is the arithmetic derivative of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A190274 Numbers n such that n' = p^2-1, with n = semiprime = p*q, n' is the arithmetic derivative of n. Also: semiprimes of the form p*(p^2-p-1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A190274 Numbers n such that n' = p^2-1, with n = semiprime = pq, n' is the arithmetic derivative of n. Also: semiprimes of the form p(p^2-p-1).