From 1a230bd61e640bb70290402fac95d03da870a0b3 Mon Sep 17 00:00:00 2001
From: Omar <59146036+omarmahamid@users.noreply.github.com>
Date: Fri, 29 Apr 2022 16:30:05 +0300
Subject: [PATCH] Add New Prime Check (#3036)

---
 .../com/thealgorithms/maths/PrimeCheck.java   | 47 ++++++++++++++++++-
 1 file changed, 45 insertions(+), 2 deletions(-)

diff --git a/src/main/java/com/thealgorithms/maths/PrimeCheck.java b/src/main/java/com/thealgorithms/maths/PrimeCheck.java
index f2587f70..fd06cb7f 100644
--- a/src/main/java/com/thealgorithms/maths/PrimeCheck.java
+++ b/src/main/java/com/thealgorithms/maths/PrimeCheck.java
@@ -10,9 +10,15 @@ public class PrimeCheck {
         System.out.print("Enter a number: ");
         int n = scanner.nextInt();
         if (isPrime(n)) {
-            System.out.println(n + " is a prime number");
+            System.out.println("algo1 verify that " + n + " is a prime number");
         } else {
-            System.out.println(n + " is not a prime number");
+            System.out.println("algo1 verify that " + n + " is not a prime number");
+        }
+
+        if (fermatPrimeChecking(n, 20)) {
+            System.out.println("algo2 verify that " + n + " is a prime number");
+        } else {
+            System.out.println("algo2 verify that " + n + " is not a prime number");
         }
         scanner.close();
     }
@@ -38,4 +44,41 @@ public class PrimeCheck {
         }
         return true;
     }
+
+    /**
+     * *
+     * Checks if a number is prime or not
+     *
+     * @param n the number
+     * @return {@code true} if {@code n} is prime
+     */
+    public static boolean fermatPrimeChecking(int n, int iteration){
+      long a;
+      int up = n - 2, down = 2;
+      for(int i=0;i<iteration;i++){
+        a = (long)Math.floor(Math.random()*(up - down + 1) + down);
+        if(modPow(a,n-1,n) != 1){
+          return false;
+        }
+      }
+      return true;
+    }
+
+
+    /**
+     * *
+     * @param a basis
+     * @param b exponent
+     * @param c modulo
+     * @return (a^b) mod c
+     */
+    private static long modPow(long a, long b, long c){
+        long res = 1;
+        for (int i = 0; i < b; i++)
+        {
+            res *= a;
+            res %= c;
+        }
+        return res % c;
+    }
 }