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.

A293435 a(n) is the number of the proper divisors of n that are Fibonacci numbers (A000045).

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Oct 09 2017

Keywords

Examples

			For n = 55, its proper divisors are [1, 5, 11], of which only two, namely 1 and 5 are in A000045, thus a(55) = 2.
		

Crossrefs

Programs

  • Mathematica
    With[{s = Fibonacci@ Range[2, 40]}, Table[DivisorSum[n, 1 &, And[MemberQ[s, #], # != n] &], {n, 105}]] (* Michael De Vlieger, Oct 09 2017 *)
  • PARI
    A010056(n) = { my(k=n^2); k+=(k+1)<<2; (issquare(k) || (n>0 && issquare(k-8))) }; \\ This function from Charles R Greathouse IV, Jul 30 2012
    A293435(n) = sumdiv(n,d,(dA010056(d));

Formula

a(n) = Sum_{d|n, dA010056(d).
a(n) = A005086(n) - A010056(n).
G.f.: Sum_{k>=2} x^(2*Fibonacci(k)) / (1 - x^Fibonacci(k)). - Ilya Gutkovskiy, Apr 14 2021
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A079586 - 1 = 2.359885... . - Amiram Eldar, Jul 05 2025

A293434 a(n) is the sum of the proper divisors of n that are Jacobsthal numbers (A001045).

Original entry on oeis.org

0, 1, 1, 1, 1, 4, 1, 1, 4, 6, 1, 4, 1, 1, 9, 1, 1, 4, 1, 6, 4, 12, 1, 4, 6, 1, 4, 1, 1, 9, 1, 1, 15, 1, 6, 4, 1, 1, 4, 6, 1, 25, 1, 12, 9, 1, 1, 4, 1, 6, 4, 1, 1, 4, 17, 1, 4, 1, 1, 9, 1, 1, 25, 1, 6, 15, 1, 1, 4, 6, 1, 4, 1, 1, 9, 1, 12, 4, 1, 6, 4, 1, 1, 25, 6, 44, 4, 12, 1, 9, 1, 1, 4, 1, 6, 4, 1, 1, 15, 6, 1, 4, 1, 1, 30
Offset: 1

Views

Author

Antti Karttunen, Oct 09 2017

Keywords

Examples

			For n = 15, whose proper divisors are [1, 3, 5], all of them are in A001045, thus a(15) = 1 + 3 + 5 = 9.
For n = 21, whose proper divisors are [1, 3, 7], both 1 and 3 are in A001045, thus a(21) = 1 + 3 = 4.
For n = 21845, whose proper divisors are [1, 5, 17, 85, 257, 1285, 4369], only 1, 5, 85 are in A001045, thus a(21845) = 1 + 5 + 85 = 91.
		

Crossrefs

Programs

  • Mathematica
    With[{s = LinearRecurrence[{1, 2}, {0, 1}, 24]}, Table[DivisorSum[n, # &, And[MemberQ[s, #], # != n] &], {n, 105}]] (* Michael De Vlieger, Oct 09 2017 *)
  • PARI
    A147612aux(n,i) = if(!(n%2),n,A147612aux((n+i)/2,-i));
    A147612(n) = 0^(A147612aux(n,1)*A147612aux(n,-1));
    A293434(n) = sumdiv(n,d,(dA147612(d)*d);
    
  • Python
    from sympy import divisors
    def A293434(n): return sum(d for d in divisors(n,generator=True) if d(m-3).bit_length()) # Chai Wah Wu, Apr 18 2025

Formula

a(n) = Sum_{d|n, dA147612(d)*d.
a(n) = A293432(n) - (A147612(n)*n).
Showing 1-2 of 2 results.