Review: A066889

Klaus Brockhaus klaus-brockhaus at t-online.de
Sat Jun 1 19:27:53 CEST 2002



David Wilson wrote:
> 
> Change name to shorter:
> 
> %N A066889 nextprime(m)-m, where m = PROD(k = 1 to n, fib(n)).

Rather

%N A066889 nextprime(m+2)-m, where m = PROD(k = 1 to n, fib(k)).
%O A066889 1,1

PARI gives

? for(n=1,56,m=prod(k=1,n,fibonacci(k));print1(nextprime(m+2)-m,","))
2,2,3,5,7,11,17,17,37,23,47,37,29,19,47,59,19,37,71,59,31,67,239,101,739,409,43,
367,167,251,73,71,419,1567,107,83,223,191,227,449,97,173,103,523,79,137,223,1163
,661,103,103,541,227,2383,433,71,

which is A066889, whereas

? for(n=1,56,m=prod(k=1,n,fibonacci(k));print1(nextprime(m)-m,","))
1,1,0,1,1,1,1,1,37,23,47,37,29,19,47,59,19,37,71,59,31,1,239,101,739,409,43,1,16
7,251,73,71,419,1567,107,83,223,191,227,449,97,173,103,523,79,137,223,1163,661,1
03,103,541,227,2383,433,71,

differs from A066889 in a(1),...,a(8),a(22),a(28),....





More information about the SeqFan mailing list