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.

Previous Showing 21-22 of 22 results.

A247884 Number of positive integers < 10^n divisible by their first digit.

Original entry on oeis.org

9, 41, 327, 3158, 31450, 314349, 3143320, 31433005, 314329833, 3143298089, 31432980631, 314329806030, 3143298060001, 31432980599686, 314329805996514, 3143298059964770, 31432980599647312, 314329805996472711, 3143298059964726682, 31432980599647266367
Offset: 1

Views

Author

Derek Orr, Sep 25 2014

Keywords

Comments

a(n)/10^n seems to converge to a number around .3143...
a(n)/10^n converges to 7129/22680. - Hiroaki Yamanouchi, Sep 26 2014

Crossrefs

Programs

  • PARI
    a(n)=c=0;for(k=1,10^n-1,d=digits(k);if(k%d[1]==0,c++));c
    n=1;while(n<10,print1(a(n),", ");n++)
    
  • Python
    count = 9 # Start with the first 9 digits
    print(1, 9)
    n = 2
    while n < 101:
        for a in range(1, 10):
            count += 10**(n-1)//a
            if 10**(n-1) % a != 0:
                count += 1
        print(n, count)
        n += 1
    # David Consiglio, Jr., Sep 26 2014

Formula

G.f.: x*(9 - 67*x + 24*x^2 + 14*x^3 - 56*x^4 + 21*x^5 + 7*x^6 + 5*x^7)/((1 - x)^2*(1 + x)*(1 - 10*x)*(1 - x + x^2)). - Robert Israel, Mar 10 2025

Extensions

a(9)-a(20) from Hiroaki Yamanouchi, Sep 26 2014

A275591 a(n) = n^2 + 9*n + 1.

Original entry on oeis.org

1, 11, 23, 37, 53, 71, 91, 113, 137, 163, 191, 221, 253, 287, 323, 361, 401, 443, 487, 533, 581, 631, 683, 737, 793, 851, 911, 973, 1037, 1103, 1171, 1241, 1313, 1387, 1463, 1541, 1621, 1703, 1787, 1873, 1961, 2051, 2143, 2237, 2333, 2431, 2531, 2633, 2737
Offset: 0

Views

Author

Miquel Cerda, Aug 02 2016

Keywords

Comments

Also, nonnegative integers m such that 4*m + 77 is a square. The negative values of m are -7, -13, -17, -19.
The product of two consecutive terms belongs to the sequence. In fact: a(k)*a(k+1) = a(k*(k+1)+9*k+1).

Crossrefs

Cf. A028569.
Subsequence of A007775.

Programs

Formula

O.g.f.: (1 + 8*x - 7*x^2)/(1 - x)^3. - Colin Barker, Aug 03 2016
E.g.f.: (1 + 10*x + x^2)*exp(x).
a(n) = a(-n-9) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Colin Barker, Aug 03 2016
a(n) = A048058(n-1) + A008592(n-1) for n>0.
a(n) = 1 + A028569(n). - Omar E. Pol, Aug 02 2016
a(n) + a(-n) = (n-1)^2 + (n+1)^2.
Sum_{i>=0} 1/a(i) = 9736/29393 + tan(sqrt(77)*Pi/2)*Pi/sqrt(77) = 1.301517...

Extensions

Edited and extended by Bruno Berselli, Aug 05 2016
Previous Showing 21-22 of 22 results.