mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2012-01-02, 11:04   #45
kar_bon
 
kar_bon's Avatar
 
Mar 2006
Germany

2·3·11·43 Posts
Default

Quote:
Originally Posted by firejuggler View Post
Would you detect any error, please say so,it is all done manually.
Here're some:

ends in '43': append to "1114" this: 1114->560->928->962->634

doubled:
483168->29351
474750->82
765380->82
11386->5696

Last fiddled with by kar_bon on 2012-01-02 at 11:06 Reason: 11386 added
kar_bon is offline   Reply With Quote
Old 2012-01-02, 13:54   #46
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

44758 Posts
Default

here is a lua script that will get it into form.

it detect 2 long amicable cycle, perfect number, "one away' merge
find any lua interpreter and it will create a "graph.txt" wherever the script is.

you can control what you want :=>

aliquot sequence from 3 to 50001, with iteration value up to 10000000 take about 1 or 2 minute to write and take 1800ko
for sequences up to 500001, with iteration value up to 100000000 ,about 5minutes and 14800 ko


adit-> speeding the script , replace

for i=1,n-1 do
by
for i=1,math.floor(math.sqrt(n)) do
Attached Files
File Type: txt lua-aliquot.txt (3.3 KB, 71 views)

Last fiddled with by firejuggler on 2012-01-02 at 14:10
firejuggler is online now   Reply With Quote
Old 2012-01-02, 14:07   #47
kar_bon
 
kar_bon's Avatar
 
Mar 2006
Germany

1011000101102 Posts
Default

Quote:
Originally Posted by firejuggler View Post
here is a lua script that will get it into form.
Not yet run it but line 29 should be 'pres3="n"' not 'pre3="n"', I think.

PS: Run it, but...

1. 50000 is too high: dotty nor gvedit will display them (no response for a long time)
2. The calculations are not correct: try to n=100
The result will give: 80->20->7 but should like this: 80->106->56->64->63->41

Last fiddled with by kar_bon on 2012-01-02 at 15:04
kar_bon is offline   Reply With Quote
Old 2012-01-02, 16:38   #48
kar_bon
 
kar_bon's Avatar
 
Mar 2006
Germany

2·3·11·43 Posts
Default

I've played a bit with the source and here're some new things:

- Calculation of SumOfDivisors changed
- the prime-endings were not correct: two different displays for example "41" and "41p"
-> changed printing node: "41" with "label=41p"
- function "aliqout" with parameter

Issue to solve:
- some edges/arrows are doubled

Example source is set for n=1000, so dotty will display it.
Attached Files
File Type: txt lua-aliquot.txt (3.4 KB, 189 views)
kar_bon is offline   Reply With Quote
Old 2012-01-03, 10:32   #49
schickel
 
schickel's Avatar
 
"Frank <^>"
Dec 2004
CDP Janesville

2×1,061 Posts
Default

Wolfgang has posted scans online of some beautiful hand-drawn graphs he did of aliquot trees in the 90s. There is one drawing linked on the English side of his site and 4 more linked on the German side.....

Last fiddled with by schickel on 2012-01-03 at 10:34 Reason: keyboard made a typo....
schickel is offline   Reply With Quote
Old 2012-01-03, 11:26   #50
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

5×11×43 Posts
Default

took care of the issue of kar bon
just add a line

Code:
if sudiv<range and j<range then break    end
now to take care of the 'circle' issue...
firejuggler is online now   Reply With Quote
Old 2012-01-03, 12:10   #51
MichelMarcus
 
Nov 2011
Saint Maur, France

2·52 Posts
Default

See http://www.strw.leidenuniv.nl/~mathar/public/mathar20100305.pdf
at http://www.strw.leidenuniv.nl/~mathar/public/index.html
MichelMarcus is offline   Reply With Quote
Old 2012-01-03, 12:47   #52
kar_bon
 
kar_bon's Avatar
 
Mar 2006
Germany

2·3·11·43 Posts
Default

How about this one:

- if sigma > range
-> calculate sigma until < range (or ends in prime) and count the steps
-> display the edge/arrow as dotted with number of steps as comment

Not all is done correctly but have no time for next hours.

As example see 126 (range = 166):

shown as "126" ......6......> 136 (so steps 126 -> [186 -> 198 -> 270 -> 450 -> 759 -> 393 ->] 135 (in "[]") are only counted)
Attached Files
File Type: txt lua-aliquot.txt (4.0 KB, 185 views)
kar_bon is offline   Reply With Quote
Old 2012-01-03, 13:37   #53
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

5×11×43 Posts
Default

hmm nice.
the case where it doesn't come back under the range isn't handled but is easilly handled with

Code:
    while j>range do
                    tmp=j            -- hold j
                    j=divisorSum(j)
                    pres4 = pres4 + 1    -- one step more
                    if j==1 then        -- if ends in prime (sigma=1) -> set back hold value
                        j=tmp
                        io.write(j.."[label=\""..j.."p\",shape=plaintext];\n")
                        break
                    end
                    if j>1000000 then break end
                end
now to resolve another problem... why does some even value are 'isolated'?
firejuggler is online now   Reply With Quote
Old 2012-01-03, 13:42   #54
kar_bon
 
kar_bon's Avatar
 
Mar 2006
Germany

2·3·11·43 Posts
Default

Quote:
Originally Posted by firejuggler View Post
now to resolve another problem... why does some even value are 'isolated'?
Because the sequence is bigger than this:

Code:
   if sudiv>10000000 then
      break
   end
Look for example n=138.

Until later.
kar_bon is offline   Reply With Quote
Old 2012-01-03, 14:59   #55
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

5×11×43 Posts
Default

yup i know about those, but what about 104,102,96,90,88,84,80,78,72,66,60,54,48,42,36,30,28,24,18?
those have no arrow going to them of going from them.
firejuggler is online now   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
aliquot escape firejuggler Aliquot Sequences 26 2012-01-19 08:15
Special project #3b - Project 400 schickel Aliquot Sequences 307 2011-10-28 01:29
New σ for Aliquot JohnFullspeed Aliquot Sequences 18 2011-08-20 21:11
Special project #3a - Project 300 schickel Aliquot Sequences 29 2011-08-12 17:45
serious bug in aliquot.ub Andi47 Aliquot Sequences 3 2009-03-08 10:18

All times are UTC. The time now is 20:59.

Mon Aug 3 20:59:58 UTC 2020 up 17 days, 16:46, 0 users, load averages: 1.37, 1.41, 1.44

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, 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.