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.
is(n)=ispseudoprime((48^n+47^n)/95) \\ Charles R Greathouse IV, Jun 13 2017
is(n)=ispseudoprime((138^n+137^n)/275) \\ Charles R Greathouse IV, Jun 13 2017
is(n)=ispseudoprime((140^n+139^n)/279) \\ Charles R Greathouse IV, Jun 13 2017
a(10) = 4 because (5^29 + 4^29)/9 = 2149818248341 is prime and (2^29 + 1^29)/3, (3^29 + 2^29)/5 and (4^29 + 3^29)/7 are all composite.
Table[p = Prime[n]; k = 1; While[q = ((b+1)^n+b^n)/(2*b+1); ! PrimeQ[q], k++]; k, {n, 200}]
f[n_] := Block[{b = 1, p = Prime@ n}, While[! PrimeQ[((b +1)^p + b^p)/(2b +1)], b++]; b]; Array[f, 70, 2] (* Robert G. Wilson v, Jun 13 2018 *)
for(n=2, 200, b=0; until(isprime((((b+1)^prime(n)+b^prime(n))/(2*b+1))), b++); print1(b,", ")) \\ corrected by Eric Chen, Jun 06 2018
[n: n in [1..10000] |IsPrime((9^n+7^n)/16)]
select(n->isprime((9^n+7^n)/16),[seq(n,n=1..10000,2)]); # Muniru A Asiru, Mar 27 2018
Select[Range[1, 10000], PrimeQ[(9^n+7^n)/16] &]
forprime(n=3, 10000, if(isprime((9^n+7^n)/16), print1(n, ", ")))
Comments