fact2p52-1.

16.06.00  17:18  pm.


YES.

BEST NEEDS basic64


 progm--    fin.soc.acorncpusr.donWn994/2.pi.facthcf2/. byte/zeta4

don mcdonald,    18.07.99  22:25  16.06.00  15:01


Examples : 1E15-1,  1E18/999 999,  2^37-1,  2^31-1

2^33-9 prime,  6763*10627*29947 Maple IsPrime

1E10-1, 1E11-1,  1E12-1, 1E14-1, 2^32+1 Euler, 13^10+1

1E13-1, 1E16-1,  100 895 598 169 Mersenne

10 662 526 601,  15 527 402 881. 

===== Below. (n+1)^p-n^p, generalised mersenne #

5^19-4^19= 830339* 8.min

6^13 prime, 3^23 =47*  ,  8 (157*) 9 (6553*) 11 (937*) 12^13 (157*) ,

 10^13- (2081)divide by zero?

13^13 = 443*442 171 334 567, 14^13=490 839 66666 1891 = 

 1 198 679* 409 483 829, 15^13 (79*547*677*39,394,031) = (13n+1)

16^13-15^13 max = 53*403391*1196 18227

mults of  2

mults of  3

mults of  5

mults of  7


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?9955953199793821

9955953199793821 = 9955953199793820

$10/20 S.N. #

too big.


proc Hhcf(x = 9955953199793820


 HCF = 30

9955953199793820  Factor = 30 * 331865106659794

2 * 15

3 * 5

5 prime.           centisec = 8


 HCF = 2

331865106659794  Factor = 2 * 165932553329897

2 * 1

1 prime.           centisec = 12


 HCF = 401

165932553329897  Factor = 401 * 413796891097

401 prime.         centisec = 83


 HCF = 4951

413796891097  Factor = 4951 * 83578447

4951 prime.        centisec = 1166

83578447 prime.    centisec = 1578


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?2^51-1

*******2^51-1 = 2251799813685247


proc Hhcf(x = 2251799813685247


 HCF = 7

2251799813685247  Factor = 7 * 321685687669321

7 prime.           centisec = 8


 HCF = 103

321685687669321  Factor = 103 * 3123162016207

103 prime.         centisec = 21


 HCF = 2143

3123162016207  Factor = 2143 * 1457378449

2143 prime.        centisec = 492

11119 * 131071		****************  what about the other one **

131071 prime.      centisec = 998


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?2^51+1

****	2^51+1 = 2251799813685249


proc Hhcf(x = 2251799813685249


 HCF = 3

2251799813685249  Factor = 3 * 750599937895083

3 * 1

1 prime.           centisec = 8


 HCF = 3

750599937895083  Factor = 3 * 250199979298361

3 * 1

1 prime.           centisec = 11


 HCF = 307

250199979298361  Factor = 307 * 814983645923

307 prime.         centisec = 59


 HCF = 2857

814983645923  Factor = 2857 * 285258539

2857 prime.        centisec = 671

6529 * 43691   ************ 3 factors over 2800.******

		what about the other one.

43691 prime.       centisec = 967


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?2^52-1

2^52-1 = 4503599627370495                	**********


proc Hhcf(x = 4503599627370495


 HCF = 15

4503599627370495  Factor = 15 * 300239975158033

3 * 5

5 prime.           centisec = 8


 HCF = 53

300239975158033  Factor = 53 * 5664905191661

53 prime.          centisec = 15


 HCF = 157

5664905191661  Factor = 157 * 36082198673

157 prime.         centisec = 31


 HCF = 1613

36082198673  Factor = 1613 * 22369621

1613 prime.        centisec = 48

2731 * 8191

8191 prime.       ****** 2^13-1. ** centisec = 172


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?2^52+1

2^52+1 = 4503599627370497


proc Hhcf(x = 4503599627370497


 HCF = 17

4503599627370497  Factor = 17 * 264917625139441

17 prime.          centisec = 8

  x  > 2 ^31.  TESTING slow e.g. minutes.

TRY FOR j = 246723  TO SQRT = 16276290.2757182661  STEP 2.  .. | centisec

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  303451  4145cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  361201  7184cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  418951  10222cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  476701  13265cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  534451  16309cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  592201  19355cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  649951  22401cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  707701  25447cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  765451  28493cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  823201  31393cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,**858001


factors ?   x1 =  j*x1 =   858001 *  308761441

   centiseconds = 33178


 You have 10 seconds to press a key.

858001 prime.      centisec = 34220

308761441 prime.   centisec = 35014

858001 * 308761441                 	**  very interesting ** 6 mins.**



Unknown or missing variable  at  erl ...          91


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?503646201388795

503646201388795 = 503646201388795


proc Hhcf(x = 503646201388795


 HCF = 5

503646201388795  Factor = 5 * 100729240277759

5 prime.           centisec = 7

  x  > 2 ^31.  TESTING slow e.g. minutes.

TRY FOR j = 93811  TO SQRT = 10036395.7812433336  STEP 2.  .. | centisec

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  151201  4182cs. / 

,,,,,,,,,,,,,**165463


factors ?   x1 =  j*x1 =   165463 *  608771993

   centiseconds = 4932


 You have 10 seconds to press a key.

165463 prime.      centisec = 5951

608771993 prime.   centisec = 7067

165463 * 608771993	pin.*********


Unknown or missing variable  at  erl ...          91


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?1E15-1

1E15-1 = 999999999999999             ***************


proc Hhcf(x = 999999999999999


 HCF = 3

999999999999999  Factor = 3 * 333333333333333

3 * 1

1 prime.           centisec = 7


 HCF = 3441

333333333333333  Factor = 3441 * 96871064613

3 * 1147

31 * 37

37 prime.          centisec = 12


 HCF = 123

96871064613  Factor = 123 * 787569631

3 * 41

41 prime.          centisec = 15

271 * 2906161

2906161 prime.     centisec = 94


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?1E16+1

1E16+1 = 10000000000000000                          ****too big.


proc Hhcf(x = 10000000000000000


 HCF = 10

10000000000000000  Factor = 10 * 1000000000000000

2 * 5

5 prime.           centisec = 8


 HCF = 1250

1000000000000000  Factor = 1250 * 800000000000

2 * 625

[...]


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?1E15+1

1E15+1 = 1000000000000001


proc Hhcf(x = 1000000000000001


 HCF = 1001

1000000000000001  Factor = 1001 * 999000999001

7 * 143

11 * 13

13 prime.          centisec = 8


 HCF = 211

999000999001  Factor = 211 * 4734601891

211 prime.         centisec = 36


 HCF = 241

4734601891  Factor = 241 * 19645651

241 prime.         centisec = 45

2161 * 9091

9091 prime.        centisec = 144


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?1E15+3

1E15+3 = 1000000000000003	********* Is prime ??


proc Hhcf(x = 1000000000000003

  x  > 2 ^31.  TESTING slow e.g. minutes.

TRY FOR j = 931323  TO SQRT = 31622776.6016838402  STEP 2.  .. | centisec

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  988051  4168cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1045801  7220cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1103551  10273cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1161301  13326cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1219051  16379cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1276801  19432cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1334551  22485cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1392301  25538cs. / 

,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,  1450051  28591cs. / 

,,,,,,,,

Escape  at  erl ...          82


factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?1E16-1

1E16-1 = 10000000000000000



factor 5 byte integer, max 2^51  gets most factors .TEST bugs

enter no. / expression  . Qq/ 0 <CR> = quit  ?

progm facthcf2/.   e n d.    CLOSEs * RAM SPOOL

