cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A156942 Odd abundant numbers whose abundance is odd.

Original entry on oeis.org

11025, 99225, 245025, 275625, 342225, 540225, 893025, 1334025, 1863225, 2205225, 2480625, 3080025, 3186225, 3980025, 4601025, 4862025, 5832225, 6125625, 6890625, 7868025, 8037225, 8555625, 9272025, 9828225, 10595025, 10989225
Offset: 1

Views

Author

Robert G. Wilson v, Feb 18 2009

Keywords

Comments

Number of terms <10^n: 0, 0, 0, 0, 2, 7, 24, 83, 250, 792, 2484, 7988, 25383, 80082, ..., . Not all are a multiple of 25, i.e.; 81162081 = 9009^2 = (9*7*11*13)^2. See A156943.
Any term must be an odd square. Square roots are in A174830.
Indeed, the sum of divisors of any number isn't odd unless it's a square or twice a square (A028982), and to get the abundance, twice the number is subtracted, so the parity remains the same. - M. F. Hasler, Jan 26 2020
Question: Is this a subsequence of A379503? (Is A379504(a(n)) > 0 for all n? See A379951). The first 15000 terms are all included there. - Amiram Eldar and Antti Karttunen, Jan 06 2025
Question 2: Is A379505(a(n)) > 1 for all n, especially if there are no quasiperfect numbers (numbers k such that sigma(k) = 2k+1)? - Antti Karttunen, Jan 06 2025
From Amiram Eldar, Jan 16 2025: (Start)
The least term that is not divisible by 5 is a(75) = 81162081.
The least term that is not divisible by 3 is a(296889) = 1382511906801025.
The least term that is coprime to 15 is 15285071557677427358507559514565648611799881. (End)

Crossrefs

Subsequences: A156943, A325311 (thus also A379490), A347890, A379949 (terms that are primitive abundant).

Programs

  • Mathematica
    fQ[n_] := Block[{ds = DivisorSigma[1, n] - 2 n}, ds > 0 && OddQ@ ds]; Select[ Range[1, 12006223, 2], fQ @# &]
  • PARI
    is(n)=my(s=sigma(n)); n%2 && s>2*n && (s-2*n)%2 \\ Charles R Greathouse IV, Feb 21 2017

Formula

a(n) = A174830(n)^2. - M. F. Hasler, Jan 26 2020

Extensions

Edited by Robert G. Wilson v at the suggestion of T. D. Noe, Mar 30 2010

A379949 Primitive abundant numbers (A091191) that are odd squares.

Original entry on oeis.org

342225, 1029447225, 1757705625, 2177622225, 14787776025, 18114198921, 32871503025, 45018230625, 150897287025, 245485566225, 296006724225, 705373218225, 1126920249225, 1329226832241, 1358425215225, 1545732725625, 1555265892609, 1783322538921, 2811755495241, 4627123655625, 5248080775161, 6140855705625, 7683069267225
Offset: 1

Views

Author

Antti Karttunen, Jan 07 2025

Keywords

Comments

Question: Does A379504(.) obtain generally smaller values for the terms of this subsequence of A156942 than for its non-primitive terms? (See A379951, with A379951(5) = 5969, where A156942(5) = 342225, the first term of this sequence). Is A103977(.) = 1 for all terms, i.e., is this a subsequence of A379503?

Crossrefs

Cf. A103977, A379504, A379950 (square roots).
Intersection of A016754 and A091191.
Intersection of A006038 and A156942.
Subsequences of the following sequences: A306796 (odd terms, but only if there are no odd perfect numbers), A363176, A379503 (conjectured).

Programs

  • PARI
    is_A379949(n) = if(!(n%2) || !issquare(n) || sigma(n)<=2*n, 0, fordiv(n, d, if(d>1 && sigma(n/d, -1)>2, return(0))); (1));
    
  • PARI
    is1(k) = {my(f = factor(k)); for(i = 1, #f~, f[i, 2] *= 2); if(sigma(f, -1) <= 2, return(0)); for(i = 1, #f~, f[i, 2] -= 1; if(sigma(f, -1) > 2, return(0)); f[i, 2] += 1); 1;}
    list(lim) = forstep(k = 1, lim, 2, if(is1(k), print1(k^2, ", "))); \\ Amiram Eldar, Mar 12 2025

Formula

a(n) = A379950(n)^2.
Showing 1-2 of 2 results.