A258882
Primitive weird numbers of the form 2^k*p*q with k > 0 and where p < q are odd primes.
Original entry on oeis.org
70, 836, 7192, 7912, 9272, 10792, 17272, 73616, 83312, 113072, 519712, 539744, 555616, 682592, 786208, 1188256, 1229152, 1901728, 2081824, 2189024, 3963968, 4128448, 4145216, 4486208, 4559552, 4632896, 4960448, 5440192, 5568448, 6460864, 6621632, 7354304, 7470272, 8000704, 8134208
Offset: 1
a(1) = A002975(1) = 70 = 2*5*7.
a(2) = A002975(2) = 836 = 2^2*11*19.
A002975(3) = 4030 = 2*5*13*31 is not in this sequence since it is not of the required form.
The same is true for A002975(4) = 5830.
a(3) = A002975(5) = 7192 = 2^3*29*31, etc.
A002975(179) = 2319548096 = 2^6 * 137^2 * 1931 is the first term of A002975 with only two odd prime divisors, but not of the required form. - _M. F. Hasler_, Nov 20 2018
- S. Kravitz, A search for large weird numbers. J. Recreational Math. 9 (1976), 82-85 (1977). Zbl 0365.10003
- Douglas E. Iannucci and Robert G. Wilson v, Table of n, a(n) for n = 1..15384, updated Dec 06 2015; corrected by _M. F. Hasler_, Jul 16 2016
- R. Bagula et al., A very big weird number, Number Theory group on LinkedIn (web.archive.org snapshot; page no longer available). Dec. 2013
- Central Washington University, CWU Math Students Break World Record for Largest Weird Number [alternate article]
- Douglas E. Iannucci, On primitive weird numbers of the form 2^k*p*q, arXiv:1504.02761 [math.NT], 2015.
- Giuseppe Melfi, On the conditional infiniteness of primitive weird numbers, Journal of Number Theory, Vol. 147, Feb 2015, pp 508-514.
- Eric Weisstein's World of Mathematics, Weird Number.
- Wikipedia, Weird number
Cf.
A320875 (more general case of Karavitz' formula).
-
(* copy the terms from A002975, assign them equal to 'lst' and then *) fQ[n_] := Block[{m = n}, While[ Mod[m, 2] == 0, m /= 2]; PrimeOmega@ m == 2]; Select[lst, fQ]
-
select(t->factor(t)[, 2][^1]=[1, 1]~, A002975) \\ Assuming that A002975 is defined as set or vector. - M. F. Hasler, Jul 11 2016
A242998
Number of integers k such that R = (2^k*Q - Q - 1)/(Q + 1 - 2^k) is a prime number, when Q = A000668(n) is the n-th Mersenne prime.
Original entry on oeis.org
0, 1, 1, 2, 1, 1, 2, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1
For given p=A000043(n), the following k's yield a prime R:
p : k's (and resulting primes R, Q=2^p-1 and/or weird W=2^(k-1)*Q*R)
2 : -
3 : 2 (R=5, Q=7, W=70)
5 : 4 (R=29, Q=31, W=7192)
7 : 4 (R=17, Q=127, W=17272), 5 (R=41, Q=127, W=83312)
13 : 11 (R=2729, Q=8191, W=22889716736)
17 : 13 (R=8737, Q=131071, W=4690605371392)
19 : 16 (R=74897, W=1286718208049152), 17 (R=174761, W=6004730783793152)
31 : 16 (R=65537, W=2^15*(2^31-1)*R), 29 (R=715827881, W=2^28*(2^31-1)*R)
61 : 57 (R=153722867280912929, W=2^56*(2^61-1)*R)
89 : 78 (R=302379100949042568368129, W=2^77*(2^89-1)*R)
107 through 86243 : none.
107 through 3021377: none. _Robert Price_, Sep 05 2019
The present sequence lists the number of k's in each line.
- S. Kravitz, A search for large weird numbers. J. Recreational Math. 9(1976), 82-85 (1977). Zbl 0365.10003
- E. Weisstein, Weird numbers, on MathWorld - a Wolfram web ressource.
See also
A320875 for more general solutions to R = (MQ-1)/(Q-M) = prime.
-
A000043 = {2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607,
1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937,
21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433,
1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011,
24036583, 25964951, 30402457, 32582657, 37156667, 42643801,
43112609};
lst = {};
For[i = 1, i <= 28, i++,
p = A000043[[i]];
kc = 0;
For[k = 1, k < p, k++,
r = 2^k - 1 + (2^k - 2)/(2^(p - k) - 1);
If[! IntegerQ[r], Continue[]];
If[PrimeQ[r], kc++]];
AppendTo[lst, kc]];
lst (* Robert Price, Sep 05 2019 *)
-
A242998(n,p=A000043[n])={sum(k=p\2+1, p-1, Mod(2, 2^(p-k)-1)^k==2 && ispseudoprime(2^k-1+(2^k-2)/(2^(p-k)-1)))}
A242025
Primes of the form R = 2^k-1+(2^k-2)/(2^(p-k)-1), where p are Mersenne prime exponents listed in A000043.
Original entry on oeis.org
5, 17, 29, 41, 2729, 8737, 65537, 74897, 174761, 715827881, 153722867280912929, 302379100949042568368129
Offset: 1
For given p = A000043(n), the following k yield a prime R and an associated (primitive) weird number W = 2^(k-1)*(2^p-1)*R in A258882 c A002975 c A006037:
For p = 2, no k yields a prime R = 2^k-1+(2^k-2)/(2^(p-k)-1).
For p = 3, k = 2 yields R = 5 and the (smallest) weird number W = 70 = A006037(1).
For p = 5, k = 4 yields R = 29 = a(3) and W = 7192 = A258882(3).
For p = 7, k = 4 yields R = 17 = a(2) and W = 17272 = A258882(7),
and k = 5 yields R = 41 = a(4) and W = 83312 = A258882(9).
For p = 13, k = 11 yields R = 2729 = a(5) and W = 22889716736 = A258882(288)
For p = 17, k = 13 yields R = 8737 = a(6) and W = 4690605371392 = A258882(1203).
For p = 19, k = 16 yields R = 74897 = a(8), W = 1286718208049152 = A258882(7154),
and k = 17 yields R = 174761 = a(9), W = 6004730783793152 = A258882(11466).
For p = 31, k = 16 yields R = 65537 = a(7) (smaller than both R's for p = 19),
and k = 29 yields R = 715827881 = a(10).
For p = 61, only k = 57 yields a prime R = 153722867280912929 = a(11).
For p = 89, only k = 78 yields a prime R = 302379100949042568368129 = a(12).
For p = 107 through p = 86243, no k yields a prime R.
For p = 107 through p = 3021377, no k yields a prime R. - _Robert Price_, Sep 04 2019
Cf.
A258882 (PWN of the form 2^k*p*q).
Cf.
A320875 for a more general pattern leading to more and larger PWN.
-
A000043 = {2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607,
1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937,
21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433,
1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011,
24036583, 25964951, 30402457, 32582657, 37156667, 42643801,
43112609};
lst = {};
For[i = 1, i <= Length[A000043], i++,
p = A000043[[i]];
For[k = 1, k < p, k++,
r = 2^k - 1 + (2^k - 2)/(2^(p - k) - 1);
If[! IntegerQ[r], Continue[]];
If[PrimeQ[r], AppendTo[lst, r]]]];
Union[lst] (* Robert Price, Sep 04 2019 *)
Showing 1-3 of 3 results.
Comments