Sophie's dyslexic sister's primes? In[60]:= Table[PowerMod[k, # - 1, #], {k, 22}] & /@ ((2*# - 1)*# & /@ Select[Prime[Range[105]], PrimeQ[2 # - 1] &]) Out[60]= {{1, 2, 3, 4, 5, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 4}, {1, 4, 9, 1, 10, 6, 4, 4, 6, 10, 1, 9, 4, 1, 0, 1, 4, 9, 1, 10, 6, 4}, {1, 64, 1, 1, 64, 64, 77, 64, 1, 1, 64, 1, 78, 14, 64, 1, 1, 64, 64, 64, 77, 1}, {1, 628, 1, 1, 628, 628, 1, 628, 1, 1, 1, 1, 628, 628, 628, 1, 628, 628, 665, 628, 1, 628}, {1, 1768, 1, 1, 1, 1768, 1768, 1768, 1, 1768, 1768, 1, 1, 1, 1, 1, 1768, 1768, 1, 1, 1768, 1}, {1, 1, 1, 1, 2554, 1, 2554, 1, 1, 2554, 2554, 1, 2554, 2554, 2554, 1, 2554, 1, 1, 2554, 2554, 2554}, {1, 12088, 1, 1, 12088, 12088, 12088, 12088, 1, 1, 1, 1, 1, 1, 12088, 1, 1, 12088, 1, 12088, 12088, 12088}, {1, 1, 1, 1, 18334, 1, 1, 1, 1, 18334, 18334, 1, 18334, 1, 18334, 1, 18334, 1, 18334, 18334, 1, 18334}, {1, 37948, 1, 1, 37948, 37948, 1, 37948, 1, 1, 37948, 1, 1, 37948, 37948, 1, 37948, 37948, 1, 37948, 1, 1}, {1, 1, 1, 1, 48514, 1, 48514, 1, 1, 48514, 1, 1, 1, 48514, 48514, 1, 48514, 1, 1, 48514, 48514, 1}, {1, 78208, 1, 1, 78208, 78208, 78208, 78208, 1, 1, 1, 1, 78208, 1, 78208, 1, 78208, 78208, 1, 78208, 78208, 78208}, {1, 87988, 1, 1, 1, 87988, 1, 87988, 1, 87988, 1, 1, 87988, 87988, 1, 1, 1, 87988, 87988, 1, 1, 87988}, {1, 1, 1, 1, 103738, 1, 1, 1, 1, 103738, 103738, 1, 103738, 1, 103738, 1, 1, 1, 1, 103738, 1, 103738}, {1, 145528, 1, 1, 1, 145528, 1, 145528, 1, 145528, 145528, 1, 145528, 145528, 1, 1, 145528, 145528, 1, 1, 1, 1}, {1, 186964, 1, 1, 186964, 186964, 1, 186964, 1, 1, 186964, 1, 186964, 186964, 186964, 1, 1, 186964, 1, 186964, 1, 1}, {1, 217468, 1, 1, 1, 217468, 217468, 217468, 1, 217468, 1, 1, 217468, 1, 1, 1, 1, 217468, 217468, 1, 217468, 217468}, {1, 1, 1, 1, 225454, 1, 1, 1, 1, 225454, 225454, 1, 1, 1, 225454, 1, 225454, 1, 225454, 225454, 1, 225454}, {1, 267544, 1, 1, 267544, 267544, 267544, 267544, 1, 1, 267544, 1, 267544, 1, 267544, 1, 1, 267544, 1, 267544, 267544, 1}, {1, 285388, 1, 1, 285388, 285388, 1, 285388, 1, 1, 1, 1, 1, 285388, 285388, 1, 1, 285388, 1, 285388, 1, 285388}, {1, 383248, 1, 1, 383248, 383248, 1, 383248, 1, 1, 383248, 1, 383248, 383248, 383248, 1, 383248, 383248, 383248, 383248, 1, 1}, {1, 495508, 1, 1, 495508, 495508, 495508, 495508, 1, 1, 495508, 1, 1, 1, 495508, 1, 495508, 495508, 1, 495508, 495508, 1}, {1, 595684, 1, 1, 595684, 595684, 1, 595684, 1, 1, 1, 1, 1, 595684, 595684, 1, 595684, 595684, 595684, 595684, 1, 595684}} Interesting source of bit sequences, anyway. --rwg On Mon, Jun 2, 2014 at 5:05 PM, Bill Gosper <billgosper@gmail.com> wrote:
In[49]:= Table[PowerMod[k, 497502, 497503], {k, 69}]
Out[49]= {1, 495508, 1, 1, 495508, 495508, 495508, 495508, 1, 1, 495508, 1, 1, 1, 495508, 1, 495508, 495508, 1, 495508, 495508, 1, 1, 495508, 1, 495508, 1, 495508, 495508, 1, 1, 495508, 495508, 1, 1, 1, 495508, 495508, 1, 1, 495508, 1, 495508, 495508, 495508, 495508, 495508, 1, 1, 495508, 495508, 1, 1, 495508, 1, 1, 1, 1, 1, 495508, 495508, 495508, 495508, 1, 495508, 1, 1, 495508, 1} ? --rwg