This site is supported by donations to The OEIS Foundation.
Calendar for Sequence of the Day in July
Calendar for Sequence of the Day in June * Calendar for Sequence of the Day in August
Template:Sequence of the Day for July 1 A070308: Numbers such that the sum of each base 10 digit squared adds up to the sum of the nontrivial divisors.
If we take the base 10 digits of these "Canada perfect numbers," square each of them, and add them up, we find that they equal the sum of their nontrivial divisors (excluding 1 and the number itself). For example, with 125, we have . The divisors of 125 are 1, 5, 25, and 125, and we see that 5 and 25 ad up to 30. Prime numbers are excluded from consideration since their nontrivial divisors add up to 0. These numbers were defined by mathematicians at the University of Manitoba to celebrate Canada's 125th anniversary.
|
Template:Sequence of the Day for July 2 A024702: , where is a prime congruent to –1 or +1 (mod 6).
For all numbers and congruent to –1 or +1 (mod 6) we have is divisible by 24, so this applies for the unit 1 or any primes and other than 2 or 3. Iterate through all primes other than 2 and 3 and let be 1. This is a subsequence of the generalized pentagonal numbers, A001318.
|
Template:Sequence of the Day for July 3 A020639 lpf(n): least prime dividing n (a(1) = 1)
I'm listing the elements without the initial 1. I'm a purist, and given that 1 has no prime factors, I feel a(1) is undefined and should not appear in the sequence. I'm sure I (David W. Wilson) did not include it when I authored the sequence. Why do I find this sequence SOTD-worthy? It is because of a question I posed back in 1997. I noticed that in this sequence, between every pair of 2's, there was an element > 2, easy enough to prove. Similarly, I could show that between every pair of 3's there is an element > 3. I found the same to be true of 5, 7, 11, and 13. I began to think it might true for every prime number, but I could not find a general proof. So I posed my question the seqfan list. Apparently it was interesting, because it sent John Conway, Johan de Jong, Derek Smith, and Manjul Bhargava to the blackboard find a counterexample. After two hours, they found one:
Later, Fred Helenius found that the smallest prime that generates a counterexample is 71, for which
He also found the earliest known counterexample,
I was rather amazed that you had to look so far into the sequence.
|
Template:Sequence of the Day for July 4 A123456: Sequence name
Paragraph or two of info.
|
Template:Sequence of the Day for July 5 A177854: Smallest prime of rank .
This is a re-imagination of the Erdős-Selfridge classification of primes. 2 is rank 0 by definition, and for other primes the rank is the lesser of the maximum rank of the primes dividing and the maximum rank of the primes dividing . There are many interesting questions that could be asked about this sequence. What are , , and ? The sequence is infinite: how closely can its growth be bounded? (Look at the logarithmic graph, it seems to follow quite close to a regular pattern) How does the closely-related sequence A169818 behave around random large integers? (This last question is related to the performance of modern algorithms!) Trivial bounds: If and are the smallest primes of rank , then the smallest prime of rank is at least and at most where is Linnik's constant. These give rise to bounds like and , respectively.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Template:Sequence of the Day for July 6 A036236: Least inverse of A015910: smallest integer such that or 0 if no such exists.
Fermat's little theorem says that for odd prime , and on some random day long ago I got to wondering what values might take for other . I did some experimenting and found that that for small , took on many small values, but 1 and 3 remained elusive. I asked about them on the seqfan list, and found that is provably insoluble, while D. H. Lehmer had found the impressively large minimal solution = 4700063497 for , an impressively large solution for such a simple identity. Joe Crump has done some amazing work with this sequence, see 2^n mod n = c.
|
Template:Sequence of the Day for July 7 A000032: Lucas numbers (beginning with 2):
John Blythe Dobson connected some dots regarding Lucas and Fibonacci numbers: The comments submitted by Miklos Kristof from Mar 19, 2007 for the Fibonacci numbers (A000045) contain four important identities which have close analogues in the Lucas numbers (A000204): For and odd , . For and even , . For and odd , . For and even , .
|
Template:Sequence of the Day for July 8 A181780: Numbers which are Fermat pseudoprimes to some base , with .
For example 15 is Fermat pseudoprime to every base of the form and with . The sequence shows that the distribution of the Fermat pseudoprimes is the opposite of rare, if it is not determined on a fixed base . There are 429 Fermat pseudoprimes below 1000, but only 168 prime numbers. Fermat pseudoprimes to all bases , with are called Carmichael numbers.
|
Template:Sequence of the Day for July 9 A156695: Odd numbers which are not of the form , with , and prime.
Can you find a member (> 5) which is not divisible by 255? There should be infinitely many, but none have been found!
|
Template:Sequence of the Day for July 10
Simon Plouffe published a helpful list of odd zeta values from ζ(9) up to ζ(2051) with a precision of 1000 decimal digits — an older version of this list also covered ζ(3), ζ(5) and ζ(7) and is now available in separate files with more digits. You could hardwire these values in projects (e.g. RxShell (ZETA and RX.MATH)) where computing small odd
|
Template:Sequence of the Day for July 11 A004022: Repunit primes in base 10
These primes are of the form . The exponent must itself also be prime. Those three repunit primes have 2, 19 and 23 digits, respectively. The next repunit prime has 317 digits!
|
Template:Sequence of the Day for July 12 A097478: Atomic numbers of halogens in the periodic table
This sequence demonstrates that there is a mathematical structure to physical matter. We can compute this sequence mathematically by taking the formula for magic numbers, which happen to be the atomic numbers of the noble gases and then subtracting 1. The outermost shell of a noble gas atom has as many electrons as it can possibly have. Halogen atoms, on the other hand, can admit precisely one more electron to their outer shell (specifically in their outer p subshell.) Note that as you go down the periodic table the elements tend to acquire more metallic properties, which happens to be the case of astatine (85), which has somewhat more metallic properties than iodine (53), although astatine is still considered a halogen.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Template:Sequence of the Day for July 13 A118372: -perfect numbers
Also called Granville numbers after the mathematician Andrew Granville. In October 2010, William Rex Marshall offered a much clearer explanation† of these numbers as well as a few more terms to list in the sequence (which includes the familiar 2-perfect numbers of A000396).
|
Template:Sequence of the Day for July 14 A106108: Rowland's sequence: , then .
The difference between two consecutive terms is a noncomposite number (either 1 or a prime number.)
|
Template:Sequence of the Day for July 15 A001783: -phi-torial, or phi-torial of is the product of all positive integers up to and coprime to . Here means that and are relatively prime and is the GCD of and
The phi-torial of is a divisor of the factorial of . We take the positive integers below , cull out those that have prime factors in common with and then multiply the residual together. For example, with = 8, we cull out 2, 4, 6, 8, and multiply 1, 3, 5, 7, giving 105. Of course when is prime this works out to .
|
Template:Sequence of the Day for July 16 A004215: Numbers that are the sum of 4 but no fewer nonzero squares.
From Waring's problem we know that all positive integers can be expressed as a sum of at most four squares. In fact, most numbers can be expressed as the sum of fewer than four squares, even though the perfect squares thin out quadratically.
|
Template:Sequence of the Day for July 17 A083216: Fibonacci-like primefree sequence
Like the sequence of Fibonacci numbers, this sequence has two initial terms set, and the rest are given by the familiar recurrence relation . The amazing thing about this sequence is that there are no prime numbers among any of its terms even though the initial terms are coprime. And whereas the search for prime numbers seeks larger and larger primes, the search for primefree sequences seeks smaller initial terms. When in 1964, Ronald Graham (of Graham's number fame) proved that this kind of sequence is possible, his example had initial terms with more than thirty digits each (in base 10). Donald Knuth found a 17-digit pair in 1990, and later that same year Herbert Wilf found this sequence with slightly smaller initial terms. The current record is a 12-digit pair found by John W. Nicol in 1999.
|
Template:Sequence of the Day for July 18 A002372: Number of decompositions of into ordered sums of two odd primes
If the Goldbach conjecture is true, then for all . Note that this sequence does not require the primes to be distinct to be counted, and counts solutions of different primes twice (for example, for 8, both 3 + 5 and 5 + 3 are counted). Also note in the b-file that above about 3000, all values of this sequence are over a hundred.
|
Template:Sequence of the Day for July 19
If instead of the first seven terms I had given the first five, you might have guessed that the next five terms are 13, 17, 19, 23, 29. That is, unless you know by heart that the sequence of prime numbers has a much lower A-number in the OEIS, namely A000040. Rainer Rosenthal "deliberately contrived" this sequence "to begin with [the] first five primes" in order to illustrate the "absurdity of many 'guess the next term' puzzles." I think this sequence also makes two other important points. First, it shows the futility of searching for an arithmetic formula that will give the prime numbers in order, since such a level of complication is required for a formula that gives just the first five primes in order. Second, it shows the usefulness of the OEIS in finding formulae when one only knows a few terms of the sequence.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Template:Sequence of the Day for July 20 A007955: Product of divisors of where is the Iverson bracket.
In Only Problems, Not Solutions!, Florentin Smarandache lists this sequence (page 24), describes it, and urges the reader to study it, saying nothing else. It is easy enough to see that each prime number will occur once and only once in this sequence, namely, at . T. D. Noe goes much further than that and actually proves that each term of this sequence occurs just once. Of course not every positive integer occurs in this sequence (e.g., squares of primes). The product of divisors of is sometimes called the divisorial of .
|
Template:Sequence of the Day for July 21 A175607 Largest number such that the greatest prime factor of is
It surprised me, but after a little reflection, it makes a lot of sense: for any prime , there is a short (that is to say, finite) list of numbers such that has as its largest prime factor (Cf. Largest prime dividing n^2 - 1, n ≥ 2.) For example, for , we find that . No larger number satisfies . Even-indexed powers of two are squares of a smaller power of two, and therefore is not an integer if is even. The reason why higher odd-indexed powers of two don't work for this purpose is, as so many math textbooks say, left as an exercise for the Reader. According to Artur Jasinski, every prime has a corresponding with the property described here.
|
Template:Sequence of the Day for July 22 A038770 Numbers divisible by at least one of their base 10 digits.
I choose to start the listing at 20 because 1-digit numbers are trivially members of this sequence, while 2-digit numbers with a 1 among their digits are also trivially in this sequence. The choice of base 10 here encourages multiples of 5 but discourages odd multiples of 3 or 7. Regardless of the base, most prime numbers are excluded if in that base they are more than one digit and none of their digits in that base are 1.
|
Template:Sequence of the Day for July 23 A000945: The Euclid-Mullin sequence: , is the smallest prime factor of
The first four terms coincide with the Sylvester sequence (where each term is the product of all preceding terms plus 1) (A000058), an (thus strictly increasing) infinite coprime sequence whose terms may be prime or composite.
|
Template:Sequence of the Day for July 24 A115368: Decimal expansion of first zero of the Bessel function
The most famous zeros of a function are of course the nontrivial zeros of the Riemann zeta function. But other functions also have interesting zeros which mathematicians have studied, the Bessel function being one of them. (Mathematicians have also studied the zeros of the Bessel functions with other integer values of .) See Weisstein, Eric W., Bessel Function of the First Kind, from MathWorld—A Wolfram Web Resource. [http://mathworld.wolfram.com/BesselFunctionoftheFirstKind.html]
|
Template:Sequence of the Day for July 25 A000058: Sylvester's sequence, which gives the denominators for the greedy Egyptian representation of 1: (Of course the numerators are given by A000012.) The denominators are given by the quadratic recurrence and by the formula (which shows that it is an infinite coprime sequence) where for we get 1 + (empty product, i.e. 1) = 2.
|
Template:Sequence of the Day for July 26 A034872: Central column of Lozanić's triangle
Note that although the values of odd-numbered columns are only listed once, they appear twice in the triangle.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Template:Sequence of the Day for July 27 A008892: Aliquot sequence starting at 276.
|
Template:Sequence of the Day for July 28
In The Book of Numbers, John Conway and Richard Guy coyly describe this sequence as "not very different from the powers of 2."
|
Template:Sequence of the Day for July 29 A029578: An obvious mixture of two sequences
Can you tell what the two sequences being mixed are?
|
Template:Sequence of the Day for July 30 A105417: Numbers that are pandigital in Roman numerals, using each of the symbols I, V, X, L, C, D and M at least once.
We consider the largest member of this sequence to be MMMDCCCLXXXVIII, since none of the many systems to cope with greater numbers became standard prior to Roman numerals being displaced by modern numerals. There are 192 = 3 * 43 pandigital Roman numerals
Note that the Roman numeral representation of numbers is not admissible in the OEIS, thus the decimal representation of those numbers is used in A105417.
|
Template:Sequence of the Day for July 31 A173419: Length of shortest computation yielding using addition, subtraction and multiplication.
So you start with and then through successive , , or operations for some , you get to . The fun is of course getting there in as few steps as possible. Prime are not necessarily more difficult to get to than their neightbors: 29, for example, requires five operations, just as 28 and 30 do.
|