mersenneforum.org  

Go Back   mersenneforum.org > Math Stuff > Computer Science & Computational Number Theory > PARI/GP

Reply
 
Thread Tools
Old 2010-11-21, 23:10   #1618
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
The former: A power check would be necessary, in order to avoid wasted time.
Do you want it to do trial division by small primes, some rho and SQUFOF, and a probable-prime test too, in case you forgot to do those? How about stopping ECM at some point and switching to a general-purpose factoring algorithm like SIQS?



Quote:
Originally Posted by 3.14159 View Post
The latter: Please show me a link which shows the factors of RSA1024.
So just because no one has made the factors of a particular composite public, you feel safe assuming no one can find the factors?

And of course when you're using an RSA key you're looking to keep it secret in the future, not the past.

But hey, maybe your security just isn't valuable enough to protect.
CRGreathouse is offline   Reply With Quote
Old 2010-11-21, 23:43   #1619
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

168010 Posts
Default

sigma(71608085354600097239391307^4) took an hour to factor via SIQS. The disappointing part is, that the cofactor's factors were a p7, p18, and p77. It's a surprise that ECM missed a p7 and a p18.

Last fiddled with by 3.14159 on 2010-11-21 at 23:43
3.14159 is offline   Reply With Quote
Old 2010-11-22, 00:03   #1620
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3·1,993 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
sigma(71608085354600097239391307^4) took an hour to factor via SIQS. The disappointing part is, that the cofactor's factors were a p7, p18, and p77. It's a surprise that ECM missed a p7 and a p18.
Pari took 2.15 seconds on my machine to factor that number.

Last fiddled with by CRGreathouse on 2010-11-22 at 00:04
CRGreathouse is offline   Reply With Quote
Old 2010-11-22, 01:09   #1621
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

17×251 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
sigma(71608085354600097239391307^4) took an hour to factor via SIQS. The disappointing part is, that the cofactor's factors were a p7, p18, and p77. It's a surprise that ECM missed a p7 and a p18.
You'd have to do a very VERY small amount of ECM, or be impossibly unlucky, to miss a p7 and p18. I think most likely, this number was never ECM'd.
Mini-Geek is offline   Reply With Quote
Old 2010-11-22, 01:26   #1622
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

69016 Posts
Default

For me, it was the latter, for YAFU's ECM. PARI's ECM factored it with ease.

And, about QS; How does the comp determine the amount of relations that are necessary to factor a number?

An example would be the number 2802104782336128351749706435932730996816588058828391669. SIQS reports the following, "2312 relations needed".

I've recently taken on learning something about how it works. I know the very basic stuff about it, but not much else.

Last fiddled with by 3.14159 on 2010-11-22 at 01:32
3.14159 is offline   Reply With Quote
Old 2010-11-22, 01:37   #1623
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

102538 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
For me, it was the latter, for YAFU's ECM. PARI's ECM factored it with ease.
AFAIK, YAFU doesn't ECM at all before running SIQS if you tell it to run with the siqs() command. If you use the factor() command, then it runs ECM and should've found the factors.
Mini-Geek is offline   Reply With Quote
Old 2010-11-22, 02:03   #1624
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

24·3·5·7 Posts
Default

Quote:
Originally Posted by Mini-Geek View Post
AFAIK, YAFU doesn't ECM at all before running SIQS if you tell it to run with the siqs() command. If you use the factor() command, then it runs ECM and should've found the factors.
Indeed. siqs(n) goes directly to the quadratic sieve. In factor(n), SIQS is pretty much the "last resort" method, should trial division, Fermat's, Rho, P+1(?), and ECM all fail to factor the number completely. Though, it's still puzzling that it missed a p7 and a p18. I would have expected it to at least be present as a composite, the product of the p7 and p18.

(?) = "Does it use P+1"?

Last fiddled with by 3.14159 on 2010-11-22 at 02:06
3.14159 is offline   Reply With Quote
Old 2010-11-22, 02:41   #1625
bsquared
 
bsquared's Avatar
 
"Ben"
Feb 2007

2·3·587 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
Indeed. siqs(n) goes directly to the quadratic sieve. In factor(n), SIQS is pretty much the "last resort" method, should trial division, Fermat's, Rho, P+1(?), and ECM all fail to factor the number completely. Though, it's still puzzling that it missed a p7 and a p18. I would have expected it to at least be present as a composite, the product of the p7 and p18.
I'm not sure what you were expecting if you just ran SIQS. The only thing SIQS does prior to sieving is check for small factors (it does trial division up to 10k, IIRC) and perfect squares. After that, it will not find any factors until sieving, linalg, and sqrt are complete.

Quote:
Originally Posted by 3.14159 View Post

(?) = "Does it use P+1"?
yes, factor() uses P+1 with three different random bases. default B1= 20,000.
bsquared is offline   Reply With Quote
Old 2010-11-22, 02:43   #1626
bsquared
 
bsquared's Avatar
 
"Ben"
Feb 2007

352210 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
Pari took 2.15 seconds on my machine to factor that number.
I know this is the pari/gp thread, but here's a gratuitous yafu plug:

Code:
11/21/10 20:31:08 v1.20 @ AVA-343686, System/Build Info:
Using GMP-ECM 6.3, Powered by GMP 5.0.1
detected AMD Phenom(tm) II X4 955 Processor
detected L1 = 65536 bytes, L2 = 6291456 bytes, CL = 64 bytes
measured cpu frequency ~= 3209.987220

===============================================================
======= Welcome to YAFU (Yet Another Factoring Utility) =======
=======             bbuhrow@gmail.com                   =======
=======     Type help at any time, or quit to quit      =======
===============================================================
cached 78504 primes. pmax = 1000099


>> factoring 262934907404708576862758446841091336882604808443805210895581843776020374194863215112535
94913093697297001

div: primes less than 10000
fmt: 1000000 iterations
rho: x^2 + 1, starting 1000 iterations on C102
rho: x^2 + 3, starting 1000 iterations on C102
rho: x^2 + 2, starting 1000 iterations on C102
pp1: starting B1 = 20K, B2 = 1M on C102
pm1: starting B1 = 100K, B2 = 5M on C95
Total factoring time = 0.3290 seconds


***factors found***

P3 = 211
PRP7 = 2487811
PRP18 = 887712475607927971
PRP77 = 56425586867064181439411544449942718884745465011098050341509361224221109755611
looks like rho found the p7, and P-1 found the P18.

Last fiddled with by bsquared on 2010-11-22 at 02:44
bsquared is offline   Reply With Quote
Old 2010-11-22, 02:48   #1627
bsquared
 
bsquared's Avatar
 
"Ben"
Feb 2007

2·3·587 Posts
Default

errr, sorry. P+1 found the p7.

ran out of edit time...
bsquared is offline   Reply With Quote
Old 2010-11-22, 03:08   #1628
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

24×3×5×7 Posts
Default

Quote:
Originally Posted by bsquared View Post
errr, sorry. P+1 found the p7.

ran out of edit time...
... Isn't editing time one hour? It;s only been 23 minutes since you made the post.

I noted the source of the problem. Another difference between 32/64 bit programs.

Last fiddled with by 3.14159 on 2010-11-22 at 03:13
3.14159 is offline   Reply With Quote
Reply



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why do I sometimes see all the <> formatting commands when I quote or edit? cheesehead Forum Feedback 3 2013-05-25 12:56
Passing commands to PARI on Windows James Heinrich Software 2 2012-05-13 19:19
Ubiquity commands Mini-Geek Aliquot Sequences 1 2009-09-22 19:33
64-bit Pari? CRGreathouse Software 2 2009-03-13 04:22
Are these commands correct? jasong Linux 2 2007-10-18 23:40

All times are UTC. The time now is 23:06.


Fri Aug 6 23:07:00 UTC 2021 up 14 days, 17:35, 1 user, load averages: 4.19, 3.93, 3.92

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.