This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A356794 #10 Aug 29 2022 10:30:05 %S A356794 15,35,39,45,51,55,75,87,91,95,105,111,115,117,119,123,135,143,153, %T A356794 155,159,165,175,183,187,195,203,215,219,225,235,245,247,255,259,261, %U A356794 267,273,275,285,287,291,295,299,303,315,319,323,327,333,335,339,345,351 %N A356794 Odd numbers that have at least one prime factor congruent to 1 (mod 4) and at least one prime factor congruent to 3 (mod 4). %C A356794 This sequence consists of the odd positive integers that are in neither A004613 nor A004614. %e A356794 15 is a term: its prime factors are 3 and 5; 3 == 3 (mod 4), 5 == 1 (mod 4). %t A356794 Select[Range[1, 351, 2], Length[Union[Mod[FactorInteger[#][[;; , 1]], 4]]] == 2 &] (* _Amiram Eldar_, Aug 28 2022 *) %o A356794 (PARI) isok(m) = if (m % 2, my(vf=factor(m)[,1]); (#select(x->((x%4)==1), vf) >= 1) && (#select(x->((x%4)==3), vf) >= 1)); \\ _Michel Marcus_, Aug 29 2022 %Y A356794 Cf. A004613, A004614. %K A356794 nonn %O A356794 1,1 %A A356794 _Jon E. Schoenfield_, Aug 27 2022