User Portlet
| Discussions |
|---|
| try: Select[{#, pQ[# - 1] == Boole[PrimeQ[#]]} & /@ Range[1, 1000, 2], (! Last[#]) &] you will get all primes plus Fermat's little friends in the form of: PowerMod[2, (n - 1), n ] == 1 eg: 341, 561, 645, 1105, 1387, 1729,... |
| There are a lot of smart people here, I'm sure someone has found a way to integrate the OEIS database so I don't have to leave Mathematica! If not I might be willing to implement that if there are others that would find it useful. ... |