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

2021-03-05, 03:26   #947
EdH

"Ed Hall"
Dec 2009

13×192 Posts

Quote:
 Originally Posted by RichD Base 38 can be added at the next update. At least one merge and a cycle can be found in this table.
I was updating my set of sequences and I found the cycle for 38^9 (1184/1210). I also see yoyo has finished 439^46 into the same cycle (1184/1210). In fact, they merge at the previous term (2362) to 1184.

2021-03-05, 04:08   #948
Happy5214

"Alexander"
Nov 2008
The Alamo City

3·52·11 Posts

Quote:
 Originally Posted by henryzz Assuming no factors below 1e20 then the maximum abundance contribution of the C353 can be upper bounded by ((1e20+1)/1e20)^(floor(353/20))= (1+1e-20)^17= 1.00000000000000000017 or there about.
I ran a t35 and 428 curves at B1=3M on the remaining cofactor (a C27 had already been found since that post), and according to http://www.wraithx.net/math/ecmprobs/ecmprobs.html, the probability of a uncaught factor 30 digits or less is about 1 in 2 million. It's not a proof unless/until we factor that number, but that term is most likely deficient and the associated conjecture about the 38^12n series is probably false.

Last fiddled with by Happy5214 on 2021-03-05 at 04:08

2021-03-05, 07:44   #949
henryzz
Just call me Henry

"David"
Sep 2007
Liverpool (GMT/BST)

176816 Posts

Quote:
 Originally Posted by Happy5214 I ran a t35 and 428 curves at B1=3M on the remaining cofactor (a C27 had already been found since that post), and according to http://www.wraithx.net/math/ecmprobs/ecmprobs.html, the probability of a uncaught factor 30 digits or less is about 1 in 2 million. It's not a proof unless/until we factor that number, but that term is most likely deficient and the associated conjecture about the 38^12n series is probably false.
The p27 factor just strengthens the argument. Choosing 20 digits probably wasn't sensible given that the trivially provable 6 digits indicates it is deficient (1+10^-6)^floor(353/6) is approximately 1.000058. These techniques are regularly used by odd perfect number searchers as part of proofs. I think we can call this a proof of contradiction now.

2021-03-05, 14:35   #950
EdH

"Ed Hall"
Dec 2009

469310 Posts

Quote:
 Originally Posted by RichD In keeping with the spirit of investigating base 2*p, next up is base 46. Preliminary results show a similar phenomena. Additionally, 46^6n shows an abundance of (3 * 5 * 7) along with 46^12n (3 * 5 * 7 * 13) to advance the sequence.
Base 46 shows some merges:
Code:
46^9:i231 merges with 11408:i15
46^21:i452 merges with 9852:i11
46^29:i913 merges with 67260:i2

2021-03-05, 17:49   #951
garambois

"Garambois Jean-Luc"
Oct 2011
France

16368 Posts

Quote:
 Originally Posted by EdH I was updating my set of sequences and I found the cycle for 38^9 (1184/1210). I also see yoyo has finished 439^46 into the same cycle (1184/1210). In fact, they merge at the previous term (2362) to 1184.

Quote:
 Originally Posted by EdH Base 46 shows some merges: Code: 46^9:i231 merges with 11408:i15 46^21:i452 merges with 9852:i11 46^29:i913 merges with 67260:i2

OK, many thanks !
We will visualize all of this in the next update ...

Last fiddled with by garambois on 2021-03-05 at 17:53 Reason: Addition of the second quote

 2021-03-06, 12:49 #952 garambois     "Garambois Jean-Luc" Oct 2011 France 2·463 Posts Two new conjectures and general method to find increasing sequences from index 1 for all bases (with bases and exponents of the same parity) After reading all your exchanges on the abundances about base 38, I had some new ideas. It helped me find several thousand sequences for base 6 with even exponents with abundant index 1 term. I remind you that until today, we had nothing for base 6. I therefore set out below two new conjectures : (138) et (139). I explain my approach for each of the two conjectures. This lead me to propose a general method which makes it possible to find for each base sequences with exponents of the same parity as the base with terms of index 1 which are abundant. Conjecture (138) : Base 6 sequences starting with 6^((2^3 * 3^2 * 5 * 7) * k) are increasing at least from index 1 to 2. This is proved for all k integers such that 0 < k < 5001. It remains to prove it for all integer values ​​of k, but this is too difficult for me. Indeed, s(6^(2^3 * 3^2 * 5 * 7)) = d * c, with d = 5^2 * 7^2 * 11 * 13 * 19 * 29 * 31 * 37 * 41 * 43 * 61 * 71 * 73 * 127 * 181 * 211 * 281 * 337 * 421 * 631. But, d is abundant, because s(d) / d > 1.2. And using a very small program, it is very easy to verify for all k from 1 to 5000 that d is a divisor of s(6^((2^3 * 3^2 * 5 * 7) * k)). Here is the little program in python : Code: k=1 d=5^2*7^2*11*13*19*29*31*37*41*43*61*71*73*127*181*211*281*337*421*631 print('d = 5^2*7^2*11*13*19*29*31*37*41*43*61*71*73*127*181*211*281*337*421*631') print('s(d)/d =',(sigma(d)-d)*1.0/d) while True: n=6^((2^3*3^2*5*7)*k) m=sigma(n)-n if m%d!=0: print ('***** ',k) if k%1000==0: print (k) k+=1 Here is the result after a few hours : Code: sage: load("conj_6.sage") d=5^2*7^2*11*13*19*29*31*37*41*43*61*71*73*127*181*211*281*337*421*631 s(d)/d = 1.20464241970184 1000 2000 3000 4000 5000 On the other hand, I don't know how to prove that d divides the terms of index 1 for all k. Now let me introduce you a general method to find increasing sequences from index 1 for all bases with bases and exponents of the same parity. You can proceed with any base as with base 6, as in the example above. The problem is to find a good "starting exponent" i = 2^3 * 3^2 * 5 * 7 in our example. This is done by trial and error, but i think the search for a good starting exponent i can be automated. It must be chosen so that we can find a compound factor d which is abundant and which then divides the terms of index 1 for all the sequences for k from 1 to 10. We can finally verify using of a program similar to the one shown above that d divides well all the multiple exponents of i for k from 1 to 1000 for example. Then, we can formulate a solid conjecture, like the conjecture (138). For conjecture 137, for the sequences 2^(12 * k), we have d = 3^2 * 5 * 7 * 13 and s(d) / d > 1.13, so d is abundant. For the sequences 2^(40 * k), we have d = 3 * 5^2 * 11 * 17 * 31 * 41 and s(d) / d > 1.01. I did not check for the other starting exponents of conjecture 137 (90, 140, 210, 220 and 330) if the corresponding d factors are indeed abundant. Strictly speaking, if the factor d corresponding to one of these exponent was not abundant, it would be necessary to remove this exponent with its multiples from the conjecture. To find this starting exponent i, we noticed that for almost all bases, the more the exponent breaks down into small prime numbers, the more chance we have of finding an abundant value of d. But, we don't know how to prove it for even bases. It's just an observation that seems to be true. Thus, even for base 3 and in the general case, for base b, we could take a starting exponent of this type : i = 3 * 5 * 7 * 11 * 13 * 17 * 19 ... and we conjecture that after several tries, we would find a abundant d that would divide the b^(i * k) for the first values of k. But the most difficult will certainly then be to prove that d divides s (b^(i * k)) for any integer k ! And also to find the i as small as possible. This general method should work in theory. But in practice, I have not been able to find an exponent i and an abundant compound factor d for the base 3. The numbers grow to astronomical size very quickly ! But we feel that we have enough good reasons to state the following more general conjecture. Conjecture (139) : There exists for each base b a starting exponent i, such that for any integer k, the sequences b^(i * k) are increasing from index 1 during at least one iteration. One last little idea To finish, I would like to present to you an idea which should allow us to demonstrate that the conjecture for the base 38 was false, with the d = 3 * 5 * 7 * 13 which is deficient. I wrote a little program. I calculate for as many k as possible the numbers c = s(38^(12 * k)) / d = s(38^(12 * k)) / (3 * 5 * 7 * 13). And then I do a primality test on c. If c is a prime number for a single value of k, the conjecture is invalidated. Here is the little program : Code: k=1 while True: n=38^(12*k) m=sigma(n)-n c=m/(3*5*7*13) if is_pseudoprime(c): print ('***** ',k) if k%100==0: print (k) k+=1 Unfortunately, after several hours no c was prime even though I had tested all k from 1 to 2300. But in principle, this test should work, to my mind, if we let the program run longer ... Here is the result of the program : Code: sage: load("conj_38_12k.sage") 100 200 300 ... ... 2000 2100 2200 2300 Last fiddled with by Uncwilly on 2021-03-06 at 22:40 Reason: Per OP request, see post #955
2021-03-06, 16:18   #953
henryzz
Just call me Henry

"David"
Sep 2007
Liverpool (GMT/BST)

10111011010002 Posts

Quote:
 Originally Posted by garambois This is proved for all k integers such that 0 < k < 5001. It remains to prove it for all integer values ​​of k, but this is too difficult for me. Indeed, s(6^(2^3 * 3^2 * 5 * 7)) = d * c, with d = 5^2 * 7^2 * 11 * 13 * 19 * 29 * 31 * 37 * 41 * 43 * 61 * 71 * 73 * 127 * 181 * 211 * 281 * 337 * 421 * 631. But, d is abundant, because s(d) / d > 1.2. And using a very small program, it is very easy to verify for all k from 1 to 5000 that d is a divisor of s(6^((2^3 * 3^2 * 5 * 7) * k)). On the other hand, I don't know how to prove that d divides the terms of index 1 for all k.
I think the trick would be to prove that:
2^(2^3 * 3^2 * 5 * 7 * k + 1) - 1 == 1 mod d
3^(2^3 * 3^2 * 5 * 7 * k + 1) - 1) / 2 == 1 mod d
6^(2^3 * 3^2 * 5 * 7 * k) == 1 mod d

If this is true then s(6^(2^3 * 3^2 * 5 * 7 * k)) = (2^(2^3 * 3^2 * 5 * 7 * k + 1) - 1)*3^(2^3 * 3^2 * 5 * 7 * k + 1) - 1) / 2 - 6^(2^3 * 3^2 * 5 * 7 * k) = 1*1-1 mod d

Quote:
 Conjecture (139) : There exists for each base b a starting exponent i, such that for any integer k, the sequences b^(i * k) are increasing from index 1 during at least one iteration. One last little idea To finish, I would like to present to you an idea which should allow us to demonstrate that the conjecture for the base 38 was false, with the i = 3 * 5 * 7 * 13 which is deficient. I wrote a little program. I calculate for as many k as possible the numbers c = 38^(12 * k) / i = 38^(12 * k) / (3 * 5 * 7 * 13). And then I do a primality test on c. If c is a prime number for a single value of k, the conjecture is invalidated.
Surely you mean c = s(38^(12 * k)) / i

I suspect that the way to attack this for individual b at least would be to find an i that has a d that is abundant.

As far as I can tell it is usually quite easy. Finding that for 22^(180 * k), d= 3^3*7*13*19*31*5*37*61 which is abundant was quite easy.

d=3^2*5*7*13 is abundant.
2^(180 * k + 1) - 1 == 1 mod d
p^(180 * k) == 1 mod d for a lot of p. I haven't worked out the logic behind this yet. http://factordb.com/index.php?query=...at=1&sent=Show
(2*p)^(180 * k) == 1 mod d if gcd(d,p)==1 as indicated by http://factordb.com/index.php?query=...at=1&sent=Show

I would be very interested if someone could come up with a case where (2*p)^(180*k) does not have an abundant divisor for some p.

2021-03-06, 18:19   #954
warachwe

Aug 2020

19 Posts

Quote:
 Originally Posted by henryzz I think the trick would be to prove that: 2^(2^3 * 3^2 * 5 * 7 * k + 1) - 1 == 1 mod d 3^(2^3 * 3^2 * 5 * 7 * k + 1) - 1) / 2 == 1 mod d 6^(2^3 * 3^2 * 5 * 7 * k) == 1 mod d
For p prime, and and a such that gcd(a,p)=1, we must have a^(p-1) == 1 mod p.
So if p-1 divides 2^3 * 3^2 * 5 * 7, it means 2^(2^3 * 3^2 * 5 * 7 * k) == 1^k == 1 mod p, hence 2^(2^3 * 3^2 * 5 * 7 * k+1)-1 == 1^k*2-1 == 1 mod p.
Same thing for 3 and 6.

That take care of p=5+,7+,11,13,19,29,31,37,41,43,61,71,73,127,181,211,281,421,631
(+ This only prove that 5,7 divide s(6^(2^3 * 3^2 * 5 * 7 * k)), not 5^2,7^2)
For p=337, we have (p-1)/2 divides 2^3 * 3^2 * 5 * 7. Possible value for a^((p-1)/2) is +-1 mod p. I guess that mean 2^((337-1)/2),3^((337-1)/2),6^((337-1)/2) are happened to be 1 mod 337.

To prove that 5^2,7^2 divide s(6^(2^3 * 3^2 * 5 * 7 * k)), we use the fact that if a == 1 mod p, then a^p == 1 mod p^2.
(5-1) divides 2^3 * 3^2 * 7, so a^(2^3 * 3^2 * 5 * 7) == 1 mod 5^2 for (a,p)=1
(7-1) divides 2^3 * 3^2 * 5, so a^(2^3 * 3^2 * 5 * 7) == 1 mod 7^2 for (a,p)=1

2021-03-06, 18:23   #955
garambois

"Garambois Jean-Luc"
Oct 2011
France

2×463 Posts

Quote:
 Originally Posted by henryzz Surely you mean c = s(38^(12 * k)) / i

Many thanks henryzz !
Too much fatigue after all this work !
Please, would it be possible for someone (Edwin maybe?) to make the following correction in my post # 952 :
Quote:
 One last little idea To finish, I would like to present to you an idea which should allow us to demonstrate that the conjecture for the base 38 was false, with the i = 3 * 5 * 7 * 13 which is deficient. I wrote a little program. I calculate for as many k as possible the numbers c = 38^(12 * k) / i = 38^(12 * k) / (3 * 5 * 7 * 13).
Should be :
Quote:
 One last little idea To finish, I would like to present to you an idea which should allow us to demonstrate that the conjecture for the base 38 was false, with the d = 3 * 5 * 7 * 13 which is deficient. I wrote a little program. I calculate for as many k as possible the numbers c = s(38^(12 * k)) / d = s(38^(12 * k)) / (3 * 5 * 7 * 13).

2021-03-06, 18:51   #956
garambois

"Garambois Jean-Luc"
Oct 2011
France

2×463 Posts

Quote:
 Originally Posted by henryzz I suspect that the way to attack this for individual b at least would be to find an i that has a d that is abundant.

Yes, that is exactly what I am saying above the statement of the conjecture (139).
The small paragraph below titled "One last little idea" was just intended to attempt to provide a method to invalidate a conjecture that would use a deficient d in its statement.
I did not have fun looking for a conjecture true with a larger i and an abundant d for base 38.
But it would indeed not be difficult ...
And in my opinion, this is easy for any base in the form 2*p.
The difficulty is much greater for an odd base like 3, I have tried and we have to take values of i which give starting numbers of sequences of several million digits.

Quote:
 Originally Posted by henryzz I would be very interested if someone could come up with a case where (2*p)^(180*k) does not have an abundant divisor for some p.

It seems difficult to me !
I'll think about it...

Last fiddled with by garambois on 2021-03-06 at 18:56 Reason: Addition of the second quote and my answer

 2021-03-06, 19:21 #957 garambois     "Garambois Jean-Luc" Oct 2011 France 39E16 Posts @henryzz and warachwe : It is difficult for me to understand your reasoning, I am not used to this theoretical reasoning ! But does this mean that conjecture (138) is fully proven ? I will try to understand all this, because these mechanisms interest me a lot, like everything you taught us in the summer of 2020. I often use your explanations to understand phenomena on aliquot sequences.

 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 17:55.

Tue Aug 9 17:55:36 UTC 2022 up 33 days, 12:42, 1 user, load averages: 1.13, 1.25, 1.21

Copyright ©2000 - 2022, 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.

≠ ± ∓ ÷ × · − √ ‰ ⊗ ⊕ ⊖ ⊘ ⊙ ≤ ≥ ≦ ≧ ≨ ≩ ≺ ≻ ≼ ≽ ⊏ ⊐ ⊑ ⊒ ² ³ °
∠ ∟ ° ≅ ~ ‖ ⟂ ⫛
≡ ≜ ≈ ∝ ∞ ≪ ≫ ⌊⌋ ⌈⌉ ∘ ∏ ∐ ∑ ∧ ∨ ∩ ∪ ⨀ ⊕ ⊗ 𝖕 𝖖 𝖗 ⊲ ⊳
∅ ∖ ∁ ↦ ↣ ∩ ∪ ⊆ ⊂ ⊄ ⊊ ⊇ ⊃ ⊅ ⊋ ⊖ ∈ ∉ ∋ ∌ ℕ ℤ ℚ ℝ ℂ ℵ ℶ ℷ ℸ 𝓟
¬ ∨ ∧ ⊕ → ← ⇒ ⇐ ⇔ ∀ ∃ ∄ ∴ ∵ ⊤ ⊥ ⊢ ⊨ ⫤ ⊣ … ⋯ ⋮ ⋰ ⋱
∫ ∬ ∭ ∮ ∯ ∰ ∇ ∆ δ ∂ ℱ ℒ ℓ
𝛢𝛼 𝛣𝛽 𝛤𝛾 𝛥𝛿 𝛦𝜀𝜖 𝛧𝜁 𝛨𝜂 𝛩𝜃𝜗 𝛪𝜄 𝛫𝜅 𝛬𝜆 𝛭𝜇 𝛮𝜈 𝛯𝜉 𝛰𝜊 𝛱𝜋 𝛲𝜌 𝛴𝜎𝜍 𝛵𝜏 𝛶𝜐 𝛷𝜙𝜑 𝛸𝜒 𝛹𝜓 𝛺𝜔