MathGroup Archive 2001

[Date Index] [Thread Index] [Author Index]

Search the Archive

Two factors of (10^71-1)/9 = R71

  • To: mathgroup at smc.vnet.net
  • Subject: [mg30125] Two factors of (10^71-1)/9 = R71
  • From: seidovzf at yahoo.com (Zakir F. Seidov)
  • Date: Sun, 29 Jul 2001 21:26:17 -0400 (EDT)
  • Organization: The Math Forum
  • References: <6kajko$2bn@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

The number William Mopppett  wrote:
           ================ 
(10^71 - 1)/9 
is R71, repunit, number with "all ones".

And it has two factors:
R71 = 
241573142393627673576957439049
*45994811347886846310221728895223034301839

See, e.g.,  
http://www.ping.be/~ping6758/repunits.htm

But suprisingly enough,
when,  before looking for "repunit"s in 37.com,
I asked my PC to work for saturday, 
its Mathematica session was:

 $Version
 4.0 for Microsoft Windows (December 5, 1999)

<< NumberTheory`FactorIntegerECM`

Timing[f = FactorIntegerECM[(2^128 + 1)/f] ]
{793.95 Second, 59649589127497217} (* good, but...*)

Timing[FactorIntegerECM[(10^71 - 1)/9 ]]
{112642. Second, \
111111111111111111111111111111111111111111111111111111111111111\
11111111} (* ??!! *)

That is Mathematica (or my PC?) considers R71 as prime!

And this is a real suprise 
and hence the question to Mathematica experts:
how it can be?

Zakir F. Seidov
seidovzf at yahoo.com

%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%
Subject: [mg30125]      RE: factor
Author:       Ersek_Ted%PAX1A at mr.nawcad.navy.mil
Organization: Steven M. Christensen and Associates, Inc and
MathTensor, Inc.


William Mopppett  wrote:
|
|Can anyone give me a factor of (10^71 - 1)/9 Mathematica 3.0 says it
is
|not prime email
|wmoppet at nsw.bigpond.net.au
|
|

  The attempt below is probably the most effective way to do this. My
90
Mhz Pentium worked on it for 14 hours and still didn't find a factor.

The documentation makes it clear you will may not get an answer in a
short  amount of time.  It also says it is designed to find factors up
to about 18  digits in 3 hours on a "workstation" (what ever that is).

Well it seems all the prime factors of this number have 18 to 35
digits,
and  it has two, three, or four prime factors.  The number of
potential
prime  factors to consider is roughly,
(  10^35/Log[10^35] - 10^18/Log[10^18] ) = ( 1.2 * 10^33  )     !!


In[1]:=
<<NumberTheory`FactorIntegerECM`

In[2]:=
FactorIntegerECM[(10^71 - 1)/9 ]

(*  still waiting  *)
___________________________________
Perhaps you can do better by specifying a value for one or more of the

options.

In[3]:=
Options[FactorIntegerECM]

Out[3]=
{FactorSize -> Automatic,
  CurveNumber -> Automatic,
  CurveCountLimit -> 10000}

__________________________________

Ted Ersek



  • Prev by Date: Great Michael Trott & six-fold integral
  • Next by Date: Re: ListCorrelate[] ??
  • Previous by thread: Great Michael Trott & six-fold integral
  • Next by thread: Re: Two factors of (10^71-1)/9 = R71