mersenneforum.org #of divisors Series
 User Name Remember Me? Password
 Register FAQ Search Today's Posts Mark Forums Read

2022-06-13, 00:49   #12
oreotheory

"NT"
May 2022
U.S.

19 Posts

Quote:
 Looks like 578 is the first number that goes to infinity!
That's really interesting! I had imagined that every number would end in a loop (of course who knows maybe this one does eventually). The post of @henryzz using powers of prime seemed to at least superficially confirm this impression.
(Though I like @Villiam Furik's idea of making the behavior of primes less predictable, I think that sequence isn't as exciting because the end term always ends up being a factor or multiple of the original number, whereas the original one I proposed is often random and unpredictable. (Furik Conjecture for primes: p*2/1=2p))
I'd be curious to see a tree for the original series, to see what numbers connect to which. Also very curious which numbers go to infinity. By the way, @Stargate38's 10^32-1 is a multiple of 578 (unless google's calculator is fooling me...)

2022-06-13, 01:35   #13
Batalov

"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

33×367 Posts

Quote:
 Originally Posted by oreotheory By the way, @Stargate38's 10^32-1 is a multiple of 578 (unless google's calculator is fooling me...)
Apparently it does.
How can an odd number (1032-1 is odd) be a multiple of an even number (578 is even)?

2022-06-13, 20:53   #14
Viliam Furik

"Viliam Furík"
Jul 2018
Martin, Slovakia

52×31 Posts

Quote:
 Originally Posted by oreotheory (Though I like @Villiam Furik's idea of making the behavior of primes less predictable, I think that sequence isn't as exciting because the end term always ends up being a factor or multiple of the original number, whereas the original one I proposed is often random and unpredictable. (Furik Conjecture for primes: p*2/1=2p))
The conjecture originally was that these prime loops (p -> 2p -> 2p) are the only loops and that all numbers terminate in such a loop. @henryzz pointed out more formats (actually a pattern for infinitely many patterns of infinitely many loops), extending the conjecture.

I will try to whip up a proof of this, as it seems to me to be obvious (that for every starting number it will eventually get into a loop of that format (or some another unknown format), and none will go forever). If I get somewhere, I will post an update on this.

2022-06-13, 21:05   #15
oreotheory

"NT"
May 2022
U.S.

238 Posts

Quote:
 Originally Posted by Viliam Furik The conjecture originally was that these prime loops (p -> 2p -> 2p) are the only loops and that all numbers terminate in such a loop. @henryzz pointed out more formats (actually a pattern for infinitely many patterns of infinitely many loops), extending the conjecture. I will try to whip up a proof of this, as it seems to me to be obvious (that for every starting number it will eventually get into a loop of that format (or some another unknown format), and none will go forever). If I get somewhere, I will post an update on this.
Oh, got it now. It's been true for every number I've tried in the pattern. Also not sure how to prove it. Squares and primes seem to be the building blocks, but different squares and primes will act differently...

 2022-06-14, 10:34 #16 bur     Aug 2020 79*6581e-4;3*2539e-3 58610 Posts I knew there'd be a simple in-house solution... thanks.
 2022-06-15, 15:50 #17 bur     Aug 2020 79*6581e-4;3*2539e-3 2×293 Posts I somehow messed up my post, it was supposed to quote Batalov's elegant tau(n) = #divisors(n). Some variations: Instead of tau, one could use omega, bigomega or even sigma. Or instead of division/multiplication, use subtraction/addition. Which to use could be determined by e.g. µ. To prevent a quick stop at µ = 0 one could do a(n) = a(n-1) + [c+d*µ] * tau(a(n-1)). The choice of c,d influences the growth or shrink of the series. Or a(n) = a(n-1) + (-1)^omega(a(n-1)) * phi(a(n-1)). I toyed with that function for a while. Most even numbers with additional odd prime factors grow to infinity, containing ever larger powers of 2. Odd numbers are more interesting. The script below computes the sequence and prints the factors. Again first parameter is n0, second parameter the number of terms to calculate. If omitted, it runs until a prime is reached and a(n) = 1. omegaPhi(n) = { return(n + (-1)^omega(n) * eulerphi(n)) } omPhiSeq(n,j) = { i = 0; until( n == 1 || i == j, print(i,": ",n," = ",factor(n)); n = omegaPhi(n); i++ ) } Sorry for hijacking the thread, but I thought it was a similar idea. :) Last fiddled with by bur on 2022-06-15 at 16:06 Reason: Fixed script
 2022-06-15, 19:05 #18 oreotheory     "NT" May 2022 U.S. 19 Posts The idea is very cool. Forgive me, but I have zero knowledge of coding. Could you explain to me (like I'm four) where in the code I type in the input number? I've been trying to do so on this online pari executor with no luck: https://www.tutorialspoint.com/execute_pari_online.php
2022-06-19, 20:47   #19
Viliam Furik

"Viliam Furík"
Jul 2018
Martin, Slovakia

52×31 Posts

Quote:
 Originally Posted by Batalov Apparently it does. How can an odd number (1032-1 is odd) be a multiple of an even number (578 is even)?
It must be a rational multiple

 2022-06-20, 13:14 #20 bur     Aug 2020 79*6581e-4;3*2539e-3 10010010102 Posts Just copy&paste the code into GP. Then you can run it via omPhiSeq(n,j) where n is the first value of the sequence and j (which is optional) is the number of terms you want to calculate. For example omPhiSeq(63,3) calculated the first 3 terms, whereas omPhiSeq(63) calculates terms until a prime is reached. Here, we'd enter a loop since for prime p: eulerphi(p) is p-1, omega(p) is 1, hence: p - 1 * (p - 1) = 1. Followed by 1 + 1 * 1 = 2, followed by 2 - 1 * 1 = 1. Thus, reaching a prime is a natural stop for the sequence. If you want to modify something, this is the line you want to edit: return(n + (-1)^omega(n) * eulerphi(n)). Substitute omega by bigomega or sigma or whatever. You can also a constant, such as return(n + (-1)^omega(n) * eulerphi(n) + 3) which'll add 3 to each term. No idea what'll happen then. ;) Very briefly spoken, name(a) means that you run command name with the parameter a. edit: By GP I meant the PARI/GP software if you're not familiar with that. Either download the standalone or try the browser version. Last fiddled with by bur on 2022-06-20 at 13:17
2022-06-21, 16:12   #21
Viliam Furik

"Viliam Furík"
Jul 2018
Martin, Slovakia

52×31 Posts

Quote:
 Originally Posted by Viliam Furik 9 -> 18 -> 3 -> 6 -> 6
I've noticed a mistake I made. Oddly enough, it fixed itself. Under my gcd rule, 9 goes to 3, not 18 (while doing those first few numbers in my head, I probably forgot to add 1 to the power 2 in 3^2), but 18 does indeed go to 3.

 2022-06-28, 14:46 #22 Viliam Furik     "Viliam Furík" Jul 2018 Martin, Slovakia 52·31 Posts While searching for information about the number of factors of a number and the properties of a function of the number of divisors, I've found something which I knew about, but forgot. https://en.wikipedia.org/wiki/Refactorable_number - Curtis Cooper researched them. Those are precisely the numbers for which the first-definition sequence decreases.

 Thread Tools

 Similar Threads Thread Thread Starter Forum Replies Last Post MattcAnderson MattcAnderson 0 2021-06-16 07:09 MattcAnderson MattcAnderson 1 2021-06-15 17:36 MattcAnderson MattcAnderson 3 2021-06-14 16:24 firejuggler Prime Sierpinski Project 2 2012-01-10 17:14 Citrix Math 10 2006-02-08 04:09

All times are UTC. The time now is 13:26.

Thu Aug 18 13:26:11 UTC 2022 up 10:54, 0 users, load averages: 1.45, 1.34, 1.38

Powered by vBulletin® Version 3.8.11
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.

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