A128283 Numbers of the form m = p1 * p2 where for each d|m we have (d+m/d)/2 prime and p1 < p2 both prime.
21, 33, 57, 85, 93, 133, 145, 177, 205, 213, 217, 253, 393, 445, 553, 565, 633, 697, 793, 817, 865, 913, 933, 973, 1137, 1285, 1345, 1417, 1437, 1465, 1477, 1513, 1537, 1717, 1765, 1837, 1857, 1893, 2101, 2173, 2245, 2305, 2517, 2577, 2581, 2605, 2641, 2653, 2733, 2761
Offset: 1
Keywords
Examples
85=5 * 17, (5 * 17+1)/2=43, (5+17)/2=11 are both primes and 85 is in the sequence. From _Hartmut F. W. Hoft_, Jan 09 2021: (Start) 9=3*3 is not in the sequence even though (1+9)/2 and (3+3)/2 are primes, see also A340482. a(33) = 1537 = 29*53 is the first number for which the symmetric representation of sigma consists of three regions ( 769, 82, 769 ) with 5 units of width 2 straddling the diagonal in the central region; (1537+1)/2 = 769 and (29+53)/2 = 41 are primes. (End)
Links
- David A. Corneth, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
ppQ[s_, k_] := Last[Transpose[FactorInteger[s]]]==Table[1, k] dQ[s_] := Module[{d=Divisors[s]}, AllTrue[Map[(d[[#]]+d[[-#]])/2&, Range[Length[d]/2]], PrimeQ]] goodL[{m_, n_}, k_] := Module[{i=m, list={}}, While[i<=n, If[ppQ[i, k] && dQ[i], AppendTo[list, i]]; i+=2]; list]/;OddQ[m] a128283[n_] := goodL[{1, n}, 2] a128283[2761] (* Hartmut F. W. Hoft, Jan 09 2021 *)
Extensions
Added "distinct" for clarification since 9 satisfies the divisor property. See also A340482. - Hartmut F. W. Hoft, Jan 09 2021
New name from David A. Corneth, Jan 09 2021
Comments