[seqfan] Re: A005850 Primes p such that the NSW number A002315(p) is prime

[Charles Greathouse]

I concur -- the values I get from the definition are

2, 3, 23, 29, ...

using

forprime(p=2, 1e4, w=3+quadgen(32); if(ispseudoprime(imag((1+w)*w^p)),
print1(p", ")))

Charles Greathouse
Analyst/Programmer
Case Western Reserve University

On Fri, Oct 19, 2012 at 7:02 AM, Jean-François Alcover
<jf.alcover at gmail.com> wrote:
> Dear SeqFans,
>
> Using the given definition I can't correctly recompute
> A005850<https://oeis.org/A005850>
>
> Could someone please give me an example ?
> Thanks in advance
>
> jfa
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/