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-3 of 3 results.

A070843 Largest n-digit number with exactly n distinct prime divisors. There are no further terms.

Original entry on oeis.org

9, 99, 996, 9996, 99990, 999570, 9999990, 99981420, 999068070, 9592993410
Offset: 1

Views

Author

Amarnath Murthy, May 12 2002

Keywords

Comments

Primes counted without multiplicity. - Harvey P. Dale, Sep 30 2021

Examples

			a(4) = 9996 with the prime divisors 2,3,7 and 17.
		

Crossrefs

Cf. A070842.

Programs

  • Mathematica
    Do[k = 10^n - 1; While[Length[FactorInteger[k]] != n, k-- ]; Print[k], {n, 1, 10}] (* Ryan Propper, Nov 05 2005 *)
    lnpd[n_]:=Module[{k=10^n-1},While[PrimeNu[k]!=n,k--];k]; Array[lnpd,10] (* Harvey P. Dale, Sep 30 2021 *)

Extensions

More terms from Rick L. Shepherd, May 14 2002
a(10) from Ryan Propper, Nov 05 2005
Edited by Ray Chandler, Feb 15 2009
Definition clarified by Harvey P. Dale, Sep 30 2021

A342108 Smallest positive integer m with n digits and such that omega(m) = bigomega(m) = n.

Original entry on oeis.org

2, 10, 102, 1110, 10010, 101010, 1009470, 11741730, 223092870, 6469693230
Offset: 1

Views

Author

Bernard Schott, Feb 28 2021

Keywords

Comments

Equivalently: smallest n-digit squarefree number with n distinct prime factors.
Differs from A036336 where length(m) = bigomega(m) = n, when length(m) is the number of digits of m (A055642) and the n prime factors of m are counted with multiplicity (A001222).
Differs from A070842 where length(m) = omega(m) = n, when length(m) is the number of digits of m (A055642) and omega(m) is the number of distinct prime factors dividing m (A001221).
The first index for which these three sequences give three distinct terms is 4:
-> a(4) = 1110 = 2 * 3 * 5 * 37 , with length(1110) = omega(1110) = bigomega(1110) = 4.
-> A036336(4) = 1012 = 2 * 2 * 11 * 23 with length(1012) = bigomega(1012) = 4 > omega(1012) = 3.
-> A070842(4) = 1020 = 2 * 2 * 3 * 5 * 17 with length(1020) = omega(1020) = 4 < bigomega(1020) = 5.
As these terms are the smallest n-digit numbers in A167050 that is finite, this sequence is also finite with 10 terms, as for A070842.

Examples

			10010 = 2*5*7*11*13 is the smallest 5-digit number such that omega(10010) = bigomega(10010) = 5, hence a(5) = 10010.
		

Crossrefs

Subsequence of A167050.

Programs

  • Mathematica
    a={};For[n=1,n<=10,n++,For[m=10^(n-1),m<10^n,m++,If[PrimeOmega[m]==PrimeNu[m]==n,AppendTo[a,m];Break[]]]];a (* Stefano Spezia, Mar 04 2021 *)

Formula

A036336(n) <= A070842(n) <= a(n).

A364468 Number of primitive n-bead necklaces (turning over is allowed) comprising elements of two flavors where complements and scalings are equivalent.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 3, 8, 12, 20, 35, 62, 106, 189, 343, 603, 1130, 2055, 3860, 7154, 13562, 25463, 48607, 92204, 176646, 337587, 649151, 1246819, 2404519, 4636389, 8963497, 17334800, 33585928, 65107935, 126385919, 245492221, 477345359, 928772649, 1808662015, 3524337599, 6872457828, 13409202675, 26179870365
Offset: 0

Views

Author

Richard B. Canty, Jul 25 2023

Keywords

Comments

A complement necklace is one where the flavor of each element is inverted ("010" is equivalent to "101"). A scaled necklace is one where each element in the sequence is repeated by the same integer scalar ("010" is equivalent to "001100", "000111000", etc.).

Examples

			For a(4) = 1, there is one solution: "1110". The other primitive sequence "1100" can be reduced to "10", which no longer uses 4 elements.
For a(6) = 3, there are three solutions: "111110", "111010", and "110010". The other primitive sequences "111100" and "111000" can be reduced to "110" and "10", respectively, which no longer use 6 elements.
		

Crossrefs

Programs

  • PARI
    a11(n) = if( n<1, n==0, 2^(n\2) / 2 + sumdiv(n, k, eulerphi(2*k) * 2^(n/k)) / (4*n));
    a46(n) = {
      my(s=0);
      fordiv (n, d,
        s+=moebius(d)*a11(n/d));
      s};
    a364468(n) = {
      my(s=a46(n));
      fordiv (n, k,
        s-=if(k!=1&&k!=n, a364468(k), 0));
      s};
    for (k=1,42, print1(a364468(k),", "))  \\ Hugo Pfoertner, Jul 26 2023

Formula

a(n) = A000046(n) - Sum_{k = nontrivial divisors of n} a(k). (nontrivial divisors, d: 1 < d < n.)
Showing 1-3 of 3 results.