From efc4b04903e8fcaa074cba27034628a0eaa9e22e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Marcin=20Wo=C5=BAniak?= <y0rune@aol.com>
Date: Sat, 19 Dec 2020 22:24:14 +0100
Subject: [PATCH] Added module3
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: Marcin Woźniak <y0rune@aol.com>
---
 8-lab/szyfrowanie-krzywe.rb | 157 +++++++++++++++++++++++++++++++++---
 1 file changed, 145 insertions(+), 12 deletions(-)
 mode change 100644 => 100755 8-lab/szyfrowanie-krzywe.rb

diff --git a/8-lab/szyfrowanie-krzywe.rb b/8-lab/szyfrowanie-krzywe.rb
old mode 100644
new mode 100755
index c0aeec4..3307932
--- a/8-lab/szyfrowanie-krzywe.rb
+++ b/8-lab/szyfrowanie-krzywe.rb
@@ -9,23 +9,156 @@
 ###################################
 
 load '../module.rb'
-load '../6/krzywa.rb'
+load '../module2.rb'
 
-def genKeysElGamalKrzywaEliptyczna
-  p = generate(300)
+puts `clear`
 
-  if primalityTest(p)
+def wielokrotnoscPunktu2(a,b,p,n,x,y)
+  counter = n - 1
+  punktQ = []
+  punktP = [x,y]
+  temp = [x,y]
+  while counter > 0
+    punktQ = sumaPunktow(a,b,p,temp[0],temp[1],punktP[0],punktP[1])
+    temp = punktQ
+    counter = counter - 1
+  end
+  return punktQ
+end
 
-    krzywa = generatorKrzywej(p)
-    a = krzywa[0].to_i
-    b = krzywa[1].to_i
+def wielokrotnoscPunktu(a,b,p,n,x,y)
+  punktQ = [x,y]
+  punktR = ["e","e"]
 
-    while 
+  while n > 0
+    if n % 2 == 1
+      punktR = sumaPunktow(a,b,p,punktR[0],punktR[1],punktQ[0],punktQ[1])
+      n = n - 1
+    end
+    punktQ = sumaPunktow(a,b,p,punktQ[0],punktQ[1],punktQ[0],punktQ[1])
+    n = n / 2
+  end
+  return punktR
+end
 
-    # if x < ord = #E(Fp) 
-    # x = SecureRandom.random_number(1..ord)
-    # end
+def liczenieOrd(a,b,p)
+  threads = []
+  punktyNaKrzywej = Array.new
+  unikalnePunktyNaKrzywej = Array.new
 
-    
+  threads << Thread.new {
+    counter = 10
+    while counter > 0
+      punktyNaKrzywej <<  punktNaKrzywej(a,b,p)[0]
+      counter = counter -1
+    end
+  }
+
+  threads << Thread.new {
+    counter = 10
+    while counter > 0
+      punktyNaKrzywej <<  punktNaKrzywej(a,b,p)[0]
+      counter = counter -1
+    end
+  }
+
+  threads << Thread.new {
+    counter = 10
+    while counter > 0
+      punktyNaKrzywej <<  punktNaKrzywej(a,b,p)[0]
+      counter = counter -1
+    end
+  }
+
+  ThreadsWait.all_waits(*threads)
+
+  punktyNaKrzywej.sort.each { |e|
+    unikalnePunktyNaKrzywej.push(e) if e != unikalnePunktyNaKrzywej[-1]
+  }
+
+  return unikalnePunktyNaKrzywej.count
+end
+
+def generowanieKluczyElGamalKrzywaEliptyczna(k)
+  while true
+  p = generate(k)
+    if (primalityTest(p)) && (p % 4 == 3)
+      krzywa = generatorKrzywej(p)
+      a = krzywa[0].to_i
+      b = krzywa[1].to_i
+
+      punktyNaKrzywej = Array.new
+
+      punktP = punktNaKrzywej(a,b,p)
+      ord = liczenieOrd(a,b,p)
+
+      while true
+        x = SecureRandom.random_number(1..ord)
+        if x < ord
+          punktQ = wielokrotnoscPunktu(a,b,p,x,punktP[0],punktP[1])
+          pubKey = [a,b,p,punktP[0],punktP[1],punktQ[0],punktQ[1]]
+          privKey = [a,b,p,punktP[0],punktP[1],punktQ[0],punktQ[1],x]
+          return a,b,p,punktP[0],punktP[1],punktQ[0],punktQ[1],x
+        end
+      end
+    end
   end
 end
+
+def algorytmKodowania(a,b,p,m,n,u)
+  if (m < n) && (p > n*u)
+    for i in 1..u
+      x = (m * u % p) + (i % p)
+      fx = rownanieKrzywej(a,b,p,x)
+      if remSqEuler(fx,p)
+        y = betterExponentiation(fx,((p+1)/4),p)
+      end
+    end
+  else
+    puts "Nieprawidołowe dane"
+  end
+  return [x,y]
+end
+
+def algorytmDeKodowania(a,b,p,x,y,n,u)
+  m = (x - 1) / u
+  return m
+end
+
+def szyfrowanieElGamala(m,n,u,a,b,p,px,py,qx,qy)
+  c = algorytmKodowania(a,b,p,m,n,u)
+  puts "Wiadomość na prostej jako punkt #{c.inspect}"
+  return c
+end
+
+def main
+  k = 30
+  m = 73
+  n = m + SecureRandom.random_number(0..1000000)
+  u = SecureRandom.random_number(30..50)
+  a,b,p,px,py,qx,qy,x = generowanieKluczyElGamalKrzywaEliptyczna(k)
+
+  pubKey = [a,b,p,px,py,qx,qy]
+  privKey = [a,b,p,px,py,qx,qy,x]
+
+  pm = szyfrowanieElGamala(m,n,u,a,b,p,px,py,qx,qy)
+
+  # BOB
+  y = SecureRandom.random_number(0..liczenieOrd(a,b,p))
+  c1 = wielokrotnoscPunktu(a,b,p,y,px,py)
+  yq = wielokrotnoscPunktu(a,b,p,y,qx,qy)
+  c2 = sumaPunktow(a,b,p,pm[0],pm[1],yq[0],yq[1])
+  puts "Ciphers: #{c1} + #{c2}"
+
+  # ALICE
+  xc1 = wielokrotnoscPunktu(a,b,p,x,c1[0],c1[1])
+  pmd = sumaPunktow(a,b,p,c2[0],c2[1],xc1[0],-xc1[1])
+  d = algorytmDeKodowania(a,b,p,pmd[0],pmd[1],n,u)
+  puts "Dekodowana wiadomość #{d.inspect}"
+end
+
+#puts generowanieKluczyElGamalKrzywaEliptyczna(2048).inspect
+#puts wielokrotnoscPunktu(8,10,19,3,15,3).inspect
+#puts algorytmSzyfrowania(8,10,19,29102901920190,29102901920199,50).inspect
+#puts algorytmDeSzyfrowania(8,10,19,12,17,50)
+puts main.inspect