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 16 results. Next

A262783 Terms of A173143 (partial sums of squarefree numbers) that are square.

Original entry on oeis.org

1, 144, 6084, 33856, 37636, 6310144, 10150596, 14668900, 20007729, 3292234884, 39600602001, 92437853419755481, 2849708967414258609
Offset: 1

Views

Author

Michel Marcus, Oct 01 2015

Keywords

Comments

Intersection of A000290 and A173143.
a(12) > 10^15. - Jon E. Schoenfield, Oct 11 2015
a(14) > 3*10^21. The first 13 terms are the sums of the first 1, 21, 141, 331, 354, 4555, 5785, 6949, 8110, 104070, 360930, 551460217, 3061889301 squarefree numbers. - Giovanni Resta, Aug 19 2018

Examples

			The sum of the first 14 squarefree integers is 1 + 2 + 3 + 5 + 6 + 7 + 10 + 11 + 13 + 14 + 15 + 17 + 19 + 21 = 144, a square.
		

Crossrefs

Programs

  • PARI
    lista(nn) = {s = 0; i = 0; for (n=1, nn, if (issquarefree(n), i++; s += n; if (issquare(s), print1(s, ", "));););}

Extensions

a(12)-a(13) from Giovanni Resta, Aug 19 2018

A005117 Squarefree numbers: numbers that are not divisible by a square greater than 1.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103, 105, 106, 107, 109, 110, 111, 113
Offset: 1

Views

Author

Keywords

Comments

1 together with the numbers that are products of distinct primes.
Also smallest sequence with the property that a(m)*a(k) is never a square for k != m. - Ulrich Schimke (ulrschimke(AT)aol.com), Dec 12 2001
Numbers k such that there is only one Abelian group with k elements, the cyclic group of order k (the numbers such that A000688(k) = 1). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 25 2001
Numbers k such that A007913(k) > phi(k). - Benoit Cloitre, Apr 10 2002
a(n) is the smallest m with exactly n squarefree numbers <= m. - Amarnath Murthy, May 21 2002
k is squarefree <=> k divides prime(k)# where prime(k)# = product of first k prime numbers. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 30 2004
Numbers k such that omega(k) = Omega(k) = A072047(k). - Lekraj Beedassy, Jul 11 2006
The LCM of any finite subset is in this sequence. - Lekraj Beedassy, Jul 11 2006
This sequence and the Beatty Pi^2/6 sequence (A059535) are "incestuous": the first 20000 terms are bounded within (-9, 14). - Ed Pegg Jr, Jul 22 2008
Let us introduce a function D(n) = sigma_0(n)/2^(alpha(1) + ... + alpha(r)), sigma_0(n) number of divisors of n (A000005), prime factorization of n = p(1)^alpha(1) * ... * p(r)^alpha(r), alpha(1) + ... + alpha(r) is sequence (A001222). Function D(n) splits the set of positive integers into subsets, according to the value of D(n). Squarefree numbers (A005117) has D(n)=1, other numbers are "deviated" from the squarefree ideal and have 0 < D(n) < 1. For D(n)=1/2 we have A048109, for D(n)=3/4 we have A060687. - Ctibor O. Zizka, Sep 21 2008
Numbers k such that gcd(k,k')=1 where k' is the arithmetic derivative (A003415) of k. - Giorgio Balzarotti, Apr 23 2011
Numbers k such that A007913(k) = core(k) = k. - Franz Vrabec, Aug 27 2011
Numbers k such that sqrt(k) cannot be simplified. - Sean Loughran, Sep 04 2011
Indices m where A057918(m)=0, i.e., positive integers m for which there are no integers k in {1,2,...,m-1} such that k*m is a square. - John W. Layman, Sep 08 2011
It appears that these are numbers j such that Product_{k=1..j} (prime(k) mod j) = 0 (see Maple code). - Gary Detlefs, Dec 07 2011. - This is the same claim as Mohammed Bouayoun's Mar 30 2004 comment above. To see why it holds: Primorial numbers, A002110, a subsequence of this sequence, are never divisible by any nonsquarefree number, A013929, and on the other hand, the index of the greatest prime dividing any n is less than n. Cf. A243291. - Antti Karttunen, Jun 03 2014
Conjecture: For each n=2,3,... there are infinitely many integers b > a(n) such that Sum_{k=1..n} a(k)*b^(k-1) is prime, and the smallest such an integer b does not exceed (n+3)*(n+4). - Zhi-Wei Sun, Mar 26 2013
The probability that a random natural number belongs to the sequence is 6/Pi^2, A059956 (see Cesàro reference). - Giorgio Balzarotti, Nov 21 2013
Booker, Hiary, & Keating give a subexponential algorithm for testing membership in this sequence without factoring. - Charles R Greathouse IV, Jan 29 2014
Because in the factorizations into prime numbers these a(n) (n >= 2) have exponents which are either 0 or 1 one could call the a(n) 'numbers with a fermionic prime number decomposition'. The levels are the prime numbers prime(j), j >= 1, and the occupation numbers (exponents) e(j) are 0 or 1 (like in Pauli's exclusion principle). A 'fermionic state' is then denoted by a sequence with entries 0 or 1, where, except for the zero sequence, trailing zeros are omitted. The zero sequence stands for a(1) = 1. For example a(5) = 6 = 2^1*3^1 is denoted by the 'fermionic state' [1, 1], a(7) = 10 by [1, 0, 1]. Compare with 'fermionic partitions' counted in A000009. - Wolfdieter Lang, May 14 2014
From Vladimir Shevelev, Nov 20 2014: (Start)
The following is an Eratosthenes-type sieve for squarefree numbers. For integers > 1:
1) Remove even numbers, except for 2; the minimal non-removed number is 3.
2) Replace multiples of 3 removed in step 1, and remove multiples of 3 except for 3 itself; the minimal non-removed number is 5.
3) Replace multiples of 5 removed as a result of steps 1 and 2, and remove multiples of 5 except for 5 itself; the minimal non-removed number is 6.
4) Replace multiples of 6 removed as a result of steps 1, 2 and 3 and remove multiples of 6 except for 6 itself; the minimal non-removed number is 7.
5) Repeat using the last minimal non-removed number to sieve from the recovered multiples of previous steps.
Proof. We use induction. Suppose that as a result of the algorithm, we have found all squarefree numbers less than n and no other numbers. If n is squarefree, then the number of its proper divisors d > 1 is even (it is 2^k - 2, where k is the number of its prime divisors), and, by the algorithm, it remains in the sequence. Otherwise, n is removed, since the number of its squarefree divisors > 1 is odd (it is 2^k-1).
(End)
The lexicographically least sequence of integers > 1 such that each entry has an even number of proper divisors occurring in the sequence (that's the sieve restated). - Glen Whitney, Aug 30 2015
0 is nonsquarefree because it is divisible by any square. - Jon Perry, Nov 22 2014, edited by M. F. Hasler, Aug 13 2015
The Heinz numbers of partitions with distinct parts. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} prime(j) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] the Heinz number is 2*2*3*7*29 = 2436. The number 30 (= 2*3*5) is in the sequence because it is the Heinz number of the partition [1,2,3]. - Emeric Deutsch, May 21 2015
It is possible for 2 consecutive terms to be even; for example a(258)=422 and a(259)=426. - Thomas Ordowski, Jul 21 2015. [These form a subsequence of A077395 since their product is divisible by 4. - M. F. Hasler, Aug 13 2015]
There are never more than 3 consecutive terms. Runs of 3 terms start at 1, 5, 13, 21, 29, 33, ... (A007675). - Ivan Neretin, Nov 07 2015
a(n) = product of row n in A265668. - Reinhard Zumkeller, Dec 13 2015
Numbers without excess, i.e., numbers k such that A001221(k) = A001222(k). - Juri-Stepan Gerasimov, Sep 05 2016
Numbers k such that b^(phi(k)+1) == b (mod k) for every integer b. - Thomas Ordowski, Oct 09 2016
Boreico shows that the set of square roots of the terms of this sequence is linearly independent over the rationals. - Jason Kimberley, Nov 25 2016 (reference found by Michael Coons).
Numbers k such that A008836(k) = A008683(k). - Enrique Pérez Herrero, Apr 04 2018
The prime zeta function P(s) "has singular points along the real axis for s=1/k where k runs through all positive integers without a square factor". See Wolfram link. - Maleval Francis, Jun 23 2018
Numbers k such that A007947(k) = k. - Kyle Wyonch, Jan 15 2021
The Schnirelmann density of the squarefree numbers is 53/88 (Rogers, 1964). - Amiram Eldar, Mar 12 2021
Comment from Isaac Saffold, Dec 21 2021: (Start)
Numbers k such that all groups of order k have a trivial Frattini subgroup [Dummit and Foote].
Let the group G have order n. If n is squarefree and n > 1, then G is solvable, and thus by Hall's Theorem contains a subgroup H_p of index p for all p | n. Each H_p is maximal in G by order considerations, and the intersection of all the H_p's is trivial. Thus G's Frattini subgroup Phi(G), being the intersection of G's maximal subgroups, must be trivial. If n is not squarefree, the cyclic group of order n has a nontrivial Frattini subgroup. (End)
Numbers for which the squarefree divisors (A206778) and the unitary divisors (A077610) are the same; moreover they are also the set of divisors (A027750). - Bernard Schott, Nov 04 2022
0 = A008683(a(n)) - A008836(a(n)) = A001615(a(n)) - A000203(a(n)). - Torlach Rush, Feb 08 2023
From Robert D. Rosales, May 20 2024: (Start)
Numbers n such that mu(n) != 0, where mu(n) is the Möbius function (A008683).
Solutions to the equation Sum_{d|n} mu(d)*sigma(d) = mu(n)*n, where sigma(n) is the sum of divisors function (A000203). (End)
a(n) is the smallest root of x = 1 + Sum_{k=1..n-1} floor(sqrt(x/a(k))) greater than a(n-1). - Yifan Xie, Jul 10 2024
Number k such that A001414(k) = A008472(k). - Torlach Rush, Apr 14 2025
To elaborate on the formula from Greathouse (2018), the maximum of a(n) - floor(n*Pi^2/6 + sqrt(n)/17) is 10 at indices n = 48715, 48716, 48721, and 48760. The maximum is 11, at the same indices, if floor is taken individually for the two addends and the square root. If the value is rounded instead, the maximum is 9 at 10 indices between 48714 and 48765. - M. F. Hasler, Aug 08 2025

References

  • Jean-Marie De Koninck, Ces nombres qui nous fascinent, Entry 165, p. 53, Ellipses, Paris, 2008.
  • David S. Dummit and Richard M. Foote, Abstract algebra. Vol. 1999. Englewood Cliffs, NJ: David S.Prentice Hall, 1991.
  • Ivan M. Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers. 2nd ed., Wiley, NY, 1966, p. 251.
  • Michael Pohst and Hans J. Zassenhaus, Algorithmic Algebraic Number Theory, Cambridge Univ. Press, page 432.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Complement of A013929. Subsequence of A072774 and A209061.
Characteristic function: A008966 (mu(n)^2, where mu = A008683).
Subsequences: A000040, A002110, A235488.
Subsequences: numbers j such that j*a(k) is squarefree where k > 1: A056911 (k = 2), A261034 (k = 3), A274546 (k = 5), A276378 (k = 6).

Programs

  • Haskell
    a005117 n = a005117_list !! (n-1)
    a005117_list = filter ((== 1) . a008966) [1..]
    -- Reinhard Zumkeller, Aug 15 2011, May 10 2011
    
  • Magma
    [ n : n in [1..1000] | IsSquarefree(n) ];
    
  • Maple
    with(numtheory); a := [ ]; for n from 1 to 200 do if issqrfree(n) then a := [ op(a), n ]; fi; od:
    t:= n-> product(ithprime(k),k=1..n): for n from 1 to 113 do if(t(n) mod n = 0) then print(n) fi od; # Gary Detlefs, Dec 07 2011
    A005117 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if numtheory[issqrfree](a) then return a; end if; end do: end if; end proc:  # R. J. Mathar, Jan 09 2013
  • Mathematica
    Select[ Range[ 113], SquareFreeQ] (* Robert G. Wilson v, Jan 31 2005 *)
    Select[Range[150], Max[Last /@ FactorInteger[ # ]] < 2 &] (* Joseph Biberstine (jrbibers(AT)indiana.edu), Dec 26 2006 *)
    NextSquareFree[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sf = n + sgn; While[c < Abs[k], While[ ! SquareFreeQ@ sf, If[sgn < 0, sf--, sf++]]; If[ sgn < 0, sf--, sf++]; c++]; sf + If[ sgn < 0, 1, -1]]; NestList[ NextSquareFree, 1, 70] (* Robert G. Wilson v, Apr 18 2014 *)
    Select[Range[250], MoebiusMu[#] != 0 &] (* Robert D. Rosales, May 20 2024 *)
  • PARI
    bnd = 1000; L = vector(bnd); j = 1; for (i=1,bnd, if(issquarefree(i),L[j]=i; j=j+1)); L
    
  • PARI
    {a(n)= local(m,c); if(n<=1,n==1, c=1; m=1; while( cMichael Somos, Apr 29 2005 */
    
  • PARI
    list(n)=my(v=vectorsmall(n,i,1),u,j); forprime(p=2,sqrtint(n), forstep(i=p^2, n, p^2, v[i]=0)); u=vector(sum(i=1,n,v[i])); for(i=1,n,if(v[i],u[j++]=i)); u \\ Charles R Greathouse IV, Jun 08 2012
    
  • PARI
    for(n=1, 113, if(core(n)==n, print1(n, ", "))); \\ Arkadiusz Wesolowski, Aug 02 2016
    
  • PARI
    S(n) = my(s); forsquarefree(k=1,sqrtint(n),s+=n\k[1]^2*moebius(k)); s;
    a(n) = my(min=1, max=231, k=0, sc=0); if(n >= 144, min=floor(zeta(2)*n - 5*sqrt(n)); max=ceil(zeta(2)*n + 5*sqrt(n))); while(min <= max, k=(min+max)\2; sc=S(k); if(abs(sc-n) <= sqrtint(n), break); if(sc > n, max=k-1, if(sc < n, min=k+1, break))); while(!issquarefree(k), k-=1); while(sc != n, my(j=1); if(sc > n, j = -1); k += j; sc += j; while(!issquarefree(k), k += j)); k; \\ Daniel Suteu, Jul 07 2022
    
  • PARI
    first(n)=my(v=vector(n),i); forsquarefree(k=1,if(n<268293,(33*n+30)\20,(n*Pi^2/6+0.058377*sqrt(n))\1), if(i++>n, return(v)); v[i]=k[1]); v \\ Charles R Greathouse IV, Jan 10 2023
    
  • PARI
    A5117=[1..3]; A005117(n)={if(n>#A5117, my(N=#A5117); A5117=Vec(A5117, max(n+999, N*5\4)); iferr(forsquarefree(k=A5117[N]+1, #A5117*Pi^2\6+sqrtint(#A5117)\17+11, A5117[N++]=k[1]),E,)); A5117[n]} \\ M. F. Hasler, Aug 08 2025
    
  • Python
    from sympy.ntheory.factor_ import core
    def ok(n): return core(n, 2) == n
    print(list(filter(ok, range(1, 114)))) # Michael S. Branicky, Jul 31 2021
    
  • Python
    from itertools import count, islice
    from sympy import factorint
    def A005117_gen(startvalue=1): # generator of terms >= startvalue
        return filter(lambda n:all(x == 1 for x in factorint(n).values()),count(max(startvalue,1)))
    A005117_list = list(islice(A005117_gen(),20)) # Chai Wah Wu, May 09 2022
    
  • Python
    from math import isqrt
    from sympy import mobius
    def A005117(n):
        def f(x): return n+x-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 22 2024

Formula

Limit_{n->oo} a(n)/n = Pi^2/6 (see A013661). - Benoit Cloitre, May 23 2002
Equals A039956 UNION A056911. - R. J. Mathar, May 16 2008
A122840(a(n)) <= 1; A010888(a(n)) < 9. - Reinhard Zumkeller, Mar 30 2010
a(n) = A055229(A062838(n)) and a(n) > A055229(m) for m < A062838(n). - Reinhard Zumkeller, Apr 09 2010
A008477(a(n)) = 1. - Reinhard Zumkeller, Feb 17 2012
A055653(a(n)) = a(n); A055654(a(n)) = 0. - Reinhard Zumkeller, Mar 11 2012
A008966(a(n)) = 1. - Reinhard Zumkeller, May 26 2012
Sum_{n>=1} 1/a(n)^s = zeta(s)/zeta(2*s). - Enrique Pérez Herrero, Jul 07 2012
A056170(a(n)) = 0. - Reinhard Zumkeller, Dec 29 2012
A013928(a(n)+1) = n. - Antti Karttunen, Jun 03 2014
A046660(a(n)) = 0. - Reinhard Zumkeller, Nov 29 2015
Equals {1} UNION A000040 UNION A006881 UNION A007304 UNION A046386 UNION A046387 UNION A067885 UNION A123321 UNION A123322 UNION A115343 ... - R. J. Mathar, Nov 05 2016
|a(n) - n*Pi^2/6| < 0.058377*sqrt(n) for n >= 268293; this result can be derived from Cohen, Dress, & El Marraki, see links. - Charles R Greathouse IV, Jan 18 2018
From Amiram Eldar, Jul 07 2021: (Start)
Sum_{n>=1} (-1)^(a(n)+1)/a(n)^2 = 9/Pi^2.
Sum_{k=1..n} 1/a(k) ~ (6/Pi^2) * log(n).
Sum_{k=1..n} (-1)^(a(k)+1)/a(k) ~ (2/Pi^2) * log(n).
(all from Scott, 2006) (End)

A375707 First differences minus 1 of nonsquarefree numbers.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 16 2024

Keywords

Comments

Also the number of squarefree numbers between the nonsquarefree numbers A013929(n) and A013929(n+1).
Delete all 0's to get A120992.
The image is {0,1,2,3}.
Add 1 to all terms for A078147.

Examples

			The runs of squarefree numbers begin:
  (5,6,7)
  ()
  (10,11)
  (13,14,15)
  (17)
  (19)
  (21,22,23)
  ()
  (26)
  ()
  (29,30,31)
  (33,34,35)
		

Crossrefs

Positions of 0, 1, 2, 3 are A375709, A375710, A375711, A375712. This is a set partition of the positive integers into four blocks.
For runs of squarefree numbers:
- length: A120992, anti A373127
- min: A072284, anti A373408
- max: A373415, anti A007674
- sum: A373413, anti A373411
For runs of nonsquarefree numbers:
- length: A053797, anti A373409
- min: A053806, anti A373410
- max: A376164, anti A068781
- sum: A373414, anti A373412
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
A046933 counts composite numbers between consecutive primes.
A073784 counts primes between consecutive composite numbers.
A093555 counts non-prime-powers between consecutive prime-powers.

Programs

  • Mathematica
    Differences[Select[Range[100],!SquareFreeQ[#]&]]-1
  • PARI
    lista(nmax) = {my(prev = 4); for (n = 5, nmax, if(!issquarefree(n), print1(n - prev - 1, ", "); prev = n));} \\ Amiram Eldar, Sep 17 2024

Formula

Asymptotic mean: lim_{n->oo} (1/n) Sum_{k=1..n} a(k) = 6/(Pi^2-6) = 1.550546... . - Amiram Eldar, Sep 17 2024

A373411 Sum of the n-th maximal antirun of squarefree numbers differing by more than one.

Original entry on oeis.org

1, 2, 8, 6, 17, 24, 14, 72, 22, 78, 30, 64, 34, 72, 38, 80, 42, 89, 263, 58, 120, 127, 66, 136, 70, 144, 151, 78, 161, 168, 86, 360, 94, 293, 102, 208, 106, 216, 110, 224, 114, 233, 241, 379, 130, 264, 271, 138, 280, 142, 288, 600, 312, 158, 648, 166, 510, 351
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2024

Keywords

Comments

The length of this antirun is given by A373127.
An antirun of a sequence (in this case A005117) is an interval of positions at which consecutive terms differ by more than one.

Examples

			Row-sums of:
   1
   2
   3  5
   6
   7 10
  11 13
  14
  15 17 19 21
  22
  23 26 29
  30
  31 33
  34
  35 37
  38
  39 41
  42
  43 46
  47 51 53 55 57
		

Crossrefs

The partial sums are a subset of A173143.
Functional neighbors: A007674, A373127 (firsts A373128, sorted firsts A373200), A373404, A373405, A373408, A373412, A373413.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.

Programs

  • Mathematica
    Total/@Split[Select[Range[100],SquareFreeQ],#1+1!=#2&]//Most

A373413 Sum of the n-th maximal run of squarefree numbers.

Original entry on oeis.org

6, 18, 21, 42, 17, 19, 66, 26, 90, 102, 114, 126, 93, 51, 53, 55, 174, 123, 198, 210, 147, 234, 165, 258, 89, 91, 282, 97, 306, 318, 330, 342, 237, 245, 127, 390, 267, 414, 426, 291, 149, 151, 309, 474, 161, 163, 498, 170, 347, 534, 546, 558, 381, 582, 197
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2024

Keywords

Comments

The length of this run is given by A120992.
A run of a sequence (in this case A005117) is an interval of positions at which consecutive terms differ by one.

Examples

			Row-sums of:
   1   2   3
   5   6   7
  10  11
  13  14  15
  17
  19
  21  22  23
  26
  29  30  31
  33  34  35
  37  38  39
  41  42  43
  46  47
  51
  53
  55
  57  58  59
		

Crossrefs

The partial sums are a subset of A173143.
Functional neighbors: A054265, A072284, A120992, A373406, A373411, A373414, A373415.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.

Programs

  • Mathematica
    Total/@Split[Select[Range[100],SquareFreeQ],#1+1==#2&]//Most

A013932 Integers that are squarefree and also the sum of first k squarefrees for some k.

Original entry on oeis.org

1, 3, 6, 11, 17, 34, 58, 87, 123, 166, 215, 274, 305, 407, 482, 521, 562, 647, 791, 899, 1073, 1261, 1327, 1394, 1463, 1533, 1677, 1751, 1906, 1985, 2067, 2235, 2321, 2497, 2681, 2870, 2967, 3170, 3273, 3378, 3810, 3921, 4034, 4381, 4622, 4745, 5001, 5131, 5262
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    Select[Accumulate[Select[Range[150],SquareFreeQ]],SquareFreeQ] (* Harvey P. Dale, Jul 27 2011 *)
  • PARI
    lista(nn) = {my(s=0); for(k=1, nn, if(issquarefree(k)==1, s+=k; if(issquarefree(s)==1, print1(s, ", ")))); } \\ Jinyuan Wang, Feb 26 2020

Formula

a(n) = A005117(A020643(n)). - Amiram Eldar, Aug 19 2025

Extensions

More terms from Jinyuan Wang, Feb 26 2020

A329472 Partial sums of numbers that are not squarefree.

Original entry on oeis.org

4, 12, 21, 33, 49, 67, 87, 111, 136, 163, 191, 223, 259, 299, 343, 388, 436, 485, 535, 587, 641, 697, 757, 820, 884, 952, 1024, 1099, 1175, 1255, 1336, 1420, 1508, 1598, 1690, 1786, 1884, 1983, 2083, 2187, 2295, 2407, 2523, 2640, 2760, 2881, 3005, 3130, 3256
Offset: 1

Views

Author

Metin Sariyar, Nov 13 2019

Keywords

Examples

			49 is a term because 4, 8, 9, 12, 16 are the first five nonsquarefree numbers and 4 + 8 + 9 + 12 + 16 = 49.
		

Crossrefs

Partial sums of A013929.

Programs

  • Mathematica
    Accumulate[Select[Range[1000], !SquareFreeQ[#]&]]

Formula

a(n) ~ (Pi^2/(2*Pi^2 - 12)) * n^2. - Amiram Eldar, Oct 21 2020

A322608 Values of k such that (product of squarefree numbers <= k) / (sum of squarefree numbers <= k) is an integer.

Original entry on oeis.org

1, 3, 11, 14, 17, 21, 23, 33, 34, 37, 46, 47, 55, 58, 59, 61, 62, 67, 69, 73, 82, 83, 87, 94, 95, 97, 101, 106, 107, 109, 114, 115, 119, 123, 127, 133, 134, 141, 146, 151, 157, 158, 159, 161, 165, 166, 173, 181, 187, 197, 202, 203, 210, 218, 219, 223, 226, 230
Offset: 1

Views

Author

Paolo P. Lava, Dec 20 2018

Keywords

Examples

			3 is in the sequence because (1*2*3)/(1+2+3) = 1.
11 is in the sequence because (1*2*3*5*6*7*10*11)/(1+2+3+5+6+7+10+11) = 138600/45 = 3080.
		

Crossrefs

Cf. A005117, A051838, A111059, A116536, A141092, A173143, A322607 (values of the quotient), A347690 (numbers of terms in the numerators).

Programs

  • Maple
    with(numtheory): P:=proc(q) local a,b,c,n; a:=1; b:=0; c:=[];
    for n from 1 to q do if issqrfree(n) then a:=a*n; b:=b+n;
    if frac(a/b)=0 then c:=[op(c),n];
    fi; fi; od; op(c); end: P(60);
  • Mathematica
    seq = {}; sum = 0; prod = 1; Do[If[SquareFreeQ[n], sum += n; prod *= n; If[Divisible[prod, sum], AppendTo[seq, n]]], {n, 1, 230}]; seq (* Amiram Eldar, Mar 05 2021 *)

Extensions

Definition corrected by N. J. A. Sloane, Sep 19 2021 at the suggestion of Harvey P. Dale.

A322607 Numbers that can be expressed as the ratio between the product and the sum of consecutive squarefree numbers starting from 1.

Original entry on oeis.org

1, 3080, 350350, 61850250, 17823180375, 6871260396000, 88909822914869880000, 2746644314348614680000, 2980109081068246927800000, 9638057975990853416623724908800000, 424217819372970387341691005411520000, 51912228216508515627667235880347808000000, 152157812632066726080765311397008321568000000
Offset: 1

Views

Author

Paolo P. Lava, Dec 20 2018

Keywords

Examples

			1 is a term because 1/1 = (1*2*3)/(1+2+3) = 1.
3080 is a term because (1*2*3*5*6*7*10*11)/(1+2+3+5+6+7+10+11) = 138600/45 = 3080.
		

Crossrefs

Programs

  • Maple
    with(numtheory): P:=proc(q) local a,b,c,n; a:=1; b:=0; c:=[];
    for n from 1 to q do if issqrfree(n) then a:=a*n; b:=b+n;
    if frac(a/b)=0 then if n>1 then c:=[op(c),a/b];
    fi; fi; fi; od; op(c); end: P(60);

Formula

a(n) = A111059(A322608(n+1))/A173143(A322608(n+1)).

A336025 Numbers m providing record values for the proportion of nonsquarefree integers in the interval [1, m].

Original entry on oeis.org

4, 9, 20, 25, 27, 28, 172, 176
Offset: 1

Views

Author

Javier Múgica, Jul 05 2020

Keywords

Comments

Also, numbers providing record low values for the proportion of squarefree integers.
The proportion of nonsquarefree integers approaches 1-6/Pi^2. For low values of m the proportion in [1, m] tends to be lower, since squares appear late. But values of m for which the ratio in the interval [1, m] is larger than the limit value do exist. The first such one is 28. Therefore this sequence is finite and it can be proved that 176 is indeed its last term. The proportion of nonsquarefree numbers in [1, 176] is 70/176 = 0.397727272... and that of squarefree ones is 0.6022727...
Rogers (1964) proved that if A(x) is the number of squarefree numbers not exceeding x then A(x)/x >= A(176)/176 = 53/88. Therefore 176 is the last term of this sequence. - Amiram Eldar, Mar 12 2021

Examples

			Up to m = 9 there are 3 numbers which are divisible by some square: 4, 8 and 9, for a proportion of 3/9 = 1/3. No interval [1, k] for k < 9 has a ratio as high, so 9 is in the sequence.
		

Crossrefs

Programs

  • Mathematica
    s = {}; c = 0; rm = 0; Do[If[!SquareFreeQ[n], c++]; If[(r = c/n) > rm, rm = r; AppendTo[s, n]], {n, 1, 10^3}]; s (* Amiram Eldar, Mar 12 2021 *)
    Module[{nn=200},DeleteDuplicates[{#[[1]],#[[2]]/#[[1]]}&/@Thread[{Range[nn],Accumulate[ Table[If[SquareFreeQ[n],0,1],{n,nn}]]}],GreaterEqual[#1[[2]],#2[[2]]]&]][[;;,1]]//Rest (* Harvey P. Dale, Jun 15 2024 *)
Showing 1-10 of 16 results. Next