|
%I #17 Dec 13 2022 08:24:50
%S 1,3,7,11,19,25,29,37,43,51,57,65,75,79,83,99,105,115,125,133,139,149,
%T 153,175,179,187,197,207,213,221,231,247,251,269,279,287,299,307,319,
%U 329,339,349,357,369,373,387,397,407,423,441,455,465,471,485,491,501
%N a(n) = prime(2*n) - 2*n.
%F a(n) = A031215(n) - A005843(n). - _Michel Marcus_, Feb 01 2019 [Corrected by _Georg Fischer_, Dec 13 2022]
%F a(n) = A014689(A005843(n)). - _Michel Marcus_, Dec 13 2022
%e a(1) = prime(2*1) - 1*2 = prime(2) - 2 = 3 - 2 = 1.
%e a(4) = prime(2*4) - 4*2 = prime(8) - 8 = 19 - 8 = 11.
%p p:= ithprime: a:= seq(p(2*n)-2*n,n=1..100); # _D. S. McNeil_, Mar 21 2009
%o (PARI) a(n) = prime(2*n) - 2*n; \\ _Michel Marcus_, Dec 13 2022
%Y Cf. A005843, A031215, A014689.
%K nonn
%O 1,2
%A _Juri-Stepan Gerasimov_, Jul 31 2008
%E a(32), a(40) corrected and extended by _D. S. McNeil_, Mar 21 2009
|
