|
REFERENCES
|
Prime Obsession, John Derbyshire, Joseph Henry Press, April 2004, p 16.
|
|
PROG
|
(PARI) primenum(n, k, typ) = \\ k=mult, typ=1 num, 2 denom. output prime num or denom.
{ local(a, b, x, tmp, v);
a=1; b=1; for(x=1, n, tmp=b; b=a+b; a=k*tmp+a; if(typ==1, v=a, v=b); if(isprime(v), print1(v", "); ) );
print(); print(a/b+.) }
|