20120102, 11:04  #45  
Mar 2006
Germany
2·3·5^{2}·19 Posts 
Quote:
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 20120102 at 11:06 Reason: 11386 added 

20120102, 13:54  #46 
Apr 2010
Over the rainbow
2·1,217 Posts 
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,n1 do by for i=1,math.floor(math.sqrt(n)) do Last fiddled with by firejuggler on 20120102 at 14:10 
20120102, 14:07  #47 
Mar 2006
Germany
2·3·5^{2}·19 Posts 
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 20120102 at 15:04 
20120102, 16:38  #48 
Mar 2006
Germany
2·3·5^{2}·19 Posts 
I've played a bit with the source and here're some new things:
 Calculation of SumOfDivisors changed  the primeendings 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. 
20120103, 10:32  #49 
"Frank <^>"
Dec 2004
CDP Janesville
2×1,049 Posts 
Wolfgang has posted scans online of some beautiful handdrawn 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 20120103 at 10:34 Reason: keyboard made a typo.... 
20120103, 11:26  #50 
Apr 2010
Over the rainbow
4602_{8} Posts 
took care of the issue of kar bon
just add a line Code:
if sudiv<range and j<range then break end 
20120103, 12:10  #51 
Nov 2011
Saint Maur, France
2×5^{2} Posts 

20120103, 12:47  #52 
Mar 2006
Germany
2×3×5^{2}×19 Posts 
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) 
20120103, 13:37  #53 
Apr 2010
Over the rainbow
100110000010_{2} Posts 
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 
20120103, 13:42  #54 
Mar 2006
Germany
2×3×5^{2}×19 Posts 

20120103, 14:59  #55 
Apr 2010
Over the rainbow
2434_{10} Posts 
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. 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
aliquot escape  firejuggler  Aliquot Sequences  26  20120119 08:15 
Special project #3b  Project 400  schickel  Aliquot Sequences  307  20111028 01:29 
New σ for Aliquot  JohnFullspeed  Aliquot Sequences  18  20110820 21:11 
Special project #3a  Project 300  schickel  Aliquot Sequences  29  20110812 17:45 
serious bug in aliquot.ub  Andi47  Aliquot Sequences  3  20090308 10:18 