mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2021-05-13, 14:31   #56
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

22×7×137 Posts
Default

Quote:
Originally Posted by garambois View Post
Thanks Edwin for all your hard work !
I have some questions, remarks, answers, information...

1) Would it be possible to show us an example of an entry to achieve an "Advanced Search" ?
Using your previous example:
Code:
[n%2==0, a==0, b, c, d, e==0, f, g, h, i, j, 1.7<k<2.3, l, m]
and translating for the change of elements:
Code:
$ ./seqinfo2 regina_file_1e5 
Data available for sequences 2 through 100000 
Sequence endings - prime: 85544, cycle: 2225, open: 12230 
Enter sequence (##/a/h/p/p##/q/u): a 
Advanced routine 
Sequence filter a (mod b) (a b): 2 
Sequence range (begin end):  
All, Prime, Open, Cylce (a/p/o/c): o 
Range for number of relative minimums (min max):  
Range for number of relative maximums (min max): 0 
Range for number of parity changes (min max):  
Range for number of consecutive even abundant terms (min max):  
Range for number of consecutive even deficient terms (min max):  
Range for number of consecutive odd deficient terms (min max):  
Range for number of consecutive odd abundant terms (min max):  
Range for number of downdriver extracts (min max):  
Range for smallest quotient of two consecutive terms (min max):  
Range for largest quotient of two consecutive terms (min max):  
Range for arith mean of all quotients of 2 consequtive terms (min max): 1.7 2.3 
Range for geo mean of all quotients of 2 consequtive terms (min max):  
Range for arith mean of # of digits of all the mins (min max):  
Range for arith mean of # of digits of all the maxes (min max):  
List filtered sequences for 0 (mod 2) from 2 through 100000? (y/n/c/f): y
. . .
17360     17544     17994     18006     18018     18138     18150     18490
18702     18786     18798     18888     19040     19410     19500     19560
19656     20204     20382     20640     21096     21198     21210     21234
21246     21258     21384     21582     21858     21870     22132     22392     
. . .
842 sequences found. 
 Enter sequence (##/a/h/p/p##/q/u):
Of course, this was only with a regina_file that was limited to 1e5 and I showed only 32 of 842 sequences, to include 19560 which you referenced earlier.

I entered "a" to call up the Advanced Routine. Then "2" for 0 (mod 2), "o" for open sequences, "0" for maximums, which became the range 0-0, and "1.7 2.3" for the arithmetic mean range. In this case I was inclusive of the values, whereas the original wasn't. I suppose "1.700001 2.29999" could have been used.
Quote:
Originally Posted by garambois View Post
2) After less than 3 weeks of computation, the new program came to more than 4,000,000. So I think it should reach 20 M sometime this summer. But then, the RAM needed to run the program should become too large for almost all computers, unless I am wrong ! For sequences up to 20 M, the file will be over 2 GB. And in python, in my experience, it is impossible to process such a file by storing it in RAM : it takes much more memory than the file size. But maybe it's different in C, I hope so !

3) Just for information, I am thinking of letting the program run for several months and this new program could process all the sequences up to 50 M, or even more... if my old computer does not break down.
I've placed a limit of 20M in the code at present, but I have no idea what the memory usage would actually be. Or, what happens when there's not enough RAM. Theoretically, I think it should move into swap space, but then what?
Quote:
Originally Posted by garambois View Post
4)

There can be no other negative value than B in the table, that is for sure.
Great! Thanks.
EdH is offline   Reply With Quote
Old 2021-05-14, 08:32   #57
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

647 Posts
Default

Edwin, this new regina file analysis tool is very powerful !
Thank you very much for all this work !

I still have some suggestions to improve it a bit, if in the next times you work on this program again :

1) In the advanced search, would it be possible to add a query "Range for number of iterations (min max) :" ?

2) Unless I'm mistaken, there are only two types of searches I can't do with this wonderful program (The equivalent of ##p, but for open end sequences and for sequences that end in cycles) :
I) Search for the number or list of sequences that merge with, for example, 276 or another open end sequence.
II) Search for the number or list of sequences that arrive on a perfect number, and on another cycle. This is much more difficult, even impossible, because to enter for example a cycle of length 4, there are 4 possible entry doors and there are many different cycles 4 !

Maybe I will find other things to add to the program when I start to analyze the data very seriously in a few months, when the regina file will have 20 M sequences or more !
But for the moment, this is all I see...

Congratulations, once again for the ease of use of this program and for the possibilities it offers !


garambois is online now   Reply With Quote
Old 2021-05-14, 13:37   #58
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

22×7×137 Posts
Default

Quote:
Originally Posted by garambois View Post
Edwin, this new regina file analysis tool is very powerful !
Thank you very much for all this work !

I still have some suggestions to improve it a bit, if in the next times you work on this program again :

1) In the advanced search, would it be possible to add a query "Range for number of iterations (min max) :" ?

2) Unless I'm mistaken, there are only two types of searches I can't do with this wonderful program (The equivalent of ##p, but for open end sequences and for sequences that end in cycles) :
I) Search for the number or list of sequences that merge with, for example, 276 or another open end sequence.
II) Search for the number or list of sequences that arrive on a perfect number, and on another cycle. This is much more difficult, even impossible, because to enter for example a cycle of length 4, there are 4 possible entry doors and there are many different cycles 4 !

Maybe I will find other things to add to the program when I start to analyze the data very seriously in a few months, when the regina file will have 20 M sequences or more !
But for the moment, this is all I see...

Congratulations, once again for the ease of use of this program and for the possibilities it offers !


Those that merge with other sequences are shown in the basic search, but I'm guessing you would like to also filter those (and the others) through the Advanced Function?:
Code:
Enter sequence (##/a/h/p/p##/q/u): 70508 
70508 is open ended. It merges with 276. 
Enter sequence (##/a/h/p/p##/q/u): 276 
276 is open ended. List any sequences that merge with 276? (y/n/c/f): y 
. . . 
77764     78718     79678     82604     82810     84094     85772     86158 
87566     88210     90010     90772     91684     92158     93798     93932 
94718     96098     96878     96912     97534     98008     99058      
167 sequences found. 
Enter sequence (##/a/h/p/p##/q/u): 3142 
3142 ends with the perfect number 496, 
Display all sequences that end with 496? (y/n/c/f): y 
496       608       650       652       790       1294      1574      1778 
2162      2582      3142      5158       
12 sequences found.
To perform cycle searches, you need to have wget and an Internet connection to reach factordb:
Code:
Enter sequence (##/a/h/p/p##/q/u): 83076 
83076 ends in a cycle. Display cycle? (starts at entry point) (y/n/f): y 
376736 
381028 
285778 
152990 
122410 
97946 
48976 
45946 
22976 
22744 
19916 
17716 
14316 
19116 
31704 
47616 
83328 
177792 
295488 
629072 
589786 
294896 
358336 
418904 
366556 
274924 
275444 
243760 
Display all sequences that end in this cycle? (y/n/c/f): y 
. . . 
66312     67078     67128     67894     69542     70846     71612     75536 
77506     77576     81898     82768     83076     83114     83328     84000 
85438     88692     91886     94214     97946      
85 sequences found.
But, again, this does not include the Advance filtering, "yet."

Last fiddled with by EdH on 2021-05-14 at 13:42 Reason: fixt spiiln errrer
EdH is offline   Reply With Quote
Old 2021-05-14, 16:07   #59
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

12078 Posts
Default

Quote:
Originally Posted by EdH View Post
but I'm guessing you would like to also filter those (and the others) through the Advanced Function ?

YES !

Yes, but no problem if it's too complicated, especially when it comes to sequences that arrive into a cycle !
The possibilities of the program are really already fabulous now...

And as said in my previous message, I think it would be very useful to add also the filter on the number of iterations of the sequence "Range for number of iterations (min max) :".
garambois is online now   Reply With Quote
Old 2021-05-14, 16:20   #60
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

22×7×137 Posts
Default

Quote:
Originally Posted by garambois View Post
. . .
And as said in my previous message, I think it would be very useful to add also the filter on the number of iterations of the sequence "Range for number of iterations (min max) :".
How strange. . . I added that and thought I posted it. I must have skipped hitting "Submit."

Anyway, the basics were in the earlier source. A little copy/paste/editing and voilĂ . . .
Attached Files
File Type: cpp seqinfo2.cpp (36.8 KB, 12 views)
EdH is offline   Reply With Quote
Old 2021-05-14, 16:48   #61
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

647 Posts
Default

Quote:
Originally Posted by EdH View Post
How strange. . . I added that and thought I posted it. I must have skipped hitting "Submit."

Anyway, the basics were in the earlier source. A little copy/paste/editing and voilĂ . . .

OK, many thanks.
I will do some more testing in the next few days with a larger regina file.
I'll let you know if I run into any problems, but everything should work fine.
garambois is online now   Reply With Quote
Old 2021-05-14, 17:56   #62
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

22·7·137 Posts
Default

Quote:
Originally Posted by garambois View Post
OK, many thanks.
I will do some more testing in the next few days with a larger regina file.
I'll let you know if I run into any problems, but everything should work fine.
Sounds Good! Could you take a look at how much RAM is being used? If you run System Monitor, the Processes tab should show you the resting use of memory for seqinfo2. For the _1e5 file I have, my system is using 1.5 GB.
EdH is offline   Reply With Quote
Old 2021-05-14, 19:30   #63
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

647 Posts
Default

Quote:
Originally Posted by EdH View Post
Sounds Good! Could you take a look at how much RAM is being used? If you run System Monitor, the Processes tab should show you the resting use of memory for seqinfo2. For the _1e5 file I have, my system is using 1.5 GB.

This is truly amazing !
The program uses 1.5 GB if I use a 0.1M (1e5) regina file and it uses 1.7GB if I use a 3.5M (35e5) regina file !
I don't understand how this is possible ?
garambois is online now   Reply With Quote
Old 2021-05-14, 21:44   #64
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

22·7·137 Posts
Default

Quote:
Originally Posted by garambois View Post
This is truly amazing !
The program uses 1.5 GB if I use a 0.1M (1e5) regina file and it uses 1.7GB if I use a 3.5M (35e5) regina file !
I don't understand how this is possible ?
I "think" it's because I'm initializing a huge structure (actually two), which allocates memory and then that memory is filled with the data. So, there is probably 1.5 GB set aside from the outset, and then it enlarges it as needed.
EdH is offline   Reply With Quote
Old 2021-05-15, 10:31   #65
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

10100001112 Posts
Default

OK, thanks for the explanation.
Let's wait and see in a few weeks when we're at 10M.
garambois is online now   Reply With Quote
Old 2021-05-18, 20:13   #66
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

73748 Posts
Default

I have added the capability to perform an Advanced Search on the results of a basic search. This doesn't cover everything that was mentioned, but it might present a way of doing most of the other things. Basically, the program now asks if you want to perform the Advanced Routine on the previous results.

I've provided an interesting example, where in the case of the regina_file only containing up to 1e5, there is only one even sequence that ends with the prime 269. I wonder if there are any others within a larger set.

Code:
$ ./seqinfo2 regina_file_1e5 
Data available for sequences 2 through 100000 
Sequence endings - prime: 85544, cycle: 2225, open: 12230 
Enter sequence (##/a/h/p/p##/q/u): 45987 
45987 terminates with prime 269. 
List all sequences that terminate with 269? (y/n/c/f): y 
269       595       1315      2093      2827      2888      4109      4267
6267      6357      6931      8451      9149      9307      9809      10435
12933     13133     13987     14333     15333     15931     16867     16913
17453     17469     17947     18789     19059     19733     20609     21305
22049     22891     25341     25709     28347     28673     29693     29813
32021     33627     36099     37133     38033     38165     39773     39955
40107     41451     42485     45675     45711     45987     47249     47587
47907     49015     49275     50573     51045     52109     52145     52347
54585     55289     57609     59013     59101     59309     59827     60653
61587     62609     63045     64469     65635     66555     66893     68695
69045     69547     70053     70253     71635     72833     72989     73853
76605     77393     78209     78453     79891     80813     81113     82459
83083     83309     85035     85073     86007     86609     87053     87809
87953     88013     89921     91533     92851     96031     97589     98635      
112 sequences found. 
Perform Advanced Filtering on these results? (y/n): y 
Advanced Routine using previous results 
Sequence filter a (mod b) (a b): 2 
Sequence range (begin end):  
All, Prime, Open, Cylce (a/p/o/c):  
Range for number of iterations (min max):  
Range for number of relative minimums (min max):  
Range for number of relative maximums (min max):  
Range for number of parity changes (min max):  
Range for number of consecutive even abundant terms (min max):  
Range for number of consecutive even deficient terms (min max):  
Range for number of consecutive odd deficient terms (min max):  
Range for number of consecutive odd abundant terms (min max):  
Range for number of downdriver extracts (min max):  
Range for smallest quotient of two consecutive terms (min max):  
Range for largest quotient of two consecutive terms (min max):  
Range for arith mean of all quotients of 2 consequtive terms (min max):  
Range for geo mean of all quotients of 2 consequtive terms (min max):  
Range for arith mean of # of digits of all the mins (min max):  
Range for arith mean of # of digits of all the maxes (min max):  
List filtered sequences for 0 (mod 2) from 2 through 100000? (y/n/c/f):  
2888       
1 sequence found. 
Enter sequence (##/a/h/p/p##/q/u): q
Attached Files
File Type: cpp seqinfo2.cpp (38.6 KB, 8 views)

Last fiddled with by EdH on 2021-05-18 at 20:18 Reason: Helps to attach attachment. . .
EdH is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Aliquot sequence reservations schickel Aliquot Sequences 3533 2021-08-04 04:41
A new tool to identify aliquot sequence margins and acquisitions garambois Aliquot Sequences 24 2021-02-25 23:31
Extending an aliquot sequence backwards arbooker Aliquot Sequences 5 2020-11-07 15:58
Another Aliquot Sequence site schickel Aliquot Sequences 67 2012-01-20 17:53
Useful aliquot-sequence links 10metreh Aliquot Sequences 2 2009-07-31 17:43

All times are UTC. The time now is 19:50.


Wed Aug 4 19:50:33 UTC 2021 up 12 days, 14:19, 0 users, load averages: 2.64, 2.68, 2.63

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