mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Aliquot Sequences (https://www.mersenneforum.org/forumdisplay.php?f=90)
-   -   Aliquot sequences that start on the integer powers n^i (https://www.mersenneforum.org/showthread.php?t=23612)

garambois 2021-03-16 17:50

[QUOTE=warachwe;573819]How high can you test with base 30?
It might work for exponent 2^4*3^2*5*7*11=55440, or even some lower exponents (15120, 27720,30240, etc).
If not, 2^4*3^3*5*7*11=166320 should work.[/QUOTE]


OK, so I misunderstood what you were saying in your post #984.
I understood that you had calculated that the smallest exponent that could be suitable for base 30 was :
12 * (23 #) = 12 * 2 * 3 * 5 * 7 * 11 * 13 * 17 * 19 * 23 = 2677114440.
And 30^2677114440 has over 3.9 billion digits, if I'm not mistaken !
I think I can test up to 30^(2 * 10^6) which has over 2.9 million digits.
Maybe even for exponents > 2 * 10^6, I don't know exactly, I have to be careful not to overload the RAM !
But I have never tried something so huge.
I will try to run the test for base 30 for the 5 exponents you suggest in the order : 15120, 27720, 30240, 55440, 166320.
But I have no idea how long it will take, or even if it will be successful.
I'll keep you informed...

@VBCurtis :
Do not forget that here, the complete factorization does not interest us.
It is just for example for s(n) = s(30^15120) to find all its prime factors less than 10^4 in order to have a factor of s(n) which is abundant in order to finally prove the abundance of s(n).

warachwe 2021-03-16 19:23

[QUOTE=VBCurtis;573823]Factoring algorithms on general-form numbers such as these can be reasonably solved up to 180 digits or so, with 200 digits possible via concerted effort (and a few CPU-years of computation).

We can split off small factors up to 50-60 digits fairly easily, so a number of roughly 240 digits has a reasonable chance of a full factorization (by finding small factors summing to 50-70 digits, and cracking the rest with a full NFS algorithm).[/QUOTE]

Thank you for the information. However, as @garambois said, we don't need whole factorisation.


@garambois The 12 * (23 #) is just an exponent that I'm quite sure would be abundant, but very likely that there are lowers exponent than work.
15120, even if it work, might not be the smallest one either, just the smallest one I can come up with that have reasonable chance (depend on factor of primes factor >200 that I didn't check) to be abundant.

Happy5214 2021-03-16 20:06

30^55440 is a winner! I (ab)used FactorDB as a glorified trial-factorer, feeding it a bunch of small primes. I did check the smaller ones first, and they were not shown to be abundant.

garambois 2021-03-17 21:50

[QUOTE=Happy5214;573867]30^55440 is a winner! I (ab)used FactorDB as a glorified trial-factorer, feeding it a bunch of small primes. I did check the smaller ones first, and they were not shown to be abundant.[/QUOTE]


I'm sorry, but I don't understand your message !
What do you mean ?
Have you determined that s(30^55440) is abundant or not ?
Because when you say "30^55440 is a winner!", I would tend to understand that s(30^55440) is abundant, which seems to me to contradict the end of your message.
I'm sorry, but I am unable to grasp some of the subtleties of the English language, I work with machine translation a lot ...

:smile:

I started running my program 28 hours ago and so far have obtained the following results :
s(30^15120) not abundant if we consider all its prime factors < 10^4
s(30^27720) not abundant if we consider all its prime factors < 10^4
s(30^30240) not abundant if we consider all its prime factors < 10^4
s(30^55440) the calculation is not yet complete ! I hope to have the result tomorrow.

garambois 2021-03-17 21:54

[QUOTE=warachwe;573865]
@garambois The 12 * (23 #) is just an exponent that I'm quite sure would be abundant, but very likely that there are lowers exponent than work.
15120, even if it work, might not be the smallest one either, just the smallest one I can come up with that have reasonable chance (depend on factor of primes factor >200 that I didn't check) to be abundant.[/QUOTE]

Ok, thanks !

Happy5214 2021-03-18 09:56

[QUOTE=garambois;573974]I'm sorry, but I don't understand your message !
What do you mean ?
Have you determined that s(30^55440) is abundant or not ?
Because when you say "30^55440 is a winner!", I would tend to understand that s(30^55440) is abundant, which seems to me to contradict the end of your message.
I'm sorry, but I am unable to grasp some of the subtleties of the English language, I work with machine translation a lot ...

:smile:

I started running my program 28 hours ago and so far have obtained the following results :
s(30^15120) not abundant if we consider all its prime factors < 10^4
s(30^27720) not abundant if we consider all its prime factors < 10^4
s(30^30240) not abundant if we consider all its prime factors < 10^4
s(30^55440) the calculation is not yet complete ! I hope to have the result tomorrow.[/QUOTE]

Sorry, the original post was clearer, but I edited my message too much for my own good. The factors I found for s(30^55440) and posted to FactorDB (linked [URL="http://factordb.com/sequences.php?se=1&aq=30%5E55440&action=last20&fr=0&to=100"]here[/URL]) are enough to show that the number is abundant. I meant to say the ones smaller than 55440 are not abundant, as your results showed.

Happy5214 2021-03-18 18:12

22^43 terminates at the prime 240349.

garambois 2021-03-18 18:44

[QUOTE=Happy5214;574011]Sorry, the original post was clearer, but I edited my message too much for my own good. The factors I found for s(30^55440) and posted to FactorDB (linked [URL="http://factordb.com/sequences.php?se=1&aq=30%5E55440&action=last20&fr=0&to=100"]here[/URL]) are enough to show that the number is abundant. I meant to say the ones smaller than 55440 are not abundant, as your results showed.[/QUOTE]


OK, thank you very much Happy !
When I saw your result on FactorDB, I really thought that my program must have a problem and that it was abnormally slow !
And I tried to understand why my program was so slow compared to FactorDB and [U]I found why[/U] !
It was a stupid instruction that did an unnecessary primality test on very large numbers.
And so, by removing this unnecessary instruction, I multiplied the speed of the program by more than 1000.

Result :
In a few minutes, I tested all the even exponents for i from 1 to 100,000 for base 30.
And indeed, the exponent i = 55440 = 2^4*3^2*5*7*11 is the first one that is suitable : there is no smaller one.
Congratulations for warachwe who had foreseen it !
And so, we have s(30^55440) which is abundant.
But we also have s(30^(55440*2)) and s(30^(55440*3)) which are abundant.
I have not tested the k>3.

And the next i that is suitable is : i = 65520 = 2^4*3^2*5*7*13.
So, we have s(30^65520), s(30^(65520*2)) and s(30^(65520*3)) which are abundant.
[URL="http://factordb.com/sequences.php?se=1&aq=30%5E65520&action=last20&fr=0&to=100"]See on FactorDB the factors less than 10^4 for s(30^65520)[/URL].
The next is i = 75600 = 2^4*3^3*5^2*7, as well as its double and its triple.
[URL="http://factordb.com/sequences.php?se=1&aq=30%5E75600&action=last20&fr=0&to=100"]See on FactorDB the factors less than 10^4 for s(30^75600)[/URL].
And then next is : 80640 = 2^8*3^2*5*7, but this time, its double and its triple are not suitable !
The next is i = 90720 = 2^5*3^4*5*7, as well as its double and its triple.

This new program allows me to restart the search for base 3.
In less than an hour, I have already tested all the odd exponents up to 100,000.
Let's wait and see.

I'll keep you posted if I find anything interesting ...

Happy5214 2021-03-18 21:08

Good work! Hopefully the faster code proves fruitful for base 3.

[QUOTE=garambois;574070] The next is i = 75600 = 2^4*3^3*5^2*7, as well as [B]his[/B] double and [B]his[/B] tripl. [sic][/QUOTE]

I didn't know numbers were men. :wink: I'm assuming that's a translation error. (The proper English possessive pronoun for an abstract concept like a number is the neuter "its".)

garambois 2021-03-19 17:56

[QUOTE=Happy5214;574082]Good work! Hopefully the faster code proves fruitful for base 3.
[/QUOTE]

I hope too !
For base 3, in almost 24 hours, the program tested all odd exponents up to 800,000 : still nothing !
I have also just launched the execution of calculations for base 5.


[QUOTE=Happy5214;574082]
I didn't know numbers were men. :wink: I'm assuming that's a translation error. (The proper English possessive pronoun for an abstract concept like a number is the neuter "its".)[/QUOTE]

Sorry, I didn't notice the machine translator error.
But honestly, without the machine translation, I don't think I could communicate with all of you !
There must still be a lot of other errors in my posts that you don't point out to me, because my English is very poor.
If an administrator wants to make the correction in post #998, no problem : I understand that it hurts the eyes of English speakers !

:hello:

Happy5214 2021-03-20 12:05

[QUOTE=garambois;574137]I hope too !
For base 3, in almost 24 hours, the program tested all odd exponents up to 800,000 : still nothing !
I have also just launched the execution of calculations for base 5.[/QUOTE]

Please don't tell me you're trying to do anything efficient in Python. This is a job for C. You'd probably get a significant speed boost with a compiled C or C++ program with GMP compared to Python.

[QUOTE=garambois;574137]Sorry, I didn't notice the machine translator error.
But honestly, without the machine translation, I don't think I could communicate with all of you !
There must still be a lot of other errors in my posts that you don't point out to me, because my English is very poor.
If an administrator wants to make the correction in post #998, no problem : I understand that it hurts the eyes of English speakers !

:hello:[/QUOTE]

That's a really bad job by the translator. Are numbers considered masculine nouns in French, because I don't believe Romance languages inherited the neuter from Latin?

English is a hard language to master. Many native English speakers haven't done it. :smile: It's a shame that such a complicated language has become the global lingua franca (thanks to the global hegemonic status of the British Empire and the US over the past two centuries). I would have preferred Esperanto, which was designed for the purpose.


All times are UTC. The time now is 22:27.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.