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.

A351554 Numbers k such that there are no odd prime factors p of sigma(k) such that p does not divide A003961(k) and the valuation(k, p) is different from valuation(sigma(k), p), where A003961 is fully multiplicative with a(p) = nextprime(p), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 2, 3, 6, 7, 10, 14, 15, 20, 21, 22, 24, 27, 28, 30, 31, 33, 34, 40, 42, 46, 54, 57, 60, 62, 66, 69, 70, 84, 87, 91, 93, 94, 102, 105, 106, 110, 114, 120, 127, 130, 138, 140, 141, 142, 154, 160, 168, 170, 174, 177, 182, 186, 189, 190, 195, 198, 210, 214, 216, 217, 220, 224, 230, 231, 237, 238, 254, 260, 264, 270, 273
Offset: 1

Views

Author

Antti Karttunen, Feb 16 2022

Keywords

Comments

Numbers k for which A351555(k) = 0. This is a necessary condition for the terms of A349169 and of A349745, therefore they are subsequences of this sequence.
All six known 3-perfect numbers (A005820) are included in this sequence.
All 65 known 5-multiperfects (A046060) are included in this sequence.
Moreover, all multiperfect numbers (A007691) seem to be in this sequence.
From Antti Karttunen, Aug 27 2025: (Start)
Multiperfect number m is included in this sequence only if its abundancy sigma(m)/m has only such odd prime factors p that prevprime(p) [A151799] divides m for each p. E.g., all 65 known 5-multiperfects are multiples of 3, and all known terms of A005820 and A046061 are even.
This sequence contains natural numbers k such that the odd primes in the prime factorization of sigma(k) have the same valuation there as in k, except that the primes in A003961(k) [or equally in A003961(A007947(k))] stand for "don't care primes", that are "masked off" from the comparison.
(End)

Crossrefs

Positions of zeros in A351555.
Subsequences: A000396, A351553 (even terms), A386430 (odd terms), A351551, A349169, A349745, A387160 (terms of the form prime * m^2), also these, at least all the currently (Feb 2022) known terms: A005820, A007691, A046060.

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A351555(n) = { my(s=sigma(n),f=factor(s),u=A003961(n)); sum(k=1,#f~,if((f[k,1]%2) && 0!=(u%f[k,1]), (valuation(n,f[k,1])!=f[k,2]), 0)); };
    isA351554(n) = (0==A351555(n));
    
  • PARI
    isA351554(n) = { my(sh=A351546(n),f=factor(sh)); for(i=1,#f~, if((f[i,1]%2)&&valuation(n,f[i,1])!=f[i,2],return(0))); (1); }; \\ Uses also program given in A351546.

Extensions

Definition corrected by Antti Karttunen, Aug 22 2025

A386430 Odd numbers k such that there are no prime factors p of sigma(k) such that p does not divide A003961(k) and the valuation(k, p) is different from valuation(sigma(k), p), where A003961 is fully multiplicative with a(p) = nextprime(p), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 7, 15, 21, 27, 31, 33, 57, 69, 87, 91, 93, 105, 127, 141, 177, 189, 195, 217, 231, 237, 273, 285, 301, 381, 399, 447, 465, 483, 495, 513, 567, 573, 597, 609, 627, 651, 717, 775, 819, 837, 861, 889, 903, 987, 1023, 1029, 1149, 1185, 1239, 1311, 1365, 1419, 1431, 1437, 1455, 1497, 1561, 1653, 1659, 1687, 1743
Offset: 1

Views

Author

Antti Karttunen, Aug 22 2025

Keywords

Comments

Conjecture: After the initial 1, and apart from any hypothetical odd perfect numbers, all other terms are in A248150, i.e., sigma(k) == 0 (mod 4). This would imply (with the same caveat), that this sequence has no common terms with A228058 and no squares larger than one. This is true at least for the first 709203 terms (terms in range [1..2^34]).
Terms k such that A162642(k) = 1 are rare: 3, 7, 27, 31, 127, 567, 775, 8191, 27783, 131071, 524287, 2147483647, ... (odd terms of A387160).

Examples

			a(386548) = 5919068925 = 3^4 * 5^2 * 7^2 * 11^2 * 17 * 29. sigma(5919068925) = 15355618740 = 2^2 * 3^4 * 5 * 7 * 11^2 * 19^2 * 31. The "don't care primes" is given by A003961(A007947(5919068925))) = 2947945 = 5*7*11*13*19*31, thus only odd prime factor that matters here is 3, which in case has the same exponent (4) in both n = 5919068925 and sigma(n). In a way, this number is very close to satisfying Euler's criterion for odd perfect numbers (A228058), except that it has two unitary prime factors of the form 4k+1, instead of just one, apart from the square factor. Both n/17 and n/29 are in A228058.
		

Crossrefs

Odd terms of A351554.
Cf. A349169 (subsequence).

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A351555(n) = { my(s=sigma(n),f=factor(s),u=A003961(n)); sum(k=1,#f~,if((f[k,1]%2) && 0!=(u%f[k,1]), (valuation(n,f[k,1])!=f[k,2]), 0)); };
    isA386430(n) = ((n%2) && (0==A351555(n)));

Formula

{k | k odd, A351555(k) = 0}.
Showing 1-2 of 2 results.