A179213 -id:A179213 - cp's OEIS frontend

A008966 a(n) = 1 if n is squarefree, otherwise 0.

1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0

Offset: 1

N. J. A. Sloane

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3, 1).

The infinite lower triangular matrix with A008966 on the main diagonal and the rest zeros is the square of triangle A143255. - Gary W. Adamson, Aug 02 2008

Daniel Forgues, Table of n, a(n) for n = 1..100000
Eric Weisstein's World of Mathematics, Moebius Function
Index entries for characteristic functions
Index entries for sequences computed from exponents in factorization of n

Cf. A005117, A008836 (Dirichlet inverse), A013928 (partial sums).
Cf. A179211, A179212, A179213, A179214, A179215.
Cf. A020639, A073576, A087188.
Cf. A124010, A212793, A085357, A156552.
Parity of A002033.
Cf. A047999, A036987
Cf. A082020 (Dgf at s=2), A157289 (Dgf at s=3), A157290 (Dgf at s=4).

a008966 = abs . a008683
-- Reinhard Zumkeller, Dec 13 2015, Dec 15 2014, May 27 2012, Jan 25 2012

Magma
```
[ Abs(MoebiusMu(n)) : n in [1..100]];
```

A008966 := proc(n) if numtheory[issqrfree](n) then 1 ; else 0 ; end if; end proc: # R. J. Mathar, Mar 14 2011

A008966[n_] := Abs[MoebiusMu[n]]; Table[A008966[n], {n, 100}] (* Enrique Pérez Herrero, Apr 15 2010 *)
Table[If[SquareFreeQ[n],1,0],{n,100}] (* or *) Boole[SquareFreeQ/@ Range[ 100]] (* Harvey P. Dale, Feb 28 2015 *)

MuPAD
```
func(abs(numlib::moebius(n)), n):
```
PARI
```
a(n)=if(n<1,0,direuler(p=2,n,1+X))[n]
```

a(n)=issquarefree(n) \\ Michel Marcus, Feb 22 2015

from sympy import factorint
def A008966(n): return int(max(factorint(n).values(),default=1)==1) # Chai Wah Wu, Apr 05 2023

Dirichlet g.f.: zeta(s)/zeta(2s).
a(n) = abs(mu(n)), where mu is the Moebius function (A008683).
a(n) = 0^(bigomega(n) - omega(n)), where bigomega(n) and omega(n) are the numbers of prime factors of n with and without repetition (A001222, A001221, A046660). - Reinhard Zumkeller, Apr 05 2003
Multiplicative with p^e -> 0^(e - 1), p prime and e > 0. - Reinhard Zumkeller, Jul 15 2003
a(n) = 0^(A046951(n) - 1). - Reinhard Zumkeller, May 20 2007
a(n) = 1 - A107078(n). - Reinhard Zumkeller, Oct 03 2008
a(n) = floor(rad(n)/n), where rad() is A007947. - Enrique Pérez Herrero, Nov 13 2009
A175046(n) = a(n)*A073311(n). - Reinhard Zumkeller, Apr 05 2010
a(n) = floor(A000005(n^2)/A007425(n)). - Enrique Pérez Herrero, Apr 15 2010
a(A005117(n)) = 1; a(A013929(n)) = 0; a(n) = A013928(n + 1) - A013928(n). - Reinhard Zumkeller, Jul 05 2010
a(n) * A112526(n) = A063524(n). - Reinhard Zumkeller, Sep 16 2011
a(n) = mu(n) * lambda(n) = A008836(n) * A008683(n). - Enrique Pérez Herrero, Nov 29 2013
a(n) = Sum_{d|n} 2^omega(d)*mu(n/d). - Geoffrey Critzer, Feb 22 2015
a(n) = A085357(A156552(n)). - Antti Karttunen, Mar 06 2017
Limit_{n->oo} (1/n)*Sum_{j=1..n} a(j) = 6/Pi^2. - Andres Cicuttin, Aug 13 2017
a(1) = 1; a(n) = -Sum_{d|n, d < n} (-1)^bigomega(n/d) * a(d). - Ilya Gutkovskiy, Mar 10 2021

Deleted an unclear comment. - N. J. A. Sloane, May 30 2021

A045943 Triangular matchstick numbers: a(n) = 3n(n+1)/2.

Original entry on oeis.org

0, 3, 9, 18, 30, 45, 63, 84, 108, 135, 165, 198, 234, 273, 315, 360, 408, 459, 513, 570, 630, 693, 759, 828, 900, 975, 1053, 1134, 1218, 1305, 1395, 1488, 1584, 1683, 1785, 1890, 1998, 2109, 2223, 2340, 2460, 2583, 2709, 2838, 2970, 3105, 3243, 3384, 3528

Offset: 0

R. K. Guy

Also, 3 times triangular numbers, a(n) = 3*A000217(n).
In the 24-bit RGB color cube, the number of color-lattice-points in r+g+b = n planes at n < 256 equals the triangular numbers. For n = 256, ..., 765 the number of legitimate color partitions is less than A000217(n) because {r,g,b} components cannot exceed 255. For n = 256, ..., 511, the number of non-color partitions are computable with A045943(n-255), while for n = 512, ..., 765, the number of color points in r+g+b planes equals A000217(765-n). - Labos Elemer, Jun 20 2005
If a 3-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is also the smallest number that may be written both as the sum of n-1 consecutive positive integers and n consecutive positive integers. - Claudio Meller, Oct 08 2010
For n >= 3, a(n) equals 4^(2+n)*Pi^(1 - n) times the coefficient of zeta(3) in the following integral with upper bound Pi/4 and lower bound 0: int x^(n+1) tan x dx. - John M. Campbell, Jul 17 2011
The difference a(n)-a(n-1) = 3*n, for n >= 1. - Stephen Balaban, Jul 25 2011 [Comment clarified by N. J. A. Sloane, Aug 01 2024]
Sequence found by reading the line from 0, in the direction 0, 3, ..., and the same line from 0, in the direction 0, 9, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. This is one of the orthogonal axes of the spiral; the other is A032528. - Omar E. Pol, Sep 08 2011
A005449(a(n)) = A000332(3n + 3) = C(3n + 3, 4), a second pentagonal number of triangular matchstick number index number. Additionally, a(n) - 2n is a pentagonal number (A000326). - Raphie Frank, Dec 31 2012
Sum of the numbers from n to 2n. - Wesley Ivan Hurt, Nov 24 2015
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 5376 or 17920 or 20160. - Philippe A.J.G. Chevalier, Dec 28 2015
Also the number of 4-cycles in the (n+4)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Number of terms less than 10^k, k=0,1,2,3,...: 1, 3, 8, 26, 82, 258, 816, 2582, 8165, 25820, 81650, 258199, 816497, 2581989, 8164966, ... - Muniru A Asiru, Jan 24 2018
Numbers of the form 3*m*(2*m + 1) for m = 0, -1, 1, -2, 2, -3, 3, ... - Bruno Berselli, Feb 26 2018
Partial sums of A008585. - Omar E. Pol, Jun 20 2018
Column 1 of A273464. (Number of ways to select a unit lozenge inside an isosceles triangle of side length n; all vertices on a hexagonal lattice.) - R. J. Mathar, Jul 10 2019
Total number of pips in the n-th suit of a double-n domino set. - Ivan N. Ianakiev, Aug 23 2020

			From _Stephen Balaban_, Jul 25 2011: (Start)
T(n), the triangular numbers = number of nodes,
a(n-1) = number of edges in the T(n) graph:
       o    (T(1) = 1, a(0) = 0)
       o
      / \   (T(2) = 3, a(1) = 3)
     o - o
       o
      / \
     o - o  (T(3) = 6, a(2) = 9)
    / \ / \
   o - o - o
... [Corrected by _N. J. A. Sloane_, Aug 01 2024] (End)

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 543.

Vincenzo Librandi, Table of n, a(n) for n = 0..2000
Abderrahim Arabi, Hacène Belbachir, and Jean-Philippe Dubernard, Enumeration of parallelogram polycubes, arXiv:2105.00971 [cs.DM], 2021.
Francesco Brenti and Paolo Sentinelli, Wachs permutations, Bruhat order and weak order, arXiv:2212.04932 [math.CO], 2022.
Sela Fried, Counting r X s rectangles in nondecreasing and Smirnov words, arXiv:2406.18923 [math.CO], 2024. See p. 5.
Jose Manuel Garcia Calcines, Luis Javier Hernandez Paricio, and Maria Teresa Rivas Rodriguez, Semi-simplicial combinatorics of cyclinders and subdivisions, arXiv:2307.13749 [math.CO], 2023. See p. 29.
T. Aaron Gulliver, Sums of Powers of Integers Divisible by Three, Int. J. Contemp. Math. Sciences, Vol. 7, 2012, no. 38, 1895 - 1901.
Alfred Hoehn, Illustration of initial terms of A000326, A005449, A045943, A115067
Milan Janjic, Two Enumerative Functions
Milan Janjic and Boris Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013.
Milan Janjic and Boris Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014) # 14.3.5.
E. Lábos, On the number of RGB-colors we can distinguish. Partition Spectra. Lecture at 7th Hungarian Conference on Biometry and Biomathematics. Budapest. Jul 06, 2005. Applied Ecology and Environmental Research 4(2): 159-169, 2006.
R. J. Mathar, Lozenge tilings of the equilateral triangle, arXiv:1909.06336 [math.CO], 2019.
Eric Weisstein's World of Mathematics, Graph Cycle.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A005448, A002378, A046092, A051162, A126804, A001318, A032528.
3 times n-gonal numbers: A033428, A062741, A094159, A152773, A152751, A152759, A152767, A153783, A153448, A153875.
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
A diagonal of A010027.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A115067, A008585, A005843, A001477, A000217.
Cf. A027480 (partial sums).
Cf. A002378 (3-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
This sequence: Sum_{k = n..2*n} k.
Cf. A304993:   Sum_{k = n..2*n} k*(k+1)/2.
Cf. A050409:   Sum_{k = n..2*n} k^2.
Similar sequences are listed in A316466.

List([0..10^4], n -> 3*n*(n+1)/2); # Muniru A Asiru, Jan 24 2018

a n = sum [x | x <- [n..2*n]] -- Peter Kagey, Jul 27 2015

[3*n*(n+1)/2: n in [0..50]]; // Vincenzo Librandi, May 02 2011

seq(3*binomial(n+1,2), n=0..49); # Zerinvary Lajos, Nov 24 2006

Table[3 n (n + 1)/2, {n, 0, 50}] (* Vladimir Joseph Stephan Orlovsky, Oct 31 2008 *)
3 Accumulate@Range[0, 48] (* Arkadiusz Wesolowski, Oct 29 2012 *)
CoefficientList[Series[-3 x/(x - 1)^3, {x, 0, 47}], x] (* Robert G. Wilson v, Jan 29 2015 *)
LinearRecurrence[{3, -3, 1}, {0, 3, 9}, 50] (* Jean-François Alcover, Dec 12 2016 *)

a(n)=3*binomial(n+1,2) \\ Charles R Greathouse IV, Jun 16 2011

(3 to 150 by 3).scanLeft(0)( + ) // Alonso del Arte, Sep 12 2019

a(n) is the sum of n+1 integers starting from n, i.e., 1+2, 2+3+4, 3+4+5+6, 4+5+6+7+8, etc. - Jon Perry, Jan 15 2004
a(n) = A126890(n+1,n-1) for n>1. - Reinhard Zumkeller, Dec 30 2006
a(n) + A145919(3*n+3) = 0. - Matthew Vandermast, Oct 28 2008
a(n) = A000217(2*n) - A000217(n-1); A179213(n) <= a(n). - Reinhard Zumkeller, Jul 05 2010
a(n) = a(n-1)+3*n, n>0. - Vincenzo Librandi, Nov 18 2010
G.f.: 3*x/(1-x)^3. - Bruno Berselli, Jan 21 2011
a(n) = A005448(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A001477(n)+A000290(n)+A000217(n). - J. M. Bergot, Dec 08 2012
a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2. - Wesley Ivan Hurt, Nov 24 2015
a(n) = A027480(n)-A027480(n-1). - Peter M. Chema, Jan 18 2017.
2*a(n)+1 = A003215(n). - Miquel Cerda, Jan 22 2018
a(n) = T(2*n) - T(n-1), where T(n) = A000217(n). In general, T(k)*T(n) = Sum_{i=0..k-1} (-1)^i*T((k-i)*(n-i)). - Charlie Marion, Dec 06 2020
E.g.f.: 3*exp(x)*x*(2 + x)/2. - Stefano Spezia, May 19 2021
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(2*log(2)-1)/3. (End)
Product_{n>=1} (1 - 1/a(n)) = -(3/(2*Pi))*cos(sqrt(11/3)*Pi/2). - Amiram Eldar, Feb 21 2023

A073837 Sum of primes p satisfying n <= p <= 2n.

Original entry on oeis.org

2, 5, 8, 12, 12, 18, 31, 24, 41, 60, 60, 72, 72, 59, 88, 119, 119, 102, 139, 120, 161, 204, 204, 228, 228, 228, 281, 281, 281, 311, 372, 341, 341, 408, 408, 479, 552, 515, 515, 594, 594, 636, 636, 593, 682, 682, 682, 635, 732, 732, 833, 936, 936, 990, 1099, 1099

Offset: 1

Amarnath Murthy, Aug 12 2002

nonn

a(n) = A034387(2*n) - A034387(n-1); a(n) <= A179213(n). [Reinhard Zumkeller, Jul 05 2010]

			a(7) = 31 = 7+11+13 (sum of primes between 7 and 14).

T. D. Noe, Table of n, a(n) for n = 1..10000

Cf. A073838.

for n from 1 to 150 do l := 0:for j from n to 2*n do if isprime(j) then l := l+j:fi:od:a[n] := l:od:seq(a[j],j=1..150);

Table[Total[Select[Range[n, 2 n], PrimeQ]], {n, 56}] (* Jayanta Basu, Aug 12 2013 *)

PARI
```
a(n)=sum(i=n,2*n,i*isprime(i))
```

More terms from Sascha Kurz and Benoit Cloitre, Aug 14 2002

A179211 Number of squarefree numbers between n and 2*n (inclusive).

Original entry on oeis.org

2, 2, 3, 3, 4, 4, 5, 5, 6, 7, 8, 8, 9, 8, 9, 9, 11, 11, 13, 13, 15, 15, 15, 15, 15, 16, 16, 17, 19, 19, 20, 19, 21, 21, 22, 22, 24, 23, 24, 24, 25, 25, 26, 26, 27, 28, 29, 29, 30, 30, 32, 32, 34, 34, 36, 36, 38, 38, 38, 38, 39, 39, 38, 39, 41, 41, 42, 41, 43, 43, 44, 44, 46, 45, 45

Offset: 1

Reinhard Zumkeller, Jul 05 2010

nonn

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Cf. A005117, A008966, A013928, A035250, A061398, A179212, A179213, A179214.

a := n -> nops(select(issqrfree, [$n..(2*n)])):
seq(a(n), n=1..75); # Peter Luschny, Mar 02 2017

a[n_] := Select[Range[n, 2n], SquareFreeQ] // Length;
Array[a, 75] (* Jean-François Alcover, Jun 22 2018 *)

f(n)=my(s); forfactored(k=1,sqrtint(n), s += n\k[1]^2*moebius(k)); s
a(n)=f(2*n)-f(n-1) \\ Charles R Greathouse IV, Nov 05 2017

a(n) = Sum_{k=n..2*n} A008966(k).
a(n) > A035250(n) for n>2;
A179212(n) = a(n+1) - a(n);
a(n) = A013928(2*n+1) - A013928(n).
a(n) ~ (6/Pi^2) * n. - Amiram Eldar, Mar 03 2021

A066779 Sum of squarefree numbers <= n.

Original entry on oeis.org

1, 3, 6, 6, 11, 17, 24, 24, 24, 34, 45, 45, 58, 72, 87, 87, 104, 104, 123, 123, 144, 166, 189, 189, 189, 215, 215, 215, 244, 274, 305, 305, 338, 372, 407, 407, 444, 482, 521, 521, 562, 604, 647, 647, 647, 693, 740, 740, 740, 740, 791, 791, 844, 844, 899, 899

Offset: 1

Benoit Cloitre, Jan 18 2002

D. Suryanarayana, The number and sum of k-free integers <= x which are prime to n, Indian J. Math., Vol. 11 (1969), pp. 131-139.

Harry J. Smith, Table of n, a(n) for n = 1..1000

Cf. A008966, A013928, A034387, A179213, A179215.

Table[ n*Boole[ SquareFreeQ[n] ], {n, 1, 56}] // Accumulate (* Jean-François Alcover, Jun 18 2013 *)

s=0; for (n=1, 1000, write("b066779.txt", n, " ", s+=moebius(n)^2*n) ) \\ Harry J. Smith, Mar 24 2010

a(n)=sum(d=1,sqrtint(n),moebius(d)*d^2*binomial(n\d^2+1,2)) \\ Charles R Greathouse IV, Apr 26 2012

a(n)=my(s,k2); forsquarefree(k=1,sqrtint(n), k2=k[1]^2; s+= k2*binomial(n\k2+1,2)*moebius(k)); s \\ Charles R Greathouse IV, Jan 08 2018

from sympy.ntheory.factor_  import core
def a(n): return sum ([i for i in range(1, n + 1) if core(i) == i]) # Indranil Ghosh, Apr 16 2017

a(n) = Sum_{i=1..n} mu(i)^2*i.
a(n) = Sum_{k=1..n} k*A008966(k). - Reinhard Zumkeller, Jul 05 2010
a(n) = Sum_{d=1..sqrt(n)} mu(d)*d^2*floor(n/d^2)*floor(n/d^2+1)/2. - Charles R Greathouse IV, Apr 26 2012
G.f.: Sum_{k>=1} mu(k)^2*k*x^k/(1 - x). - Ilya Gutkovskiy, Apr 16 2017
a(n) ~ (3/Pi^2) * n^2 + O(n^(3/2)) (Suryanarayana, 1969). - Amiram Eldar, Mar 07 2021

cp's OEIS Frontend

A008966 a(n) = 1 if n is squarefree, otherwise 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

MuPAD

PARI

PARI

Python

Formula

Extensions

A045943 Triangular matchstick numbers: a(n) = 3*n*(n+1)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Scala

Formula

A073837 Sum of primes p satisfying n <= p <= 2n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A179211 Number of squarefree numbers between n and 2*n (inclusive).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A066779 Sum of squarefree numbers <= n.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Python

Formula

A179214 Product of squarefree numbers between n and 2*n (inclusive).

Views

Author

Keywords

Comments

A045943 Triangular matchstick numbers: a(n) = 3n(n+1)/2.