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-10 of 15 results. Next

A013929 Numbers that are not squarefree. Numbers that are divisible by a square greater than 1. The complement of A005117.

Original entry on oeis.org

4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, 52, 54, 56, 60, 63, 64, 68, 72, 75, 76, 80, 81, 84, 88, 90, 92, 96, 98, 99, 100, 104, 108, 112, 116, 117, 120, 121, 124, 125, 126, 128, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160
Offset: 1

Views

Author

Keywords

Comments

Sometimes misnamed squareful numbers, but officially those are given by A001694.
This is different from the sequence of numbers k such that A007913(k) < phi(k). The two sequences differ at the values: 420, 660, 780, 840, 1320, 1560, 4620, 5460, 7140, ..., which is essentially A070237. - Ant King, Dec 16 2005
Numbers k such that Sum_{d|k} (d/phi(d))*mu(k/d) = 0. - Benoit Cloitre, Apr 28 2002
Also, k with at least one x < k such that A007913(x) = A007913(k). - Benoit Cloitre, Apr 28 2002
Numbers k for which there exists a partition into two parts p and q such that p + q = k and p*q is a multiple of k. - Amarnath Murthy, May 30 2003
Numbers k such that there is a solution 0 < x < k to x^2 == 0 (mod k). - Franz Vrabec, Aug 13 2005
Numbers k such that moebius(k) = 0.
a(n) = k such that phi(k)/k = phi(m)/m for some m < k. - Artur Jasinski, Nov 05 2008
Appears to be numbers such that when a column with index equal to a(n) in A051731 is deleted, there is no impact on the result in the first column of A054525. - Mats Granvik, Feb 06 2009
Numbers k such that the number of prime divisors of (k+1) is less than the number of nonprime divisors of (k+1). - Juri-Stepan Gerasimov, Nov 10 2009
Orders for which at least one non-cyclic finite abelian group exists: A000688(a(n)) > 1. This follows from the fact that not all exponents in the prime factorization of a(n) are 1 (moebius(a(n)) = 0). The number of such groups of order a(n) is A192005(n) = A000688(a(n)) - 1. - Wolfdieter Lang, Jul 29 2011
Subsequence of A193166; A192280(a(n)) = 0. - Reinhard Zumkeller, Aug 26 2011
It appears that terms are the numbers m such that Product_{k=1..m} (prime(k) mod m) <> 0. See Maple code. - Gary Detlefs, Dec 07 2011
A008477(a(n)) > 1. - Reinhard Zumkeller, Feb 17 2012
A057918(a(n)) > 0. - Reinhard Zumkeller, Mar 27 2012
A056170(a(n)) > 0. - Reinhard Zumkeller, Dec 29 2012
Numbers k such that A001221(k) != A001222(k). - Felix Fröhlich, Aug 13 2014
Numbers k such that A001222(k) > A001221(k), since in this case at least one prime factor of k occurs more than once, which implies that k is divisible by at least one perfect square > 1. - Carlos Eduardo Olivieri, Aug 02 2015
Lexicographically least sequence such that each term has a positive even number of proper divisors not occurring in the sequence, cf. the sieve characterization of A005117. - Glen Whitney, Aug 30 2015
There are arbitrarily long runs of consecutive terms. Record runs start at 4, 8, 48, 242, ... (A045882). - Ivan Neretin, Nov 07 2015
A number k is a term if 0 < min(A000010(k) + A023900(k), A000010(k) - A023900(k)). - Torlach Rush, Feb 22 2018
Every squareful number > 1 is nonsquarefree, but the converse is false and the nonsquarefree numbers that are not squareful (see first comment) are in A332785. - Bernard Schott, Apr 11 2021
Integers m where at least one k < m exists such that m divides k^m. - Richard R. Forberg, Jul 31 2021
Consider the Diophantine equation S(x,y) = (x+y) + (x-y) + (x*y) + (x/y) = z, when x and y are both positive integers with y | x. Then, there is a solution (x,y) iff z is a term of this sequence; in this case, if x = K*y, then z = S(K*y,y) = K*(y+1)^2 (see A351381, link and references Perelman); example: S(12,4) = 75 = a(28). The number of solutions for S(x,y) = a(n) is A353282(n). - Bernard Schott, Mar 29 2022
For each positive integer m, the number of unitary divisors of m = the number of squarefree divisors of m (see A034444); but only for the terms of this sequence does the set of unitary divisors differ from the set of squarefree divisors. Example: the set of unitary divisors of 20 is {1, 4, 5, 20}, while the set of squarefree divisors of 20 is {1, 2, 5, 10}. - Bernard Schott, Oct 15 2022

Examples

			For the terms up to 20, we compute the squares of primes up to floor(sqrt(20)) = 4. Those squares are 4 and 9. For every such square s, put the terms s*k^2 for k = 1 to floor(20 / s). This gives after sorting and removing duplicates the list 4, 8, 9, 12, 16, 18, 20. - _David A. Corneth_, Oct 25 2017
		

References

  • I. Perelman, L'Algèbre récréative, Deux nombres et quatre opérations, Editions en langues étrangères, Moscou, 1959, pp. 101-102.
  • Ya. I. Perelman, Algebra can be fun, Two numbers and four operations, Mir Publishers Moscow, 1979, pp. 131-132.

Crossrefs

Complement of A005117. Subsequences: A130897, A190641, A332785.
Partitions into: A114374, A256012.

Programs

  • Haskell
    a013929 n = a013929_list !! (n-1)
    a013929_list = filter ((== 0) . a008966) [1..]
    -- Reinhard Zumkeller, Apr 22 2012
    
  • Magma
    [ n : n in [1..1000] | not IsSquarefree(n) ];
    
  • Maple
    a := n -> `if`(numtheory[mobius](n)=0,n,NULL); seq(a(i),i=1..160); # Peter Luschny, May 04 2009
    t:= n-> product(ithprime(k),k=1..n): for n from 1 to 160 do (if t(n) mod n <>0) then print(n) fi od; # Gary Detlefs, Dec 07 2011
    with(NumberTheory): isQuadrateful := n -> irem(Radical(n), n) <> 0:
    select(isQuadrateful, [`$`(1..160)]);  # Peter Luschny, Jul 12 2022
  • Mathematica
    Union[ Flatten[ Table[ n i^2, {i, 2, 20}, {n, 1, 400/i^2} ] ] ]
    Select[ Range[2, 160], (Union[Last /@ FactorInteger[ # ]][[ -1]] > 1) == True &] (* Robert G. Wilson v, Oct 11 2005 *)
    Cases[Range[160], n_ /; !SquareFreeQ[n]] (* Jean-François Alcover, Mar 21 2011 *)
    Select[Range@160, ! SquareFreeQ[#] &] (* Robert G. Wilson v, Jul 21 2012 *)
    Select[Range@160, PrimeOmega[#] > PrimeNu[#] &] (* Carlos Eduardo Olivieri, Aug 02 2015 *)
    Select[Range[200], MoebiusMu[#] == 0 &] (* Alonso del Arte, Nov 07 2015 *)
  • PARI
    {a(n)= local(m,c); if(n<=1,4*(n==1), c=1; m=4; while( cMichael Somos, Apr 29 2005 */
    
  • PARI
    for(n=1, 1e3, if(omega(n)!=bigomega(n), print1(n, ", "))) \\ Felix Fröhlich, Aug 13 2014
    
  • PARI
    upto(n)=my(res = List()); forprime(p = 2, sqrtint(n), for(k = 1, n \ p^2, listput(res, k * p^2))); listsort(res, 1); res \\ David A. Corneth, Oct 25 2017
    
  • Python
    from sympy.ntheory.factor_ import core
    def ok(n): return core(n, 2) != n
    print(list(filter(ok, range(1, 161)))) # Michael S. Branicky, Apr 08 2021
    
  • Python
    from math import isqrt
    from sympy import mobius
    def A013929(n):
        def f(x): return n+sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
        m, k = n, f(n)
        while m != k:
            m, k = k, f(k)
        return m # Chai Wah Wu, Jul 20 2024

Formula

A008966(a(n)) = 0. - Reinhard Zumkeller, Apr 22 2012
Sum_{n>=1} 1/a(n)^s = (zeta(s)*(zeta(2*s)-1))/zeta(2*s). - Enrique Pérez Herrero, Jul 07 2012
a(n) ~ n/k, where k = 1 - 1/zeta(2) = 1 - 6/Pi^2 = A229099. - Charles R Greathouse IV, Sep 13 2013
A001222(a(n)) > A001221(a(n)). - Carlos Eduardo Olivieri, Aug 02 2015
phi(a(n)) > A003958(a(n)). - Juri-Stepan Gerasimov, Apr 09 2019

Extensions

More terms from Erich Friedman
More terms from Franz Vrabec, Aug 13 2005

A087188 Number of partitions of n into distinct squarefree parts.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 3, 3, 4, 4, 5, 6, 6, 8, 9, 10, 13, 14, 16, 18, 20, 24, 27, 30, 35, 37, 42, 47, 51, 59, 64, 72, 81, 88, 98, 109, 120, 134, 147, 163, 179, 195, 216, 236, 258, 284, 310, 339, 371, 403, 441, 480, 523, 572, 621, 675, 734, 796, 865, 937, 1014, 1100, 1189
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 24 2003

Keywords

Examples

			n=9: 5+3+1 = 6+2+1 = 6+3 = 7+2: a(9)=4;
n=10: 5+3+2 = 6+3+1 = 7+2+1 = 7+3 = 10: a(10)=5.
		

Crossrefs

Programs

  • Haskell
    a087188 = p a005117_list where
       p _      0 = 1
       p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
    -- Reinhard Zumkeller, Jun 01 2015
    
  • Maple
    with(numtheory):
    b:= proc(n, i) option remember;
          `if`(i*(i+1)/2 b(n$2):
    seq(a(n), n=0..100);  # Alois P. Heinz, Jun 02 2015
  • Mathematica
    b[n_, i_] := b[n, i] = If[i*(i+1)/2 < n, 0, If[n == 0, 1, b[n, i-1] + If[i <= n && SquareFreeQ[i], b[n-i, i-1], 0]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Jun 24 2015, after Alois P. Heinz *)
    nmax = 100; CoefficientList[Series[Exp[Sum[(-1)^(j + 1)/j * Sum[Abs[MoebiusMu[k]] * x^(j*k), {k, 1, Floor[nmax/j] + 1}], {j, 1, nmax}]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Mar 31 2018 *)
  • PARI
    ok(v)=for(i=2,#v, if(v[i]==v[i-1] || !issquarefree(v[i]), return(0))); #v==0 || issquarefree(v[1])
    a(n)=my(s,u); forpart(v=n, if(ok(v), s++)); s \\ Charles R Greathouse IV, Nov 05 2017

Formula

O.g.f.: product_{i=1,2,...infinity} [1+x^A005117(i)]. - R. J. Mathar, May 16 2008
a(n) ~ exp(sqrt(2*n)) / (2^(1/4) * sqrt(Pi) * n^(3/4)). - Vaclav Kotesovec, Mar 24 2018

Extensions

Offset changed and a(0)=1 prepended by Reinhard Zumkeller, Jun 01 2015

A114374 Number of partitions of n into parts that are not squarefree.

Original entry on oeis.org

1, 0, 0, 0, 1, 0, 0, 0, 2, 1, 0, 0, 3, 1, 0, 0, 5, 2, 2, 0, 7, 3, 2, 0, 11, 6, 4, 3, 15, 8, 6, 3, 22, 13, 11, 6, 34, 18, 15, 9, 46, 27, 24, 17, 64, 43, 33, 23, 89, 60, 51, 37, 124, 84, 78, 51, 166, 119, 109, 78, 226, 168, 152, 118, 300, 228, 215, 166, 404, 313, 300, 230, 546, 421, 409
Offset: 0

Views

Author

Reinhard Zumkeller, Feb 09 2006

Keywords

Comments

a(A078135(n)) = 0; a(A078137(n)) > 0.

Examples

			a(12) = #{2*2*3, 2*2*2 + 2*2, 2*2 + 2*2 + 2*2} = 3;
a(13) = #{3*3 + 2*2} = 1.
		

Crossrefs

Programs

  • Haskell
    a114374 = p a013929_list where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Jun 01 2015
  • Maple
    with(numtheory):
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+
          `if`(i>n or issqrfree(i), 0, b(n-i, i))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..100);  # Alois P. Heinz, Jun 03 2015
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n || SquareFreeQ[i], 0, b[n-i, i]]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Jun 30 2015, after Alois P. Heinz *)

Formula

a(n) = A000041(n) - A073576(n) - A117395(n). - Reinhard Zumkeller, Mar 11 2006
G.f.: Product_{k>=1} (1 - mu(k)^2*x^k)/(1 - x^k), where mu(k) is the Moebius function (A008683). - Ilya Gutkovskiy, Dec 30 2016

Extensions

Offset changed and a(0)=1 prepended by Reinhard Zumkeller, Jun 01 2015

A379301 Positive integers whose prime indices include a unique composite number.

Original entry on oeis.org

7, 13, 14, 19, 21, 23, 26, 28, 29, 35, 37, 38, 39, 42, 43, 46, 47, 52, 53, 56, 57, 58, 61, 63, 65, 69, 70, 71, 73, 74, 76, 77, 78, 79, 84, 86, 87, 89, 92, 94, 95, 97, 101, 103, 104, 105, 106, 107, 111, 112, 113, 114, 115, 116, 117, 119, 122, 126, 129, 130, 131
Offset: 1

Views

Author

Gus Wiseman, Dec 25 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 70 are {1,3,4}, so 70 is in the sequence.
The prime indices of 98 are {1,4,4}, so 98 is not in the sequence.
		

Crossrefs

For no composite parts we have A302540, counted by A034891 (strict A036497).
For all composite parts we have A320629, counted by A023895 (strict A204389).
For a unique prime part we have A331915, counted by A379304 (strict A379305).
Positions of one in A379300.
Partitions of this type are counted by A379302 (strict A379303).
A000040 lists the prime numbers, differences A001223.
A002808 lists the composite numbers, nonprimes A018252, differences A073783 or A065310.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A066247 is the characteristic function for the composite numbers.
A377033 gives k-th differences of composite numbers, see A073445, A377034-A377037.
Other counts of prime indices:
- A087436 postpositive, see A038550.
- A330944 nonprime, see A002095, A096258, A320628, A330945.
- A379306 squarefree, see A302478, A379308, A379309, A379316.
- A379310 nonsquarefree, see A114374, A256012, A379307.
- A379311 old prime, see A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Length[Select[prix[#],CompositeQ]]==1&]

A379308 Number of integer partitions of n with a unique squarefree part.

Original entry on oeis.org

0, 1, 1, 1, 0, 2, 2, 2, 0, 3, 5, 5, 1, 6, 9, 9, 2, 10, 14, 18, 6, 18, 24, 30, 11, 28, 39, 47, 24, 48, 63, 76, 41, 74, 95, 118, 65, 120, 149, 181, 107, 181, 221, 266, 169, 266, 335, 398, 262, 394, 487, 578, 391, 578, 697, 844, 592, 834, 997, 1198, 867
Offset: 0

Views

Author

Gus Wiseman, Dec 26 2024

Keywords

Examples

			The a(1) = 1 through a(11) = 5 partitions:
  (1)  (2)  (3)  .  (5)    (6)    (7)    .  (5,4)    (10)     (11)
                    (4,1)  (4,2)  (4,3)     (8,1)    (6,4)    (7,4)
                                            (4,4,1)  (8,2)    (8,3)
                                                     (9,1)    (9,2)
                                                     (4,4,2)  (4,4,3)
		

Crossrefs

If all parts are squarefree we have A073576 (strict A087188), ranks A302478.
If no parts are squarefree we have A114374 (strict A256012), ranks A379307.
For composite instead of squarefree we have A379302 (strict A379303), ranks A379301.
For prime instead of squarefree we have A379304, (strict A379305), ranks A331915.
The strict case is A379309.
For old prime instead of squarefree we have A379314, (strict A379315), ranks A379312.
Ranked by A379316, positions of 1 in A379306.
A000041 counts integer partitions, strict A000009.
A005117 lists the squarefree numbers, differences A076259.
A013929 lists the nonsquarefree numbers, differences A078147.
A377038 gives k-th differences of squarefree numbers.
A379310 counts nonsquarefree prime indices.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Count[#,_?SquareFreeQ]==1&]],{n,0,30}]

A379309 Number of strict integer partitions of n with a unique squarefree part.

Original entry on oeis.org

0, 1, 1, 1, 0, 2, 2, 2, 0, 2, 4, 4, 1, 4, 7, 7, 2, 6, 8, 11, 4, 9, 13, 17, 7, 13, 20, 22, 13, 20, 29, 33, 21, 29, 40, 47, 27, 41, 56, 64, 42, 59, 77, 85, 60, 74, 104, 115, 83, 101, 141, 155, 113, 138, 179, 206, 156, 183, 236, 272, 212, 239, 309, 343, 282, 315
Offset: 0

Views

Author

Gus Wiseman, Dec 27 2024

Keywords

Examples

			The a(9) = 2 through a(15) = 7 partitions:
  (5,4)  (10)   (11)   (9,3)  (13)     (14)     (15)
  (8,1)  (6,4)  (7,4)         (8,5)    (8,6)    (8,7)
         (8,2)  (8,3)         (12,1)   (9,5)    (9,6)
         (9,1)  (9,2)         (8,4,1)  (10,4)   (11,4)
                                       (12,2)   (12,3)
                                       (8,4,2)  (8,4,3)
                                       (9,4,1)  (9,4,2)
		

Crossrefs

If all parts are squarefree we have A087188, non-strict A073576 (ranks A302478).
If no parts are squarefree we have A256012, non-strict A114374 (ranks A379307).
For composite instead of squarefree we have A379303, non-strict A379302 (ranks A379301).
For prime instead of squarefree we have A379305, non-strict A379304 (ranks A331915).
The non-strict version is A379308, ranks A379316.
For old prime instead of squarefree we have A379315, non-strict A379314 (ranks A379312).
Ranked by A379316 /\ A005117 = squarefree positions of 1 in A379306.
A000041 counts integer partitions, strict A000009.
A005117 lists the squarefree numbers, differences A076259.
A013929 lists the nonsquarefree numbers, differences A078147.
A377038 gives k-th differences of squarefree numbers.
A379310 counts nonsquarefree prime indices.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Count[#,_?SquareFreeQ]==1&]],{n,0,30}]
  • PARI
    lista(nn) = my(r=1, s=0); for(k=1, nn, if(issquarefree(k), s+=x^k, r*=1+x^k)); concat(0, Vec(r*s+O(x^(1+nn)))); \\ Jinyuan Wang, Feb 21 2025

Extensions

More terms from Jinyuan Wang, Feb 21 2025

A379312 Positive integers whose prime indices include a unique 1 or prime number.

Original entry on oeis.org

2, 3, 5, 11, 14, 17, 21, 26, 31, 35, 38, 39, 41, 46, 57, 58, 59, 65, 67, 69, 74, 77, 83, 86, 87, 94, 95, 98, 106, 109, 111, 115, 119, 122, 127, 129, 141, 142, 143, 145, 146, 147, 157, 158, 159, 178, 179, 182, 183, 185, 191, 194, 202, 206, 209, 211, 213, 214
Offset: 1

Views

Author

Gus Wiseman, Dec 28 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The terms together with their prime indices begin:
    2: {1}
    3: {2}
    5: {3}
   11: {5}
   14: {1,4}
   17: {7}
   21: {2,4}
   26: {1,6}
   31: {11}
   35: {3,4}
   38: {1,8}
   39: {2,6}
   41: {13}
   46: {1,9}
   57: {2,8}
   58: {1,10}
   59: {17}
   65: {3,6}
   67: {19}
   69: {2,9}
   74: {1,12}
   77: {4,5}
		

Crossrefs

These "old" primes are listed by A008578.
For no composite parts we have A302540, counted by A034891 (strict A036497).
For all composite parts we have A320629, counted by A023895 (strict A204389).
For a unique prime part we have A331915, counted by A379304 (strict A379305).
Positions of ones in A379311, see A379313.
Partitions of this type are counted by A379314, strict A379315.
A000040 lists the prime numbers, differences A001223.
A002808 lists the composite numbers, nonprimes A018252, differences A073783 or A065310.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A080339 is the characteristic function for the old prime numbers.
A376682 gives k-th differences of old prime numbers, see A030016, A075526.
Other counts of prime indices:
- A330944 nonprime, see A002095, A096258, A320628, A330945.
- A379306 squarefree, see A302478, A379308, A379309, A379316.
- A379310 nonsquarefree, see A114374, A256012, A379307.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],Length[Select[prix[#],#==1||PrimeQ[#]&]]==1&]

A379316 Positive integers whose prime indices include a unique squarefree number.

Original entry on oeis.org

2, 3, 5, 11, 13, 14, 17, 21, 29, 31, 35, 38, 41, 43, 46, 47, 57, 59, 67, 69, 73, 74, 77, 79, 83, 91, 95, 98, 101, 106, 109, 111, 113, 115, 119, 122, 127, 137, 139, 142, 147, 149, 157, 159, 163, 167, 178, 179, 181, 183, 185, 191, 194, 199, 203, 206, 209, 211
Offset: 1

Views

Author

Gus Wiseman, Dec 29 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The terms together with their prime indices begin:
    2: {1}
    3: {2}
    5: {3}
   11: {5}
   13: {6}
   14: {1,4}
   17: {7}
   21: {2,4}
   29: {10}
   31: {11}
   35: {3,4}
   38: {1,8}
   41: {13}
   43: {14}
   46: {1,9}
		

Crossrefs

For all squarefree parts we have A302478, zeros of A379310.
Positions of 1 in A379306.
For no squarefree parts we have A379307, counted by A114374, strict A256012.
Partitions of this type are counted by A379308, strict A379309.
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
Other counts of prime indices:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379311 prime or 1, see A204389, A320629, A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],Length[Select[prix[#],SquareFreeQ]]==1&]

A379300 Number of prime indices of n that are composite.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 2, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 25 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 39 are {2,6}, so a(39) = 1.
The prime indices of 70 are {1,3,4}, so a(70) = 1.
The prime indices of 98 are {1,4,4}, so a(98) = 2.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 3.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 3.
		

Crossrefs

Positions of first appearances are A000420.
Positions of zero are A302540, counted by A034891 (strict A036497).
Positions of one are A379301, counted by A379302 (strict A379303).
A000040 lists the prime numbers, differences A001223.
A002808 lists the composite numbers, nonprimes A018252, differences A073783 or A065310.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A066247 is the characteristic function for the composite numbers.
A377033 gives k-th differences of composite numbers, see A073445, A377034-A377037.
Other counts of prime indices:
- A087436 postpositive, see A038550.
- A330944 nonprime, see A002095, A096258, A320628, A330945.
- A379306 squarefree, see A302478, A379308, A379309, A379316.
- A379310 nonsquarefree, see A114374, A256012, A379307.
- A379311 old prime, see A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[prix[n],CompositeQ]],{n,100}]

Formula

Totally additive with a(prime(k)) = A066247(k).

A379311 Number of prime indices of n that are 1 or prime.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 27 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 39 are {2,6}, so a(39) = 1.
The prime indices of 70 are {1,3,4}, so a(70) = 2.
The prime indices of 98 are {1,4,4}, so a(98) = 1.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 1.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 2.
		

Crossrefs

Positions of first appearances are A000079.
These "old" primes are listed by A008578.
Positions of zero are A320629, counted by A023895 (strict A204389).
Positions of one are A379312, counted by A379314 (strict A379315).
Positions of nonzero terms are A379313.
A000040 lists the prime numbers, differences A001223.
A002808 lists the composite numbers, nonprimes A018252, differences A073783 or A065310.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A080339 is the characteristic function for the old prime numbers.
A376682 gives k-th differences of old prime numbers, see A030016, A075526, A173390, A376683, A376855.
Other counts of prime indices:
- A330944 nonprime, see A002095, A096258, A320628, A330945.
- A379306 squarefree, see A302478, A379308, A379309, A379316.
- A379310 nonsquarefree, see A114374, A256012, A379307.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[prix[n],#==1||PrimeQ[#]&]],{n,100}]

Formula

Totally additive with a(prime(k)) = A080339(k).
Showing 1-10 of 15 results. Next