cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A284597 a(n) is the least number that begins a run of exactly n consecutive numbers with a nondecreasing number of divisors, or -1 if no such number exists.

Original entry on oeis.org

46, 5, 43, 1, 1613, 241, 17011, 12853, 234613, 376741, 78312721, 125938261, 4019167441, 16586155153, 35237422882, 1296230533473, 42301168491121, 61118966262061
Offset: 1

Views

Author

Fred Schneider, Mar 29 2017

Keywords

Comments

The words "begins" and "exactly" in the definition are crucial. The initial values of tau (number of divisors function, A000005) can be partitioned into nondecreasing runs as follows: {1, 2, 2, 3}, {2, 4}, {2, 4}, {3, 4}, {2, 6}, {2, 4, 4, 5}, {2, 6}, {2, 6}, {4, 4}, {2, 8}, {3, 4, 4, 6}, {2, 8}, {2, 6}, {4, 4, 4, 9}, {2, 4, 4, 8}, {2, 8}, {2, 6, 6}, {4}, {2, 10}, ... From this we can see that a(1) = 46 (the first singleton), a(2)=5 (the first pair), a(3)=43 (the first triple), a(4)=1, etc. - Bill McEachen and Giovanni Resta, Apr 26 2017. (see also A303577 and A303578 - N. J. A. Sloane, Apr 29 2018)
Initial values computed with a brute force C++ program.
It seems very likely that one can always find a(n) and that we never need to take a(n) = -1. But this is at present only a conjecture. - N. J. A. Sloane, May 04 2017
Conjecture follows from Dickson's conjecture (see link). - Robert Israel, Mar 30 2020
If a(n) > 1, then A013632(a(n)) >= n. Might be useful to help speed up brute force search. - Chai Wah Wu, May 04 2017
The analog sequence for sigma (sum of divisors) instead of tau (number of divisors) is A285893 (see also A028965). - M. F. Hasler, May 06 2017
a(n) > 3.37*10^14 for n > 18. - Robert Gerbicz, May 14 2017

Examples

			241 = 241^1 => 2 divisors
242 = 2^1 * 11^2 => 6 divisors
243 = 3^5 => 6 divisors
244 = 2^2 * 61^1 => 6 divisors
245 = 5^1 * 7^2 => 6 divisors
246 = 2^1 * 3^1 * 41^1 => 8 divisors
247 = 13^1 * 19^1 => 4 divisors
So, 247 breaks the chain. 241 is the lowest number that is the beginning of exactly 6 consecutive numbers with a nondecreasing number of divisors. So it is the 6th term in the sequence.
Note also that a(5) is not 242, even though tau evaluated at 242, 243,..., 246 gives 5 nondecreasing values, because here we deal with full runs and 242 belongs to the run of 6 values starting at 241.
		

Crossrefs

Programs

  • Mathematica
    Function[s, {46}~Join~Map[Function[r, Select[s, Last@ # == r &][[1, 1]]], Range[2, Max[s[[All, -1]] ] ]]]@ Map[{#[[1, 1]], Length@ # + 1} &, DeleteCases[SplitBy[#, #[[-1]] >= 0 &], k_ /; k[[1, -1]] < 0]] &@ MapIndexed[{First@ #2, #1} &, Differences@ Array[DivisorSigma[0, #] &, 10^6]] (* Michael De Vlieger, May 06 2017 *)
  • PARI
    genit()={for(n=1,20,q=0;ibgn=1;for(m=ibgn,9E99,mark1=q;q=numdiv(m);if(mark1==0,summ=0;dun=0;mark2=m);if(q>=mark1,summ+=1,dun=1);if(dun>0&&summ==n,print(n," ",mark2);break);if(dun>0&&summ!=n,q=0;m-=1)));} \\ Bill McEachen, Apr 25 2017
    
  • PARI
    A284597=vector(19);apply(scan(N,s=1,t=numdiv(s))=for(k=s+1,N,t>(t=numdiv(k))||next;k-s>#A284597||A284597[k-s]||printf(" a(%d)=%d,",k-s,s)||A284597[k-s]=s;s=k);done,[10^6]) \\ Finds a(1..10) in ~ 1 sec, but would take 100 times longer to get one more term with scan(10^8). You may extend the search using scan(END,START). - M. F. Hasler, May 06 2017
  • Python
    from sympy import divisor_count
    def A284597(n):
        count, starti, s, i = 0,1,0,1
        while True:
            d = divisor_count(i)
            if d < s:
                if count == n:
                    return starti
                starti = i
                count = 0
            s = d
            i += 1
            count += 1 # Chai Wah Wu, May 04 2017
    

Extensions

a(1), a(2), a(4) corrected by Bill McEachen and Giovanni Resta, Apr 26 2017
a(17)-a(18) from Robert Gerbicz, May 14 2017

A339910 After 1, numbers k > 1 such that k has fewer prime divisors than k-1, when they are counted with multiplicity.

Original entry on oeis.org

1, 5, 7, 9, 11, 13, 17, 19, 21, 23, 25, 29, 31, 33, 37, 41, 43, 46, 47, 49, 51, 53, 55, 57, 59, 61, 65, 67, 69, 71, 73, 77, 79, 81, 82, 83, 85, 89, 91, 93, 97, 101, 103, 105, 106, 107, 109, 111, 113, 115, 118, 121, 127, 129, 131, 133, 137, 139, 141, 145, 149, 151, 153, 155, 157, 161, 163, 166, 167, 169, 173, 177, 179
Offset: 1

Views

Author

Antti Karttunen, Dec 22 2020

Keywords

Comments

List of starts of nondecreasing runs of values of A001222 (= bigomega, number of prime factors, when counted with multiplicity).
Differs from A303578 (a similar sequence computed for tau, the number of divisors) for the first time at n=44, where a(44) = 105, a term missing from A303578.

Crossrefs

Cf. A001222.
One more than A322839.
Cf. A339911, A339912 for subsequences.
Cf. also A303578.

Programs

  • Mathematica
    1 + {0}~Join~Position[Sign@ Differences[Array[PrimeOmega, 180]], -1][[All, 1]] (* Michael De Vlieger, Dec 22 2020 *)
  • PARI
    isA339910(n) = ((1==n)||(bigomega(n)
    				

A303577 Break up the list of values of the divisor function d(k) into nondecreasing runs; sequence gives lengths of successive runs.

Original entry on oeis.org

4, 2, 2, 2, 2, 4, 2, 2, 2, 2, 4, 2, 2, 4, 4, 2, 3, 1, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 4, 2, 2, 1, 1, 2, 4, 2, 2, 4, 4, 2, 3, 1, 2, 2, 2, 2, 3, 3, 4, 2, 2, 2, 2, 4, 2, 2, 4, 4, 2, 2, 2, 2, 4, 2, 3, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 4, 2, 2, 4, 4, 2, 2, 4, 4, 2, 2, 1, 1, 2, 4, 2, 2, 2, 2, 6
Offset: 1

Views

Author

N. J. A. Sloane, Apr 29 2018

Keywords

Examples

			The initial values of d(k) = A000005(k) for k = 1,2,3,... are
1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 6, 2, 4, 4, 5, 2, 6, 2, 6, 4, 4, 2, 8, 3, 4, 4, 6, 2, 8, 2, 6, 4, 4, 4, 9, 2, 4, 4, 8, 2, 8, 2, 6, 6, 4, 2, 10, 3, 6, 4, 6, 2, 8, 4, 8, 4, 4, 2, 12, 2, 4, 6, 7, 4, 8, 2, 6, 4, 8, 2, 12, 2, 4, 6, 6, 4, 8, 2, ...
Breaking this up into nondecreasing runs we get:
[1, 2, 2, 3], [2, 4], [2, 4], [3, 4], [2, 6], [2, 4, 4, 5], [2, 6], [2, 6], [4, 4], [2, 8], [3, 4, 4, 6], [2, 8], [2, 6], [4, 4, 4, 9], [2, 4, 4, 8], [2, 8], [2, 6, 6], [4], [2, 10], [3, 6], [4, 6], [2, 8], [4, 8], [4, 4], [2, 12], [2, 4, 6, 7], ...
whose successive lengths are
4,2,2,2,2,4,2,2,2,2,4,2,2,4,4,2,3,1,2,...
		

Crossrefs

Cf. A000005.
A303578(m) gives value of n that starts the m-th run.
A284597(m) is the smallest number that starts a run of length m.

Extensions

More terms from Seiichi Manyama, Apr 29 2018

A340157 Numbers m such that numbers m, m + 1, m + 2 and m + 3 have k, 2k, 3k and 4k divisors respectively.

Original entry on oeis.org

421, 3013, 5029, 5223, 5245, 5893, 6487, 10533, 11911, 14677, 17173, 23077, 23573, 24613, 25141, 25213, 27637, 27973, 28357, 30661, 32407, 34117, 37477, 38282, 39751, 43495, 45973, 47365, 48423, 50821, 50965, 53413, 53989, 54421, 55141, 56103, 57877, 58165
Offset: 1

Views

Author

Jaroslav Krizek, Dec 29 2020

Keywords

Comments

Numbers m such that tau(m) = tau(m + 1)/2 = tau(m + 2)/3 = tau(m + 3)/4, where tau(k) = the number of divisors of k (A000005).
Quadruplets of [tau(a(n)), tau(a(n) + 1), tau(a(n) + 2), tau(a(n) + 3)] = [tau(a(n)), 2*tau(a(n)), 3*tau(a(n)), 4*tau(a(n))]: [2, 4, 6, 8], [4, 8, 12, 16], [4, 8, 12, 16], [4, 8, 12, 16], [4, 8, 12, 16], [4, 8, 12, 16], [4, 8, 12, 16], ...
Corresponding values of tau(a(n)): 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, ...
Subsequence of A063446 and A339778. Supersequence of A340158.
Prime terms (primes p such that p, p + 1, p + 2 and p + 3 have 2, 4, 6 and 8 divisors respectively): 421, 30661, 50821, 54421, 130021, 195541, 423781, 635461, 1003381, 1577941, 1597381, 1883941, ...

Examples

			tau(421) = 2, tau(422) = 4, tau(423) = 6, tau(424) = 8.
		

Crossrefs

Programs

  • Magma
    [m: m in [1..10^5] | #Divisors(m) eq #Divisors(m + 1)/2 and #Divisors(m) eq #Divisors(m + 2)/3 and #Divisors(m) eq #Divisors(m + 3)/4]
    
  • Mathematica
    Select[Range[60000], Equal @@ (DivisorSigma[0, # + {0, 1, 2, 3}]/{1, 2, 3, 4}) &] (* Amiram Eldar, Dec 30 2020 *)
  • PARI
    isok(m, n=4) = {my(k=numdiv(m)); for (i=1, n-1, if (numdiv(m+i) != (i+1)*k, return (0));); return(1);} \\ Michel Marcus, Dec 30 2020
Showing 1-4 of 4 results.