|
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
All terms are odd. Of the 63 initial terms of A349169, only term 13923 occurs also in this sequence. The first common term with A332458 is 161257. - Antti Karttunen, Mar 10 2024
|
|
LINKS
|
|
|
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 *)
|
|
PROG
|
(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
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|