Hallo,
auf dem Forum "Matheplanet" läuft gerade eine Art Wettbewerb, bei dem es um möglichst schnelle Berechnung von langen Zahlen geht. Die Berechnung hat eigentlich wenig Sinn, es geht nur um Geschwindigkeit.
25 Millionen mal werden zwei lange Zahlen multipliziert, dann unter eine Grenze reduziert und das Ergebnis mit 2^128 div und mod geteilt.
Nachdem ich mit Delphi 5 im Vergleich zu anderen Programmiersprachen nicht ganz so gut war, habe ich es jetzt auf Console umgestellt und es wurde etwas schneller. Es nutzt die Langzahlroutinen von Wolfgang Ehrhardt.
Delphi-Quellcode:
program matheplanet_projekt2;
{$APPTYPE CONSOLE}
uses SysUtils,
mp_base, mp_types;
var
a, b, soll, sub, produkt, filter: mp_int;
Time1: TDateTime;
nn, nad, nsu: integer;
begin
Writeln('hyperG-Problem');
Writeln;
mp_init6(a, b, soll, sub, produkt, filter);
nad:=0;
nsu:=0;
mp_read_decimal_astr(a,'170141183460469231731687303715884105757');
mp_read_decimal_astr(b,'170141183460469231731687303715884105703');
mp_read_decimal_astr(soll,'57896044618658097711785492504343953926634992332820282019728792003956564819968');
mp_read_decimal_astr(sub,'11111111111111111111111111111111111111111111111111111111111111111111111111111');
mp_read_decimal_astr(filter,'340282366920938463463374607431768211455');
time1 := now;
for nn:=1 to 25000000 do begin
mp_mul(a,b,produkt);
while mp_cmp_mag(produkt,soll)<0 do begin
mp_add(produkt,sub,produkt);
inc(nad);
end;
repeat
mp_sub(produkt,sub,produkt);
inc(nsu);
until mp_cmp_mag(soll,produkt)>0;
mp_shr(produkt,128,a);
mp_and(produkt,filter,b);
end;
Writeln(FormatDateTime('SS:ZZZ', now - Time1));
Writeln(mp_adecimal(produkt));
Writeln(nad);
Writeln(nsu);
Writeln;
Writeln('Programmende mit ENTER');
mp_clear6(a, b, soll, sub, produkt, filter);
readln;
end.
Das Ergebnis ist
13:325
47306219819553867530681020915755484915944251665601 642944764215683054098939729
83884810
25000000
Da ich gern die "Ehre von Delphi" verteidigen würde, suche ich noch nach Ideen zur Beschleunigung. Bei mir läuft die Routine in etwa 13,3 Sekunden durch.
In der ZIP-Datei sind auch die nötigen Quelltexte von Wolfgang Erhardts Routinen.
Danke für evtl. Ideen
Steffen