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

A344380 Complement of A344378 in A172186.

Original entry on oeis.org

6, 14, 38, 42, 57, 65, 70, 93, 106, 114, 118, 138, 154, 158, 182, 186, 190, 205, 210, 217, 218, 222, 266, 277, 281, 285, 309, 326, 334, 366, 381, 390, 393, 394, 397, 398, 401, 406, 434, 457, 469, 473, 478, 493, 498, 505, 518, 542, 561, 570, 581, 606, 614, 618
Offset: 1

Views

Author

René Gy, May 16 2021

Keywords

Comments

Terms belong to A172186 but not to A344378. Even though a(n)*(a(n)+1)*(2*a(n)+1) is squarefree, Sum_{j=1..a(n)} j^(2k) always has a prime divisor which is smaller than 2*a(n)+3, whatever k. For the integers m such that m*(m+1)*(2*m+1) is nonsquarefree, Sum_{j=1..m} j^(2k) always has a prime divisor which is smaller than 2*m+3, whatever k, because it is divisible by any prime p such that p^2 divides m*(m+1)*(2*m+1).

Examples

			14 belongs to the sequence, because it is squarefree, and Sum_{j=1..14} j^(2k) is always divisible by 29 when 14 does not divide k, and when 14 divides k, it is divisible by 13 or by 7.
		

Crossrefs

Extensions

More terms added and incorrect Mathematica program removed by Jinyuan Wang, Mar 07 2025

A007674 Numbers m such that m and m+1 are squarefree.

Original entry on oeis.org

1, 2, 5, 6, 10, 13, 14, 21, 22, 29, 30, 33, 34, 37, 38, 41, 42, 46, 57, 58, 61, 65, 66, 69, 70, 73, 77, 78, 82, 85, 86, 93, 94, 101, 102, 105, 106, 109, 110, 113, 114, 118, 122, 129, 130, 133, 137, 138, 141, 142, 145
Offset: 1

Views

Author

Keywords

Comments

m and m+1 squarefree implies that m*(m+1) is a squarefree oblong number and that m*(m+1)/2 is a squarefree triangular number. - Daniel Forgues, Aug 18 2012
Numbers m such that A002378(m) is squarefree. - Thomas Ordowski, Sep 01 2015

References

  • P. R. Halmos, Problems for Mathematicians Young and Old. Math. Assoc. America, 1991, p. 28.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    ff = {}; gg = {}; Do[kk = FactorInteger[n]; tak = False; Do[If[kk[[m]][[2]] > 1, tak = True], {m, 1, Length[kk]}]; If[tak == False, jj = FactorInteger[n + 1]; tak1 = False; Do[If[jj[[m]][[2]] > 1, tak1 = True], {m, 1, Length[jj]}]; If[tak1 == False, AppendTo[ff, n]]], {n, 1, 500}]; ff (* Artur Jasinski, Jan 28 2010 *)
    Select[Range[400],SquareFreeQ[#(#+1)]&] (* Vladimir Joseph Stephan Orlovsky, Mar 30 2011 *)
  • PARI
    list(lim)=my(v=vectorsmall(lim\1,i,1),u=List()); for(n=2, sqrt(lim), forstep(i=n^2,lim,n^2, v[i]=v[i-1]=0)); for(i=1,lim, if(v[i], listput(u,i))); v=0; Vec(u) \\ Charles R Greathouse IV, Aug 10 2011

Formula

A008966(a(n))*A008966(a(n)+1) = 1. - Reinhard Zumkeller, Dec 03 2009
a(n) ~ k*n, where k = 1/A065474. This result is originally due to Carlitz; for the (current) best error term, see Reuss. - Charles R Greathouse IV, Aug 10 2011, expanded Sep 18 2019

Extensions

Initial 1 added at the suggestion of Zak Seidov, Sep 19 2007

A179017 Odd numbers c such that c*(c^2 - 1)/4 is squarefree.

Original entry on oeis.org

3, 5, 11, 13, 21, 29, 43, 59, 61, 67, 69, 77, 83, 85, 93, 115, 123, 131, 133, 139, 141, 155, 157, 165, 173, 187, 203, 205, 211, 213, 219, 221, 227, 229, 237, 259, 267, 277, 283, 285, 291, 309, 317, 331, 347, 355, 357, 365, 371, 373, 381, 389, 403, 411, 419, 421
Offset: 1

Views

Author

Artur Jasinski, Jun 24 2010

Keywords

Comments

Original title was: "Numbers c such that (c^2-1)c is square free and gcd(c-1,c,c+1)=1", but (c^2-1)c is never squarefree for odd c, and gcd(n,n+1) is always = 1. - M. F. Hasler, Nov 03 2013
These numbers c with distribution a+b=c such that a=(c-1)/2 (see A172186) and b=(c+1)/2 (see A179019) have minimal possible values with function L(a,b,c) = log(c)/log(N(a,b,c)) = log(c)/log((c^2-1)c/4).
This function is minimal orbital in hypothesis (a,b,c).
There are no numbers or distributions which have value L less than log(c)/log((c^2-1)*c/4).
Equivalently, odd squarefree numbers c such that (c^2 - 1)/4 is also squarefree. - Jon E. Schoenfield, Feb 13 2023
The asymptotic density of this sequence is Product_{p prime} (1 - 3/p^2) = A206256 = 0.125486980905... (Tsang, 1985). - Amiram Eldar, Feb 26 2024

Crossrefs

Programs

  • Mathematica
    aa = {}; Do[If[(GCD[x, (x - 1)/2] == 1) && (GCD[x, (x + 1)/2] == 1) && (GCD[(x - 1)/2, (x + 1)/2] == 1), If[SquareFreeQ[(x^2 - 1) x/4], AppendTo[aa, x]]], {x, 2, 1000}]; aa
  • PARI
    forstep(n=3,421,2,issquarefree(n*(n^2-1)/4)&&print1(n",")) \\ M. F. Hasler, Nov 03 2013
    
  • PARI
    is(n)=n%2 && issquarefree(n) && issquarefree(n^2\4) \\ Charles R Greathouse IV, Mar 11 2014

Formula

a(n) = 2*A172186(n) + 1. - Bernard Schott, Mar 06 2023

Extensions

Edited by M. F. Hasler, Nov 03 2013

A172187 Numbers k such that k and k+1 are squarefree but 2*k+1 is not.

Original entry on oeis.org

13, 22, 37, 58, 73, 85, 94, 122, 130, 137, 157, 166, 181, 193, 202, 229, 237, 238, 253, 262, 265, 301, 302, 310, 318, 346, 373, 382, 409, 418, 433, 437, 445, 454, 462, 465, 481, 514, 517, 526, 537, 541, 553, 562, 589, 598, 634, 661, 662, 670, 697, 706, 733
Offset: 1

Views

Author

Artur Jasinski, Jan 28 2010

Keywords

Crossrefs

Complement of A007674 and A172186.

Programs

  • Mathematica
    Select[Range[750], And @@ SquareFreeQ /@ {#, # + 1} && !SquareFreeQ[2# + 1] &] (* Amiram Eldar, Mar 22 2020 *)

A119356 Numbers k such that A000330(k) is squarefree.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 10, 11, 14, 17, 18, 19, 20, 21, 22, 28, 29, 30, 33, 34, 35, 36, 38, 41, 42, 43, 44, 45, 46, 51, 52, 57, 58, 59, 61, 65, 66, 68, 69, 70, 76, 77, 78, 82, 83, 85, 86, 89, 90, 91, 92, 93, 101, 102, 105, 106, 109, 110, 113, 114, 115, 116, 117, 118, 123, 126
Offset: 1

Views

Author

Zak Seidov, May 16 2006

Keywords

Comments

The numbers of terms that do not exceed 10^k, for k = 1, 2, ..., are 8, 53, 504, 5029, 50187, 501925, 5019527, 50194688, 501948054, 5019478733, ... . Conjecture: The asymptotic density of this sequence is 4 * Product_{p prime} (1 - 3/p^2) = 4 * A206256 = 0.50194792... . - Amiram Eldar, Sep 24 2024

Examples

			10 is a term because 10*11*(2*10+1)/6 = 5*7*11 is squarefree.
		

Crossrefs

Cf. A000330, A005117, A172186 (subsequence), A206256,

Programs

  • Maple
    filter:= n -> numtheory:-issqrfree(n*(n+1)*(2*n+1)/6):
    select(filter, [$1..200]); # Robert Israel, Aug 04 2020
  • Mathematica
    Select[Range[126],SquareFreeQ[#(#+1)(2#+1)/6]&] (* James C. McMahon, Sep 15 2024 *)
  • PARI
    lista(nn) = {for (n=1, nn, if (issquarefree(n*(n+1)*(2*n+1)/6), print1(n, ", ")););} \\ Michel Marcus, May 18 2013

A179019 a(n) = (A179017(n)+1)/2.

Original entry on oeis.org

2, 3, 6, 7, 11, 15, 22, 30, 31, 34, 35, 39, 42, 43, 47, 58, 62, 66, 67, 70, 71, 78, 79, 83, 87, 94, 102, 103, 106, 107, 110, 111, 114, 115, 119, 130, 134, 139, 142, 143, 146, 155, 159, 166, 174, 178, 179, 183, 186, 187, 191, 195, 202, 206, 210, 211, 214, 215, 218, 219
Offset: 1

Views

Author

Artur Jasinski, Jun 24 2010

Keywords

Comments

For numbers a and c, see A172186 and A179017. Numbers b are this sequence.
These numbers c, with distribution a+b=c such that a=(c-1)/2 and b=(c+1)/2, have minimal possible values with function L(a,b,c) = log(c)/log(N[a,b,c]) = log(c)/log((c^2-1)c/4).
There exist no numbers or distributions for which L < log(c)/log((c^2-1)c/4). - Artur Jasinski

Crossrefs

Programs

  • Mathematica
    aa = {}; Do[If[(GCD[x, (x - 1)/2] == 1) && (GCD[x, (x + 1)/2] == 1) && (GCD[(x - 1)/2, (x + 1)/2] == 1), If[SquareFreeQ[(x^2 - 1) x/4], AppendTo[aa, (x + 1)/2]]], {x, 2, 1000}]; aa

Formula

a(n) = A179017(n) - A172186(n). - Hugo Pfoertner, Mar 22 2020

A344378 Positive integers m for which there exists a positive even integer 2k such that Sum_{j=1..m} j^(2k) has no prime divisor smaller than 2*m + 3.

Original entry on oeis.org

1, 2, 5, 10, 21, 29, 30, 33, 34, 41, 46, 61, 66, 69, 77, 78, 82, 86, 101, 102, 105, 109, 110, 113, 129, 133, 141, 142, 145, 165, 173, 177, 178, 185, 194, 201, 209, 213, 214, 221, 226, 230, 246, 254, 257, 258, 273, 282, 286, 290, 298, 313, 317, 321, 322, 329, 330
Offset: 1

Views

Author

René Gy, May 16 2021

Keywords

Comments

a(n)*(a(n)+1)*(2a(n)+1) must be squarefree, so this is a subsequence of A172186. It is the complement of A344380 in A172186.
Let L= LCM_j[(p_j-1)/2], where p_j run through the set of the prime divisors of a(n)*(a(n)+1)*(2a(n)+1). For a given member a(n) any admissible k must be a multiple of L, and for any prime p smaller than a(n) such that (p-1)/2 divides L, it holds that p does not divide a(n)-Floor[a(n)/p]. But the converse is not true: 397 is squarefree and satisfies the former condition, but Sum_{j=1..397} j^(2k) is always divisible either by 17 or by 73. 397 is the smallest "false positive" with the above test. Other "false positives" are rather scarce: 397,469,478,561,885,1002,1554,1658,1702,1977,... - René Gy, Apr 15 2025

Examples

			2 belongs to the sequence since 1 + 2^(2*2) = 17 is a prime number which is larger than 2*2 + 1 = 5.
5 belongs to the sequence because 1 + 2^20 + 3^20 + 4^20 + 5^20 = 96470431101379 = 137*704163730667 has no prime divisor smaller than 2*5 + 3 = 13.
		

Crossrefs

Extensions

Incorrect Mathematica program removed by Jinyuan Wang, Mar 18 2025
Showing 1-7 of 7 results.