mrd_prime 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: d17143e5a145a9c98fc99c910d98bfb27a9fe6a097998ce57685e14a6a0e28db
+  data.tar.gz: 8f3ad79e530c346e598298c94605a39f8c4460e7eef05e5a6bc3c86465fba9ee
+SHA512:
+  metadata.gz: 6f18db85331ff59361692fb557ce7811fdd443e07a86eb20493f295d0d64ab25e7d4cfcfa78641a6fc70f638f69fe72d0dd185cccefdb1556e2b5f63ed2b2fd9
+  data.tar.gz: f77022bba5faa89d1778316cf06fe5ffb720088bcd6b23990af4deed6a71db6092f85caaf466e7e27dd90da53effa651d321f3953374c3d9b039f7de875fcff4

data/.gitignore

@@ -0,0 +1,23 @@
+/.bundle/
+/.yardoc
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+
+# rspec failure tracking
+.rspec_status
+
+.vscode/*
+!.vscode/settings.json
+!.vscode/tasks.json
+!.vscode/launch.json
+!.vscode/extensions.json
+*.code-workspace
+
+# Local History for Visual Studio Code
+.history/
+
+Gemfile.lock

data/.rspec

@@ -0,0 +1,3 @@
+--format documentation
+--color
+--require spec_helper

data/.travis.yml

@@ -0,0 +1,6 @@
+---
+language: ruby
+cache: bundler
+rvm:
+  - 2.7.0
+before_install: gem install bundler -v 2.1.4

data/Gemfile

@@ -0,0 +1,7 @@
+source "https://rubygems.org"
+
+# Specify your gem's dependencies in mrd_prime.gemspec
+gemspec
+
+gem "rake", "~> 12.0"
+gem "rspec", "~> 3.0"

data/README.md

@@ -0,0 +1,62 @@
+# MrdPrime
+
+Miller-Rabin Prime Number Determination.
+
+A small number(~2^81) makes a definitive decision.
+A large number will fallback to `OpenSSL::BN::prime?`.
+
+### About the amount of calculation
+
+In order to reduce the amount of calculation, a decision is made using the minimum amount of calculation.
+For example, just testing on 2 will determine up to 2047.
+Testing on 2 and 3, gives a decision up to 1373653.
+
+	https://en.wikipedia.org/wiki/Miller%E2%80%93Rabin_primality_test#Testing_against_small_sets_of_bases
+	"if n < 3,317,044,064,679,887,385,961,981,
+	it is enough to test a = 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, and 41.
+	and the result will be deterministic."
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```
+gem 'mrd_prime'
+```
+
+And then execute:
+
+```
+$ bundle install
+```
+
+Or install it yourself as:
+
+```
+$ gem install mrd_prime
+```
+
+## Usage
+
+```
+require 'mrd_prime'
+
+small_number = (2**79 + 23)
+
+puts( small_number.mrd_prime? )  #=> true
+```
+
+If you want to use the name `prime?`, you can alias it.
+
+```rb
+class Integer
+	alias_method :prime?, :mrd_prime?
+end
+
+puts( (2**79 + 23).prime? )  #=> true
+```
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/Matsuyanagi/mrd_prime .
+

data/Rakefile

@@ -0,0 +1,6 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "mrd_prime"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start(__FILE__)

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/lib/mrd_prime.rb

@@ -0,0 +1,79 @@
+require "mrd_prime/version"
+
+require 'openssl'
+
+class Integer
+  # Returns true if +self+ is a prime number, else returns false.
+  # 
+  # https://en.wikipedia.org/wiki/Miller%E2%80%93Rabin_primality_test
+  # https://oeis.org/A014233/list
+  #
+  # "if n < 3,317,044,064,679,887,385,961,981,
+  # it is enough to test a = 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, and 41.
+  # and the result will be deterministic."
+  #
+  def mrd_prime?
+
+    return false if self < 2
+    return true  if self < 4
+    return false if self.even?
+
+    prime_and_max = {
+      2  =>  8321,      # 2047.  Strong pseudoprimes to base 2, [ 2047(=23x89), 3277(=29x113), 4033(=37x109), 4681(=31x151), 8321(=53x157) ] ~ 4681 is divided by a small prime number. 
+      3  =>  1373653,
+      5  =>  25326001,
+      7  =>  3215031751,
+      11  =>  2152302898747,
+      13  =>  3474749660383,
+      17  =>  341550071728321,
+      19  =>  341550071728321,
+      23  =>  3825123056546413051,
+      29  =>  3825123056546413051,
+      31  =>  3825123056546413051,
+      37  =>  318665857834031151167461,
+      41  =>  3317044064679887385961981,
+    }
+    
+    last_p = 1
+    prime_and_max.each do |p,m|
+      return true if self == p
+      return false if self % p == 0
+      last_p = p
+    end
+    return true if self < last_p * last_p
+
+    p_1 = self - 1
+
+    d = p_1
+    while d.even?
+      d >>= 1
+    end
+
+    prime_and_max.each do |a,m|
+      x = a.pow( d, self )
+
+      if x == 1
+        return true if self < m
+        next
+      end
+
+      td = d
+      while td != p_1 && x != p_1
+        x = x.pow( 2, self )
+        td <<= 1
+      end
+
+      if td == p_1
+        return false
+      else
+        return true if x == p_1 && self < m
+      end
+
+    end
+
+    # OpenSSL fallback
+    # return true
+    return OpenSSL::BN.new( self ).prime?
+  end
+end
+

data/lib/mrd_prime/version.rb

@@ -0,0 +1,3 @@
+module MrdPrime
+  VERSION = "0.1.0"
+end

data/mrd_prime.gemspec

@@ -0,0 +1,27 @@
+require_relative 'lib/mrd_prime/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "mrd_prime"
+  spec.version       = MrdPrime::VERSION
+  spec.authors       = ["Matsuyanagi"]
+  spec.email         = ["1475814+Matsuyanagi@users.noreply.github.com"]
+
+  spec.summary       = %q{A definitive Miller-Rabin prime determination test for small numbers.}
+  spec.description   = %q{A definitive Miller-Rabin prime determination test for small numbers (~2^81). For large numbers, fall back to OpenSSL::BN.}
+  spec.homepage      = "https://github.com/Matsuyanagi/mrd_prime"
+  spec.required_ruby_version = Gem::Requirement.new(">= 2.3.0")
+  spec.license       = "MIT"
+
+  spec.metadata["homepage_uri"] = spec.homepage
+  spec.metadata["source_code_uri"] = "https://github.com/Matsuyanagi/mrd_prime"
+  spec.metadata["changelog_uri"] = "https://github.com/Matsuyanagi/mrd_prime"
+
+  # Specify which files should be added to the gem when it is released.
+  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
+  spec.files         = Dir.chdir(File.expand_path('..', __FILE__)) do
+    `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  end
+
+  spec.require_paths = ["lib"]
+  spec.add_dependency 'openssl', '~> 2'
+end

metadata

@@ -0,0 +1,72 @@
+--- !ruby/object:Gem::Specification
+name: mrd_prime
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- Matsuyanagi
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2020-08-07 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: openssl
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2'
+description: A definitive Miller-Rabin prime determination test for small numbers
+  (~2^81). For large numbers, fall back to OpenSSL::BN.
+email:
+- 1475814+Matsuyanagi@users.noreply.github.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- lib/mrd_prime.rb
+- lib/mrd_prime/version.rb
+- mrd_prime.gemspec
+homepage: https://github.com/Matsuyanagi/mrd_prime
+licenses:
+- MIT
+metadata:
+  homepage_uri: https://github.com/Matsuyanagi/mrd_prime
+  source_code_uri: https://github.com/Matsuyanagi/mrd_prime
+  changelog_uri: https://github.com/Matsuyanagi/mrd_prime
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 2.3.0
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.1.4
+signing_key: 
+specification_version: 4
+summary: A definitive Miller-Rabin prime determination test for small numbers.
+test_files: []