Aliquot genealogy project

EdH · 2011-12-30, 05:34

Is there a "family tree" of sorts that shows all the merges for sequences under 1M in the db?

schickel · 2011-12-30, 06:56

Quote:

Originally Posted by EdH

Is there a "family tree" of sorts that shows all the merges for sequences under 1M in the db?

No. The closest thing would be Wolfgang's C30C90 table. If you look in there, you will find a list of open sequences, with the first c9 and first c30 encountered during their lifetime. Some merges and terminations are listed, but at the top it says:

Quote:

Originally Posted by Wolfgang

All numbers you will not find in the following table are key-numbers of side-sequences. Have a look into the matrix of ALIQUOT. You will find them there.

Last updated 2009, so not real current.

The matrix listed is one that is maintained by a TurboPascal program that is listed in the download section. Said matrix is a listing of sequence starting numbers and a status which is indicated by a number: 0 for open, the prime a sequence terminates at, or the member of an aliquot cycle that a sequence merges with. Not so useful if you want to fiind out what merges with what, from either side, but you could get a count of how many sequences terminate at which prime. Major drawback: uses "1" for the starting line of a sequence file, rather than "0"....another drawback: the record file has not been updated since 1994.

science_man_88 · 2011-12-30, 12:54

Quote:

Originally Posted by schickel

No. The closest thing would be Wolfgang's C30C90 table. If you look in there, you will find a list of open sequences, with the first c9 and first c30 encountered during their lifetime. Some merges and terminations are listed, but at the top it says:Last updated 2009, so not real current.

The matrix listed is one that is maintained by a TurboPascal program that is listed in the download section. Said matrix is a listing of sequence starting numbers and a status which is indicated by a number: 0 for open, the prime a sequence terminates at, or the member of an aliquot cycle that a sequence merges with. Not so useful if you want to fiind out what merges with what, from either side, but you could get a count of how many sequences terminate at which prime. Major drawback: uses "1" for the starting line of a sequence file, rather than "0"....another drawback: the record file has not been updated since 1994.

yeah I could play with my ali and aligen scripts ( they work backwards) again but I doubt anything interesting comes out of it.

EdH · 2011-12-30, 14:44

@schickel: Thanks. I looked at the c30c90 file a while ago, but not the other reference, yet. You have posted an "AllSeq.zip" file in the past, so there must be a method to query just the last line from the db for all the sequences. I can get full .elf files, but all my last line attempts end up downloading an .html page. How can I get just the last lines, either individually or all at once (preferred)? For that matter, is there a reference that shows what queries are available for the db?

If I can get all the last lines, I'm pretty sure I can write something that can build a merge listing. Of course, then I need to figure out the best method for displaying merges...

schickel · 2011-12-30, 19:59

Quote:

Originally Posted by EdH

@schickel: Thanks. I looked at the c30c90 file a while ago, but not the other reference, yet.

Ah, OK. I wasn't sure how much of the background info you were aware of, since it has gotten harder to keep track of who joined when and who knows how much....

Quote:

You have posted an "AllSeq.zip" file in the past, so there must be a method to query just the last line from the db for all the sequences. I can get full .elf files, but all my last line attempts end up downloading an .html page.

I use a 5-lb sledgehammer to do the job. I have a list of the open sequences and just download the .elf file for each one.

I don't do the whole range but once a week or so, since it is so query intensive. I keep meaning to sit down and do a "last line" tool Real Soon Now.

Quote:

How can I get just the last lines, either individually or all at once (preferred)?

Unfortunately, if you use the "show last element" radio button, the last line is shown as an html page. The last line is there, but you have to parse it out of the html....

Quote:

For that matter, is there a reference that shows what queries are available for the db?

Not really. Some of the queries available were discussed over in the FactorDb thread, but I would say that the only one who knows what all is available is Syd....

Quote:

If I can get all the last lines, I'm pretty sure I can write something that can build a merge listing. Of course, then I need to figure out the best method for displaying merges...

I have been thinking along those lines, too. So far it seems to me there is no really easy way to do things, since if you want to make it comprehensive, you almost have to have 1.000.000 entries (less the primes since they're easy), even if most of them are "S_n merges with x after y lines" or "S_nterminates in x after y lines" with y<10.....

science_man_88 · 2011-12-30, 20:16

Quote:

Originally Posted by schickel

even if most of them are "S_n merges with x after y lines" or "S_nterminates in x after y lines" with y<10.....

this good enough for you:

Code:

aligen5(w,s)=for(z=w,s,for(x=1,#ali(z),print1(ali(z)[x]":i2=");for(y=1,#ali(ali(z)[x]),print1(ali(ali(z)[x])[y]":i3=")));print())

? it was a slight change to my aligen script.

EdH · 2011-12-31, 05:08

Quote:

Originally Posted by schickel

Ah, OK. I wasn't sure how much of the background info you were aware of, since it has gotten harder to keep track of who joined when and who knows how much....

Wow, I guess I've been around for two years now, according to my stats... But I should still be considered as someone who doesn't know much.

Quote:

Originally Posted by schickel

I use a 5-lb sledgehammer to do the job. I have a list of the open sequences and just download the .elf file for each one.

The more I look at this, the larger a task it seems. To download 1M .elfs, or even just last lines, is a huge endeavor, itself. Even If I d/l 1k per day, it would take 3 years to complete...

Quote:

Originally Posted by schickel

I don't do the whole range but once a week or so, since it is so query intensive. I keep meaning to sit down and do a "last line" tool Real Soon Now.

I'm going to have to think a bit more about this. Things like maybe I should eliminate those sequences that end in primes or cycles.

Quote:

Originally Posted by schickel

Unfortunately, if you use the "show last element" radio button, the last line is shown as an html page. The last line is there, but you have to parse it out of the html....

It would be easier to parse the .elf files, although they would be larger.

Quote:

Originally Posted by schickel

Not really. Some of the queries available were discussed over in the FactorDb thread, but I would say that the only one who knows what all is available is Syd....

I kind of wondered if there was a .php method that queries available functions.

Quote:

Originally Posted by schickel

I have been thinking along those lines, too. So far it seems to me there is no really easy way to do things, since if you want to make it comprehensive, you almost have to have 1.000.000 entries (less the primes since they're easy), even if most of them are "S_n merges with x after y lines" or "S_nterminates in x after y lines" with y<10.....

Yeah, I was thinking along the lines of a list:

Code:

seq1:line = seq12:line
seq1:line = seq23:line
seq34:line = seq44:line
...

But, let's address two sequences that end in the same cycle. Did they merge?

This is seeming more complex, the further I dig...

schickel · 2011-12-31, 06:23

Quote:

Originally Posted by EdH

Wow, I guess I've been around for two years now, according to my stats... But I should still be considered as someone who doesn't know much.

The more I look at this, the larger a task it seems. To download 1M .elfs, or even just last lines, is a huge endeavor, itself. Even If I d/l 1k per day, it would take 3 years to complete...

Yeah, but the trick is, I just do the open ones, and there are only ~~9244~~ 9243 open ones. I have it broken into 3 parts because I have a memory leak somewhere that kills the task after ~6000 downloads. It takes about 4 hours to do each 1/3 of the list. (Most of the time is spent reformatting the files a little; if you pad the line numbers out to 10 spaces, you can peg the size of a line easier.)

Quote:

I'm going to have to think a bit more about this. Things like maybe I should eliminate those sequences that end in primes or cycles.

Eliminate from the download? Yes, you can, since they're done once they merge....

Quote:

It would be easier to parse the .elf files, although they would be larger.

That's true, though the download time would be lots less just for the last line.

Quote:

I kind of wondered if there was a .php method that queries available functions.

Nothing was ever brought up about that. It looked like Syd was adding stuff as people brought things up.

Quote:

Yeah, I was thinking along the lines of a list:

Code:

seq1:line = seq12:line
seq1:line = seq23:line
seq34:line = seq44:line
...

But, let's address two sequences that end in the same cycle. Did they merge?

I would look at where they merged. If they merge before the resulting confluence hits the cycle, I would say that was a merge, since there was some activity before the merge.

This actually came up on one of our terminations. I marked it as a merge because it hit a sequence and headed up before the last downdriver run to the termination. Wolfgang marked it down as a temination since it actually terminated rather than ran into another open sequence.

Quote:

This is seeming more complex, the further I dig...

Welcome to our world....

------------------------

The ultimate solution is to draw a tree where each number is a node with an arrow pointing to its $\sigma(n)-n$ value and has arrows pointing from numbers that point to it.

After this post, I'll look back; someone was here wondering what software to use to visualize aliquot cycles in this manner, and I think they posted a paper, with some of the aliquot trees in it. Unfortunately, I forgot where I found it.....:face palm:

schickel · 2011-12-31, 06:26

Here's the thread I was thinking of. Geez.....it was only in April that this came up.

Now I'll go and see if I can find that paper.

schickel · 2011-12-31, 06:53

Aha! Here we go. Richard Mathar at Leiden University; here is an index of his papers, and here's a deep link to the paper I was thinking of....

[Edit: I wonder if this is the first citation of the FactoDB in a paper.....]

schickel · 2011-12-31, 08:19

Quote:

Originally Posted by schickel

After this post, I'll look back; someone was here wondering what software to use to visualize aliquot cycles in this manner, and I think they posted a paper, with some of the aliquot trees in it. Unfortunately, I forgot where I found it.....:face palm:

OK, so it looks like it wasn't the same person....I wonder if Wini left a working email on file....

2011-12-30, 05:34	#1
EdH "Ed Hall" Dec 2009 Adirondack Mtns 2³×677 Posts	Aliquot genealogy project Is there a "family tree" of sorts that shows all the merges for sequences under 1M in the db? Last fiddled with by schickel on 2011-12-31 at 18:29 Reason: Change thread title (I myspelled) [& change icon]

2011-12-31, 06:53	#10
schickel "Frank <^>" Dec 2004 CDP Janesville 2·1,061 Posts	Aha! Here we go. Richard Mathar at Leiden University; here is an index of his papers, and here's a deep link to the paper I was thinking of.... [Edit: I wonder if this is the first citation of the FactoDB in a paper.....] Last fiddled with by schickel on 2011-12-31 at 06:56 Reason: Adding edit

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

2011-12-30, 14:44	#4
EdH "Ed Hall" Dec 2009 Adirondack Mtns 2³×677 Posts	@schickel: Thanks. I looked at the c30c90 file a while ago, but not the other reference, yet. You have posted an "AllSeq.zip" file in the past, so there must be a method to query just the last line from the db for all the sequences. I can get full .elf files, but all my last line attempts end up downloading an .html page. How can I get just the last lines, either individually or all at once (preferred)? For that matter, is there a reference that shows what queries are available for the db? If I can get all the last lines, I'm pretty sure I can write something that can build a merge listing. Of course, then I need to figure out the best method for displaying merges...

2011-12-31, 06:26	#9
schickel "Frank <^>" Dec 2004 CDP Janesville 2·1,061 Posts	Here's the thread I was thinking of. Geez.....it was only in April that this came up. Now I'll go and see if I can find that paper.