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

A226120 Decimal expansion of Sum_{n>=1} n^3/(exp(2*Pi*n/7)-1).

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 9, 0, 1, 6, 1, 7, 6, 7, 8, 8, 8, 6, 6, 2, 6, 7, 5, 5, 8, 4, 3, 5, 9, 3, 0, 5, 8, 5, 5, 4, 4, 5, 3, 3, 3, 4, 8, 0, 2, 5, 4, 8, 9, 7, 8, 4, 3, 4, 0, 6, 1, 0, 9, 9, 4, 3, 8, 7, 3, 7, 8, 5, 0, 6, 7, 1, 4, 8, 0, 1, 7, 9, 1, 6, 2, 7, 1, 3, 6, 6, 2, 1
Offset: 2

Views

Author

Jean-François Alcover, May 27 2013

Keywords

Comments

An almost-integer discovered by Simon Plouffe. The computed sum equals 10 within 15 digits.

Examples

			10.00000000000000019016176788866267558435930585544533348025489784340610994387...
		

Crossrefs

Cf. A060295 (famous almost-integer: Ramanujan's constant), A226121 (another surprising almost-integer by Simon Plouffe), A007775, A089034.

Programs

  • Mathematica
    NSum[n^3/(Exp[2*Pi*n/7] - 1), {n, 1, Infinity}, NSumTerms -> 220,    WorkingPrecision -> 100] // RealDigits[#, 10, 100] & // First

Extensions

Offset corrected by Rick L. Shepherd, Jan 01 2014

A306174 a(n) = (prime(n)^6 - 1)/504.

Original entry on oeis.org

3515, 9577, 47892, 93345, 293722, 1180205, 1760920, 5090727, 9424810, 12542387, 21387332, 43976907, 83691535, 102222965, 179480917, 254167230, 300266322, 482316380, 648691217, 986073990, 1652722232, 2106190775, 2369151382, 2977639587, 3327579585, 4130856652
Offset: 5

Views

Author

Jianing Song, Jul 03 2018

Keywords

Comments

Note that 504 = 7*8*9. For odd prime we have p^2 == 1 (mod 8). By Fermat's theorem and Euler's totient theorem we have p^6 == 1 (mod 7) for p != 7 and p^6 == 1 (mod 9) for p != 3. So 504 divides p^6 - 1 for p != 2, 3, 7.
There are no primes in this sequence except for a term not shown here, which is a(3) = (5^6 - 1)/504 = 31. Furthermore, omega(a(n)) = A001221(a(n)) >= 4 for n >= 7, and is exactly 4 for n = 7, 8, 9, 10, 13, 14, 16, 17, 23, ...
The set of prime factors of this sequence include all primes. First, a(7) is divisible by 2, and a(8) is divisible by 3 and 7. And also, for any prime q != 2, 3, 7, by Dirichlet's theorem on arithmetic progressions, there exists a prime p of the form k*q + d with d^6 == 1 (mod p), 0 < d < q. Since gcd(q,504) = 1, we have p^6 == d^6 == 1 (mod 504*q), so q is divisible by (p^6 - 1)/504.
Note that a(3)=31 for prime 5 is also an integer. - Michel Marcus, Jul 05 2018

Examples

			a(5) = (11^6 - 1)/504 = 3515, a(6) = (13^6 - 1)/504 = 9577, a(7) = (17^6 - 1)/504 = 47892, ...
		

Crossrefs

Programs

  • Magma
    [(NthPrime(n)^6 - 1) div 504: n in [5..40]]; // Vincenzo Librandi, Jul 13 2018
  • Mathematica
    Table[(Prime[n]^6 - 1) / 504, {n, 5, 40}] (* Vincenzo Librandi, Jul 13 2018 *)
  • PARI
    a(n)=(prime(n)^6 - 1)/504
    
Showing 1-2 of 2 results.