A191218 -id:A191218 - cp's OEIS frontend

A228058 Odd numbers of the form p^(1+4k) * r^2, where p is prime of the form 1+4m, r > 1, and gcd(p,r) = 1. (Euler's criteria for odd perfect numbers).

45, 117, 153, 245, 261, 325, 333, 369, 405, 425, 477, 549, 605, 637, 657, 725, 801, 833, 845, 873, 909, 925, 981, 1017, 1025, 1053, 1233, 1325, 1341, 1377, 1413, 1421, 1445, 1525, 1557, 1573, 1629, 1737, 1773, 1805, 1813, 1825, 2009, 2057, 2061, 2097, 2169

Offset: 1

T. D. Noe, Aug 13 2013

It has been proved that if an odd perfect number exists, it belongs to this sequence. The first term of the form p^5 * n^2 is 28125 = 5^5 * 3^2, occurring in position 520.
Sequence A228059 lists the subsequence of these numbers that are closer to being perfect than smaller numbers. - T. D. Noe, Aug 15 2013
Sequence A326137 lists terms with at least five distinct prime factors. See further comments there. - Antti Karttunen, Jun 13 2019

T. D. Noe, Table of n, a(n) for n = 1..10000
Charles Greathouse and Eric W. Weisstein, MathWorld: Odd perfect number
Oliver Knill, The oldest open problem in mathematics, Handout for NEU Math Circle, December 2, 2007
P. P. Nielsen, Odd Perfect Numbers Have At Least Nine Distinct Prime Factors, arXiv:math/0602485 [math.NT], 2006.
Wikipedia, Perfect number: Odd perfect numbers
Index entries for sequences where any odd perfect numbers must occur

Subsequence of A191218, and also of A228056 and A228057 (simpler versions of this sequence).
For various subsequences with additional conditions, see A228059, A325376, A325380, A325822, A326137 (with omega(n)>=5), A324898 (conjectured, subsequence if it does not contain any prime powers), A354362, A386425 (conjectured), A386427 (nondeficient terms), A386428 (powerful terms), A386429 U A351574.
Cf. A027748, A124010, A005408, A324647, A325319, A325320, A325375, A325377, A325378, A325379, A325819, A325823, A325824.

import Data.List (partition)
a228058 n = a228058_list !! (n-1)
a228058_list = filter f [1, 3 ..] where
   f x = length us == 1 && not (null vs) &&
         fst (head us) `mod` 4 == 1 && snd (head us) `mod` 4 == 1
         where (us,vs) = partition (odd . snd) $
                         zip (a027748_row x) (a124010_row x)
-- Reinhard Zumkeller, Aug 14 2013

nn = 100; n = 1; t = {}; While[Length[t] < nn, n = n + 2; {p, e} = Transpose[FactorInteger[n]]; od = Select[e, OddQ]; If[Length[e] > 1 && Length[od] == 1 && Mod[od[[1]], 4] == 1 && Mod[p[[Position[e, od[[1]]][[1,1]]]], 4] == 1, AppendTo[t, n]]]; t (* T. D. Noe, Aug 15 2013 *)

up_to = 1000;
isA228058(n) = if(!(n%2)||(omega(n)<2),0,my(f=factor(n),y=0); for(i=1,#f~,if(1==(f[i,2]%4), if((1==y)||(1!=(f[i,1]%4)),return(0),y=1), if(f[i,2]%2, return(0)))); (y));
A228058list(up_to) = { my(v=vector(up_to), k=0, n=0); while(kA228058(n), k++; v[k] = n)); (v); };
v228058 = A228058list(up_to);
A228058(n) = v228058[n]; \\ Antti Karttunen, Apr 22 2019

From Antti Karttunen, Apr 22 2019 & Jun 03 2019: (Start)
A325313(a(n)) = -A325319(n).
A325314(a(n)) = -A325320(n).
A001065(a(n)) = A325377(n).
A033879(a(n)) = A325379(n).
A034460(a(n)) = A325823(n).
A325814(a(n)) = A325824(n).
A324213(a(n)) = A325819(n).
(End)

Note in parentheses added to the definition by Antti Karttunen, Jun 03 2019

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 15, 105, 3003, 3465, 13923, 45045, 264537, 459459, 745875, 1541475, 5221125, 8729721, 10790325, 14171625, 29288025, 34563375, 57034575, 71430975, 99201375, 109643625, 144729585, 205016175, 255835125, 295708875, 356080725, 399242025, 419159475, 449323875, 928602675, 939495375, 1083656925, 1941623775, 1962350685, 2083228875

Offset: 1

Antti Karttunen, Nov 10 2021

Numbers k such that A348990(k) [= k/gcd(k, A003961(k))] is equal to A348992(k), which is the odd part of A349162(k), thus all terms must be odd, as A348990 preserves the parity of its argument.
Equally, numbers k for which gcd(A064987(k), A191002(k)) is equal to A000265(gcd(A064987(k), A341529(k))).
Also odd numbers k for which A348993(k) = A319627(k).
Odd terms of A336702 are given by the intersection of this sequence and A349174.
Conjectures:
  (1) After 1, all terms are multiples of 3. (Why?)
  (2) After 1, all terms are in A104210, in other words, for all n > 1, gcd(a(n), A003961(a(n))) > 1. Note that if we encountered a term k with gcd(k, A003961(k)) = 1, then we would have discovered an odd multiperfect number.
  (3) Apart from 1, 15, 105, 3003, 13923, 264537, all other terms are abundant. [These apparently are also the only terms that are not Zumkeller, A083207. Note added Dec 05 2024]
  (4) After 1, all terms are in A248150. (Cf. also A386430).
  (5) After 1, all terms are in A348748.
  (6) Apart from 1, there are no common terms with A349753.
Note: If any of the last four conjectures could be proved, it would refute the existence of odd perfect numbers at once. Note that it seems that gcd(sigma(k), A003961(k)) < k, for all k except these four: 1, 2, 20, 160.
Questions:
  (1) For any term x here, can 2*x be in A349745? (Partial answer: at least x should be in A191218 and should not be a multiple of 3). Would this then imply that x is an odd perfect number? (Which could explain the points (1) and (4) in above, assuming the nonexistence of opn's).

Cf. A000203, A003961, A007949, A064989, A083207, A104210, A161942, A191002, A191218, A228058, A319627, A322361, A326134, A329963, A342671, A348748, A348990, A348992, A348993, A349162, A349164, A349174.
Cf. also A349745, A349746, A349753.
Subsequence of A349749, and of A386430.

Select[Range[10^6], #1/GCD[#1, #3] == #2/(2^IntegerExponent[#2, 2]*GCD[#2, #3]) & @@ {#, DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &] (* Michael De Vlieger, Nov 11 2021 *)

A000265(n) = (n >> valuation(n, 2));
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA349169(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == A000265(s)*gcd(n,u)); }; \\ (Program simplified Nov 30 2021)

For all n >= 1, A007949(A000203(a(n))) = A007949(a(n)). [sigma preserves the 3-adic valuation of the terms of this sequence] - Antti Karttunen, Nov 29 2021

Name changed and comment section rewritten by Antti Karttunen, Nov 29 2021

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 120, 216, 672, 2464, 22176, 228480, 523776, 640640, 837760, 5581440, 5765760, 7539840, 12999168, 19603584, 33860736, 38342304, 71344000, 95472000, 102136320, 197308800, 220093440, 345080736, 459818240, 807009280, 975576960, 1476304896, 1510831360, 1773584640

Offset: 1

Antti Karttunen, Nov 29 2021

nonn

Numbers k for which k * A342671(k) = A000203(k) * A322361(k).
Numbers k such that gcd(A064987(k), A191002(k)) = gcd(A064987(k), A341529(k)).
Obviously, all odd terms in this sequence must be squares.
All the terms k of A005820 that satisfy A007949(k) < A007814(k) [i.e., whose 3-adic valuation is strictly less than their 2-adic valuation] are also terms of this sequence. Incidentally, the first six known terms of A005820 satisfy this condition, while on the other hand, any hypothetical odd 3-perfect number would be excluded from this sequence. Also, as a corollary, any hypothetical 3-perfect numbers of the form 4u+2 must not be multiples of 3 if they are to appear here. Similarly for any k which occurs in A349169, for 2*k to occur in this sequence, it shouldn't be a multiple of 3 and k should also be a term of A191218. See question 2 and its partial answer in A349169.
From Antti Karttunen, Feb 13-20 2022: (Start)
Question: Are all terms/2 (A351548) abundant, from n > 1 onward?
Note that of the 65 known 5-multiperfect numbers, all others except these three 1245087725796543283200, 1940351499647188992000, 4010059765937523916800 are also included in this sequence. The three exceptions are distinguished by the fact that their 3 and 5-adic valuations are equal. In 62 others the former is larger.
If k satisfying the condition were of the form 4u+2, then it should be one of the terms of A191218 doubled as only then both k and sigma(k) are of the form 4u+2, with equal 2-adic valuations for both. More precisely, one of the terms of A351538.
(End)

Amiram Eldar, Table of n, a(n) for n = 1..52 (terms below 10^11)
Antti Karttunen, Ratio A342671(n)/A322361(n) plotted with OEIS Plot2 tool
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000203, A003961, A005101, A005820, A007814, A007949, A046060, A064987, A191002, A191218, A248150, A322361, A341529, A342671, A347870, A351534, A351536, A351538.
Cf. also A349169, A349746, A351458, A351549 for other variants.
Subsequence of A351554 and also of its subsequence A351551.
Cf. A351459 (subsequence, intersection with A351458), A351548 (terms halved).

f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := NextPrime[p]^e; q[1] = True; q[n_] := n * GCD[(s = Times @@ f1 @@@ (f = FactorInteger[n])), (r = Times @@ f2 @@@ f)] == s*GCD[n, r]; Select[Range[10^6], q] (* Amiram Eldar, Nov 29 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA349745(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == (s*gcd(n,u))); };

For all n >= 1, A007814(A000203(a(n))) = A007814(a(n)). [sigma preserves the 2-adic valuation of the terms of this sequence]

A351458 Numbers k for which k * gcd(sigma(k), A276086(k)) is equal to sigma(k) * gcd(k, A276086(k)), where A276086 is the primorial base exp-function, and sigma gives the sum of divisors of its argument.

Original entry on oeis.org

1, 10, 56, 9196, 9504, 56160, 121176, 239096, 354892, 411264, 555520, 716040, 804384, 904704, 1063348, 1387386, 1444352, 1454112, 1884800, 2708640, 3317248, 3548920, 4009824, 4634784, 6179712, 6795360, 7285248, 14511744, 16328466, 28377216, 29855232, 31940280, 37444736, 42711552, 49762944, 52815744

Offset: 1

Antti Karttunen, Feb 13 2022

nonn

Numbers k such that k * A324644(k) = A000203(k) * A324198(k).
Numbers k such that gcd(A064987(k), A324580(k)) = gcd(A064987(k), A351252(k)).
Numbers k such that their abundancy index [sigma(k)/k] is equal to A324644(k)/A324198(k). See A364286.
A324644 gives odd values for even numbers and for the odd squares. A324198 is odd on all arguments, therefore on odd squares the above equation reduces to odd * odd = odd * odd, and on odd nonsquares as odd * even = even * odd. It is an open question whether there are any odd terms after the initial a(1)=1.
If k is even, but not a multiple of 3, then A276086(k) is a multiple of 3, but not even (i.e., is an odd multiple of 3). If for such k also sigma(k) = 3*k, then A007949(A324644(k)) = min(A007949(sigma(k)), A007949(A276086(k))) = 1, while A007949(A324198(k)) = min(A007949(k), A007949(A276086(k))) = 0, therefore all such k's do occur in this sequence, for example, the two known terms of A005820 (3-perfect numbers) that are not multiples of three: 459818240, 51001180160, but also any hypothetical term of A005820 of the form 4u+2, where 2u+1 is not multiple of 3, and which by necessity is then also an odd perfect number.
Similarly, of the 65 known 5-multiperfect numbers (A046060), those 20 that are not multiples of five are included in this sequence. Note that all 65 are multiples of six.
It is conjectured that the intersection of this sequence with the multiperfect numbers (A007691) gives A323653, see comments in the latter.
For all even terms k of this sequence, A007814(A000203(k)) = A007814(k), sigma preserves the 2-adic valuation, and A007949(A000203(k)) >= A007949(k), i.e., does not decrease the 3-adic valuation. The condition is equivalence (=) when k is a multiple of 6. With odd terms, any hypothetical odd perfect number x would yield a one greater 2-adic valuation for sigma(x) than for x, but would satisfy the main condition of this sequence. - Corrected Feb 17 2022
If k is a nonsquare positive odd number (in A088828), then it must be a term of A191218. - Antti Karttunen, Mar 10 2024

Cf. A000203, A005820, A007691, A007814, A007949, A046060, A088828, A191218, A276086, A323653, A324198, A324644, A327938, A364286, A351459 (intersection with a similar sequence, A349745).

Cf. also A351549.

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
isA351458(n) = { my(s=sigma(n), z=A276086(n)); (n*gcd(s,z))==(s*gcd(n,z)); };

A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]); \\ Works OK with rationals also!
isA351458(n) = { my(orgn=n, s=sigma(n), abi=s/n, p=2, q=A006530(abi), d, e1, e2); while((1!=abi)&&(p<=q), d = n%p; e1 = min(d, valuation(s, p)); e2 = min(d, valuation(orgn, p)); d = e1-e2; if(valuation(abi,p)!=d, return(0), abi /= (p^d)); n = n\p; p = nextprime(1+p)); (abi==1); }; \\ (This implementation does not require the construction of largish intermediate numbers, A276086, but might still be slower and return a few false positives on the long run, so please check the results with the above program). - Antti Karttunen, Feb 19 2022

A324647 Odd numbers k such that 2k is equal to bitwise-AND of 2k and sigma(k).

Original entry on oeis.org

1116225, 1245825, 1380825, 2127825, 10046025, 16813125, 203753025, 252880425, 408553425, 415433025, 740361825, 969523425, 1369580625, 1612924425, 1763305425, 2018027025, 2048985225, 2286684225, 3341556225, 3915517725, 3985769025, 4051698525, 7085469825, 7520472225

Offset: 1

Antti Karttunen, Mar 14 2019

nonn

If this sequence has no terms common with A324649 (A324897, A324898), or no terms common with A324727, then there are no odd perfect numbers.
First 22 terms factored:
     1116225 = 3^2 * 5^2 * 11^2 * 41
     1245825 = 3^2 * 5^2 *  7^2 * 113
     1380825 = 3^2 * 5^2 * 19^2 * 17    [Here the unitary prime is not the largest]
     2127825 = 3^2 * 5^2 *  7^2 * 193
    10046025 = 3^4 * 5^2 * 11^2 * 41
    16813125 = 3^2 * 5^4 *  7^2 * 61
   203753025 = 3^2 * 5^2 *  7^2 * 18481
   252880425 = 3^2 * 5^2 *  7^2 * 22937
   408553425 = 3^2 * 5^2 *  7^2 * 37057
   415433025 = 3^2 * 5^2 *  7^4 * 769
   740361825 = 3^2 * 5^2 *  7^2 * 67153
   969523425 = 3^4 * 5^2 * 13^2 * 2833
  1369580625 = 3^2 * 5^4 *  7^2 * 4969
  1612924425 = 3^2 * 5^2 *  7^2 * 146297
  1763305425 = 3^2 * 5^2 *  7^2 * 159937
  2018027025 = 3^2 * 5^2 *  7^2 * 183041
  2048985225 = 3^2 * 5^2 *  7^2 * 185849
  2286684225 = 3^2 * 5^2 *  7^2 * 207409
  3341556225 = 3^2 * 5^2 *  7^2 * 303089
  3915517725 = 3^4 * 5^2 *  7^2 * 39461
  3985769025 = 3^4 * 5^2 *  7^2 * 40169
  4051698525 = 3^2 * 5^2 *  7^2 * 367501.
Compare the above factorizations to the various constraints listed for odd perfect numbers in the Wikipedia article. However, this is NOT a subsequence of A191218 (A228058), see below.
The first terms that do not belong to A191218 are 399736269009 = (3 * 7^2 * 11 * 17 * 23)^2 and 1013616036225 = (3^2 * 5 * 13 * 1721)^2, that both occur instead in A325311. The first terms with omega(n) <> 4 are 9315603297, 60452246925, 68923392525, and 112206463425. They factor as 3^2 * 7^2 * 11^2 * 13^2 * 1033, 3^2 * 5^2 * 7^2 * 17^2 * 18973, 3^2 * 5^2 * 13^2 * 19^2 * 5021, 3^2 * 5^2 * 7^2 * 199^2 * 257. - Giovanni Resta, Apr 21 2019
From Antti Karttunen, Jan 13 2025: (Start)
Because of the "monotonic property" of bitwise-and, this is a subsequence of nondeficient numbers (A023196).
Both odd perfect numbers, and quasiperfect numbers, if such numbers exist at all, would satisfy the condition for being included in this sequence. Furthermore, any term must be either an odd square with an odd abundancy (in A156942), which subset is given in A379490 (where quasiperfect numbers must thus reside, if they exist), or be included in A228058, i.e., satisfy the Euler's criteria for odd perfect numbers.
(End)

Giovanni Resta, Table of n, a(n) for n = 1..500
Charles Greathouse and Eric W. Weisstein, MathWorld: Odd perfect number
Wikipedia, Perfect number: Odd perfect numbers
Index entries for sequences where any odd perfect numbers must occur

Odd terms of A324652.

Cf. A191218, A228058, A318468, A324649, A324659, A324718, A324719, A324722, A324727, A324880, A324897, A324898, A325311, A379490 (square terms).

for(n=1,oo,if((n%2)&&((2*n)==bitand(2*n,sigma(n))),print1(n,", ")));

{Odd k such that 2k = A318468(k)}.

a(23)-a(24) from Giovanni Resta, Apr 21 2019

A191217 Numbers n such that sigma(n) is congruent to 2 modulo 4.

Original entry on oeis.org

5, 10, 13, 17, 20, 26, 29, 34, 37, 40, 41, 45, 52, 53, 58, 61, 68, 73, 74, 80, 82, 89, 90, 97, 101, 104, 106, 109, 113, 116, 117, 122, 136, 137, 146, 148, 149, 153, 157, 160, 164, 173, 178, 180, 181, 193, 194, 197, 202, 208, 212, 218, 226, 229, 232, 233, 234, 241, 244, 245, 257, 261, 269

Offset: 1

Luis H. Gallardo, May 26 2011

nonn

These numbers are exactly the numbers of the form 2^a * p^(4b+1) * m^2 where p is a prime number congruent to 1 modulo 4, a is a nonnegative integer, and m is a positive integer coprime to p. In particular, they are also sums of two squares: the sequence has the first 12 terms in common with A132777.

I corrected the above comment by adding the exponent (4b+1) to p, because otherwise it would miss terms like a(614) = 3125 = 5^5, a(1140) = 6250 = 2 * 5^5, a(4421) = 28125 = 5^5 * 3^2, etc. - Antti Karttunen, May 25 2022

			For n=2, a(2) = 10 since sigma(10) = 18 = 4*4 + 2 is congruent to 2 modulo 4

Ray Chandler, Table of n, a(n) for n = 1..10000

Similar to, but different from, A230779, which is a subsequence.
Cf. A191218, A228058, A332226 for other subsequences.
Cf. A353812 (characteristic function).

with(numtheory): gen := proc(b) local n,s,d; for n from 1 to b do s := sigma(n);
if modp(s,4)=2 then print(n); fi; od; end;

for(n=1,10^3,if(2==(sigma(n)%4),print1(n,", "))) /* Joerg Arndt, May 27 2011 */

A324652 Numbers k such that A318468(k) (bitwise-AND of 2k and sigma(k)) is equal to 2k.

Original entry on oeis.org

6, 12, 18, 20, 24, 28, 36, 40, 48, 56, 80, 88, 96, 100, 104, 112, 160, 176, 192, 196, 200, 204, 208, 220, 224, 260, 264, 272, 304, 320, 336, 352, 368, 384, 392, 416, 448, 464, 496, 544, 550, 580, 608, 640, 648, 650, 672, 704, 736, 768, 784, 832, 896, 928, 992, 1032, 1040, 1044, 1056, 1060, 1068, 1088, 1104, 1120, 1184, 1216

Offset: 1

Antti Karttunen, Mar 14 2019

nonn

Positions of zeros in A324658, fixed points of A324659.
Intersection with A324649 gives A324643.
Intersection with A324726 gives A000396.
In the range 1..2^32 there are only 22 odd terms. See A324647.

Cf. A000203, A191218, A318468, A324647, A324649, A324658, A324659, A324722, A324726.

Some subsequences: A000396, A324643, A324647 (the odd terms).

Select[Range[2000], 2*# == BitAnd[2*#, DivisorSigma[1, #]] &] (* Paolo Xausa, Mar 11 2024 *)

for(n=1,oo,if((n+n)==bitand(2*n,sigma(n)), print1(n, ", ")))

A348739 Numbers k for which A326042(k) > k.

Original entry on oeis.org

4, 9, 12, 16, 18, 25, 32, 36, 44, 48, 49, 64, 72, 81, 96, 99, 100, 108, 124, 144, 147, 162, 169, 176, 180, 192, 196, 225, 236, 243, 252, 256, 279, 284, 288, 300, 320, 324, 361, 372, 396, 400, 405, 432, 441, 448, 450, 468, 484, 486, 496, 507, 512, 529, 531, 567, 576, 588, 604, 612, 625, 639, 648, 675, 676, 700, 704

Offset: 1

Antti Karttunen, Nov 02 2021

nonn

Terms that occur also in A337386 are: 180, 300, 720, 900, 960, 1008, 1200, 1440, 1620, 1800, 2016, 2400, ...

Positions of negative terms in A348736.
Cf. A326182 (subsequence after its initial 1), A348738.
Cf. A000203, A003961, A064989, A161942, A191218, A326042, A337386, A348742, A348749 (corresponding odd numbers), A348942.

f1[2, e_] := 1; f1[p_, e_] := NextPrime[p, -1]^e; s1[1] = 1; s1[n_] := Times @@ f1 @@@ FactorInteger[n]; f2[p_, e_] := NextPrime[p]^e; s2[1] = 1; s2[n_] := Times @@ f2 @@@ FactorInteger[n]; Select[Range[700], s1[DivisorSigma[1, s2[#]]] > # &] (* Amiram Eldar, Nov 04 2021 *)

A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A326042(n) = A064989(sigma(A003961(n)));
isA348739(n) = (A326042(n)>n);

A348743 Odd nonsquares k for which A161942(k) >= k, where A161942 is the odd part of sigma.

Original entry on oeis.org

2205, 19845, 108045, 143325, 178605, 187425, 236925, 266805, 319725, 353925, 372645, 407925, 452025, 462825, 584325, 637245, 646425, 658125, 672525, 789525, 796005, 804825, 845325, 920205, 972405, 981225, 1007325, 1055925, 1069425, 1102725, 1113525, 1116225, 1166445, 1201725, 1245825, 1289925, 1378125, 1380825, 1442925

Offset: 1

Antti Karttunen, Nov 02 2021

The first non-multiples of 5 are a(103) = 6243237 and a(125) = 8164233.
From Antti Karttunen, Nov 28 2024: (Start)
This is not a subsequence of A228058. At least k = A000040(28)*(A002110(27)/2)^2 = 15388519572341080054329140040512468358441210638435506649120749687401476705908239675 is a number of the form 4m+3 such that A161942(k) >= k.
Another such number is A000040(28)*81*(A002110(25)/6)^2 = 1279741205456530915782536871495922949062895982530933679752838870798129159675.
Question: What is the smallest term of this sequence that is of the form 4m+3, and thus not in A386427 (in A191218 and in A228058)?
(End)

Intersection of A088828 and A348742.
Cf. A386427 (a subsequence, which agrees for a very long time).
Cf. A033630, A065205, A083206, A161942, A191218, A228058, A336702, A348741.
Cf. also A065235, A162284.

A000265(n) = (n >> valuation(n, 2));
isA348743(n) = ((n%2)&&!issquare(n)&&A000265(sigma(n))>=n); \\ Edited Nov 28 2024

Definition changed (from > to >=) to formally include also any hypothetical odd perfect numbers - Antti Karttunen, Nov 28 2024

Comment removed, because it was more related to sequence A386427. - Antti Karttunen, Aug 21 2025

A212037 The size of the set of numbers k>=0 such that all (2^n+k)2^n-1 are prime but only the last (largest) (2^n+k)2^n+1 is also an associated twin prime.

Original entry on oeis.org

1, 6, 1, 4, 2, 2, 2, 24, 6, 2, 28, 7, 16, 47, 29, 6, 41, 16, 3, 17, 32, 10, 10, 23, 14, 15, 52, 4, 13, 20, 23, 4, 84, 26, 88, 50, 20, 35, 51, 44, 41, 87, 1, 142, 13, 188, 107, 162, 91, 96, 197, 4, 148, 71, 9, 66, 97, 41, 10, 9, 152, 234, 48, 104, 144, 40, 18, 45, 52, 204, 21, 49, 51, 9, 102, 13, 31, 108, 88

Offset: 1

Pierre CAMI, Jul 14 2012

nonn

Starting at a count of zero, we consider for increasing k>=0 the pairs (2^n+k)*2^n+-1. If the smaller of these two numbers is prime, we increase the counter. If the larger of these two numbers is also prime, we admit the counter to the sequence. It is basically a measure of how many unsuccessful primality tests on the larger of the two numbers are done before it becomes a compatible twin prime.

Heuristically, the average of a(n)/n over n=1 to N tends to 1 as N increases.

			For n=2, the 6 pairs (19,21) at k=1, (23,25) at k=2, (31,33) at k=4, (43,45) at k=7, (47,49) at k=8 and (59,61) at k=11 are counted. The smaller of these must be a prime to be counted, and at k=11 also the larger (i.e., 61) becomes prime, which finishes the search.

Pierre CAMI, Table of n, a(n) for n = 1..825

Cf. A191617, A191218, A205321.

A212037 := proc(n)
    local a,k,p ;
    a := 0 ;
    for k from 0 do
        p := (2^n+k)*2^n-1 ;
        if isprime(p) then
            a := a+1 ;
        end if;
        if isprime(p) and isprime(p+2) then
            return a;
        end if;
    end do:
end proc: # R. J. Mathar, Jul 20 2012

cp's OEIS Frontend

A228058 Odd numbers of the form p^(1+4k) * r^2, where p is prime of the form 1+4m, r > 1, and gcd(p,r) = 1. (Euler's criteria for odd perfect numbers).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A351458 Numbers k for which k * gcd(sigma(k), A276086(k)) is equal to sigma(k) * gcd(k, A276086(k)), where A276086 is the primorial base exp-function, and sigma gives the sum of divisors of its argument.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

A324647 Odd numbers k such that 2*k is equal to bitwise-AND of 2*k and sigma(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

Extensions

A191217 Numbers n such that sigma(n) is congruent to 2 modulo 4.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

PARI

A324652 Numbers k such that A318468(k) (bitwise-AND of 2*k and sigma(k)) is equal to 2*k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A348739 Numbers k for which A326042(k) > k.

A324647 Odd numbers k such that 2k is equal to bitwise-AND of 2k and sigma(k).

A324652 Numbers k such that A318468(k) (bitwise-AND of 2k and sigma(k)) is equal to 2k.

A212037 The size of the set of numbers k>=0 such that all (2^n+k)2^n-1 are prime but only the last (largest) (2^n+k)2^n+1 is also an associated twin prime.