mersenneforum.org Extending an aliquot sequence backwards
 User Name Remember Me? Password
 Register FAQ Search Today's Posts Mark Forums Read

2017-11-08, 18:30   #1
arbooker

"Andrew Booker"
Mar 2013

5×17 Posts
Extending an aliquot sequence backwards

The aliquot sequence with start value 461214 is the current longest known open ended sequence starting below 1e6 (it merges with the 4788 sequence around index 6000). Right now it's just shy of 19k iterations. Just for fun, I thought I would try to extend the sequence backwards using Goldbach so as to make a longer record.

It's easy to check that s(461214)=s(670097^2). Given an odd number n > 8, we can (conjecturally) write n=p+q+1=s(pq) for distinct primes p,q. Starting from n=670097^2, I found the smallest p such that q=n-1-p is a probable prime, replaced n by pq, and repeated for 1000 iterations. The sequence of p values is attached.

I then used Primo to certify the primality of the first 700 values of q, the largest of which has 2490 digits. It's getting very slow, but I might let it run up to 1000. I'm aware that none of this serves any real purpose, but if anyone would like to join me in this quest, feel free.
Attached Files
 p.txt (5.1 KB, 285 views)

 2017-11-08, 20:01 #2 garambois     Oct 2011 1011111002 Posts I did the same kind of exercise, see on this page (http://www.aliquotes.com/remonter_suite_envers.html) my go backwards of 3, 7, 11, 13, 17 on more than 2000 iterations, as well as the go backwards of 2005020 on more than 1000 iterations (I was talking about it here #7 : http://www.mersenneforum.org/showthread.php?t=18641)
2017-11-09, 20:48   #3
arbooker

"Andrew Booker"
Mar 2013

5×17 Posts

Quote:
 Originally Posted by garambois I did the same kind of exercise
Yet more evidence that no idea is original. The only thing I can think to add is the primality certificates. It's reasonable to generate them with Primo for 1000 iterations or so.

Good idea putting the sequence on factordb; I thought it might be too large for that, but I guess not. Will factordb compute the sequence for me if I just enter the 3000-digit starting value? Also, I don't fancy uploading 1000 certificates by hand; if anyone has a script that they don't mind sharing, that would be a great help.

 2020-10-26, 13:54 #4 sweety439   Nov 2016 22×691 Posts The number 1578 merges to 56440, but it is not known whether the sequence of 56440 (currently at 166 digits) also merges to 56440, thus it is not known whether 56440 is member of sociable number cycle (but if so, then this cycle will contain 6000+ numbers ....)
2020-10-26, 14:05   #5
EdH

"Ed Hall"
Dec 2009

2·32·197 Posts

Quote:
 Originally Posted by sweety439 The number 1578 merges to 56440, but it is not known whether the sequence of 56440 (currently at 166 digits) also merges to 56440, thus it is not known whether 56440 is member of sociable number cycle (but if so, then this cycle will contain 6000+ numbers ....)
factordb would pick up a cycle. Since all the factors of all the terms of a sequence exist within the database, it would not be shown as open ended. As to whether this eventually happens, only more factoring will tell. . .

 2020-11-07, 15:58 #6 sweety439   Nov 2016 22·691 Posts so the sequences A003416 and A122726 might miss numbers (i.e. it is unknown whether there are members of sociable number cycle <= 1799281330 other than numbers listed in A003416), however, the smallest number whose aliquot sequence has not yet been fully determined (276) cannot be member of sociable number cycle, since it is an untouchable number (A005114), however, the number 56440 might be. Last fiddled with by sweety439 on 2020-11-07 at 16:01

 Similar Threads Thread Thread Starter Forum Replies Last Post schickel Aliquot Sequences 3381 2021-01-19 17:24 schickel Aliquot Sequences 67 2012-01-20 17:53 yoyo FactorDB 6 2012-01-12 20:58 10metreh Aliquot Sequences 2 2009-07-31 17:43 philmoore Math 3 2009-03-20 19:04

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

Sun Jan 24 18:29:47 UTC 2021 up 52 days, 14:41, 0 users, load averages: 2.53, 2.84, 3.48