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.

A331903 Positions of records in A331859.

Original entry on oeis.org

1, 2, 4, 5, 6, 8, 10, 12, 14, 17, 20, 23, 26, 29, 33, 36, 40, 45, 49, 53, 58, 63, 68, 74, 79, 85, 91, 97, 104, 110, 117, 124, 131, 139, 146, 154, 162, 170, 179, 187, 196, 205, 214, 224, 233, 243, 253, 263, 274, 284, 295, 306, 318, 329, 341, 353, 365, 377, 389
Offset: 1

Views

Author

Peter Kagey, Jan 31 2020

Keywords

Crossrefs

A331904 Number of occurrences of n in A331859.

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 3, 4, 5, 4, 4, 5, 5, 5, 6, 5, 6, 6, 6, 7, 6, 7, 7, 7, 8, 7, 8, 8, 8, 9, 8, 9, 9, 9, 10, 9, 10, 10, 10, 11, 10, 11, 11, 12, 11, 12, 12, 12, 12, 12, 13, 13, 13, 13, 14, 13, 14, 14, 15, 14, 15, 15, 15, 15, 16, 15
Offset: 1

Views

Author

Peter Kagey, Jan 31 2020

Keywords

Crossrefs

A332045 Numbers k such that ceiling(Pi/arctan(1/k)) = ceiling(k*Pi)+1.

Original entry on oeis.org

6, 7, 14, 21, 28, 113, 226, 339, 452, 565, 678, 791, 904, 1017, 1130, 1243, 1356, 1469, 1582, 1695, 1808, 1921, 33215, 99532, 364913, 729826, 1725033, 3450066, 5175099, 27235615, 52746197, 131002976, 471265707, 811528438, 2774848045, 4738167652, 567663097408
Offset: 1

Views

Author

Jianing Song, Feb 05 2020

Keywords

Comments

Note that ceiling(Pi/arctan(1/k)) - ceiling(k*Pi) is equal to either 0 or 1, that is, for all other k we have ceiling(Pi/arctan(1/k)) = ceiling(k*Pi).
Numbers k such that there exists some integer m such that Pi/arctan(1/k) > m > k*Pi.
Numbers k such that A331859(k^2) = A121854(k^2)+1 = A121855(k^2).
In A331859 there is a remark that A331859(100^n) = A011545(n). I'm in doubt of this, because if k = 10^n is here, then A331859(100^n) = ceiling(k*Pi), while A011545(n) = ceiling(k*Pi)-1, this equality would be violated.
Note that for k >= 3 we have 1/k < Pi/arctan(1/k)-k*Pi < (Pi/3)/k. As a result, a necessary condition for k being a term here is that there exists some m such that 0 < m/k - Pi < (Pi/3)/k^2, and a sufficient condition is that there exists some m such that 0 < m/k - Pi < 1/k^2.
Let P(n) = A002485(n), Q(n) = A002486(n), then it is known that 1/(Q(n)*(Q(n)*Q(n+1))) < |P(n)/Q(n) - Pi| < 1/(Q(n)*Q(n+1)) for n >= 2; furthermore, P(n)/Q(n) - Pi is positive for odd n and negative for even n. As a result, let n >= 3, then we have:
- If n is even, then Q(n) can never be a term.
- If n is odd, then k = Q(n)*t is a term if t <= sqrt(Q(n+1)/Q(n)), in which case ceiling(Pi/arctan(1/k)) = P(n)*t+1 and ceiling(k*Pi) = P(n)*t. The converse is not true (e.g., n = 3, t = 4745).

Examples

			Pi/arctan(1/6) = 19.0228..., 6*Pi = 18.8495..., so 6 is a term.
113*t is here for t <= 17, because ceiling(Pi/arctan(1/(113*t))) = 355*t+1 and ceiling((113*t)*Pi) = 355*t.
		

Crossrefs

Cf. A121381 (ceiling(n*Pi)), A121854 (floor(sqrt(n)*Pi)), A121855 (ceiling(sqrt(n)*Pi)), A011545 (floor(10^n*Pi)).

Programs

  • Magma
    // See Schoenfield link.
  • PARI
    default(realprecision, 10000); isA332045(n) = ceil(Pi/atan(1/n))!=ceil(n*Pi)
    

Extensions

a(27)-a(32) from Jon E. Schoenfield, Feb 12 2020
a(33)-a(36) from Giovanni Resta, Feb 12 2020
a(37) from Jon E. Schoenfield, Feb 15 2020
Showing 1-3 of 3 results.