A005835 -id:A005835 - 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

A006038 Odd primitive abundant numbers.

Original entry on oeis.org

945, 1575, 2205, 3465, 4095, 5355, 5775, 5985, 6435, 6825, 7245, 7425, 8085, 8415, 8925, 9135, 9555, 9765, 11655, 12705, 12915, 13545, 14805, 15015, 16695, 18585, 19215, 19635, 21105, 21945, 22365, 22995, 23205, 24885, 25935, 26145, 26565, 28035, 28215

Offset: 1

N. J. A. Sloane

nonn

Dickson proves that there are only a finite number of odd primitive abundant numbers having n distinct prime factors. Sequence A188342 lists the smallest such numbers. - T. D. Noe, Mar 28 2011
Sequence A188439 sorts the numbers in this sequence by the number of distinct prime factors. Eight numbers have exactly three prime factors; 576 have exactly four prime factors. - T. D. Noe, Apr 04 2011
Any multiple of an abundant number (A005101) is again an abundant number. Primitive abundant numbers (A091191) are those not of this form, i.e., without an abundant proper divisor. We don't know any odd perfect number (A000396), so the (odd) terms here have only deficient proper divisors (A071395), and their prime factors p are less than sigma(n/p)/deficiency(n/p). See A005231 (odd abundant numbers) for an explanation why all terms have at least 3 distinct prime factors, and 5 prime factors when counted with multiplicity (A001222), whence a(1) = 3^3*5*7. All known terms are semiperfect (A005835, and thus in A006036): no odd weird number (A006037) is known, but if it exists, the smallest one is in this sequence. - M. F. Hasler, Jul 28 2016
So far, a(173) = 351351 is the only known term of A122036, i.e., which can't be written as sum of its proper divisors > 1. - M. F. Hasler, Jan 26 2020

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..10000
L. E. Dickson, Finiteness of the odd perfect and primitive abundant numbers with n distinct prime factors, American Journal of Mathematics 35 (1913), pp. 413-422.
R. K. Guy, Letter to N. J. A. Sloane with attachment, Jun. 1991
Jacob Liddy, An algorithm to determine all odd primitive abundant numbers with d prime divisors, Honors Research Projects (2018), 728.
Eric Weisstein's World of Mathematics, Primitive Abundant Number

Cf. A005101, A005231. Subsequence of A091191.

Cf. A000203, A027751, A379949 (subsequence of square terms).

a006038 n = a006038_list !! (n-1)
a006038_list = filter f [1, 3 ..] where
   f x = sum pdivs > x && all (<= 0) (map (\d -> a000203 d - 2 * d) pdivs)
         where pdivs = a027751_row x
-- Reinhard Zumkeller, Jan 31 2014

isA005101 := proc(n) is(numtheory[sigma](n) > 2*n ); end proc:
isA005100 := proc(n) is(numtheory[sigma](n) < 2*n ); end proc:
isA006038 := proc(n) local d; if type(n,'odd') and isA005101(n) then for d in numtheory[divisors](n) minus {1,n} do if not isA005100(d) then return false; end if; end do: return true;else false; end if; end proc:
n := 1 ; for a from 1 by 2 do if isA006038(a) then printf("%d %d\n",n,a) ; n := n+1 ; end if; end do: # R. J. Mathar, Mar 28 2011

t = {}; n = 1; While[Length[t] < 50, n = n + 2; If[DivisorSigma[1, n] > 2 n && Intersection[t, Divisors[n]] == {}, AppendTo[t, n]]]; t (* T. D. Noe, Mar 28 2011 *)

is(n)=n%2 && sumdiv(n,d,sigma(d,-1)>2)==1 \\ Charles R Greathouse IV, Jun 10 2013

is_A006038(n)=bittest(n,0) && sigma(n)>2*n && !for(i=1,#f=factor(n)[,1],sigma(n\f[i],-1)>2&&return) \\ More than 5 times faster. - M. F. Hasler, Jul 28 2016

A181595 Abundant numbers n for which the abundance d = sigma(n) - 2*n is a proper divisor, that is, 0 < d < n and d | n.

Original entry on oeis.org

12, 18, 20, 24, 40, 56, 88, 104, 196, 224, 234, 368, 464, 650, 992, 1504, 1888, 1952, 3724, 5624, 9112, 11096, 13736, 15376, 15872, 16256, 17816, 24448, 28544, 30592, 32128, 77744, 98048, 122624, 128768, 130304, 174592, 396896, 507392, 521728, 522752, 537248

Offset: 1

Vladimir Shevelev, Nov 01 2010

nonn

Named near-perfect numbers by sequence author.
Union of this sequence and A005820 is A153501.
Every even perfect number n = 2^(p-1)*(2^p-1), p and 2^p-1 prime, of A000396 generates three entries: 2*n, 2^p*n and (2^p-1)*n.
Every number M=2^(t-1)*P, where P is a prime of the form 2^t-2^k-1, is an entry for which (2^k)|M and sigma(M)-2^k=2*M (see A181701).
Conjecture 1: For every k>=1, there exist infinitely many entries m for which (2^k)|m and sigma(m)-2^k = 2*m.
Conjecture 2. All entries are even. [Proved to be false, see below. (Ed.)]
Conjecture 3. If the suitable (according to the definition) divisor d of an entry is not a power of 2, then it is not suitable divisor for any other entry.
Conjecture 4. If a suitable divisor for an even entry is odd, then it is a Mersenne prime (A000043).
If Conjectures 3 and 4 are true, then an entry with odd suitable divisor has the form 2^(p-1)*(2^p-1)^2, where p and 2^p-1 are primes. - Vladimir Shevelev, Nov 08 2010 to Dec 16 2010
The only odd term in this sequence < 2*10^12 is 173369889. - Donovan Johnson, Feb 15 2012
173369889 remains only odd term up to 1.4*10^19. - Peter J. C. Moses, Mar 05 2012
These numbers are obviously pseudoperfect (A005835) since they are equal to the sum of all the proper divisors except the one that is the same as the abundance. - Alonso del Arte, Jul 16 2012

			The abundance of 12 is A033880(12) = 4, which is a proper divisor of 12, so 12 is in the sequence.

Donovan Johnson, Table of n, a(n) for n = 1..200
Hùng Việt Chu, Divisibility of Divisor Functions of Even Perfect Numbers, J. Int. Seq., Vol. 24 (2021), Article 21.3.4.
Yanbin Li and Qunying Liao, A class of new near-perfect numbers, J. Korean Math. Soc. 52 (2015), No. 4, pp. 751-763.
Paul Pollack and Vladimir Shevelev, On perfect and near-perfect numbers, J. Number Theory 132 (2012), pp. 3037-3046. arXiv preprint, arXiv:1011.6160 [math.NT], 2010-2012.
X.-Z. Ren, Y.-G. Chen, On near-perfect numbers with two distinct prime factors, Bulletin of the Australian Mathematical Society, No 3 (2013), available on CJO2013. doi:10.1017/S0004972713000178.
M. Tang, X. Z. Ren and M. Li, On near-perfect and deficient-perfect numbers, Colloq. Math. 133 (2013), 221-226.

Cf. A000396, A005101, A153501, A005820.

q:= n-> (t-> t>0 and tAlois P. Heinz, May 11 2023

Select[Range[550000], 0 < (d = DivisorSigma[1, #] - 2*#) < # && Divisible[#, d] &] (* Amiram Eldar, May 12 2023 *)

is_A181595(n)=my(d=sigma(n)-2*n); (d>0) && (dA181595(n)&&print1(n","))  \\ M. F. Hasler, Apr 14 2012; corrected by Michel Marcus, May 12 2023

Definition shortened, entries checked by R. J. Mathar, Nov 17 2010

A065205 Number of subsets of proper divisors of n that sum to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 1, 0, 0, 0, 5, 0, 0, 0, 1, 0, 3, 0, 0, 0, 0, 0, 7, 0, 0, 0, 3, 0, 2, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 3, 0, 2, 0, 0, 0, 34, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 31, 0, 0, 0, 0, 0, 1, 0, 6, 0, 0, 0, 25, 0, 0, 0, 1, 0, 23, 0, 0, 0, 0, 0, 21, 0, 0, 0, 2

Offset: 1

Jonathan Ayres (jonathan.ayres(AT)btinternet.com), Oct 19 2001

nonn

Deficient and weird numbers have a(n) = 0, perfect numbers and others (see A064771) have a(n) = 1.

Number of partitions of n into distinct proper divisors of n; a(A136447(n)) = 0; a(A005835(n)) > 0; a(A064771(n)) = 1. - Reinhard Zumkeller, Jan 21 2013

			a(20) = 1 because {1, 4, 5, 10} is the only subset of proper divisors of 20 that sum to 20.
a(24) = 5 because there are five different subsets we can use to sum up to 24: {1, 2, 3, 4, 6, 8}, {1, 2, 3, 6, 12}, {1, 3, 8, 12}, {2, 4, 6, 12}, {4, 8, 12}.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)

Cf. A064771, A005835.
Cf. A065218 for records.
Cf. A027751, A210442, A211110, A033630.

a065205 n = p (a027751_row n) n where
   p _      0 = 1
   p []     _ = 0
   p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
-- Reinhard Zumkeller, Jan 21 2013

a[n_] := (dd = Most[ Divisors[n] ]; cc = Array[c, Length[dd]]; Length[ {ToRules[ Reduce[ And @@ (0 <= # <= 1 &) /@ cc && dd . cc == n, cc, Integers]]}]); Table[ a[n], {n, 1, 100}] (* Jean-François Alcover, Feb 23 2012 *)

a(n,s,d)={s || (s=sigma(n)-n) || return; d||d=vecextract(divisors(n),"^-1"); while(d[#d]>n, s-=d[#d]; d=d[1..-2]); s<=n && return(s==n); if( n>d[#d], a(n-d[#d],s-d[#d],d[1..-2]), 1)+a(n,s-d[#d],d[1..-2])} \\ M. F. Hasler, May 11 2015

a(n) = A033630(n) - 1.

More terms and additional comments from Jud McCranie, Oct 21 2001

A064771 Let S(n) = set of divisors of n, excluding n; sequence gives n such that there is a unique subset of S(n) that sums to n.

Original entry on oeis.org

6, 20, 28, 78, 88, 102, 104, 114, 138, 174, 186, 222, 246, 258, 272, 282, 304, 318, 354, 366, 368, 402, 426, 438, 464, 474, 490, 496, 498, 534, 572, 582, 606, 618, 642, 650, 654, 678, 748, 762, 786, 822, 834, 860, 894, 906, 940, 942, 978, 1002, 1014, 1038

Offset: 1

Jonathan Ayres (jonathan.ayres(AT)btinternet.com), Oct 19 2001

Perfect numbers (A000396) are a proper subset of this sequence. Weird numbers (A006037) are numbers whose proper divisors sum to more than the number, but no subset sums to the number.
Odd elements are rare: the first few are 8925, 32445, 351351, 442365; there are no more below 100 million. See A065235 for more details.
A065205(a(n)) = 1. - Reinhard Zumkeller, Jan 21 2013

			Proper divisors of 20 are 1, 2, 4, 5 and 10. {1,4,5,10} is the only subset that sums to 20, so 20 is in the sequence.

Giovanni Resta, Table of n, a(n) for n = 1..10000 (terms 1..200 from T. D. Noe, terms 201..5000 from Amiram Eldar)

A005835 gives n such that some subset of S(n) sums to n. Cf. A065205.
Cf. A006037, A065205, A378448 (characteristic function).
Subsequences: A000396, A065235 (odd terms), A378519, A378530.
Cf. A027751.

a064771 n = a064771_list !! (n-1)
a064771_list = map (+ 1) $ elemIndices 1 a065205_list
-- Reinhard Zumkeller, Jan 21 2013

filter:= proc(n)
  local P,x,d;
  P:= mul(x^d+1, d = numtheory:-divisors(n) minus {n});
  coeff(P,x,n) = 1
end proc:
select(filter, [$1..2000]); # Robert Israel, Sep 25 2024

okQ[n_]:= Module[{d=Most[Divisors[n]]}, SeriesCoefficient[Series[ Product[ 1+x^i, {i, d}], {x, 0, n}], n] == 1];Select[ Range[ 1100],okQ] (* Harvey P. Dale, Dec 13 2010 *)

from sympy import divisors
def isok(n):
    dp = {0: 1}
    for d in divisors(n)[:-1]:
        u = {}
        for k in dp.keys():
            if (s := (d + k)) <= n:
                u[s] = dp.get(s, 0) + dp[k]
                if s == n and u[s] > 1:
                    return False
        for k,v in u.items():
            dp[k] = v
    return dp.get(n, 0) == 1
print([n for n in range(1, 1039) if isok(n)]) # Darío Clavijo, Sep 17 2024

More terms from Don Reble, Jud McCranie and Naohiro Nomoto, Oct 22 2001

A103977 Zumkeller deficiency of n: Let d_1 ... d_k be the divisors of n. Then a(n) = min_{ e_1 = +-1, ... e_k = +-1 } | Sum_i e_i d_i |.

Original entry on oeis.org

1, 1, 2, 1, 4, 0, 6, 1, 5, 2, 10, 0, 12, 4, 6, 1, 16, 1, 18, 0, 10, 8, 22, 0, 19, 10, 14, 0, 28, 0, 30, 1, 18, 14, 22, 1, 36, 16, 22, 0, 40, 0, 42, 4, 12, 20, 46, 0, 41, 7, 30, 6, 52, 0, 38, 0, 34, 26, 58, 0, 60, 28, 22, 1, 46, 0, 66, 10, 42, 0, 70, 1, 72, 34, 26, 12, 58, 0, 78, 0

Offset: 1

Yasutoshi Kohmoto, Jan 01 2007

nonn

Like the ordinary deficiency (A033879) obtains 0's only at perfect numbers (A000396), the Zumkeller deficiency obtains 0's only at integer-perfect numbers, A083207. See the formula section. Unlike the ordinary deficiency, this obtains only nonnegative values. See A378600 for another version. - Antti Karttunen, Dec 03 2024

			a(6) = 1 + 2 + 3 - 6 = 0.

Antti Karttunen, Table of n, a(n) for n = 1..20000 (first 10000 terms from Amiram Eldar)

Cf. A125732, A125733, A005835, A023196, A033879, A083206, A083207 (positions of 0's), A263837, A378643 (Dirichlet inverse), A378644 (Möbius transform), A378645, A378646, A378647 (an analog of A000027), A378648 (an analog of sigma), A378649 (an analog of Euler phi), A379503 (positions of 1's), A379504, A379505.
Cf. A378600 (signed variant).
Cf. also A058377, A119347.

A103977 := proc(n) local divs,a,acandid,filt,i,p,sigs ; divs := convert(numtheory[divisors](n),list) ; a := add(i,i=divs) ; for sigs from 0 to 2^nops(divs)-1 do filt := convert(sigs,base,2) ; while nops(filt) < nops(divs) do filt := [op(filt), 0] ; od ; acandid := 0 ; for p from 0 to nops(divs)-1 do if op(p+1,filt) = 0 then acandid := acandid-op(p+1,divs) ; else acandid := acandid+op(p+1,divs) ; fi ; od: acandid := abs(acandid) ; if acandid < a then a := acandid ; fi ; od: RETURN(a) ; end: seq(A103977(n),n=1..80) ; # R. J. Mathar, Nov 27 2007
# second Maple program:
a:= proc(n) option remember; local l, b; l, b:= [numtheory[divisors](n)[]],
      proc(s, i) option remember; `if`(i<1, s,
        min(b(s+l[i], i-1), b(abs(s-l[i]), i-1)))
      end: b(0, nops(l))
    end:
seq(a(n), n=1..80);  # Alois P. Heinz, Dec 05 2024

a[n_] := Module[{d = Divisors[n], c, p, m}, c = CoefficientList[Product[1 + x^i, {i, d}], x]; p = -1 + Position[c, ?(# > 0 &)] // Flatten; m = Length[p]; If[OddQ[m], If[(d = p[[(m + 1)/2]] - p[[(m - 1)/2]]) == 1, 0, d], p[[m/2 + 1]] - p[[m/2]]]]; Array[a, 100] (* _Amiram Eldar, Dec 11 2019 *)

nonzerocoefpositions(p) = { my(v=Vec(p), lista=List([])); for(i=1,#v,if(v[i], listput(lista,i))); Vec(lista); }; \\ Doesn't need to be 0-based, as we use their differences only.
A103977(n) = { my(p=1); fordiv(n, d, p *= (1 + 'x^d)); my(plist=nonzerocoefpositions(p), m = #plist, d); if(!(m%2), plist[1+(m/2)]-plist[m/2], d = plist[(m+1)/2]-plist[(m-1)/2]; if(1==d,0,d)); }; \\ Antti Karttunen, Dec 03 2024, after Mathematica-program by Amiram Eldar

If n=p (prime), then a(n)=p-1. If n=2^m, then a(n)=1. [Corrected by R. J. Mathar, Nov 27 2007]
a(n) = 0 iff n is a Zumkeller number (A083207). - Amiram Eldar, Jan 05 2020
From Antti Karttunen, Dec 03 2024: (Start)
a(n) = A033879(n) iff n is a non-abundant number (A263837).
a(n) = abs(A378600(n)).
a(n) = 2*A378647(n) - A378648(n). [Analogously to A033879(n) = 2*n - sigma(n)]
a(n) = 0 <=> A083206(n) > 0.
(End)
a(p^e) = p^e - (1+p+...+p^(e-1)) = (p^e*(p-2) + 1)/(p-1) for prime p. - Jianing Song, Dec 05 2024
a(n) = 1 <=> A379504(n) > 0. - Antti Karttunen, Jan 07 2025

More terms from R. J. Mathar, Nov 27 2007

Name "Zumkeller deficiency" coined by Antti Karttunen, Dec 03 2024

A083209 Numbers whose divisors can be partitioned in exactly one way into two disjoint sets with the same sum.

Original entry on oeis.org

6, 12, 20, 28, 56, 70, 88, 104, 176, 208, 272, 304, 368, 464, 496, 550, 650, 736, 836, 928, 992, 1184, 1312, 1376, 1504, 1696, 1888, 1952, 2752, 3008, 3230, 3392, 3770, 3776, 3904, 4030, 4288, 4510, 4544, 4672, 5056, 5170, 5312, 5696, 5830, 6208, 6464

Offset: 1

Reinhard Zumkeller, Apr 22 2003

nonn

A083206(a(n))=1; perfect numbers (A000396) are a subset; problem: are weird numbers (A006037) a subset?
The weird numbers A006037 are not a subset of this sequence. The first missing weird number is A006037(8) = 10430. - Alois P. Heinz, Oct 29 2009
All numbers of the form p*2^k are in this sequence for k>0 and odd primes p between 2^(k+1)/3 and 2^(k+1). - T. D. Noe, Jul 08 2010
"Numbers with exactly one subset of their sets of divisors such that the complement has the same sum." - This was the original name of the sequence, but strictly taken is incorrect, because there are always two subsets that satisfy this condition: the subset and its complement. - Antti Karttunen, Dec 02 2024

			n=20: 2+4+5+10 = 1+20, 20 is a term (A083206(20)=1).

T. D. Noe, Table of n, a(n) for n=1..407 (terms < 10^6)
Eric Weisstein's World of Mathematics, Perfect Number.
Eric Weisstein's World of Mathematics, Weird Number.
Reinhard Zumkeller, Illustration of initial terms
Index entries for sequences where odd perfect numbers must occur, if they exist at all

Subsequence of A083207, Zumkeller numbers.
Positions of 1's in A083206.
Cf. A005101, A005835, A064771, A337739 (terms with record number of divisors), A378449 (characteristic function), A378530 (subsequence).
Cf. also A378652, and A335143, A335199, A335202, A335219, A335217, A339980 for variants.

with(numtheory): b:= proc(n,l) option remember; local m, ll, i; m:= nops(l); if n<0 then 0 elif n=0 then 1 elif m=0 or add(i, i=l)Alois P. Heinz, Oct 29 2009

b[n_, l_] := b[n, l] = Module[{m, ll, i}, m = Length[l]; Which[n<0, 0, n == 0, 1, m == 0 || Total[l] Nothing]; b[n, ll] + b[n - l[[m]], ll]]]; a[n_] := a[n] = Module[{i, k, l, m, r}, For[k = If[n == 1, 1, a[n-1]+1], True, k++, l = Divisors[k]; {m, r} = QuotientRemainder[Total[l], 2]; If[r==0 && b[m, l]==2, Break[]]]; k]; Table[Print["a(", n, ") = ", a[n]]; a[n], {n, 1, 50}] (* Jean-François Alcover, Jan 31 2017, after Alois P. Heinz *)

isA083209 = A378449; \\ Antti Karttunen, Nov 28 2024

More terms from Alois P. Heinz, Oct 29 2009

Improved the definition, old name moved to the comments - Antti Karttunen, Dec 02 2024

A293188 Unitary pseudoperfect numbers: numbers that equal to the sum of a subset of their aliquot unitary divisors.

Original entry on oeis.org

6, 30, 42, 60, 66, 78, 90, 102, 114, 138, 150, 174, 186, 210, 222, 246, 258, 282, 294, 318, 330, 354, 366, 390, 402, 420, 426, 438, 462, 474, 498, 510, 534, 546, 570, 582, 606, 618, 630, 642, 654, 660, 678, 690, 714, 726, 750, 762, 770, 780, 786, 798, 822

Offset: 1

Amiram Eldar, Oct 01 2017

nonn

Supersequence of A002827.

The nonsquarefree terms are 60, 90, 150, 294, 420, 630, 660, 726, 750, 780, 840, ...

			150 is in the sequence since its unitary aliquot divisors are 1, 2, 3, 6, 25, 50, 75 and 150 = 25 + 50 + 75.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A002827, A005835, A034683, A064114.

udiv[n_]:=Block[{d=Divisors[n]},Select[d,GCD[#,n/#]==1&]]; a={};n=0;While[Length[a]<100,n++;d=Most[udiv[n]];c = SeriesCoefficient[ Series[ Product[1+x^d[[i]],{i,Length[d]} ],{x,0,n}], n]; If[c>0,AppendTo[a,n]]];a (* after T. D. Noe at A005835 *)

A087167 Odd numbers such that sigma(n) - 2n = 6.

Original entry on oeis.org

8925, 32445, 442365

Offset: 1

Farideh Firoozbakht, Oct 19 2003

If m is in this sequence and 5 doesn't divide m then m is an odd Weird number. There are no other terms up to 2*10^9. Jud McCranie wrote: There are no terms between 2*10^9 and 6.5*10^9.
a(4) > 10^12. - Donovan Johnson, Dec 08 2011
a(4) > 10^13. - Giovanni Resta, Mar 29 2013
a(4) > 10^22. - Wenjie Fang, Jun 16 2014
Any term x of this sequence can be combined with any term y of A141548 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

			a(1)=8925 because sigma(8925)=2*8925+6 and 8925 is the first odd number such that sigma(n)-2n=6.

R. K. Guy, "Almost Perfect, Quasi-Perfect, Pseudoperfect, Harmonic, Weird, Multiperfect and Hyperperfect Numbers." B2 in Unsolved Problems in Number Theory, 2nd ed.New York:Springer- Verlag, pp. 45-53, 1994.

S. Benkoski, Are All Weird Numbers Even?, Problem E2308, Amer. Math. Monthly, 79 (7) (1972), 774.
C. Rivera, Puzzle 233. A little twist.
Eric Weisstein's World of Mathematics, Weird Number.

Cf. A003380, A077374, A005101, A005835, A141548 (deficiency 6).

Do[If[OddQ[n] && DivisorSigma[1, n] - 2n == 6, Print[n]], {n, 2*10^9}]

is(n)=n%2 && sigma(n)==2*n+6 \\ Charles R Greathouse IV, Mar 09 2014

A065218 Consider the subsets of proper divisors of a number that sum to the number. These are numbers that set a record number of such subsets.

Original entry on oeis.org

1, 6, 12, 24, 36, 48, 60, 120, 180, 240, 360, 720, 840, 1260, 1680, 2520, 5040, 7560, 10080, 15120, 20160, 25200, 27720, 45360, 50400, 55440, 83160, 110880, 166320, 221760, 277200, 332640, 498960, 554400, 665280, 720720, 831600, 1081080, 1441440

Offset: 1

Jud McCranie, Oct 21 2001

nonn

Indices of records in A065205 and A033630. The corresponding records (number of subsets) are in A065219.

This sequence is not a subset of A002182: 831600 belongs to this sequence but not A002182.

			Proper divisors of 12 are {1, 2, 3, 4, 6}. Two subsets of this sum to 12: {2, 4, 6} and {1, 2, 3, 6} - more than any smaller number, so 12 is in the sequence.

Michael De Vlieger, Correlation of A065218 and A065219.

Cf. A065219, A064771, A063205, A005835, A002182, A004394.

With[{s = Table[-1 + SeriesCoefficient[Series[Times @@ ((1 + z^#) & /@ Divisors[n]), {z, 0, n}], n], {n, 2520}]}, FirstPosition[s, #][[1]] & /@ Union@ FoldList[Max, s]] (* Michael De Vlieger, Oct 10 2017 *)

More terms from Franklin T. Adams-Watters, Nov 27 2006
Edited and extended by Max Alekseyev, May 29 2009
Offset changed by Andrey Zabolotskiy, Oct 10 2017

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

A006038 Odd primitive abundant numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

A181595 Abundant numbers n for which the abundance d = sigma(n) - 2*n is a proper divisor, that is, 0 < d < n and d | n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A065205 Number of subsets of proper divisors of n that sum to n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A064771 Let S(n) = set of divisors of n, excluding n; sequence gives n such that there is a unique subset of S(n) that sums to n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Python

Extensions

A103977 Zumkeller deficiency of n: Let d_1 ... d_k be the divisors of n. Then a(n) = min_{ e_1 = +-1, ... e_k = +-1 } | Sum_i e_i d_i |.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica