isok.sf 445 B

1234567891011121314151617181920212223242526272829303132333435
  1. #!/usr/bin/ruby
  2. func isok(n,k) {
  3. 1..prime(n) -> all {|b| powmod(b, (k-1)/2, k).is_congruent(kronecker(b, k), k) }
  4. }
  5. var arr = [
  6. 561
  7. 1729
  8. 1729
  9. 399001
  10. 399001
  11. 1857241
  12. 1857241
  13. 6189121
  14. 14469841
  15. 14469841
  16. 14469841
  17. 86566959361
  18. 311963097601
  19. 369838909441
  20. 6389476833601
  21. 6389476833601
  22. 1606205228509922041
  23. 1606205228509922041
  24. 1606205228509922041
  25. 1606205228509922041
  26. ]
  27. for k,v in (arr.kv) {
  28. say [k, isok(k+1, v)]
  29. }