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-30 of 77 results. Next

A133623 Binomial(n+p, n) mod n where p=3.

Original entry on oeis.org

0, 0, 2, 3, 1, 0, 1, 5, 4, 6, 1, 11, 1, 8, 6, 9, 1, 16, 1, 11, 8, 12, 1, 21, 1, 14, 10, 15, 1, 26, 1, 17, 12, 18, 1, 31, 1, 20, 14, 21, 1, 36, 1, 23, 16, 24, 1, 41, 1, 26, 18, 27, 1, 46, 1, 29, 20, 30, 1, 51, 1, 32, 22, 33, 1, 56, 1, 35, 24, 36, 1, 61, 1, 38, 26, 39, 1, 66, 1, 41, 28, 42, 1
Offset: 1

Views

Author

Hieronymus Fischer, Sep 30 2007

Keywords

Comments

Let d(m)...d(2)d(1)d(0) be the base-n representation of n+p. The relation a(n)=d(1) holds, if n is a prime index. For this reason there are infinitely many terms which are equal to 1.
Appears to satisfy the recurrence: a(n) = -2*a(n-1) - a(n-2) + 2*a(n-3) + 4*a(n-4) + 2*a(n-5) - a(n-6) - 2*a(n-7) - a(n-8) for n > 14. - Chai Wah Wu, May 25 2016

Crossrefs

Programs

  • Mathematica
    Table[Mod[Binomial[n+3,n],n],{n,90}] (* Harvey P. Dale, Nov 22 2011 *)

Formula

a(n) = binomial(n+3,3) mod n.
a(n)=1 if n is a prime > 3, since binomial(n+3,n)==(1+floor(3/n))(mod n), provided n is a prime.
From Chai Wah Wu, May 26 2016: (Start)
a(n) = (n^3 + 5*n + 6)/6 mod n.
For n > 6:
if n mod 6 == 0, then a(n) = 5*n/6 + 1.
if n mod 6 is in {1, 5}, then a(n) = 1.
if n mod 6 is in {2, 4}, then a(n) = n/2 + 1.
if n mod 6 == 3, then a(n) = n/3 + 1.
(End)

A134334 Numbers which are not divisible by the number of their prime factors (counted with multiplicity).

Original entry on oeis.org

8, 9, 15, 20, 21, 25, 28, 32, 33, 35, 39, 44, 48, 49, 50, 51, 52, 54, 55, 57, 64, 65, 68, 69, 70, 72, 76, 77, 81, 85, 87, 90, 91, 92, 93, 95, 98, 108, 110, 111, 112, 115, 116, 119, 121, 123, 124, 125, 126, 128, 129, 130, 133, 135, 141, 143, 145, 148, 150, 154, 155, 159
Offset: 1

Views

Author

Hieronymus Fischer, Oct 23 2007

Keywords

Comments

The asymptotic density of this sequence is 1 (Erdős and Pomerance, 1990). - Amiram Eldar, Jul 10 2020

Examples

			a(1) = 8, since 8 = 2*2*2 has 3 prime factors and 8 is not divisible by 3.
a(3) = 15, since 15 = 3*5 has 2 prime factors and 15 is not divisible by 2.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2,200],Mod[#,PrimeOmega[#]]!=0&] (* Harvey P. Dale, May 13 2023 *)
  • PARI
    isok(n) = (n % bigomega(n)) \\ Michel Marcus, Jul 15 2013

A133622 a(n) = 1 if n is odd, a(n) = n/2+1 if n is even.

Original entry on oeis.org

1, 2, 1, 3, 1, 4, 1, 5, 1, 6, 1, 7, 1, 8, 1, 9, 1, 10, 1, 11, 1, 12, 1, 13, 1, 14, 1, 15, 1, 16, 1, 17, 1, 18, 1, 19, 1, 20, 1, 21, 1, 22, 1, 23, 1, 24, 1, 25, 1, 26, 1, 27, 1, 28, 1, 29, 1, 30, 1, 31, 1, 32, 1, 33, 1, 34, 1, 35, 1, 36, 1, 37, 1, 38, 1, 39, 1, 40, 1, 41, 1, 42, 1, 43, 1, 44, 1
Offset: 1

Views

Author

Hieronymus Fischer, Sep 30 2007

Keywords

Comments

a(n) is the count of terms a(n+1) present so far in the sequence, with a(n+1) included in the count; example: a(1) = 1 "says" that there is 1 term "2" so far in the sequence; a(2) = 2 "says" that there are 2 terms "1" so far in the sequence... etc. This comment was inspired by A039617. - Eric Angelini, Mar 03 2020

Crossrefs

Programs

  • Haskell
    import Data.List (transpose)
    a133622 n = (1 - m) * n' + 1 where (n', m) = divMod n 2
    a133622_list = concat $ transpose [[1, 1 ..], [2 ..]]
    -- Reinhard Zumkeller, Feb 20 2015
    
  • Maple
    seq([1,n][],n=2..100); # Robert Israel, May 27 2016
  • Mathematica
    Riffle[Range[2,50],1,{1,-1,2}] (* Harvey P. Dale, Jan 19 2013 *)
  • PARI
    a(n)=if(n%2,1,n/2+1) \\ Charles R Greathouse IV, Sep 02 2015

Formula

a(n)=1+(binomial(n+1,2)mod n)=1+(binomial(n+1,n-1)mod n).
a(n)=binomial(n+2,2) mod n = binomial(n+2,n) mod n for n>2.
a(n)=1+(1+(-1)^n)*n/4.
a(n)=1+(A000217(n) mod n).
a(n)=a(n-2)+1, if n is even, a(n)=a(n-2) if n is odd.
a(n)=a(n-2)+1-(n mod 2)=a(n-2)+(1+(-1)^n)/2 for n>2.
a(n)=(a(n-3)+a(n-2))/a(n-1) for n>3.
G.f.: g(x)=x(1+2x-x^2-x^3)/(1-x^2)^2.
G.f.: (Q(0)-1-x)/x^2, where Q(k)= 1 + (k+1)*x/(1 - x/(x + (k+1)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 23 2013
a(n) = 2*a(n-2)-a(n-4) for n > 4. - Chai Wah Wu, May 26 2016
E.g.f.: exp(x) - 1 + x*sinh(x)/2. - Robert Israel, May 27 2016

A133875 n modulo 5 repeated 5 times.

Original entry on oeis.org

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

Views

Author

Hieronymus Fischer, Oct 10 2007

Keywords

Comments

Periodic with length 5^2 = 25.

Crossrefs

Programs

  • Magma
    [(1 + Floor(n/5)) mod 5 : n in [0..50]]; // Wesley Ivan Hurt, Jun 06 2014
  • Maple
    A133875:=n->((1+floor(n/5)) mod 5); seq(A133875(n), n=0..100); # Wesley Ivan Hurt, Jun 06 2014
  • Mathematica
    Table[Mod[1 + Floor[n/5], 5], {n, 0, 100}] (* Wesley Ivan Hurt, Jun 06 2014 *)
    LinearRecurrence[{1,0,0,0,-1,1,0,0,0,-1,1,0,0,0,-1,1,0,0,0,-1,1},{1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4,0},120] (* Harvey P. Dale, Dec 14 2017 *)

Formula

a(n) = (1 + floor(n/5)) mod 5.
a(n) = A010874(A002266(n+5)).
a(n) = 1 + floor(n/5) - 5*floor((n+5)/25).
a(n) = (((n+5) mod 25) - (n mod 5)) / 5.
a(n) = ((n + 5 - (n mod 5)) / 5) mod 5.
a(n) = A010874((n + 5 - A010874(n))/5).
a(n) = binomial(n+5, n) mod 5 = binomial(n+5, 5) mod 5.
a(n) = +a(n-1) -a(n-5) +a(n-6) -a(n-10) +a(n-11) -a(n-15) +a(n-16) -a(n-20) +a(n-21). - R. J. Mathar, Sep 03 2011
G.f.: ( 1+2*x^5+3*x^10+4*x^15 ) / ( (1-x)*(x^20+x^15+x^10+x^5+1) ). - R. J. Mathar, Sep 03 2011

A135120 Numbers such that the digital sum base 2 and the digital sum base 3 and the digital sum base 10 all are equal.

Original entry on oeis.org

1, 21, 222, 223, 1230, 1231, 1502, 2200, 2201, 3012, 3013, 10431, 12214, 12215, 12250, 12251, 14102, 15003, 15021, 16011, 20040, 20041, 22130, 23211, 23230, 23231, 24003, 30070, 30071, 30105, 30231, 30321, 31005, 31150, 31151, 31420
Offset: 1

Views

Author

Hieronymus Fischer, Dec 24 2007

Keywords

Examples

			a(2)=21, since ds_2(21)=ds_3(21)=ds_10(21)=3.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[5000], Total[IntegerDigits[#, 2]] == Total[IntegerDigits[#, 3]] ==  Total[IntegerDigits[#, 10]] &] (* G. C. Greubel, Sep 26 2016 *)
  • PARI
    is(n)=my(t=sumdigits(n)); t==hammingweight(n) && t==sumdigits(n,3) \\ Charles R Greathouse IV, Sep 26 2016

A133885 Binomial(n+5,n) mod 5^2.

Original entry on oeis.org

1, 6, 21, 6, 1, 2, 12, 17, 12, 2, 3, 18, 13, 18, 3, 4, 24, 9, 24, 4, 5, 5, 5, 5, 5, 6, 11, 1, 11, 6, 7, 17, 22, 17, 7, 8, 23, 18, 23, 8, 9, 4, 14, 4, 9, 10, 10, 10, 10, 10, 11, 16, 6, 16, 11, 12, 22, 2, 22, 12, 13, 3, 23, 3, 13, 14, 9, 19, 9, 14, 15, 15, 15, 15, 15, 16, 21, 11, 21, 16, 17
Offset: 0

Views

Author

Hieronymus Fischer, Oct 10 2007

Keywords

Comments

Periodic with length 5^3=125.

Crossrefs

For the sequence regarding binomial(n+5, n) mod 5 see A133875.

Programs

  • Mathematica
    Table[Mod[Binomial[n+5,n],25],{n,0,90}] (* Harvey P. Dale, Jan 12 2023 *)

Formula

a(n)=binomial(n+5,5) mod 5^2.
G.f. g(x)=sum{0<=k<125, a(k)*x^k}/(1-x^125).

A135100 Numbers which divide their digital sumorial (see A131383).

Original entry on oeis.org

1, 3, 4, 15, 26, 2573, 17226, 19786, 22083, 58133, 67693, 223657, 376460, 464713, 497068, 2621204, 4553376, 6000136, 7671158, 13975944, 14074903, 52731198, 82594577
Offset: 1

Views

Author

Hieronymus Fischer, Dec 24 2007

Keywords

Examples

			a(5)=26, since 26 divides its digital sumorial, which is A131383(26)=182.
		

Crossrefs

Programs

  • PARI
    a=1;for(n=2,10^6,if(a%(n-1)==0,print1(n-1","));x=divisors(n);L=numdiv(n);a+=n; for(i=2,L-1,d=x[i];k=n;while(k%d==0,a-=d-1;k\=d))) \\ Robert Gerbicz, May 09 2008

Extensions

a(12)-a(15) from Robert Gerbicz, May 09 2008
a(16)-a(23) from Hieronymus Fischer, Jul 31 2008

A362686 Binomial(n+p, n) mod n where p=6.

Original entry on oeis.org

0, 0, 0, 2, 2, 0, 1, 3, 1, 8, 1, 0, 1, 8, 9, 5, 1, 10, 1, 10, 15, 12, 1, 15, 6, 14, 1, 8, 1, 12, 1, 9, 12, 18, 8, 10, 1, 20, 27, 19, 1, 36, 1, 12, 10, 24, 1, 45, 1, 36, 18, 14, 1, 28, 12, 15, 39, 30, 1, 48, 1, 32, 1, 17, 14, 12, 1, 18, 24, 50, 1, 19, 1, 38
Offset: 1

Views

Author

Ray Chandler, Apr 29 2023

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Mod[Binomial[n+6, n], n], {n, 90}]

Formula

a(n)=binomial(n+6,n) mod n.
For n > 1452, a(n) = 2*a(n-720) - a(n-1440).

A362687 Binomial(n+p, n) mod n where p=7.

Original entry on oeis.org

0, 0, 0, 2, 2, 0, 2, 3, 1, 8, 1, 0, 1, 10, 9, 5, 1, 10, 1, 10, 18, 12, 1, 15, 6, 14, 1, 12, 1, 12, 1, 9, 12, 18, 13, 10, 1, 20, 27, 19, 1, 0, 1, 12, 10, 24, 1, 45, 8, 36, 18, 14, 1, 28, 12, 23, 39, 30, 1, 48, 1, 32, 10, 17, 14, 12, 1, 18, 24, 60, 1, 19, 1
Offset: 1

Views

Author

Ray Chandler, Apr 29 2023

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Mod[Binomial[n+7,n],n],{n,90}]

Formula

a(n)=binomial(n+7,n) mod n.
For n > 10122, a(n) = 2*a(n-5040) - a(n-10080).

A362688 Binomial(n+p, n) mod n where p=8.

Original entry on oeis.org

0, 1, 0, 3, 2, 3, 2, 6, 1, 8, 1, 6, 1, 10, 9, 15, 1, 1, 1, 5, 18, 1, 1, 12, 6, 14, 1, 12, 1, 12, 1, 13, 12, 1, 13, 19, 1, 1, 27, 34, 1, 0, 1, 34, 10, 24, 1, 27, 8, 11, 18, 1, 1, 1, 12, 16, 39, 30, 1, 48, 1, 32, 10, 25, 14, 45, 1, 35, 24, 25, 1, 46, 1, 38, 66
Offset: 1

Views

Author

Ray Chandler, Apr 29 2023

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Mod[Binomial[n+8,n],n],{n,90}]

Formula

a(n)=binomial(n+8,n) mod n.
For n > 645240, a(n) = 2*a(n-322560) - a(n-645120).
Previous Showing 21-30 of 77 results. Next