mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2012-03-15, 14:29   #1
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2×7×461 Posts
Default Some transition probabilities

A.B.C.D means 2^A*3^B*5^C*7^D

These are from 1250 lines from the 2^2*3 project

Code:
2.0.0.0		0.890 2.0.0.0	0.044 2.0.1.0	*0.015 3.0.0.0*	*0.015 1.0.0.0*
		0.011 2.0.2.0	0.007 2.1.0.0	*0.004 7.0.1.0*	0.004 2.0.7.0
		*0.004 4.0.0.0*	*0.004 3.2.0.0*	0.004 2.0.3.0	
2.1.0.0		0.832 2.1.0.0	0.052 2.0.0.0	0.047 2.1.1.0	0.044 2.2.0.0
		0.010 2.3.0.0	0.008 2.1.2.0	0.003 2.0.1.0	0.003 2.2.1.0
		0.003 2.4.0.0	
2.1.1.0		0.608 2.1.1.0	0.278 2.1.0.0	0.089 2.1.2.0	0.013 2.1.3.0
		0.013 2.2.1.0	
2.2.0.0		0.620 2.2.0.0	0.092 2.1.0.0	0.063 2.2.1.0	0.049 2.3.0.0
		0.042 2.0.0.0	0.021 2.2.2.0	0.021 2.1.1.0	0.021 2.5.0.0
		0.014 2.2.3.0	*0.007 1.0.0.0*	0.007 2.0.1.0	0.007 2.3.1.0
		0.007 2.4.0.0	*0.007 3.0.0.0*	0.007 2.5.1.0	*0.007 4.4.0.0*
		0.007 2.6.0.0
Read this as 'if you have currently a factor 2^2*3^2 (and are about 115 digits), probability is about 0.092 that you will be 2^2*3 next iteration'. The probabilities are taken from at least 50 (about 300 for 2.1.0.0) lines. Ones with stars are where the power of two changed.

Last fiddled with by fivemack on 2012-03-15 at 14:31
fivemack is offline   Reply With Quote
Old 2012-03-15, 16:08   #2
firejuggler
 
firejuggler's Avatar
 
"Vincent"
Apr 2010
Over the rainbow

5·11·53 Posts
Default

what about those marked with an asterisk? not enough data to be confident?
firejuggler is offline   Reply With Quote
Old 2012-03-15, 16:57   #3
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2·7·461 Posts
Default

Ones with stars are where the power of two changed.
fivemack is offline   Reply With Quote
Old 2012-03-15, 18:11   #4
firejuggler
 
firejuggler's Avatar
 
"Vincent"
Apr 2010
Over the rainbow

B6316 Posts
Thumbs up

sorry i prefer my data sorted vertically

Code:
                                
driver/guide    next iteration    % chance
2^2        
    2^2       0,89
    2^2*5     0,044
    2^3       0,015
    2^2       0,015
    2^2*5^2   0,011
    2^2*3     0,007
    2^7*5     0,004
    2^2*5^7   0,004
    2^4       0,004
    2^3*3^2   0,004
    2^2*5^3   0,004
    
2^2*3        
    2^2*3     0,832
    2^2       0,052
    2^2*3*5   0,047
    2^2*3^2   0,044
    2^2*3^3   0,01
    2^2*3*5^2 0,008
    2^2*5     0,003
    2^2*3^2*5 0,003
    2^2*3^4   0,003
  
2^2*3*5        
    2^2*3*5   0,608
    2^2*3     0,278
    2^2*3*5   0,089
    2^2*3*5^3 0,013
    2^2*3^2*5 0,013
      
2^2*3^2        
    2^2*3^2   0,62
    2^2*3     0,092
    2^2*3^2*5 0,063
    2^2*3^3   0,049
    2^2       0,042
    2^2*3^2*5^2 0,021
    2^2*3*5   0,021
    2^2*3^5   0,021
    2^2*3^2*5^3 0,014
    2         0,007
    2^2*5    0,007
    2^2*3^3*5 0,007
    2^2*3^4   0,007
    2^3       0,007
    2^2*3^5*5 0,007
    2^4*3^4   0,007
    2^2*3^6   0,007
another question, no 2^2*7 driver?

Last fiddled with by firejuggler on 2012-03-15 at 18:12
firejuggler is offline   Reply With Quote
Old 2012-03-15, 23:28   #5
henryzz
Just call me Henry
 
henryzz's Avatar
 
"David"
Sep 2007
Liverpool (GMT/BST)

22×3×7×73 Posts
Default

Quote:
Originally Posted by firejuggler View Post
another question, no 2^2*7 driver?
You can't get to 2^2*7 directly through 2^2 because sigma(2^2)=7. You have to go though another power of two.
henryzz is online now   Reply With Quote
Old 2012-03-15, 23:35   #6
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2·7·461 Posts
Default

OK, here's the analysis for all the 110-to-120 digit numbers in the aliquot sequences I have run. Hope it's vaguely interesting to someone.

Code:
1.0.0.0	(345)	0.841 1.0.0.0	
		0.070 1.0.0.1	
		0.041 1.0.1.0	
		0.014 1.0.0.2	
		0.009 1.0.1.1	
		*0.006 3.0.0.0*	
		*0.006 2.0.0.0*	
		0.006 1.0.2.0	
		0.003 1.0.4.0	
		0.003 1.0.3.0	
		*0.003 2.0.2.0*	
		
1.0.1.0	(73)	0.630 1.0.1.0	
		0.247 1.0.0.0	
		0.055 1.0.2.0	
		0.041 1.0.1.1	
		0.014 1.0.1.2	
		0.014 1.0.0.1	
		
1.1.0.0	(237)	0.844 1.1.0.0	
		0.038 1.2.0.0	
		0.038 1.1.0.1	
		0.034 1.1.1.0	
		0.017 1.3.0.0	
		0.008 1.1.2.0	
		0.008 1.1.0.2	
		0.004 1.2.3.0	
		0.004 1.6.0.0	
		0.004 1.2.0.1	
		
1.2.0.0	(172)	0.773 1.2.0.0	
		0.058 1.2.1.0	
		0.041 1.1.0.0	
		0.035 1.3.0.0	
		0.035 1.2.0.1	
		0.023 1.4.0.0	
		0.006 1.1.2.0	
		*0.006 2.1.0.0*	
		0.006 1.1.1.0	
		0.006 1.2.2.0	
		0.006 1.3.2.0	
		*0.006 4.1.1.0*	
		
1.3.0.0	(52)	0.635 1.3.0.0	
		0.135 1.2.0.0	
		0.077 1.4.0.0	
		0.077 1.1.0.0	
		0.038 1.3.0.1	
		0.038 1.5.0.0	
		
2.0.0.0	(740)	0.882 2.0.0.0	
		0.046 2.0.1.0	
		0.018 2.1.0.0	
		*0.008 1.0.0.0*	
		*0.008 3.0.0.0*	
		*0.005 4.0.0.0*	
		0.005 2.0.2.0	
		0.004 2.0.3.0	
		*0.003 3.2.0.0*	
		0.003 2.2.0.0	
		0.003 2.1.1.0	
		*0.001 4.0.1.0*	
		*0.001 7.0.0.0*	
		*0.001 3.0.1.0*	
		0.001 2.3.1.0	
		*0.001 6.0.0.0*	
		*0.001 5.0.1.0*	
		*0.001 7.0.1.0*	
		0.001 2.0.7.0	
		*0.001 7.1.0.0*	
		0.001 2.3.0.0	
		*0.001 1.0.3.0*	
		
2.0.0.1	(98)	0.847 2.0.0.1	
		0.082 2.0.1.1	
		0.051 2.0.0.2	
		0.010 2.0.3.1	
		0.010 2.0.1.2	
		
2.0.1.0	(121)	0.628 2.0.1.0	
		0.306 2.0.0.0	
		0.058 2.0.2.0	
		*0.008 3.0.0.0*	
		
2.1.0.0	(558)	0.828 2.1.0.0	
		0.057 2.0.0.0	
		0.050 2.2.0.0	
		0.034 2.1.1.0	
		0.009 2.3.0.0	
		0.009 2.1.2.0	
		0.005 2.4.0.0	
		0.004 2.0.1.0	
		0.004 2.2.1.0	
		
2.1.1.0	(95)	0.600 2.1.1.0	
		0.295 2.1.0.0	
		0.084 2.1.2.0	
		0.011 2.1.3.0	
		0.011 2.2.1.0	
		
2.2.0.0	(203)	0.635 2.2.0.0	
		0.084 2.1.0.0	
		0.059 2.2.1.0	
		0.039 2.0.0.0	
		0.039 2.3.0.0	
		0.025 2.2.2.0	
		0.025 2.5.0.0	
		0.020 2.1.1.0	
		0.010 2.3.1.0	
		0.010 2.4.0.0	
		*0.005 1.0.0.0*	
		0.005 2.2.3.0	
		0.005 2.0.1.0	
		*0.005 3.2.0.0*	
		*0.005 4.1.0.0*	
		*0.005 3.3.0.0*	
		*0.005 3.0.0.0*	
		0.005 2.5.1.0	
		*0.005 4.4.0.0*	
		0.005 2.6.0.0	
		0.005 2.0.2.0	
		
2.2.1.0	(54)	0.519 2.2.1.0	
		0.315 2.2.0.0	
		0.074 2.2.2.0	
		0.037 2.1.0.0	
		0.019 2.3.1.0	
		0.019 2.1.2.0	
		0.019 2.4.0.0	
		
3.0.0.0	(385)	0.844 3.0.0.0	
		0.075 3.0.0.1	
		*0.021 2.0.0.0*	
		*0.013 4.0.0.0*	
		0.010 3.0.0.2	
		*0.010 5.0.0.0*	
		*0.008 6.0.0.0*	
		*0.005 1.0.0.0*	
		*0.003 4.0.0.1*	
		*0.003 9.0.0.0*	
		0.003 3.0.0.3	
		*0.003 4.0.0.2*	
		0.003 3.0.0.4	
		
3.0.0.1	(76)	0.526 3.0.0.1	
		0.434 3.0.0.0	
		0.039 3.0.0.2	
		
3.1.0.0	(230)	0.865 3.1.0.0	
		0.061 3.1.0.1	
		0.022 3.1.0.2	
		0.017 3.2.0.0	
		*0.009 4.1.0.0*	
		0.009 3.3.0.0	
		*0.004 5.1.0.0*	
		*0.004 7.1.0.1*	
		0.004 3.4.0.0	
		0.004 3.4.0.1	
		
3.1.1.0	(91)	0.846 3.1.1.0	
		0.066 3.1.2.0	
		0.044 3.1.1.1	
		0.033 3.1.1.2	
		0.011 3.1.3.0	
		
3.2.0.0	(106)	0.717 3.2.0.0	
		0.085 3.3.0.0	
		0.047 3.1.0.0	
		0.019 3.5.0.0	
		*0.019 2.1.0.0*	
		0.019 3.2.0.1	
		*0.009 4.2.0.1*	
		*0.009 9.2.0.1*	
		0.009 3.2.0.2	
		0.009 3.4.0.0	
		0.009 3.1.0.4	
		*0.009 5.2.0.0*	
		*0.009 4.2.0.0*	
		*0.009 2.2.0.1*	
		0.009 3.2.0.4	
		*0.009 2.4.0.0*	
		
3.3.0.0	(61)	0.672 3.3.0.0	
		0.115 3.2.0.0	
		0.049 3.3.0.1	
		0.049 3.4.0.0	
		0.049 3.1.0.0	
		0.016 3.5.0.0	
		0.016 3.1.0.1	
		0.016 3.4.0.1	
		0.016 3.1.0.3	
		
4.0.0.0	(173)	0.734 4.0.0.0	
		0.046 4.0.0.1	
		*0.029 2.0.0.0*	
		0.029 4.1.0.0	
		*0.029 3.0.0.0*	
		*0.017 7.0.0.0*	
		*0.012 5.0.1.0*	
		0.012 4.0.2.0	
		*0.012 5.0.0.0*	
		0.006 4.0.0.4	
		0.006 4.0.3.0	
		0.006 4.0.1.0	
		0.006 4.0.2.1	
		0.006 4.0.1.2	
		0.006 4.0.0.2	
		*0.006 6.0.0.0*	
		0.006 4.1.0.1	
		*0.006 3.1.0.0*	
		*0.006 5.0.0.1*	
		*0.006 3.0.0.1*	
		0.006 4.4.0.0	
		*0.006 5.2.0.0*	
		*0.006 2.1.3.0*	
		
4.1.0.0	(121)	0.802 4.1.0.0	
		0.050 4.0.0.0	
		0.050 4.1.0.1	
		0.017 4.2.0.0	
		0.017 4.1.1.0	
		*0.008 5.1.0.0*	
		0.008 4.0.1.0	
		0.008 4.1.0.2	
		0.008 4.1.0.3	
		*0.008 6.1.0.0*	
		0.008 4.4.0.0	
		0.008 4.1.1.3	
		0.008 4.3.1.0	
		
5.0.0.0	(76)	0.724 5.0.0.0	
		*0.066 3.0.0.0*	
		*0.053 4.0.0.0*	
		0.039 5.0.1.0	
		*0.026 2.0.0.0*	
		*0.026 6.0.0.0*	
		*0.013 14.0.0.0*	
		*0.013 9.0.0.0*	
		*0.013 1.0.0.0*	
		*0.013 7.0.0.0*	
		*0.013 3.0.1.0*	
		
6.0.0.0	(109)	0.734 6.0.0.0	
		0.046 6.0.0.1	
		0.037 6.0.1.0	
		*0.028 3.0.0.0*	
		*0.018 4.0.0.0*	
		*0.018 2.0.0.0*	
		*0.018 5.0.0.0*	
		0.009 6.0.2.0	
		0.009 6.0.1.1	
		0.009 6.1.0.1	
		*0.009 7.0.0.0*	
		0.009 6.0.0.3	
		*0.009 5.0.0.1*	
		*0.009 8.0.0.1*	
		0.009 6.0.3.1	
		0.009 6.2.0.0	
		*0.009 10.1.0.0*	
		0.009 6.1.0.0
fivemack is offline   Reply With Quote
Old 2012-03-15, 23:38   #7
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2·7·461 Posts
Default

And here are the figures for 151-to-160-digit numbers

Code:
1.1.0.0	(71)	0.859 1.1.0.0	
		0.056 1.1.1.0	
		0.028 1.3.0.0	
		0.028 1.1.0.1	
		0.014 1.2.0.0	
		0.014 1.1.3.0	
		
1.2.0.0	(56)	0.786 1.2.0.0	
		0.071 1.2.0.1	
		0.054 1.3.0.0	
		0.036 1.1.0.0	
		0.018 1.2.1.0	
		0.018 1.4.0.0	
		0.018 1.2.1.1	
		
2.0.0.1	(53)	0.792 2.0.0.1	
		0.075 2.0.1.1	
		0.038 2.1.0.1	
		0.038 2.0.0.2	
		0.019 2.0.3.1	
		0.019 2.0.2.1	
		0.019 2.2.0.1	
		
2.1.0.1	(76)	0.789 2.1.0.1	
		0.053 2.1.0.2	
		0.053 2.0.0.1	
		0.026 2.2.0.1	
		0.026 2.1.1.1	
		0.013 2.1.2.1	
		0.013 2.4.0.1	
		0.013 2.3.0.1	
		0.013 2.1.1.2
(these are only the prefixes that occur 50 or more times)

Last fiddled with by fivemack on 2012-03-15 at 23:39
fivemack is offline   Reply With Quote
Old 2012-03-16, 00:15   #8
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dartmouth NS

845010 Posts
Default

Quote:
Originally Posted by fivemack View Post
A.B.C.D means 2^A*3^B*5^C*7^D

These are from 1250 lines from the 2^2*3 project

Code:
2.0.0.0		0.890 2.0.0.0	0.044 2.0.1.0	*0.015 3.0.0.0*	*0.015 1.0.0.0*
		0.011 2.0.2.0	0.007 2.1.0.0	*0.004 7.0.1.0*	0.004 2.0.7.0
		*0.004 4.0.0.0*	*0.004 3.2.0.0*	0.004 2.0.3.0	
2.1.0.0		0.832 2.1.0.0	0.052 2.0.0.0	0.047 2.1.1.0	0.044 2.2.0.0
		0.010 2.3.0.0	0.008 2.1.2.0	0.003 2.0.1.0	0.003 2.2.1.0
		0.003 2.4.0.0	
2.1.1.0		0.608 2.1.1.0	0.278 2.1.0.0	0.089 2.1.2.0	0.013 2.1.3.0
		0.013 2.2.1.0	
2.2.0.0		0.620 2.2.0.0	0.092 2.1.0.0	0.063 2.2.1.0	0.049 2.3.0.0
		0.042 2.0.0.0	0.021 2.2.2.0	0.021 2.1.1.0	0.021 2.5.0.0
		0.014 2.2.3.0	*0.007 1.0.0.0*	0.007 2.0.1.0	0.007 2.3.1.0
		0.007 2.4.0.0	*0.007 3.0.0.0*	0.007 2.5.1.0	*0.007 4.4.0.0*
		0.007 2.6.0.0
Read this as 'if you have currently a factor 2^2*3^2 (and are about 115 digits), probability is about 0.092 that you will be 2^2*3 next iteration'. The probabilities are taken from at least 50 (about 300 for 2.1.0.0) lines. Ones with stars are where the power of two changed.
one thing that came up for me when you sorted it vertically is that:

1.0.0.0 the probabilities for this add up to more than 1 from the addition I did by hand and in a check with PARI.
science_man_88 is offline   Reply With Quote
Old 2012-03-16, 00:36   #9
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

23×3×52×17 Posts
Default

I just recently saw a 2^2*3^2*P -> 2 transition in the wild (967170:i729) and was quite surprised. After a small sampling simulation, I saw that it was indeed quite rare, and your data shows the same (for sure):
2.2.0.0 -> *0.007 1.0.0.0*
I concur.
Batalov is offline   Reply With Quote
Old 2012-03-16, 08:49   #10
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

144668 Posts
Default

Quote:
1.0.0.0 the probabilities for this add up to more than 1
Yes, they do; they're all rounded independently and there are eleven of them. If you want to convert them to the integer numbers of 345ths that they represent and add them up, they will get to 345/345.
fivemack is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Caught in transition? cheesehead Forum Feedback 1 2011-12-11 05:14
How will the transition to v5 work? ixfd64 PrimeNet 3 2008-10-01 01:42
server transition news ltd Prime Sierpinski Project 4 2006-04-19 20:25
Gimps awaiting a major transition tha Hardware 42 2005-08-21 13:30

All times are UTC. The time now is 11:34.


Thu Jun 8 11:34:45 UTC 2023 up 294 days, 9:03, 0 users, load averages: 0.75, 0.94, 0.96

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2023, 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.

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