Power(123,17) mod 3233 problem

Power(123,17) giver en floating point og sådan en kan du ikke tage mod på

du skal have fat på en "big integer" pakke

Du må jo nok kunne lave det på en anden måde.

StortTal := Power(123,17);
Encrypt  := StortTal - (Round(StortTal / N) * N);

Du tager StortTal dividerer med N, afrunder til nul decimaler og ganger med N igen.
Forskellen mellem dette facit og det oprindelige tal må vel være modulus.

Med den her kode:
procedure TForm1.Button1Click(Sender: TObject);
var
  Encrypt, StortTal, Decrypt, PQ : Real48;
Const
  P = 61;
  Q = 53;
  N = 3233;
  E = 17;
  D = 2753;
  Start = 123;
begin
  StortTal := Power(123,17);
  Encrypt  := StortTal - (Round(StortTal / N) * N);
  ShowMessage(FloatToStr(Encrypt));
end;

Får jeg en fejl som hedder: "EInvalidOp" og i programmet fortæller den: "Invalid Floating Point Operation"

// Jonatan

nej - fordi en floating point har ikke bits nok til at repræsentere de heltal

Pga. de er så store?

ja

3,xxxE35 er 3,xxxE35 +/- rigtigt mange milliarder d.v.s. at modulus giver ingen mening

derfor skal du have fat på en big integer pakke

jeg lavede engang den her kode:

program BigInt;

{$APPTYPE CONSOLE}

uses
  SysUtils;

function numval(c : char) : integer;

begin
  numval := ord(c) - ord('0');
end;

function strval(v : integer) : char;

begin
  strval := chr(ord('0') + v);
end;

function add(v1,v2 : string) : string;

var
  res : string;
  ix1, ix2, tmp, carry : integer;

begin
  res := '';
  carry := 0;
  ix1 := length(v1);
  ix2 := length(v2);
  while (carry > 0) or (ix1 > 0) or (ix2 > 0) do begin
      tmp := carry;
      if(ix1 > 0) then begin
          tmp := tmp + numval(v1[ix1]);
          ix1 := ix1 - 1;
      end;
      if(ix2 > 0) then begin
          tmp := tmp + numval(v2[ix2]);
          ix2 := ix2 - 1;
      end;
      res := strval(tmp mod 10) + res;
      carry := tmp div 10;
  end;
  add := res;
end;

function mulpow10(v : string; pow10 : integer) : string;

var
  res : string;
  i : integer;

begin
  res := v;
  for i := 1 to pow10 do begin
      res := res + '0';
  end;
  mulpow10 := res;
end;

function mulone(v : string; num : integer) : string;

var
  res : string;
  ix, tmp, carry : integer;

begin
  res := '';
  carry := 0;
  ix := length(v);
  while (carry > 0) or (ix > 0) do begin
      tmp := carry;
      if(ix > 0) then begin
          tmp := tmp + num * numval(v[ix]);
          ix := ix - 1;
      end;
      res := strval(tmp mod 10) + res;
      carry := tmp div 10;
  end;
  mulone := res;
end;

function mul(v1,v2 : string) : string;

var
  res : string;
  i : integer;

begin
  res := '0';
  for i := 1 to length(v2) do begin
      res := add(res,mulpow10(mulone(v1,numval(v2[i])),length(v2)-i));
  end;
  mul := res;
end;

var
  i : integer;
  fac : string;

begin
  (* functionality test *)
  writeln(add('1234','12345678'));
  writeln(mulpow10('123',2));
  writeln(mulone('123456',2));
  writeln(mul('12345','12345'));
  (* performance test *)
  fac := '1';
  for i := 1 to 500 do begin
    fac := mul(fac, inttostr(i));
    writeln(i,' ',fac);
  end;
end.

Og hvad er det, det er?

men der findes meget meget meget meget bedre big integer pakker en det !

http://www.delphiforfun.org/Programs/Library/big_integers.htm

Det er til rigtig store tal ik?

Kan du lave et eksempel på hvordan det kombineres med mit projekt, for facit giver et tal på 3 cifre

Ps. Spm: "Og hvad er det, det er?" var til hvad en Big integer pakke var :)

http://www.efg2.com/Lab/Library/Delphi/MathFunctions/Cryptography.htm

nederst "Multiple-Precision Arithmetic"

en big integer pakke er ent bibliotek som kan regne med vilkårligt store heltal
100 cifre, 1000 cifre, 10000 cifre

og du kigger i pakken efter pow og mod og sålaver du bare

BigIntMod(BigIntPower(123,17),N)

eller hvad den nu hedder

at resultatet er 0..N-1 og dermed ret lille gør jo ikke at mellemregningerne ikke
kan blive store

Ok... Mange tak for hjælp...

Mine forældre siger at jeg skal i seng nu, så kikker lige videre på det imorgen...

Men tak for hjælpen til nu :D

// Jonatan

jeg ligger et svar just in case

Jeg har prøvet at kikke på den pakke der hedder: "LongInt" på http://www.efg2.com/Lab/Library/Delphi/MathFunctions/Cryptography.htm (nederst)...

Jeg kan bare ikke finde ud af hvordan jeg bruger den :S

// Jonat

Jeg er ikke sikker på at jeg har regnet rigtigt, men måden kan måske være at bruge logaritmer, så undgår man tit overflow.

Prøv for sjovt denne variation:


procedure TForm1.Button1Click(Sender: TObject);
var  Encrypt, StortTal, Decrypt, PQ, LogStorttal, Nl: extended;
Const
  P = 61;
  Q = 53;
  N = 3233.0;
  E = 17;
  D = 2753;
  Start = 123;
begin
  StortTal := Power(123,17);
  LogStorttal := ln(storttal);
  Nl := ln(N);
  encrypt := exp((Round(logStortTal / Nl) - Nl));
  ShowMessage(FloatToStr(Encrypt));
end;

Hej arne_v...

Har nu arbejdet på det, og har fundet en løsning med denne Dll, http://www.delphiforfun.org/Programs/Library/big_integers.htm ... og det virker :)

Mange tak, du får Point...

lsc, din løsning gav 6,xxx og den skal give 855, og det gør arne_v's link.

Min løsning:
implementation

{$R *.dfm}     
uses UBigIntsV2, math;

procedure TForm1.Button1Click(Sender: TObject);
var g,e,m :Tinteger;
begin
g:=TInteger.create;
e:=TInteger.create;
m:=TInteger.create;
g.assign(IntToStr(Start));
e.assign('17');
m.assign(IntToStr(N));
g.modpow(e,m);
ShowMessage(g.ConvertToDecimalString(true));
g.free;
e.free;
m.free;

// Jonat.