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.

Previous Showing 11-13 of 13 results.

A353124 Numbers k which have a record number of non-divisors < k (i.e., A049820(k)).

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 13, 17, 19, 22, 23, 25, 27, 29, 31, 34, 35, 37, 41, 43, 46, 47, 49, 51, 53, 57, 58, 59, 61, 65, 67, 71, 73, 77, 79, 82, 83, 86, 87, 89, 93, 94, 95, 97, 101, 103, 106, 107, 109, 113, 118, 119, 121, 123, 125, 127, 131, 134, 137, 139, 142, 143
Offset: 1

Views

Author

Johan Lindgren, Apr 24 2022

Keywords

Comments

From Jon E. Schoenfield, Apr 30 2022: (Start)
This sequence includes all noncomposite numbers, the squares of all odd primes, and the cube of every odd prime p such that p^3 - 2 is composite.
It also includes every number k of the form p*q, with p and q distinct primes, such that k-2 is composite and k-1 is neither a prime nor the square of a prime.
In general, it includes every number k such that tau(k-j) > tau(k) - j for each j in 1..tau(k)-1.
Terms with larger numbers of divisors occur less frequently. The first terms with 0, 1, 2, 3, and 4 distinct prime factors are 1, 3, 22, 2110, and 17585778, respectively (each of which is squarefree). What is the first term with 5 distinct prime factors?
(End)

Crossrefs

Programs

  • Mathematica
    s = {}; fm = -1; Do[f = n - DivisorSigma[0, n]; If[f > fm, fm = f; AppendTo[s, n]], {n, 1, 120}]; s (* Amiram Eldar, Apr 25 2022 *)
  • PARI
    f(n) = n - numdiv(n); \\ A049820
    lista(nn) = {my(m=-oo, list=List(), fn); for (n=1, nn, if ((fn=f(n)) > m, listput(list, n); m = fn;);); Vec(list);} \\ Michel Marcus, Apr 25 2022

A283147 Number n such that there are no primes of the form sigma(n)/k where 1 < k < n is a (proper) nondivisor of n.

Original entry on oeis.org

1, 2, 4, 9, 12, 16, 25, 48, 64, 112, 192, 240, 289, 448, 729, 960, 1344, 1681, 1984, 2401, 3481, 4096, 5041, 6720, 7921, 10201, 12288, 15625, 17161, 27889, 28561, 28672, 29929, 39270, 53130, 61440, 65536, 71610, 82110
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Mar 01 2017

Keywords

Crossrefs

Supersequence of A023194, A019279 and A061652.

Programs

  • PARI
    is(n)=my(s=sigma(n), p=factor(s)[,1], k); for(i=1,#p, k=s/p[i]; if(kCharles R Greathouse IV, Mar 01 2017

Extensions

Corrected by Charles R Greathouse IV, Mar 01 2017

A365533 a(n) is the nim-value of the SALIQUANT game where the option is to subtract a nondivisor from 2*n.

Original entry on oeis.org

0, 1, 1, 3, 2, 4, 6, 7, 4, 7, 5, 10, 12, 10, 13, 15, 8, 13, 9, 17, 17, 16, 11, 22, 22, 19, 25, 24, 14, 22, 30, 31, 16, 33, 32, 31, 18, 28, 19, 37, 20, 38, 21, 38, 37, 34, 23, 46, 45, 37, 42, 51, 26, 40, 27, 52, 28, 43, 29, 52, 60, 61, 52, 63, 58, 49, 66, 59, 34, 52, 35, 67, 36, 55, 62
Offset: 1

Views

Author

Michel Marcus, Sep 08 2023

Keywords

Comments

a(n) = SG(2*n) where SG(n) = mex{x in opt(n)} SG(x), where mex(A) is the least nonnegative integer not appearing in A, and opt(n) is a vector of values n-k where 1 <= k <= n is not a divisor of n. For odd n, SG(n) = (n-1)/2.
Note that SG(n) represents the nim-value (also called the Sprague-Grundy number) for position n in combinatorial game theory. It indicates the winning strategy for that position when both players play optimally.

Crossrefs

Cf. A173540 (nondivisors).

Programs

  • Mathematica
    mex[l_List]:=Module[{i=0},While[MemberQ[l,i],i++];i];SG[n_Integer?Positive]:=SG[n]=Module[{p,d},If[n==1,Return[0]];d=Select[Range[1,n],Mod[n,#]!=0&];p=n-d;mex[SG[#]&/@p]];a[n_]:=Module[{r={}},Do[AppendTo[r,SG[2*i]],{i,n}];r]; a[75] (* Robert P. P. McKone, Sep 09 2023 *)
  • PARI
    opt(n) = my(list=List()); for (k=1, n, if (n % k, listput(list, n-k))); Vec(vecsort(list));
    lista(nn) = {nn *= 2; my(vsg = vector(nn, n, if (n%2, (n-1)/2))); forstep (n=2, nn, 2, my(v=row(n), list=List()); for (i=1, #v, listput(list, vsg[v[i]])); list = Vec(vecsort(list)); if (#list==0, vsg[n] = 0, for (k=1, vecmax(list)+1, if (!vecsearch(list, k), vsg[n] = k; break)));); vector(nn\2, k, vsg[2*k]);}
Previous Showing 11-13 of 13 results.