A364286 Composite numbers k for which A324644(k)/A324198(k) = 2.
33, 51, 69, 91, 99, 135, 141, 145, 153, 159, 187, 207, 213, 217, 285, 295, 303, 321, 339, 391, 411, 423, 427, 435, 445, 477, 507, 519, 573, 637, 639, 679, 681, 699, 771, 783, 799, 843, 855, 861, 885, 895, 901, 909, 933, 951, 963, 1017, 1041, 1057, 1059, 1071, 1081, 1083, 1147, 1149, 1173, 1185, 1195, 1203, 1207
Offset: 1
Keywords
Links
Crossrefs
Programs
-
Mathematica
f[x_] := Block[{m, i, n = x, p}, m = i = 1; While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m]; Select[Select[Range[1350], CompositeQ], GCD[#2, #3]/GCD[#1, #3] == 2 & @@ {#, DivisorSigma[1, #], f[#]} &] (* Michael De Vlieger, Mar 10 2024 *)
-
PARI
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; isA364286(n) = if(isprime(n), 0, my(u=A276086(n)); (gcd(sigma(n),u)==2*gcd(n,u))); \\ Antti Karttunen, Mar 10 2024
Comments