This site is supported by donations to The OEIS Foundation.

User:Bill McEachen

From OeisWiki

Jump to: navigation, search

The email shown in sequences is obsolete...use

married, 2 kids, control systems engineer, interest in prime numbers, born 1962 live in southwest Virginia USA (near Virginia Tech), BS,MSEE,PE



My Sequences

(I'm lax, I sometimes say integers for whole numbers ...)
A129912 June 2007 relates to my main PPP conjecture (primes as f{primorials,primorial products}. This is effectively my mild claim to fame, in that it works to explain the way that the primes are distributed. For clarity I will repeat my original conjecture here:

   Every prime number >2 must have an absolute distance to a sequence entry
(primorials, primorial products) that is itself prime, aside from the special
cases prime=2 and those primes immediately adjacent to a sequence entry
(primorials, primorial products). The property is required but not sufficient considers distances no larger than the candidate

Rephrasing, this merely means every odd prime number must either be adjacent to, or a prime distance away from a primorial or primorial product. (the distance will be a prime smaller than the candidate)

The usefulness of the property is evident from the Wikicommons image referenced in the sequence. The normalization is merely -1 to +1, where 0 indicates the closest PPP entry was very close, and +/- 1 indicates the closest PPP entry was as far away as possible. Note that the generation of the entries for A129912 appears related to A071562 in a slightly convoluted way. My website details more about the conjecture, including other independent work by Sokol and Potter. Note OEIS begins the sequence at 1, while I begin it at 2 (this seems to happen a lot, I have given up trying to make sense of it).

See comments on the talk page of the Primorial article here, where a poster says more that appropriate wording is that the offset must be coprime to all of the governing factors of the PPP entries.--Bill McEachen 21:39, 25 December 2010 (UTC).

Also, as a clarification, I cannot compromise in the logic about the 1st term should start at 2. By definition, the sequence treats prime factorials, which do not exist below the first prime, which is 2.

A117825 May 2006 relates to distance from the Highly Composite Numbers and the nearest prime. This relates to Goldbach's Conjecture, my HCN conjecture (see my website), and Fortune's conjecture. This has an interesting pinplot. I posted an image based on CRG data for first 20K here:svg image. Also mentioned in A228945 & A228943.

A147517 Nov 2008 relates to the symmetrical primes as a f{A002110, primorials}. Here symmetrical primes means the pair sum to twice the primorial center.

A077287 Aug 2003 concerns a mechanism that generates a lot of prime numbers with very volatile numeric spread (see the OEIS scatterplot). I derived while looking for good prime generators. Possible rifo A023057

A113972 Feb 2006 is a form I derived that yields an interesting sequence of primes, with many hi/lo transitions. The OEIS scatterplot here depicts it:[1]

A156053 Feb 2009 shows the Twin Primes (lower of pair) encountered by the Bernoulli Number denominators. It excludes duplicates, and covers through first the 15000 BN. I have an image posted here [2]. From Mr Cloitre the BN denom are produced via a(n)=prod(p=2, 2*n+1, if(isprime(p), if((2*n)%(p-1), 1, p), 1)) . Mr Karrhus and his new TP work (Mar 2013) is a bit similar, but he produces larger Twin Primes. He produces 27 TPs evaluating thru Prime ~ 14000. I will add its plot here shortly

A172069 Jan 2010 merely relates the adjacency of primes to entries of A129912. The OEIS plots here depict it:[3]. I must analyze its mean (0.5?) and (5) pair sum ratios later (skewed to the lower half).

A147853 Nov 2008 is a spinoff of A147517 and relates to Goldbach partitions of 2*primorials (A002110), providing a very good spread of Goldbach pairs

A102648 Feb 2005 is just a nice function to generate an interesting spread of values (see OEIS pinplot). It is a simple mechanism to quickly generate prime numbers from a fixed range of small integer seeds in a fairly arbitrary way. The submitted version generates values <=20 but this range is dependent on the gain (100) used. The OEIS plots here depict it: [4]. I have a rejected sequence that does it better.

A156703 Feb 2009 is a nice infinite, diverging sequence that is pretty volatile term-by-term. It is formed from ratios of the whole numbers. All but the 3rd entry appear to contain the digit 5 or 9. I posted an image here [5]. In Jan 2014 I tied this sequence to A158911, in turn connected to the Hamming Numbers. I also made a new conjecture about the sequence entries.

I acknowledge and appreciate a big Pari code assist on this from Michel Marcus, who pleasantly pointed out my inappropriate use of reals in place of integer math.--Bill McEachen 02:32, 31 December 2015 (UTC)

A005846 this relates to a classic eqn for generating primes. I have an image posted here [6]

A181616 produces volatile list of primes from a seed, with 2 distinct resulting limit lines for the log scatterplot.

A184938 this relates to bolt-tightening patterns but the reference was deleted when approved. That reference is 1

A218391 this relates to an observation by John Richardson on distinguishing composites from primes. On its face it would not help practically with primality testing, but I'm thinking.

A226095 this produces primes from the prime factorization of integers>1 (was composites). Good "multiplier" as far as size of primes ID'd vs base composite. Nice volatility, thanks to CRG once again for his code assist ...doubt it has any improved hit rate on encountering primes...thanks to TD Noe for his observations to fix my comments.

A228446 this relates to Sun's conjecture that every odd n>5 = p+x(x+1) where p is an odd prime, x an integer. Nice plot, generates some primes higher than largest pf of originating n. It has a beautiful patterned scatterplot.

I may submit a sequence based on Papalexis' data compression method (Jan 2014) (normalized R values as F(bit string length and which bits are set). I would like to submit a sequence w/ limited terms I worked up regarding maximal production of primes from quadratics.

I submitted a candidate sequence Jan 12 2014, a companion to my A156703. It produces highly volatile entries. It stems from repeating decimals of pairs of adjacent whole numbers. (tentative 235589). Entries seem to always contain a 9,8,or 6 (akin to A156703 entries having 9 or 5).

A240563. Its form can generate families of nicely spread prime numbers via concatenation.

I submitted sequence A242902 which generates primes from sets of twin prime sets. It has a "hit" rate of ~ 30% for the first 10K evaluations (prime range 3-2500000). One can compare this to primepi of 40% thru 10K (smaller primes than the former however). A244043 - Numbers n for new peaks of floor(sigma(n)/primepi(n)) A248863 - bit patterns from Pi approximations A253899 - a(n) is the least prime greater than a(n-1) that follows a gap of exactly 2*n A257495 - The number of iterations (x -> 2x+1) until a prime is found, starting with prime(n) A264030 - Primes p such that A203069(p) is also prime (A203069 is alternating-parity rearrangement of natural #'s) A264042 - Prime numbers adjacent to Catalan numbers A265436 - pair subset sorting by sum,product to guarantee monotonicity. One of my few joint submissions and the best example for me of collaboration, thanks to Mr Alcover and Marcus. Subseq linked by conjecture to A035106 (#'s of the form k*(k+1) or k*(k+2) )

I submitted A245694 which stems from a MMF post/conjecture. It may be a short sequence but seems to consist of prime squares.

I would delete my seq A087867 as useless.


Thanks to direction from CRG, I was able to provide the solution set to the generalized conjecture mentioned by Winkler in A051451. In a nutshell, the majority of solutions are the multiples of 6. The rest stem from A053319. The generalized conjecture is: There exists a class of (even) numbers such that each member lies between two pairs of Twin primes (TP1 thru TP4 in increasing order). A governing relation is each member = (TP3-TP1). An example is member 6 and (5,7,11,13). Members of the class include all multiples of 6 as well as some others related to A053319, as described there. As the multiples of 6 are infinite, this implies the Twin Primes are as well, in a backhanded way, which most mathematicians believe.

Minor observations

A few other sequences that I've interacted with:

A005097 I made connection to Goldbach partitions
A000079 I tied this to the divisors of the primorials (see lower)
A001359 I made connection to the generalized Winkler conjecture
A109852 I conjectured regarding every odd prime >4 being encountered in order
A079149 I noted the connection between entries and a sieving involving (a+b) and (a^2-b^2) for consecutive composites
A002110 I noted that the xth root of the xth primorial seems to have magnitude on the order of its # factors.
A002182 I noted a(n) often corresponds to P(n,m) = number of permutations of n things taken m at a time.
A158869 I asked whether the pattern of rightmost digits) (1,5,7,7) will continue?
A006562 I noted similarity of plot with A013916
A077287 Here is a plot for kissing number aligned to A129912 entries for correlation (I start the latter at entry=2): [7]
A000849 These entries highly correlate to A147517 as follows:
take adjacent differences of the A147517 entries, and plot the log of those values vs A000849.
The resulting linear fit is Y=Ax+B where Y is the A000849, X =A147517, (A,B)=(1.0451,0.241579).
The fit to the admittedly small dataset is better than 0.9999
For clarity, the X dataset is (5,24,160,1374,15511,...) and Y dataset (10,46,343,3248,42331,...)
The fit yields (9,48,351,3320,41804,...)
A139270 besides a(1), it appears to be a subset of A091375


I added a comment about the frequency of primes to the sequence entries. Specifically, the "hit" rate is higher than average, after applying PrimePi with adjust for odd, non multiples of 5 as the baseline. I did the analysis in Gnumeric and can supply that data.--Bill McEachen 06:37, 6 December 2010 (UTC)

I also compared Primepi vs that derived from this sequence, the ratio of which seems to be tending towards 1.5
Specifically, for xample, entry 17 is 132. Primepi(132)=32. rat1=32/132/0.4=0.61
The sum of adjacent primes thru the 17th entry is 13 (primes adjacent to each entry). this value 13/entry ptr=17 is 0.76. I was trying to see how relevant a forecaster for PrimePi the sequence is. I looked at the first 1500 entries only and here is the ratio plot:[8]
--Bill McEachen 04:08, 23 January 2011 (UTC)


I continue to find new tie-ins to OEIS for my work. I make use of the prime sets named at [9], and they all appear in OEIS, where anciliary notes are often quite useful.--Bill McEachen 06:40, 6 December 2010 (UTC)

Radix transformations

I entered a sequence (not accepted) relating this subject with a specialized format. Considering unsigned integers only, the form is xxx.yyy where yyy is the remainder

10. Now unlike hexadecimal or other bases where we use more than 10 numeric characters (like A thru F in hex), we we constrain ourselves to digits 0-9. Then, the sequence, interest xpresses the # bases each number can be xpressed in. The first decrease is at 110, which begins the coverage gap for base 11 between 11^2 and 9*11+10. The interest to me is the local minima seen in the sequence. I can provide Python code to produce the representation for any base (fails if in coverage gap).
The sequence begins 1,1,2,2,2,2,3,3,4 considering only primes radix. The sequence was not accepted into OEIS.

Named prime sets

I have a list of named prime sets: User:Bill McEachen/Named prime sets and User:Bill McEachen/OEIS_prime_sets.

Mersenne Numbers

I've submitted a sequence A186283 related to A186253. Very volatile but well-bounded at lower end. I know of the prime power theorem, but will the factor hold for all powers? I also submitted the dual case (2^n+1) as A185343. I thank CRG and TD Noe for their assistance.

Sequences of interest

I like A166448 because of its graph and concept captured. It follows a 2nd-degree polynomial very well. Thru 500 terms I got Y=3.698x^2-231.7x+7484, R^2=0.9999+

I updated the related A166449 with a b-file if you're interested. Charles R Greathouse IV 07:13, 24 February 2011 (UTC)

As CRG et al point out, there is an established power-ln relation involved.

A084317 (from concatenation of prime factorizations) its bfile is short but I did note that the entries of A002110 seem to appear (+1). So, we have 2,6,30,210,2310,30030 and the seq contains entries 3,7,31,211,2311. Will there often be such a correspondence? One ref for prime factorizations is 1. The next would be 30031, which won't be encountered in A084317. The rest wouldn't make much sense either now that I look at them ...oh well

A229159 from M. Lagneau is very interesting (Sept 2013)

Since you find this interesting I looked at it as well. It seems to me that using the n-th prime is artificial, and I wonder if there might be another interesting sequence if you take this away and use numbers or odd numbers in place of primes. Perhaps even a sequence for "numbers n such that there is a triangle with leg n, all legs integers, and integral area". (Feel free to submit any of these if you feel so inclined and don't find them in the OEIS.) But I haven't had time to look at this in sufficient detail. Charles R Greathouse IV 20:38, 25 September 2013 (UTC)
I actually emailed Mr Lagneau giving what I think proves his conjecture (must be straightfwd for me to be able to do it)--Bill McEachen 00:37, 15 November 2013 (UTC)


This was one of those strange "hits" one gets from fiddling with numbers. I was working with bit strings, considering bit transitions, # adjacent bits set from left etc. I concluded:

the final bound is (N-1 left aligned bits ) + 2^(N -#set bits-1) - A084639(ptr) with ptr=#bit transitions-1

This takes a bit of xplanation but crudely, for a(32),input to Wolfram Alpha is:

111111111111111111111111111111110 base 2 - 101010101010101010101010101010101 base 2

this gives 2863311529. We have 32-bit transitions in the right portion. We have N=33, and there are 17 set bits. This gives 2^(33-17-1)=32768. A(32)=2863311529. So, final bound ~ 5726655829. It computes to 5726623061.
I only now see CRG's Pari code for the seq, so it transforms the above "A084639(ptr)" to 2^(ptr+2)/3-if(ptr%2, 5, 4)/3--Bill McEachen 01:46, 9 November 2013 (UTC)


My fit to the curve is 1.06766 + 1.02336*ln(x-0.574)...

FWIW there should be about 1.37*sqrt(x)/log(x) members of that sequence up to x. Put another way, a(n) should be about 0.132 n^2 log n. - Charles R Greathouse IV 07:45, 3 November 2013 (UTC)
I forgot I had fit log(Y), so the fit is log(Y)=1.06766 + 1.02336*ln(x-0.574), where x is ptr, Y=a(ptr). For example, actual is (4600,4694305229). This gives (4600,9.67). Fit gives 9.70...If I treat your eq'n above even using ln, it estimates 23556648, pretty shy of actual.--Bill McEachen 00:37, 15 November 2013 (UTC)

A248863 (proposed)

this is merely a bit pattern that results from concatenation from A011545 terms. This relates to Pi. I must graph it against A004601...--Bill McEachen 02:25, 4 March 2015 (UTC)


This goes back to Oct 2008. It begins from pf=7 1,2,4,8,16. length(divisors(7))=2 length(divisors(11*7))=4 [1,7,11,77] length(divisors(13*11*7))=8 [1, 7, 11, 13, 77, 91, 143, 1001] length(divisors(17*13*11*7))=16 power of 2 should match #terms in Prod(P), p from prime(4) to prime(n) SO, use prime(8) to prime(4). divisors should = 2^5=32 Prime(20) should yield 2^17. It does (131072). It's obviously thinking about it, but any basic connection between (powers of 2) and (products of primes) is nice.


some basic functions that are often useful (see refcard for all, or use ? xxx where xxx is command):

  • print, write
  • for, while, if
  • List(), listput,listsort,setsearch
  • sumdigits, digits
  • isprime,ispseudoprime,nextprime,precprime,prime(*)
  • x%y (mod), abs, ++ (increment), -- (decrement), max, min, floor, ceil, round, truncate
  • random, length, Str (convert to string),
  • divisors, sigma, Pi, prodeuler (product over primes), numerator, denominator
  • issquare, ispower, fibonacci, sumformal

Idea for page

I think a page covering "popular sequence topics" and/or phrases would be useful. There are obvious ones; a very quick stab might include (counts as of May 2014):

  • keyword: easy (~55000)
  • keyword: base (~30000)
  • expansion (~22000)
  • "Vincenzo Librandi" (~21600)
  • series (~18600)
  • "Clark Kimberling" (~17500)(author)
  • binomial (~ 14700)
  • binary (~11000)
  • "decimal expansion" (~10500)
  • "Reinhard Zumkeller" (~10500)(author)
  • "T. D. Noe" (~10500)(author)
  • "Robert G. Wilson" (~10000)(author)
  • divisors (~9700)
  • partitions (~9500)
  • polymonial (~9500)
  • permutations (~9000)
  • Greathouse (~ 9000)(author)
  • transform (~8600)
  • lattice (~7000)
  • Fibonacci (~ 7000)
  • roots (~6800)
  • matrix (~6200)
  • "Alois P. Heinz" (~6200)
  • compositions (~6200)
  • "Jean-Fran├žois Alcover" (~5700)(author)
  • "Amarnath Murthy" (~5500)(author)
  • phi (~5100)
  • keyword:fini (~5100)
  • keyword: hard (~4800)
  • "Colin Barker" (~4800)(author)
  • "Benoit Cloitre" (~4600)(author)
  • Euler (~4200)
  • "numerators of" (~4000)
  • sorting (~4000)
  • nodes (~3900)
  • "Vladeta Jovovic" (~3900)(author)
  • "David W. Wilson" (~3900)(author)
  • "Roger Bagula" (~3800)(author)
  • "prime number" (~3800)
  • sigma (~3800)
  • convergent (~3700)
  • trees (~3500)
  • concatenation (~3400)
  • "Jonathan Vos Post" (~3300)(author)
  • "M. F. Hasler" (~3300)(author)
  • "denominators of" (~3200)
  • "Wolfdieter Lang" (~3100)(author)
  • "Michel Marcus" (~3100)(author)
  • sin( (~3000)
  • modulo (~2950)
  • "Joerg Arndt" (~2900)
  • factorial (~2800)
  • palindrome (~2750)
  • derivative (~2600)
  • asymptote (~2500)
  • "Ralf Stephan" (~2500)(author)
  • "number of ways" (~2500)
  • "triangular numbers" (~ 2400)
  • "first difference" (~2400)
  • intersections (~2300)
  • "smallest prime" (~2300)
  • primepi (~2100)
  • Pascal (~2050)
  • Ramanujan (~2000)
  • semiprime (~2000)
  • Lucas (~2000)
  • theta (~1800)
  • "Juri-Stepan Gerasimov" (~1800)(author)
  • "Twin prime" (~1625)
  • tau (~1500)
  • "Peter Luschny" (~1500)(author)
  • "Michel Lagneau" (~1400)(author)
  • "Vladimir Shevelev" (~1350)(author)
  • gap (~1300)
  • symmetry (~1250)
  • "Lekraj Beedassy" (~1200)(author)
  • "Alonso del Arte" (~1200)(author)
  • progression (~1200)
  • Bernoulli (~1100)
  • "binary representation" (~1050)
  • Chebyshev (~1000)
  • Dirichlet (~1000)
  • Fermat (~950)
  • packing (~850)
  • irreducible (~850)
  • Mersenne (~800)
  • harmonic (~780)
  • "Zhi-Wei Sun" (~721)(author)
  • "Don Reble" (~721)(author)
  • "Jonathan Sondow" (~700)(author)
  • "Gary Detlefs" (~700)(author)
  • "Jud McCranie" (~700)(author)
  • "first occurrence" (~700)
  • "characteristic polynomial" (~700)
  • Fourier (~700)
  • primorial (~600)
  • "primitive roots" (~550)
  • Erdos (~550)
  • Hardy (~500)
  • Riemann (~500)
  • Lehmer (~500)
  • chess (~400)
  • eigen (~400)
  • Honaker (~400)
  • error-corrected (~400)
  • "upper bound" (~400)
  • truncatable (~400)
  • Goldbach (~360)
  • Poincare (314)
  • "polygonal numbers" (~275)
  • lucky (~260)
  • "number of iterations" (~250)
  • Thue-Morse (~250)
  • cellular automaton rules (~ 230)
  • greedy (~220)
  • "highly composite number" (~175)
  • Galois (~ 167)
  • square-free (~115)


I made a comment stemming from an OEIS hit for some bitwise work I was doing. Completely unexpected.

Yet another surprise hit

working with some A/D frequency response stuff led me to A140777 and A165355 (I was aware of neither). I added comments to those.--Bill McEachen 23:40, 8 January 2015 (UTC)

  • also, in the course of my dallying, I looked at the specific approximation for ln(z) associated to the bilinear transformation (2 * (z − 1) / (z + 1). The crude text one is fine for z=0.7 to 1.0 but lousy for 0+ to 0.7. A cubic fit I did up (Excel) is superior and is 6.3853z313.476z2 + 10.481 * z − 3.3064. I use the std "z" associated with the discrete transfer function eqns. I now have my old A/D work online and available. I recently have come up with potential improvements for phase match.

Seq awaiting a next large term

Seq Gap


I submitted this as a variant it turns out of A000101...relates to prime gaps--Bill McEachen 15:42, 19 January 2015 (UTC)


A058043: Wolfram Alpha input: nextprime(n^2)-prime(PrimePi(n^2)), n=2 to 81

contributor spans (some may start earlier)

just something I was curious about ...

  • "Vincenzo Librandi" 2008-present
  • "Clark Kimberling" 2004-present
  • "Reinhard Zumkeller" 2002-present
  • "T. D. Noe" 2002-2014
  • "Robert G. Wilson"  ???-present
  • Greathouse 2005-preent
  • "Alois P. Heinz" 2008-present
  • "Jean-Fran├žois Alcover" 2011-present
  • "Amarnath Murthy" 2002-2010
  • "Colin Barker" 2007-present
  • "Benoit Cloitre" 2002-2014
  • "Vladeta Jovovic" 2000-2009
  • "David W. Wilson" 2000-2014
  • "Jonathan Vos Post" 2004-2014
  • "M. F. Hasler" 2006-2014
  • "Wolfdieter Lang" ???-present
  • "Michel Marcus" 2012-present
  • "Joerg Arndt" 2001-present
  • "Ralf Stephan" 2003-present
  • "Juri-Stepan Gerasimov" 2008-present
  • "Peter Luschny" 2005-present
  • "Michel Lagneau" 2010-present
  • "Vladimir Shevelev" 2007-present
  • "Lekraj Beedassy" 2000-2012
  • "Alonso del Arte" 2004-present
  • "Zhi-Wei Sun" 2008-present
  • "Don Reble" 2001-2014
  • "Jonathan Sondow" 2004-2014
  • "Gary Detlefs" 2000-2014
  • "Jud McCranie"  ???-2013

For comparison my span is 2003-present

A256620 (proposed)

relates to maximal runs of 6n+/-1 where primes are absent.--Bill McEachen 13:10, 5 April 2015 (UTC)

must have been rejected, as that sequence is something else.

Rejected sequences and comments

My true observation on the Thue-Morse sequence was rejected(A205083). The comment was: "The sequence is characterized by an expected mean of 0.5, though pair sums ratio to 1/2/2 for 0,1,2 respectively (unlike a coin flip, which yields 1/2/1)". Additionally, a simulated clone of Thue-Morse was rejected as well if I recall.

A sequence (w/ technique) to spawn nicely volatile prime sequences from other sparse sequences was rejected as being an artificial construct. The sequence spawned from A002315 is: 17,41,2,3, 91393, 811, 947, 3, 2, 127, 5, 8071607521, 9369319, 5, 46083933182810391855077841108121860076301803820136729604319921407582209931247732827597272276147556142385931557060724704331319480811439873947611787983922003 , 7,2, 47, 5, 919, 3, 489133, 2, 8287, 2 , 43, 7 ,2 ,7 , 9285087769350986448116616132878186749607968459, 19, 5 ,7 ,5 ,61, ... It works great on A053413 as well as others. The A053413 sequence spawned was 2,2,3,7,3,12413,1216552118791 ,2, 7 ,17 ,7, 6061207171, 2,7, 8793445674707408793600011419564463863171507,5,7,64081.

My comment on A151723 connecting it to 2 sequences was rejected (A197073, A022338). The 2 sequences characterize the three main "core" ON lines--Bill McEachen 13:10, 5 April 2015 (UTC)


I just saw Vladimir Shevelev's (VS)conjecture, which looked interesting. A249943 & A251621 have no Pari code. However, we do have code for A001223. So, restating VS: prime(n) = 19 + summation (A251621(i+4) ) i=9...n, n=8,9,...

or via VS's formula in A251621, prime(n) = 19 + summation (A001223(i-1) ) i=9...n, n=8,9,...

from the sequences, we see A001223(n) via forprime(p=1,50, print1(nextprime(p+1)-p, ", "))

The trick is to ensure we use no data "ahead" of prime(n). We know A001223(11)=6 from 37-31. This requires prime(12)=37. From VA, prime(12)=19+ A001223(8) + A001223(9) + A001223(10) + A001223(11)=19+4+6+2+6=37.

SO, it appears to me to require data "ahead" of what we seek, and so is not useful...

Let's try working from A098550 w/ RK's code (no explicit apriori primes required):

a(n, show=1, a=3, o=2, u=[])={n<3&&return(n); a+1==i&&print1("1, 2"); for(i=4, n, show&&print1(", "a); u=setunion(u, Set(a)); while(#u>1 && u[2]==u[1]+1, u=vecextract(u, "^1")); for(k=u[1]+1, 9e9, gcd(k, o)>1||next; setsearch(u, k)&&next; gcd(k, a)==1||next; o=a; a=k; break)); a}

A249943(n)=k makes use of A098550 thru k, where A098550 completes run to n. A251621 uses run lengths of A249943.

It's a bit complicated to follow, so we'll use an example. Akin to the example in A249943(7)=15,as A098550(15)=7 completes the coverage, so k=15, requiring thru A098550(15). However, to capture the run length requires thru A249943(9), which supplies A251621(7). But A249943(9) of course requires add'l A098550 entries, beginning a vicious cycle projecting forward...

Similarly, A249943(16), A098550(23)=13 completes the coverage, so k=13, requiring thru A098550(13). However, to capture the run length requires thru A249943(36), which supplies A251621(16).

So, there is some variable ptr> base n required from A098550, requiring data "ahead" of what we seek, and so again I conclude it is not useful...


more to come ....--Bill McEachen 14:58, 19 April 2015 (UTC)


after 100K, we need to see 19051 at some point (permutation...) the timing on the Pari script ~ terms^2.3 (rough) generally, the largest difference betw term and prime seen is 58207 (@ n=95016 which is encouraging). 19051 will require diff>~81K

Integer Notation

This just covers my approach for some of my work. I see Mr MF Hasler covers some of the same ground in his article here:M. F. Hasler/Representing large digits.

For non pure integer notation like decimal (or any radix I suppose), we have xxx.yyy where the power=0 begins to left of the decimal point and increases moving left. Similarly, the powers decrease to right of the decimal point from the inverse of power=1,2,....

I use a form www.zzz for integers of any base, with the imposed restriction of using digits 0-9 only. Power=1 begins to left of the decimal point, and any remainder <base to the right of same. One uses additive combinations as required (examples follow).

From hex, say #00e3 =227 base 10. As defined, this falls in the "gap" where a single notation is not possible for base 16, specifically values 160 through 256, since the notation goes from 9.15 to 10.0 (159 to 256). In my notation, 227 would be (9.15+4.4) (expanded 9*16^1+15 + 4*16^1+4 = 159+68=227 ).

Of course one clarification is that the base must be known/clear, as (9.15+4.4) is valid in any other radix>15 as well. Another nuance is there are multiple representations available, another could be 256-29= (10.0-1.13). This can be a plus or minus depending on ones requirements.

An advantage is that one can effectively limit the notational representation by merely increasing the base ( I gloss over some crucial intricacies here). A trival example might be moving from 9999999999999999.0 in some base, to merely just 1.0 in a much higher base. I have not, and likely could not, mathematically prove every value can be represented for every radix but I think this is the case intuitively. When using VERY large bases one can see the desire for the restriction to digits 0-9.


As often happens, I derived this sequence coming from a different direction (then found the correlation). My sequence contained all members, but in a different ordering, starting from binary value, and treating bit counts as the prime signature exponents. I have the Python code.

Cantor pair function

I submitted a concatenated sequence for this which was not favored. It nicely encapsulated the relevant data in each entry, ordered by CP function value increasing. As it stands now, one has to go to 3 separate sequences and merge them for the same information.


Another one of the interesting observations from just fiddling around. I looked at the #digits for the sequence entries, and the first 13 returned a lone hit against A092777. I know nothing of this sequence. The match stops, but I wonder if the A092777 entries form an upper bound? For example, A240563(50) has 207 digits, while A092777(50)=208. My highest sequence entry is 100, with 494 digits, while A092777(100)=505. This begs the question why such a concatenated construct adhere to such a prescribed bound? I can confirm the bounding does not work for all starting seeds -for example 3, the bounding is violated by the 24th iteration. However, the bound seems to apply after some max iteration (for 3, perhaps iteration 53).

Presuming the sequence doesn't stop, what the above means is one knows directly the # iterations (this sequence's procedure) required to get a certain size prime, as follows:

  1. iterations ~ primepi(desired size)

So, if we want a 545-digit prime, # iterations ~ 100. Of course the time becomes significant. As example, from iteration 99 to 100 on my laptop took ~ 4.5 minutes.


This link has a sequence Drexel, with Mathematica code. It claims to be the entries for sum of consecutive primes beginning from n that equal the product of n distinct primes. It has 47 entries, which grow very large. Searching for "331,89,5297,11149,2960267" returned no hit.

more phrases

"distinct primes" returns 2109 hits; "consecutive primes" 1627 hits; prime AND concatenation returns 1079 hits;

Sloane triples

concatenated sequences (A262577/ A262579/ A262581) beg the question about sparsity (surely more primes exist going forward). We know only 40% of the entries are odd so over 5000 entries, we have 2000 odds, and yet we only encounter 2 or 3 primes. It is as if our build structure is at odds with that underlying the primes. Weird.

submitted sequence early Oct 2015

I don't know if it will be accepted, but I entered my first sequence in combinatorics. It sort of relates to Euler's "sets with distinct subset sums".

followup is the sequence was not accepted. However another OEIS user spent the time to understand it, and in turn generated a spawned sequence that actually links to A183867. So, combinatorics was tied to n+floor(2*sqrt(n)), an interesting surprise.--Bill McEachen 00:43, 9 December 2015 (UTC)
ok, the new sequence (A265436) was accepted, and with yet another assist, it is now linked to A183867, A028387 and A035106. The last connection in particular allows for much simpler generation of the sequence. This is by far the best example for me of intertwined sequences and collaboration. For reference, the (4) sequence names are:

A183867 n + floor(2*sqrt(n)) A028387 n + (n+1)^2 A035106 Numbers of the form k*(k+1) or k*(k+2), k>0 A265436 the least m (1 <= m <=n) such that the set of pairs (x, y) of distinct terms from [m, n] can be ordered in such a way that the corresponding sums (x+y) and products (x*y) are monotonic --Bill McEachen 13:41, 26 December 2015 (UTC)


I took another look at this, seeing CRG's bfile. I was curious as to any limit applicable to the moving average. Overall after 10K terms is 0.457. After about 360K, we see 0.357. It likely is slowly going to 0 lacking other info but very slowly. After 1.8 million, avg = 0.326. After 3.4, 0.315 (21:17) gp > genit(500000000) Done 0.3259 iterations = 1840168

Primitive Pythagorean triangles

I looked that this (the areas) and the overlap for from the bfiles between A024406 &A129912 is: 6,30,60,180,210,2310,4620,60060,510510,10810800,116396280. My question is will the overlap always continue ? This would be a link between primorials and Pythagorean triangles essentially. I used LibreOffice Calc and Python.

Core Sequences

Just mentioning ones I certainly think are core (personal perspective): primes (A000040), Fibonacci(A000045),primorials(A002110), squares (A000290), semiprimes(A001358), Pascals (which??), A001222 (# of prime divisors), primepi (A000720), Mersenne exponents (A000043), Twin Primes (A001097), Goldbachs (which?).


This plot is obviously very linear (in log). (a,b)= ( 0.11348, -1.6691) Y=power of 10, x=n (note no log). I used the site [10].


I just saw this and will try to extend it a bit. At quick glance will every term n>4 have a 0,7 or 9 digit? The next 2 terms are 783414130, 14428317.--Bill McEachen 02:03, 31 December 2015 (UTC)


I just saw this. I checked what I knew of a conjecture made by John Sokol, and I believe it precedes the sequence entry date by 2 years. It seems to have the same conjecture. The Sokol link is available in A129912 but is here: [11]


I was led to this from my 2009 sequence. I connected it to A254528. Which of course made me think the connections between various sequences are WAY more than are recognized, only natural I suppose. I always take satisfaction from connecting sequences. I also encouraged MMarcus to add his Pari code adapted from A265xxx. --Bill McEachen 02:03, 31 December 2015 (UTC)


I am investigating a probabilistic connection to the primes on this one. Very computationally intensive, which curtails what I conclude and iterate. This may replace a new submission relates to randomly seeking prime numbers of x-digits. Effectively it provides the answer to how many multiples of primepi are required to generate all x-digit primes. The code is "smart" in that is only treats candidates ending in 1,3,7 or 9.--Bill McEachen 02:03, 31 December 2015 (UTC)


I just saw this, and am running some Pari code to investigate.--Bill McEachen 02:03, 31 December 2015 (UTC)


I commented on the range of a(n) seen vs log10(n), using a small set of powers of 10. I also added the bfile for 10K terms.


I added this and later someone connected it to A008706. This links primitive pythagorean triples, primorial products, and tilings (planar nets).


Fiddled with a few:
A002110 vs A129912 (logs, bullet3). Primorials vs PPP. A clear relationship as expected. The image is here: Media:Log10_A002110_A129912.PNG
A000045 vs A256454 (logs, bullet3) Fibonacci vs prime gap. Good correlation, mostly near-linear, roughly Y~0.16x-4.5. (y,x) are (ratio,n). The image is here: Media:Log10_A000045_A256454.PNG
A000045 vs A000041 (logs, bullet1) Fibonacci vs Partition numbers. Clear relationship, the image here: Media:Log10_A000045_A000041.PNG. I messed about some and using Pari functions, we have a VERY slowing decreasing value of sqrt(ln(A000045)) /ln(A000041). At 50K,70K,90K I see values 0.2766,0.2758,0.2752

well CRG enlightened me that it surely converges to a known value based on work of Hardy-Ramanujan and another fella...

A000108 vs A002110 (logs, bullet1) Catalan numbers vs Primorials. Superb relationship, the image here: Media:Log10_A000108_A002110.PNG. I didn't do a detailed fit but it appears Y~ ax^b where (a,b) ~ (3?,1?).
A147517 vs A000045 (logs, bullet3) (primes centered around primorials) vs Fibonacci. Roughly, Y~0.078x^1.88 but only 8 data points. (y,x) are (ratio,n)
A200474 vs A117825 (logs, bullet1) just nice looking Caldwell's conjecture vs (dist HCN to nearest prime)
A033932 vs A117825 (logs, bullet2) Factorial/prime offset vs HCN/prime offset. The pleasing image is here: Media:Log10_A033932_A117825.PNG
A033932 vs A147853 (logs, bullet2) Factorial/prime offset vs (primorial as prime average). The image is here: Media:Log10_A033932_A147853.PNG
A001043 vs A002182 (logs, bullet1) Numbers=the sum of 2 successive primes vs Highly composite numbers. The image will be here: Media:Log10_A001043_A002182.PNG
A001043 vs A129912 (logs, bullet3) Numbers=the sum of 2 successive primes vs PPP. The image will be here: pending


I just saw this looking at core sequences. Many terms intersect with A129912.


An interesting sequence that indicates 2 integers a,b always yield a prime a^2+b, with every nonnegative integer involved in a solution (conjectured). For example a(8)=9. b=16 to give 9^2+16=97 which is prime. One question is the extremes of n vs a(n) meaning how far is an integer from its base posn? For example, 100 is seen @ a(77) but I more care where n>a(n), records for. The max is 579 in the first 50K entries (@ ~ n=32408). Every prime is encountered if the density ofa^2+b matches that of the primes ie the PNT.


My 2010 comment detailed. Take a(n)=9699690. Factor(9699690) returns 8. (9699690)^(1/8)=7.47.
Taking 32589158477190044730, the 16th power = 16.58. Via WA,"factors of 32589158477190044730" we see 16 distinct prime factors.
Again, for 7858321551080267055879090 the 19th power = 20.43. It has 19 distinct pf. Etc


I conversed with [[User:Daniel Forgues | Daniel Forgues ] about my subconjecture: "every sequence will encounter a sequence value lower than itself". As the only implied condition is that we successfully evaluated (confirmed they reached 1) starts ascending from the lowest, this by logic leads to every (higher start) sequence as satisfying Collatz by default. A002450 essentially includes the odd values that would have to be encountered for any start != 2^n. Now, treating the iteration count as described totally changes input into sequences like A082984. Again as coded, average iteration for the first 100K starts is a mere 2+.

thanks to Michel Marcus for valuable assistance in a connecting 3x+1 to (2) sequences A167135&A068228. This fell out from dabbling and searching for a resultant sequence. Mr Marcus also helped with the proper wording for the comment, not my strong suit.


This is a great example to illustrate the difference between and advantage of harmonic mean over arithmetic. Treat its bfile data. The harmonic mean is essentially linear with superb fit.

A Primality Note

I just came across this link from 2012 ( ). It got me thinking as follows.

I do NOT offer this a a competitive method, merely a method.

Besides 2&3, we know prime numbers reside @ 6n+/-1. SO, we can combine the odd composite form Barker shows with 6n+/-1 as follows.

Barker's form is every odd composite stems from 2x +4y(x+y)-1 where x,y>0 and integers.

Case A (6n-1, "lower") 6n-1 NOT within 2x +4y(x+y)-1 reducing to 3n NOT within the solution set x+2y(x+y) [ meaning x+2y(x+y)-3n=0 has NO solution ]
Case B (6n+1, "upper") 6n+1 NOT within 2x +4y(x+y)-1 reducing to 3n+1 NOT within the solution set x+2y(x+y) [ meaning x+2y(x+y)-3n-1=0 has NO solution ]
I had no luck using WolframAlpha, but Dario Alpern's site can be used, noting that he uses the integer domain, not just positive integers.
Numeric examples are 35,37 and 49. The form for Alpern's site is ax2 + bxy + cy2 + dx + ey + f = 0 Note for our specific forms we only need the "hyperbolic" method discussed by Alpern.
35 (n=6) is case A, with (a,b,c,d,e,f) = (0,2,2,1,0,-18) No solution for the equality means solutions exist for the inequality. A valid solution is returned (2,2) and thus 35 is composite.
37 (n=6) is case B, with (a,b,c,d,e,f) = (0,2,2,1,0,-19) There are no solutions returned in the positive integers, meaning 37 is prime.
49 (n=8) is case B, with (a,b,c,d,e,f) = (0,2,2,1,0,-25). A valid solution is returned (1,3) and thus 49 is composite.
Others: 25 (n=4) caseB composite
55 (n=9) caseB composite
65 (n=11) caseA composite

I hope to code this at some point, likely in C++ or Pari/GP. --Bill McEachen 23:36, 26 June 2016 (UTC)

(Update1) ok, so it's coded and performs as can be seen below. It is very fast on prime numbers, while on composite factoring results are less clear. I merely took as baseline Pari's built-in function "isprime".

Candidate new algorithm vs Pari

300-hundred digit [4] 108s vs 26403s (2902...9329 last shown at )

1+10^10000 11s vs 53143s

2^23209-1 31.6s vs ( killed after 11min 10s )

(10^71 - 1)/9 11min 36sec vs 11min 40sec

3^2833 - 2^2833 1.2sec vs 2 hrs

(2^1084+1)/17*4337*13009 has run for days w/o yet completing

2^859433 -1 aka M33 27hr 22min vs TBD

--Bill McEachen 00:05, 8 July 2016 (UTC)

First Paper

My first math paper submitted for publication is still under review as of Dec 2016. As an aside, early on I received an email that to me represents much that is wrong with some in the math community. It was from a very well-known, established person (Phd). It essentially tersely said "I doubt what you claim is true" despite my having provided no details as to my method. Suffice to say the method is solid and works extremely well. It is an original deterministic primality test, though I make clear it has use for factoring otherwise. A bonus is the method can supply the Primality certificate. I am trying to avoid having to merely "dump" the paper on one of the open sites, but will as a last resort.--Bill McEachen 01:55, 21 December 2016 (UTC)

== ?? ==

Personal tools