mersenneforum.org > Math Mersenne NUmber Factoring
 Register FAQ Search Today's Posts Mark Forums Read

 2002-11-25, 12:44 #1 andi314     Nov 2002 4A16 Posts Mersenne NUmber Factoring I recently made a program with which you can factor Mersenne Numbers over 2^80000000. Anyone who is interrested can have this program for free.( You only need UBASIC) Please contact me at andi314@gmx.de Mod edit: Title changed to remove HTML formatting...
 2002-11-25, 18:15 #2 Deamiter     Sep 2002 11710 Posts How fast can your program factor to 2^66? That seems to be the benchmark for similar factoring programs...
2002-11-26, 13:23   #3
andi314

Nov 2002

2·37 Posts

Quote:
 Originally Posted by Deamiter How fast can your program factor to 2^66? That seems to be the benchmark for similar factoring programs...
Sorry bur i don't know the exact dates because I haven't tested it yet.

But for anyone who wants to join to help me finding some factors:
:arrow: YOU ARE WELCOMED!!!

 2002-11-26, 13:32 #4 smh     "Sander" Oct 2002 52.345322,5.52471 29·41 Posts Can you post the source here? (between code tags)
2002-11-26, 14:25   #5
andi314

Nov 2002

7410 Posts

Quote:
 Originally Posted by smh Can you post the source here? (between code tags)
it is a UBASIC program:
[code:1]001 p=100000 ' Change this value if you wan to have a faster screen output
002 IF exist(&quot;allfactor.txt&quot;) print &quot;Please send the allfactor.txt file to andi314@gmx.de&quot;
003 if not exist(&quot;start.ubd&quot;) print &quot;Mersenne Factorisation Program&quot;
009 if not exist(&quot;start.ubd&quot;) print &quot;Please enter an exponent and a starting factor!&quot;
010 if exist(&quot;start.ubd&quot;) open &quot;start.ubd&quot; for input as #1
011 if exist(&quot;factor.txt&quot;) doscmd &quot;start file.exe&quot;
012 if exist(&quot;start.ubd&quot;) input#1,a
013 if exist(&quot;start.ubd&quot;) input#1,b
014 if exist(&quot;start.ubd&quot;) print &quot;exponent=&quot;;a
015 if exist(&quot;start.ubd&quot;) print &quot;current factor=&quot;;b
016 if exist(&quot;start.ubd&quot;) close #1
017 if exist(&quot;start.ubd&quot;) then e=a
018 if exist(&quot;start.ubd&quot;) goto 080
019 input&quot;Exponent=&quot;;e
020 if PRMDIV(e)&lt;>e print &quot;You didn't enter a prime number!!! Please try again!&quot;
030 if PRMDIV(e)&lt;>e goto 019
040 input&quot;Start=&quot;;b
080 c=int(b/2/e)-2
090 if c&lt;=-1 then clr c
100 inc c
110 if inkey=&quot;E&quot; goto 930
120 k=2*c*e+1
126 if c@p=0 print e,c,k,(int(log(k)/log(2)*100000000)) /100000000
127 if k@3=0 or k@5=0 or k@7=0 goto 100
150 if k@8&lt;>7 and k@8&lt;>1 goto 100
160 if c@1000000=0 open &quot;start.ubd&quot; for create as #1
170 if c@1000000=0 print #1,e, k
180 if c@1000000=0 close #1
190 q=1
200 x=e
210 y=int(log(e)/log(2))+1
300 dec y
400 if x-2^y>=0 z=1
500 if x-2^y&lt;0 clr z
510 q=(q^2)@k
600 if z=1 then q=(q*2)@k
700 if z=1 x=x-2^y
800 if y>0 goto 300
900 if q>1 goto 100
930 open &quot;start.ubd&quot; for create as #1
940 print #1,e, k
950 close #1
951 if modpow(2,e,k)=1 kill &quot;start.ubd&quot;
952 if modpow(2,e,k)=1 open &quot;factor.txt&quot; for create as #1
953 if modpow(2,e,k)=1 print#1,e,k
954 if modpow(2,e,k)=1 close #1
955 if exist(&quot;factor.txt&quot;) doscmd &quot;start file.exe&quot;
956 if modpow(2,e,k)=1 print e;&quot;has a factor:&quot;;k
957 if modpow(2,e,k)&lt;>1 print e;&quot;has no factor up to&quot;;k
958 if modpow(2,e,k)=1 print &quot;Please send the allfactor.txt file to andi314@gmx.de&quot;
960 print &quot;ENDE&quot;
970 SYSTEM
[/code:1]

 Similar Threads Thread Thread Starter Forum Replies Last Post aketilander Operazione Doppi Mersennes 1 2012-11-09 21:16 kurtulmehtap Math 12 2010-05-03 14:02 kerguilloud Math 6 2005-01-20 10:16 romanesquefr Miscellaneous Math 17 2004-12-17 12:05 Unregistered Math 4 2004-09-30 21:32

All times are UTC. The time now is 21:05.

Mon Jul 4 21:05:15 UTC 2022 up 81 days, 19:06, 0 users, load averages: 1.55, 1.41, 1.35