mersenneforum.org Aliquot sequences that start on the integer powers n^i
 Register FAQ Search Today's Posts Mark Forums Read

 2021-02-15, 19:52 #793 RichD     Sep 2008 Kansas 22·829 Posts Base 43 can be added at the next update. All expected terminating sequences are terminated and all remaining composites are greater than C78. Starting initial work on base 47.
2021-02-16, 19:34   #794
EdH

"Ed Hall"
Dec 2009

1110011010012 Posts

Quote:
 Originally Posted by garambois . . . Did I understand correctly, you also have a C++ program that does this job faster ? This other program would also interest me.
Here is the source for a C++ program. I'm afraid it's probably a bit "amateurish," and it has very little internal commenting, but it should compile on linux with:
Code:
g++ alimerge3.cpp -o alimerge3
and can then be invoked by:
Code:
./alimerge3 37 1 80
to run the entire base 37 table. Here's an example (wrapped in date calls to provide timing)*:
Code:
\$ date;./alimerge3 37 1 80;date
Tue 16 Feb 2021 02:19:01 PM EST
Running base 37 from 1 through 80 . . .
37^6:i126 merges with 109920:i1275
37^8:i151 merges with 1632:i37
37^12:i1057 merges with 10824:i28
37^18:i1430 merges with 3366:i2
37^22:i519 merges with 1567300:i0
37^30:i1193 merges with 35856:i3
Tue 16 Feb 2021 02:22:27 PM EST
That comes in at about 3:26 to run the entire base 37 table. That's about 206 vs. 1211 seconds for the Bash version.

* I tried to add an elapsed time function, but for some reason the timing was extremely inaccurate with the several methods I tried, so I left it out.
Attached Files
 alimerge3.cpp.txt (5.0 KB, 15 views)

2021-02-17, 08:06   #795
garambois

"Garambois Jean-Luc"
Oct 2011
France

10001011102 Posts

Quote:
 Originally Posted by RichD Base 43 can be added at the next update. All expected terminating sequences are terminated and all remaining composites are greater than C78. Starting initial work on base 47.

OK, a lot of thanks !
The next update will be next weekend.

2021-02-17, 08:34   #796
garambois

"Garambois Jean-Luc"
Oct 2011
France

2·32·31 Posts

Quote:
 Originally Posted by EdH Here is the source for a C++ program.

OK, a lot of thanks !
I will test this program very quickly, but I can't get it to work at the moment, I don't know why yet ?

Otherwise, more generally, in three days, I will be on vacation.
I will try to do some data analysis.
Hopefully the large amount of new data accumulated allows us to formulate new small conjectures ... that is our goal !

2021-02-17, 13:54   #797
EdH

"Ed Hall"
Dec 2009

368910 Posts

Quote:
 Originally Posted by garambois OK, a lot of thanks ! I will test this program very quickly, but I can't get it to work at the moment, I don't know why yet ? Otherwise, more generally, in three days, I will be on vacation. I will try to do some data analysis. Hopefully the large amount of new data accumulated allows us to formulate new small conjectures ... that is our goal !
The source code I provided works under various linux machines here, compiling with no issues. It even compiles on my RPi machines. Let me know of any errors/warnings if any occur. I do now have a version that includes timing. If I add that one, I'll still leave the other.

The source should compile on Windows, as well, but it's been so long since I've worked with Windows, that I wouldn't be able to even try.

If there is some data in particular you'd like harvested, let me know and I can try to do that as well. I'm trying to keep a local up-to-date set of all the tables, at the ready for such work.

Have a Great Vacation!

 2021-02-17, 18:27 #798 garambois     "Garambois Jean-Luc" Oct 2011 France 55810 Posts Most of this vacation, I am going to work on the aliquot sequences. I have already started the preliminary work, I am writing and running programs. But I would need a little help ! Would someone be able to calculate the following sequences, but only up to index 2, and of course, to enter the results on factordb : (9699690 * 23)^20 (9699690 * 23 * 29 * 31)^14 (9699690 * 23 * 29 * 31 * 37)^14 Indeed, my computer is very busy and that would allow me to have the results already next weekend. I still need more sequences, but the calculations will be shorter and I can handle them. @Edwin : I'm on a track... Since I looked again at your "2-Stats.txt" file (from here), I got a new idea because of the fourth column. But I have to do some more analysis to be sure that it is worth it that I expose you on this forum what I noticed. The calculations I am requesting above relate to this work. Last fiddled with by garambois on 2021-02-17 at 18:29
 2021-02-17, 19:04 #799 EdH     "Ed Hall" Dec 2009 Adirondack Mtns 7×17×31 Posts I have queued up index 1 for both of the first two, but they will take a little while to finish. Do you need index 2 also fully factored or just reach index 2? The bulk of my machines are running a rather large job right now, so I am somewhat limited.
2021-02-17, 20:59   #800
garambois

"Garambois Jean-Luc"
Oct 2011
France

55810 Posts

Quote:
 Originally Posted by EdH I have queued up index 1 for both of the first two, but they will take a little while to finish. Do you need index 2 also fully factored or just reach index 2? The bulk of my machines are running a rather large job right now, so I am somewhat limited.

I need just to reach index 2 but not to factor. Thank you very much Edwin !

2021-02-17, 22:23   #801
EdH

"Ed Hall"
Dec 2009

71518 Posts

Quote:
 Originally Posted by garambois I need just to reach index 2 but not to factor. Thank you very much Edwin !
The first two are at index 2. The third is in work.

 2021-02-17, 22:33 #802 garambois     "Garambois Jean-Luc" Oct 2011 France 2·32·31 Posts Thank you very much Ed !
2021-02-18, 00:49   #803
EdH

"Ed Hall"
Dec 2009

7×17×31 Posts

Quote:
 Originally Posted by garambois . . . (9699690 * 23 * 29 * 31 * 37)^14 . . .
Third one is at index 2, also.

 Similar Threads Thread Thread Starter Forum Replies Last Post fivemack FactorDB 46 2021-02-21 10:46 schickel FactorDB 18 2013-06-12 16:09 garambois Aliquot Sequences 34 2012-06-10 21:53 Andi47 FactorDB 21 2011-12-29 21:11 schickel mersennewiki 0 2008-12-30 07:07

All times are UTC. The time now is 18:16.

Sun Apr 18 18:16:24 UTC 2021 up 10 days, 12:57, 1 user, load averages: 3.07, 2.76, 2.61