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 A105934 #10 Sep 08 2022 08:45:17 %S A105934 116,176,184,300,444,470,584,690,696,950 %N A105934 Positive integers n such that n^22 + 1 is semiprime (A001358). %C A105934 We have the polynomial factorization: n^22 + 1 = (n^2 + 1) * (n^20 - n^18 + n^16 - n^14 + n^12 - n^10 + n^8 - n^6 + n^4 - n^2 + 1). Hence after the initial n=1 prime, the binomial can never be prime. It can be semiprime iff n^2+1 is prime and (n^20 - n^18 + n^16 - n^14 + n^12 - n^10 + n^8 - n^6 + n^4 - n^2 + 1) is prime. %F A105934 a(n)^22 + 1 is in A001358. a(n)^2+1 is in A000040 and (a(n)^20 - a(n)^18 + a(n)^16 - a(n)^14 + a(n)^12 - a(n)^10 + a(n)^8 - a(n)^6 + a(n)^4 - a(n)^2 + 1) is in A000040. %e A105934 116^22 + 1 = 2618639792014920380336685706161496723088736257 = 13457 * 194593133091693570657404005808240820620401, %e A105934 300^22 + 1 = 3138105960900000000000000000000000000000000000000000001 = 90001 * 34867456593815624270841435095165609271008099910001, %e A105934 950^22 + 1 = 323533544973709366507562922501564025878906250000000000000000000001 = 902501 * 358485525194663902319845543109164450653136395416736380347501. %t A105934 Select[Range[1000], PrimeOmega[#^22 + 1]==2&] (* _Vincenzo Librandi_, May 24 2014 *) %o A105934 (Magma)IsSemiprime:=func< n | &+[ k[2]: k in Factorization(n) ] eq 2 >; [n: n in [2..1000] | IsSemiprime(n^22+1)]; // _Vincenzo Librandi_, Dec 21 2010 %Y A105934 Cf. A000040, A001358, A006313, A103854, A104238, A104335, A105041, A105066, A105078, A105122, A105142, A105237, A104479, A104494, A104657, A105282. %K A105934 easy,nonn %O A105934 1,1 %A A105934 _Jonathan Vos Post_, Apr 26 2005