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 A352616 #8 Mar 29 2022 19:57:38 %S A352616 7,17,50,71,449,851,881,1467,2591,5123,9521,18167,29817,34859,38027, %T A352616 39761,52415,84359,90611,92711,106721,147267,179999,204479,206081, %U A352616 342791,388961,596231,606609,847601,922607,1060307,1256879,1292831,1389671,1590983,1976015,2157227,2169971,2268449,2571911 %N A352616 Numbers k such that A155085(k) is in A037074. %C A352616 Numbers k such that k + A000203(k) is the product of a twin prime pair. %e A352616 a(3) = 50 is a term because 50 + A000203(50) = 143 = 11*13 where (11,13) is a twin prime pair. %p A352616 filter:= proc(n) local t,s; t:= n + numtheory:-sigma(n); %p A352616 if t::even or not issqr(t+1) then return false fi; %p A352616 s:= sqrt(t+1); %p A352616 isprime(s+1) and isprime(s-1) %p A352616 end proc: %p A352616 select(filter, [$1..3000000]); %Y A352616 Cf. A000203, A037074, A155085. %K A352616 nonn %O A352616 1,1 %A A352616 _J. M. Bergot_ and _Robert Israel_, Mar 23 2022