mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2011-12-31, 22:28   #23
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

202618 Posts
Default

Quote:
Originally Posted by EdH View Post
Now, that, I like!

You guys are coming up with some good stuff.

The clickable nodes looks good as does the familyecho.com, but, yes, it would entail a whole lot of work. Of course, an automated version is what is needed.

Here's what I envision - let me know what you guys think:

1. create a program to query the db for all merges - this would be a long term project
2. have the program write the graphviz code - let's try to be less manually intensive**
3. periodically run the program against all OE sequences and update the graphviz file as needed
4. or, have the program provide whatever data is needed to update familyecho or some other choice of tree type

**@firejuggler: did you automate the file or manually enter the data?
Quote:
digraph G {
graph [URL="default.html"]
a [URL="a.html"]
b [URL="b.html"]
c [URL="c.html"]

a -> b -> c
a -> c
}
the bold looks easiest to write but Pari's been my only thing for so long, (so basically no useful language) the rest depend on searching and recording it. I think the way we'd detect it by site normally is:

1) check the sequence lines above i0 for sequences to check
2) check those sequences with a last line search
3)check if it's prime or a cycle if so the value is a termination not a merge ?

the problem is I still don't understand what to do from there to write the -> part because the example they give has multiple pathways.
science_man_88 is offline   Reply With Quote
Old 2011-12-31, 23:13   #24
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

2×1,217 Posts
Default

updated my file, there should be no mistake anymore. (added cycle, A->A, added amicable , A->B->A)
Attached Files
File Type: txt graph1.txt (5.7 KB, 98 views)
firejuggler is online now   Reply With Quote
Old 2012-01-01, 01:11   #25
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

8,369 Posts
Default

Code:
aligen7(w,s)->print("digraph G {");forprime(z=w,s,print(eval(z)""eval(Str(p))"[shape=plaintext];");for(x=1,#ali(z),print(ali(z)[x]"[shape=plaintext];");print(ali(z)[x]"->"Str(z"p"));for(y=1,#ali(ali(z)[x]),print(ali(ali(z)[x])[y]"[shape=plaintext];");print(ali(ali(z)[x])[y]"->"ali(z)[x]))));print("}")
is the best I can do without total modification of the aligen script ( and yes I know it's not in order, I guess I could break it apart a bit and do a forprime loop before to get the first few lines in order.)

Last fiddled with by science_man_88 on 2012-01-01 at 01:12
science_man_88 is offline   Reply With Quote
Old 2012-01-01, 05:47   #26
firejuggler
 
firejuggler's Avatar
 
Apr 2010
Over the rainbow

243410 Posts
Default

another update, done manually, again...
added a few more, the WIP are not done yet (still have to look for the preceding term)
It seems that there is a lot (as observed for now) of sequence that end in 43
naming convention :
prime are named "Xp" in plain text
odd are in plain text
even are default
Attached Files
File Type: txt graph1.txt (14.0 KB, 113 views)
firejuggler is online now   Reply With Quote
Old 2012-01-01, 08:28   #27
schickel
 
schickel's Avatar
 
"Frank <^>"
Dec 2004
CDP Janesville

2×1,049 Posts
Default

Quote:
Originally Posted by firejuggler View Post
another update, done manually, again...
added a few more, the WIP are not done yet (still have to look for the preceding term)
It seems that there is a lot (as observed for now) of sequence that end in 43
naming convention :
prime are named "Xp" in plain text
odd are in plain text
even are default
Wow! that's starting to get intricate. It's nice how it does draw the cycles. Reading the documentation, it said you had to change some of the internal edges on the cycles, since DOT won't render them, but I guess it does some magic to make it work.

I found a couple of duplicated links in your structure: 25->45, 455->217, 329->55, 185->82, and 377->82. Nothing major, it just results in a double-arrow.

It's taking a noticable time to reder now: 30 secs or so. Also, the resulting picutre as a PNG is 16560x3227 pixels and weighs in at 1.5MB....geez, and the jpg is a whopping 2.1MB.

Guess we'll have to use the same people that print the Mersenne Prime posters to print it out.....
schickel is offline   Reply With Quote
Old 2012-01-01, 14:42   #28
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

10AB16 Posts
Default

Quote:
Originally Posted by schickel View Post
It's taking a noticable time to reder now: 30 secs or so. Also, the resulting picutre as a PNG is 16560x3227 pixels and weighs in at 1.5MB....geez, and the jpg is a whopping 2.1MB.
I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels)
Quote:
Originally Posted by schickel View Post
Guess we'll have to use the same people that print the Mersenne Prime posters to print it out.....
At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
Attached Files
File Type: zip graph1.zip (70.8 KB, 84 views)

Last fiddled with by Mini-Geek on 2012-01-01 at 14:52
Mini-Geek is online now   Reply With Quote
Old 2012-01-01, 14:57   #29
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

8,369 Posts
Default

Quote:
Originally Posted by Mini-Geek View Post
I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels)

At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
doesn't look very clickable but still a lot of relations. does this count :


Code:
(23:55)>ali(923774)
%2 = [853356555529]

Last fiddled with by science_man_88 on 2012-01-01 at 14:59
science_man_88 is offline   Reply With Quote
Old 2012-01-01, 15:35   #30
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

100000101100012 Posts
Default

Quote:
Originally Posted by science_man_88 View Post
doesn't look very clickable but still a lot of relations. does this count :


Code:
(23:55)>ali(923774)
%2 = [853356555529]
made a new script to check for and only report under 1000000.

Code:
 (n)->a=n-1;b=[];for(x=0,floor(.5*a),if(x==0,if(sigma(a^2)-a^2==n,b=concat(b,a^2)),if(sigma((x)*(a-x))-((x)*(a-x))==n,b=concat(b,(x)*(a-x)))));b=vecsort(b,,8);for(g=1,#b,if(b[g]<1000000,b[g]=b[g],b[g]=0));
b=vecsort(b,,8)
science_man_88 is offline   Reply With Quote
Old 2012-01-01, 15:43   #31
schickel
 
schickel's Avatar
 
"Frank <^>"
Dec 2004
CDP Janesville

2·1,049 Posts
Default

Quote:
Originally Posted by Mini-Geek View Post
I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels)

At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
I was using gvedit which gets installed with the Windows package. It was helpful getting started having a GUI frontend. I assume that it's doing some post processing internally.

I hadn't considered SVG, since I'm still on IE8 and it needs a plugin to render SVG files....good to know that a portable format keeps the size down.

As far as printing, I was anticipating that as more gets added, the size is going to start to scale up, but the size info on the current set is encouraging....
schickel is offline   Reply With Quote
Old 2012-01-01, 18:44   #32
henryzz
Just call me Henry
 
henryzz's Avatar
 
"David"
Sep 2007
Cambridge (GMT/BST)

22·1,433 Posts
Default

Nice graph firejuggler. If we can get this automated this would be brilliant.
Currently there is a mistake around 341/343.
henryzz is online now   Reply With Quote
Old 2012-01-01, 19:05   #33
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

836910 Posts
Default

Quote:
Originally Posted by henryzz View Post
Nice graph firejuggler. If we can get this automated this would be brilliant.
Currently there is a mistake around 341/343.
if I read it correct I see multiple errors concerning 341:

Quote:
... 1001->341 ...
Code:
(15:02)>sigma(1001)-1001
%22 = 343

Quote:
... 341->57 ...
Code:
(15:02)>sigma(341)-341
%23 = 43
Quote:
341->82
Code:
(15:02)>sigma(341)-341
%23 = 43
science_man_88 is offline   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 12:17.

Fri Oct 23 12:17:11 UTC 2020 up 43 days, 9:28, 0 users, load averages: 1.32, 1.39, 1.38

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.