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.
Select[ Prime[ Range[1, 100000] ], PrimeQ[ (15^# + 14^#)/29 ]& ]
is(n)=ispseudoprime((15^n+14^n)/29) \\ Charles R Greathouse IV, Feb 20 2017
Select[ Prime[ Range[1, 100000] ], PrimeQ[ (11^# + 10^#)/21 ]& ]
is(n)=ispseudoprime((11^n+10^n)/21) \\ Charles R Greathouse IV, May 22 2017
Select[ Prime[ Range[1, 100000] ], PrimeQ[ (13^# + 12^#)/25 ]& ]
is(n)=ispseudoprime((13^n+12^n)/25) \\ Charles R Greathouse IV, May 22 2017
Select[ Prime[ Range[1, 100000] ], PrimeQ[ (14^# + 13^#)/27 ]& ]
forprime(p=3,10^6, if(ispseudoprime((14^p + 13^p)/27), print1(p,", ") ) ); \\ Joerg Arndt, Jul 29 2013
a(10) = 53 because (10^p + 11^p)/21 is composite for all p < 53 and prime for p = 53.
lmt = 4200; f[n_] := Block[{p = 2}, While[p < lmt && !PrimeQ[((n + 1)^p + n^p)/(2n + 1)], p = NextPrime@ p]; If[p > lmt, 0, p]]; Do[Print[{n, f[n] // Timing}], {n, 1000}] (* Robert G. Wilson v, Dec 01 2014 *)
a(n)=forprime(p=3, , if(ispseudoprime((n^p+(n+1)^p)/(2*n+1)), return(p)))
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
Comments