cp's OEIS Frontend

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.

A062459 Primes of form n^2 + mu(n), where mu is A008683.

This page as a plain text file.
%I A062459 #23 Dec 15 2017 17:35:10
%S A062459 2,3,37,101,197,677,5477,8837,17957,21317,42437,44101,98597,106277,
%T A062459 148997,217157,324901,401957,454277,476101,509797,828101,864901,
%U A062459 933157,1196837,1378277,1664101,1674437,1705637,1833317,1865957,1988101
%N A062459 Primes of form n^2 + mu(n), where mu is A008683.
%C A062459 All terms except for 3 are of the form n^2 + 1 where mu(n) = 1, since n^2 or n^2-1 can't be prime except for n=2. - _Robert Israel_, Oct 09 2015
%H A062459 Harry J. Smith, <a href="/A062459/b062459.txt">Table of n, a(n) for n = 1..1000</a>
%p A062459 R:= select(t -> numtheory:-mobius(t)=1, [$3..10000]):
%p A062459 2,3,op(select(isprime, map(t-> t^2+1, R))); # _Robert Israel_, Oct 09 2015
%t A062459 Select[Table[m^2+MoebiusMu[m],{m,0,1500}],PrimeQ]  (* _Harvey P. Dale_, Feb 08 2011 *)
%o A062459 (PARI) j=[]; for(n=1,3000,x=n^2+moebius(n); if(isprime(x),j=concat(j,x))); j
%o A062459 (PARI) n=0; for (m=1, 10^5, x=m^2 + moebius(m); if (isprime(x), write("b062459.txt", n++, " ", x); if (n==1000, break))) \\ _Harry J. Smith_, Aug 08 2009
%o A062459 (PARI) list(lim)=my(v=List([2]),t); forstep(n=2,sqrtint(lim\1),4, if(isprime(t=n^2+moebius(n)), listput(v,t))); Vec(v) \\ _Charles R Greathouse IV_, Sep 22 2015
%Y A062459 Cf. A008683, A002496, A030229.
%K A062459 nonn
%O A062459 1,1
%A A062459 _Jason Earls_, Jul 26 2001