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-10 of 42 results. Next

A152652 Least prime p with digit sum A047235(n).

Original entry on oeis.org

2, 13, 17, 19, 59, 79, 389, 499, 1889, 1999, 6899, 17989, 39989, 49999, 98999, 199999, 599999, 799999, 2999999, 4999999, 9899999, 19999999, 59999999, 189997999, 389999999, 689899999, 998999999, 2999899999, 6999999989, 9899989999, 39899999999, 68899999999, 98999999999
Offset: 1

Views

Author

Giovanni Teofilatto, Dec 10 2008

Keywords

Crossrefs

Cf. A111380 (smallest prime whose digital sum is equal to the n-th composite number not congruent to 0 (modulo 3)). - Klaus Brockhaus, Dec 12 2008

Programs

  • Magma
    T:=[ n eq 1 select 2 else Self(n-1)+2*(1+n mod 2): n in [1..22] ]; S:=[]; p:=2; for k in T do while &+Intseq(p, 10) ne k do p:=NextPrime(p); end while; Append(~S,p); end for; S; // Klaus Brockhaus, Dec 13 2008
    
  • PARI
    a(n) = {n = (n-1)\2*6+3+(-1)^n ; t = ceil(n/9); leastfound = oo; while(leastfound == oo, my(p = partitions(n, [1,9], [t,t])); v = vector(#p, i, oo); for(i = 1, #p, if(fromdigits(Vec(p[i])) > leastfound, next(2)); forperm(Vec(p[i]), q, if(isprime(fromdigits(Vec(q))), leastfound = min(leastfound, fromdigits(Vec(q))); v[i] = min(v[i], fromdigits(Vec(q))); next(2); ) ) ); t++ ); leastfound }\\ David A. Corneth, Jun 13 2020

Formula

{min A000040(i): A007605(i) = A047234(n)}. - R. J. Mathar, Dec 12 2008

Extensions

Edited and extended by R. J. Mathar, Dec 12 2008
a(20)-a(22) from Klaus Brockhaus, Dec 13 2008
More terms from Jinyuan Wang, Jun 13 2020

A228124 Number of blocks in a Steiner Quadruple System of order A047235(n+1).

Original entry on oeis.org

1, 14, 30, 91, 140, 285, 385, 650, 819, 1240, 1496, 2109, 2470, 3311, 3795, 4900, 5525, 6930, 7714, 9455, 10416, 12529, 13685, 16206, 17575, 20540, 22140, 25585, 27434, 31395, 33511, 38024, 40425, 45526, 48230, 53955, 56980, 63365, 66729, 73810, 77531, 85344
Offset: 1

Views

Author

Colin Barker, Aug 11 2013

Keywords

Comments

For order v, the number of blocks is v*(v-1)*(v-2)/24.

Examples

			For n=3, A047235(n+1)=10 and the number of blocks in SQS(10) is 30.
		

Programs

  • Mathematica
    LinearRecurrence[{1,3,-3,-3,3,1,-1},{1,14,30,91,140,285,385},50] (* Harvey P. Dale, Jul 29 2015 *)
  • PARI
    Vec(x*(x^5+4*x^4+22*x^3+13*x^2+13*x+1)/((x-1)^4*(x+1)^3) + O(x^100))

Formula

a(n) = (n*(1+3*n)*(1+3*(-1)^n+6*n))/16.
a(n) = a(n-1)+3*a(n-2)-3*a(n-3)-3*a(n-4)+3*a(n-5)+a(n-6)-a(n-7).
G.f.: x*(x^5+4*x^4+22*x^3+13*x^2+13*x+1) / ((x-1)^4*(x+1)^3).

A264606 Number of self-dual negacyclic codes of length 2n over GF(9), where 2n runs through the numbers congruent to 2 or 4 mod 6 (cf. A047235).

Original entry on oeis.org

4, 4, 8, 8, 4, 64, 8, 32, 64, 4, 8, 8, 1024, 64, 8, 32, 1024, 64, 8, 8, 4, 64, 512, 32, 64, 1024, 2048, 8, 64, 64, 8, 32, 1024, 1024, 8, 512, 64, 64, 8, 8192, 64, 4, 32768, 128, 1024, 262144, 8, 8192, 1024, 64, 512, 8, 1024, 4194304, 8, 8192, 64, 64, 8, 2147483648, 64, 64, 512, 32, 1024, 262144
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2015

Keywords

Crossrefs

A020760 Decimal expansion of 1/sqrt(3).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

If the sides of a triangle form an arithmetic progression in the ratio 1:1+d:1+2d then when d=1/sqrt(3) it uniquely maximizes the area of the triangle. This triangle has approximate internal angles 25.588 degs, 42.941 degs, 111.471 degs. - Frank M Jackson, Jun 15 2011
When a cylinder is completely enclosed by a sphere, it occupies a fraction f of the sphere volume. The value of f has a trivial lower bound of 0, and an upper bound which is this constant. It is achieved iff the cylinder diameter is sqrt(2) times its height, and the sphere is circumscribed to it. A similar constant can be associated with any n-dimensional geometric shape. For 3D cuboids it is A165952. - Stanislav Sykora, Mar 07 2016
The ratio between the thickness and diameter of a dynamically fair coin having an equal probability, 1/3, of landing on each of its two faces and on its side after being tossed in the air. The calculation is based on the dynamic of rigid body (Yong and Mahadevan, 2011). See A020765 for a simplified geometrical solution. - Amiram Eldar, Sep 01 2020
The coefficient of variation (relative standard deviation) of natural numbers: Limit_{n->oo} sqrt((n-1)/(3*n+3)) = 1/sqrt(3). - Michal Paulovic, Mar 21 2023

Examples

			0.577350269189625764509148780501957455647601751270126876018602326....
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Sections 8.4.3 and 8.17, pp. 495, 531.

Crossrefs

Cf. A002194 (sqrt(3)), A010701 (1/3).

Programs

Formula

Equals 1/A002194. - Michel Marcus, Oct 12 2014
Equals cosine of the magic angle: cos(A195696). - Stanislav Sykora, Mar 07 2016
Equals square root of A010701. - Michel Marcus, Mar 07 2016
Equals 1 + Sum_{k>=0} -(4*k+1)^(-1/2) + (4*k+3)^(-1/2) + (4*k+5)^(-1/2) - (4*k+7)^(-1/2). - Gerry Martens, Nov 22 2022
Equals (1/2)*(2 - zeta(1/2,1/4) + zeta(1/2,3/4) + zeta(1/2,5/4) - zeta(1/2,7/4)). - Gerry Martens, Nov 22 2022
Has periodic continued fraction expansion [0, 1; 1, 2] (A040001). - Michal Paulovic, Mar 21 2023
Equals Product_{k>=1} (1 + (-1)^k/A047235(k)). - Amiram Eldar, Nov 22 2024
Equals tan(Pi/6) = (1/2)/A010527. - R. J. Mathar, Aug 31 2025

A257993 Least gap in the partition having Heinz number n; index of the least prime not dividing n.

Original entry on oeis.org

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

Views

Author

Emeric Deutsch, May 18 2015

Keywords

Comments

The "least gap" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436.
In the Maple program the subprogram B yields the partition with Heinz number n.
Sum of least gaps of all partitions of m = A022567(m).
From Antti Karttunen, Aug 22 2016: (Start)
Index of the least prime not dividing n. (After a formula given by Heinz.)
Least k such that A002110(k) does not divide n.
One more than the number of trailing zeros in primorial base representation of n, A049345.
(End)
The least gap is also called the mex (minimal excludant) of the partition. - Gus Wiseman, Apr 20 2021

Examples

			a(18) = 3 because the partition having Heinz number 18 = 2*3*3 is [1,2,2], having least gap equal to 3.
		

References

  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
  • Miklós Bóna, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Crossrefs

Positions of 1's are A005408.
Positions of 2's are A047235.
The number of gaps is A079067.
The version for crank is A257989.
The triangle counting partitions by this statistic is A264401.
One more than A276084.
The version for greatest difference is A286469 or A286470.
A maximal instead of minimal version is A339662.
Positions of even terms are A342050.
Positions of odd terms are A342051.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339737 counts partitions by sum and greatest gap.

Programs

  • Maple
    with(numtheory): a := proc (n) local B, q: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: for q while member(q, B(n)) = true do  end do: q end proc: seq(a(n), n = 1 .. 150);
    # second Maple program:
    a:= n-> `if`(n=1, 1, (s-> min({$1..(max(s)+1)} minus s))(
            {map(x-> numtheory[pi](x[1]), ifactors(n)[2])[]})):
    seq(a(n), n=1..100);  # Alois P. Heinz, May 09 2016
    # faster:
    A257993 := proc(n) local p, c; c := 1; p := 2;
    while n mod p = 0 do p := nextprime(p); c := c + 1 od: c end:
    seq(A257993(n), n=1..100); # Peter Luschny, Jun 04 2017
  • Mathematica
    A053669[n_] := For[p = 2, True, p = NextPrime[p], If[CoprimeQ[p, n], Return[p]]]; a[n_] := PrimePi[A053669[n]]; Array[a, 100] (* Jean-François Alcover, Nov 28 2016 *)
    Table[k = 1; While[! CoprimeQ[Prime@ k, n], k++]; k, {n, 100}] (* Michael De Vlieger, Jun 22 2017 *)
  • PARI
    a(n) = forprime(p=2,, if (n % p, return(primepi(p)))); \\ Michel Marcus, Jun 22 2017
  • Python
    from sympy import nextprime, primepi
    def a053669(n):
        p = 2
        while True:
            if n%p!=0: return p
            else: p=nextprime(p)
    def a(n): return primepi(a053669(n)) # Indranil Ghosh, May 12 2017
    
  • Scheme
    (define (A257993 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) i (loop (/ (- n d) p) (+ 1 i))))))
    ;; Antti Karttunen, Aug 22 2016
    

Formula

a(n) = A000720(A053669(n)). - Alois P. Heinz, May 18 2015
From Antti Karttunen, Aug 22-30 2016: (Start)
a(n) = 1 + A276084(n).
a(n) = A055396(A276086(n)).
A276152(n) = A002110(a(n)).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} 1/A002110(k) = 1.705230... (1 + A064648). - Amiram Eldar, Jul 23 2022
a(n) << log n/log log n. - Charles R Greathouse IV, Dec 03 2022

Extensions

A simpler description added to the name by Antti Karttunen, Aug 22 2016

A342050 Numbers k which have an odd number of trailing zeros in their primorial base representation A049345(k).

Original entry on oeis.org

2, 4, 8, 10, 14, 16, 20, 22, 26, 28, 30, 32, 34, 38, 40, 44, 46, 50, 52, 56, 58, 60, 62, 64, 68, 70, 74, 76, 80, 82, 86, 88, 90, 92, 94, 98, 100, 104, 106, 110, 112, 116, 118, 120, 122, 124, 128, 130, 134, 136, 140, 142, 146, 148, 150, 152, 154, 158, 160, 164, 166, 170, 172, 176, 178, 180, 182, 184, 188, 190, 194, 196, 200, 202, 206, 208, 212
Offset: 1

Views

Author

Amiram Eldar, Feb 26 2021

Keywords

Comments

Numbers k such that A276084(k) is odd.
All the terms are even since odd numbers have 0 trailing zeros, and 0 is not odd.
The number of terms not exceeding A002110(m) for m>=1 is A002110(m) * Sum_{k=1..m}(-1)^k/A002110(k) = 1, 2, 11, 76, 837, 10880, 184961, ...
The asymptotic density of this sequence is Sum_{k>=1} (-1)^(k+1)/A002110(k) = 0.362306... (A132120).
Also Heinz numbers of partitions with even least gap. The least gap (mex or minimal excludant) of a partition is the least positive integer that is not a part. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions. - Gus Wiseman, Apr 23 2021
Numbers k such that A000720(A053669(k)) is even. Differences from the related A353531 seem to be terms that are multiples of 210, but not all of them, for example primorial 30030 (= 143*210) is in neither sequence. Consider also A038698. - Antti Karttunen, Apr 25 2022

Examples

			2 is a term since A049345(2) = 10 has 1 trailing zero.
4 is a term since A049345(2) = 20 has 1 trailing zero.
30 is a term since A049345(2) = 1000 has 3 trailing zeros.
From _Gus Wiseman_, Apr 23 2021: (Start)
The sequence of terms together with their prime indices begins:
      2: {1}             46: {1,9}             90: {1,2,2,3}
      4: {1,1}           50: {1,3,3}           92: {1,1,9}
      8: {1,1,1}         52: {1,1,6}           94: {1,15}
     10: {1,3}           56: {1,1,1,4}         98: {1,4,4}
     14: {1,4}           58: {1,10}           100: {1,1,3,3}
     16: {1,1,1,1}       60: {1,1,2,3}        104: {1,1,1,6}
     20: {1,1,3}         62: {1,11}           106: {1,16}
     22: {1,5}           64: {1,1,1,1,1,1}    110: {1,3,5}
     26: {1,6}           68: {1,1,7}          112: {1,1,1,1,4}
     28: {1,1,4}         70: {1,3,4}          116: {1,1,10}
     30: {1,2,3}         74: {1,12}           118: {1,17}
     32: {1,1,1,1,1}     76: {1,1,8}          120: {1,1,1,2,3}
     34: {1,7}           80: {1,1,1,1,3}      122: {1,18}
     38: {1,8}           82: {1,13}           124: {1,1,11}
     40: {1,1,1,3}       86: {1,14}           128: {1,1,1,1,1,1,1}
     44: {1,1,5}         88: {1,1,1,5}        130: {1,3,6}
(End)
		

Crossrefs

Complement of A342051.
A099800 is subsequence.
Analogous sequences: A001950 (Zeckendorf representation), A036554 (binary), A145204 (ternary), A217319 (base 4), A232745 (factorial base).
The version for reversed binary expansion is A079523.
Positions of even terms in A257993.
A000070 counts partitions with a selected part.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A079067 counts gaps in prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339662 gives greatest gap in prime indices.
Differs from A353531 for the first time at n=77, where a(77) = 212, as this sequence misses A353531(77) = 210.

Programs

  • Mathematica
    seq[max_] := Module[{bases = Prime@Range[max, 1, -1], nmax}, nmax = Times @@ bases - 1; Select[Range[nmax], OddQ @ LengthWhile[Reverse @ IntegerDigits[#, MixedRadix[bases]], #1 == 0 &] &]]; seq[4]
    Select[Range[100],EvenQ[Min@@Complement[Range[PrimeNu[#]+1],PrimePi/@First/@FactorInteger[#]]]&] (* Gus Wiseman, Apr 23 2021 *)
  • PARI
    A353525(n) = { for(i=1,oo,if(n%prime(i),return((i+1)%2))); }
    isA342050(n) = A353525(n);
    k=0; n=0; while(k<77, n++; if(isA342050(n), k++; print1(n,", "))); \\ Antti Karttunen, Apr 25 2022

Extensions

More terms added (to differentiate from A353531) by Antti Karttunen, Apr 25 2022

A047270 Numbers that are congruent to {3, 5} mod 6.

Original entry on oeis.org

3, 5, 9, 11, 15, 17, 21, 23, 27, 29, 33, 35, 39, 41, 45, 47, 51, 53, 57, 59, 63, 65, 69, 71, 75, 77, 81, 83, 87, 89, 93, 95, 99, 101, 105, 107, 111, 113, 117, 119, 123, 125, 129, 131, 135, 137, 141, 143, 147, 149
Offset: 1

Views

Author

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 10 ).
This sequence is an interleaving of A016945 with A016969. - Guenther Schrack, Nov 16 2018

Crossrefs

Cf. A047235 [(6*n-(-1)^n-3)/2], A047241 [(6*n-(-1)^n-5)/2], A047238 [(6*n-(-1)^n-7)/2]. [Bruno Berselli, Jun 24 2010]
Subsequence of A186422.
From Guenther Schrack, Nov 18 2018: (Start)
Complement: A047237.
First differences: A105397(n) for n > 0.
Partial sums: A227017(n+1) for n > 0.
Elements of odd index: A016945.
Elements of even index: A016969(n-1) for n > 0. (End)

Programs

  • Mathematica
    Select[Range@ 149, MemberQ[{3, 5}, Mod[#, 6]] &] (* or *)
    Array[(6 # - (-1)^# - 1)/2 &, 50] (* or *)
    Fold[Append[#1, 6 #2 - Last@ #1 - 4] &, {3}, Range[2, 50]] (* or *)
    CoefficientList[Series[(3 + 2 x + x^2)/((1 + x) (1 - x)^2), {x, 0, 49}], x] (* Michael De Vlieger, Jan 12 2018 *)
  • PARI
    a(n) = (6*n - 1 - (-1)^n)/2 \\ David Lovler, Aug 25 2022

Formula

a(n) = sqrt(2)*sqrt((1-6*n)*(-1)^n + 18*n^2 - 6*n + 1)/2. - Paul Barry, May 11 2003
From Bruno Berselli, Jun 24 2010: (Start)
G.f.: (3+2*x+x^2)/((1+x)*(1-x)^2).
a(n) - a(n-1) - a(n-2) + a(n-3) = 0, with n > 3.
a(n) = (6*n - (-1)^n - 1)/2. (End)
a(n) = 6*n - a(n-1) - 4 with n > 1, a(1)=3. - Vincenzo Librandi, Aug 05 2010
From Guenther Schrack, Nov 17 2018: (Start)
a(n) = a(n-2) + 6 for n > 2.
a(-n) = -A047241(n+1) for n > 0.
a(n) = A109613(n-1) + 2*n for n > 0.
a(n) = 2*A001651(n) + 1.
m-element moving averages: Sum_{k=1..m} a(n-m+k)/m = A016777(n-m/2) for m = 2, 4, 6, ... and n >= m. (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(4*sqrt(3)) - log(3)/4. - Amiram Eldar, Dec 13 2021
E.g.f.: 1 + 3*x*exp(x) - cosh(x). - David Lovler, Aug 25 2022

A342051 Numbers k which have an even number of trailing zeros in their primorial base representation A049345(k).

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 11, 12, 13, 15, 17, 18, 19, 21, 23, 24, 25, 27, 29, 31, 33, 35, 36, 37, 39, 41, 42, 43, 45, 47, 48, 49, 51, 53, 54, 55, 57, 59, 61, 63, 65, 66, 67, 69, 71, 72, 73, 75, 77, 78, 79, 81, 83, 84, 85, 87, 89, 91, 93, 95, 96, 97, 99, 101, 102, 103
Offset: 1

Views

Author

Amiram Eldar, Feb 26 2021

Keywords

Comments

Numbers k such that A276084(k) is even.
The number of terms not exceeding A002110(m) for m>=1 is A002110(m) * (1 - Sum_{k=1..m}(-1)^k/A002110(k)) = 1, 4, 19, 134, 1473, 19150, 325549 ...
The asymptotic density of this sequence is Sum_{k>=0} (-1)^k/A002110(k) = 0.637693... = 1 - A132120.
Also Heinz numbers of partitions with odd least gap. The least gap (mex or minimal excludant) of a partition is the least positive integer that is not a part. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions. - Gus Wiseman, Apr 23 2021

Examples

			1 is a term since A049345(1) = 1 has 0 trailing zero.
6 is a term since A049345(6) = 100 has 2 trailing zeros.
From _Gus Wiseman_, Apr 23 2021: (Start)
The sequence of terms together with their prime indices begins:
     1: {}           25: {3,3}          51: {2,7}
     3: {2}          27: {2,2,2}        53: {16}
     5: {3}          29: {10}           54: {1,2,2,2}
     6: {1,2}        31: {11}           55: {3,5}
     7: {4}          33: {2,5}          57: {2,8}
     9: {2,2}        35: {3,4}          59: {17}
    11: {5}          36: {1,1,2,2}      61: {18}
    12: {1,1,2}      37: {12}           63: {2,2,4}
    13: {6}          39: {2,6}          65: {3,6}
    15: {2,3}        41: {13}           66: {1,2,5}
    17: {7}          42: {1,2,4}        67: {19}
    18: {1,2,2}      43: {14}           69: {2,9}
    19: {8}          45: {2,2,3}        71: {20}
    21: {2,4}        47: {15}           72: {1,1,1,2,2}
    23: {9}          48: {1,1,1,1,2}    73: {21}
    24: {1,1,1,2}    49: {4,4}          75: {2,3,3}
(End)
		

Crossrefs

Complement of A342050.
A099788 is subsequence.
Analogous sequences: A000201 (Zeckendorf representation), A003159 (binary), A007417 (ternary), A232744 (factorial base).
The version for reversed binary expansion is A121539.
Positions of odd terms in A257993.
A000070 counts partitions with a selected part.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A079067 counts gaps in prime indices.
A238709 counts partitions by sum and least difference.
A339662 gives greatest gap in prime indices.

Programs

  • Mathematica
    seq[max_] := Module[{bases = Prime@Range[max, 1, -1], nmax}, nmax = Times @@ bases - 1; Select[Range[nmax], EvenQ @ LengthWhile[Reverse @ IntegerDigits[#, MixedRadix[bases]], #1 == 0 &] &]]; seq[4]
    Select[Range[100],OddQ[Min@@Complement[Range[PrimeNu[#]+1],PrimePi/@First/@FactorInteger[#]]]&] (* Gus Wiseman, Apr 23 2021 *)

A046729 Expansion of 4*x/((1+x)*(1-6*x+x^2)).

Original entry on oeis.org

0, 4, 20, 120, 696, 4060, 23660, 137904, 803760, 4684660, 27304196, 159140520, 927538920, 5406093004, 31509019100, 183648021600, 1070379110496, 6238626641380, 36361380737780, 211929657785304, 1235216565974040, 7199369738058940, 41961001862379596, 244566641436218640
Offset: 0

Views

Author

Keywords

Comments

Related to Pythagorean triples: alternate terms of A001652 and A046090.
Even-valued legs of nearly isosceles right triangles: legs differ by 1. 0 is smaller leg of degenerate triangle with legs 0 and 1 and hypotenuse 1. - Charlie Marion, Nov 11 2003
The complete (nearly isosceles) primitive Pythagorean triple is given by {a(n), a(n)+(-1)^n, A001653(n)}. - Lekraj Beedassy, Feb 19 2004
Note also that A046092 is the even leg of this other class of nearly isosceles Pythagorean triangles {A005408(n), A046092(n), A001844(n)}, i.e., {2n+1, 2n(n+1), 2n(n+1)+1} where longer sides (viz. even leg and hypotenuse) are consecutive. - Lekraj Beedassy, Apr 22 2004
Union of even terms of A001652 and A046090. Sum of legs of primitive Pythagorean triangles is A002315(n) = 2*a(n) + (-1)^n. - Lekraj Beedassy, Apr 30 2004

Examples

			[1,0,1]*[1,2,2; 2,1,2; 2,2,3]^0 gives (degenerate) primitive Pythagorean triple [1, 0, 1], so a(0) = 0. [1,0,1]*[1,2,2; 2,1,2; 2,2,3]^7 gives primitive Pythagorean triple [137903, 137904, 195025] so a(7) = 137904.
G.f. = 4*x + 20*x^2 + 120*x^3 + 696*x^4 + 4060*x^5 + 23660*x^6 + ...
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, pp. 122-125, 1964.
  • W. Sierpiński, Pythagorean triangles, Dover Publications, Inc., Mineola, NY, 2003, p. 17. MR2002669.

Crossrefs

Programs

  • Magma
    [4*Floor(((Sqrt(2)+1)^(2*n+1)-(Sqrt(2)-1)^(2*n+1)-2*(-1)^n) / 16): n in [0..30]]; // Vincenzo Librandi, Jul 29 2019
    
  • Mathematica
    LinearRecurrence[{5,5,-1}, {0,4,20}, 25] (* Vincenzo Librandi, Jul 29 2019 *)
  • PARI
    a(n)=n%2+(real((1+quadgen(8))^(2*n+1))-1)/2
    
  • PARI
    a(n)=if(n<0,-a(-1-n),polcoeff(4*x/(1+x)/(1-6*x+x^2)+x*O(x^n),n))
    
  • SageMath
    [(lucas_number2(2*n+1,2,-1) -2*(-1)^n)/4 for n in range(41)] # G. C. Greubel, Feb 11 2023

Formula

a(n) = ((1+sqrt(2))^(2n+1) + (1-sqrt(2))^(2n+1) + 2*(-1)^(n+1))/4.
a(n) = A089499(n)*A089499(n+1).
a(n) = 4*A084158(n). - Lekraj Beedassy, Jul 16 2004
a(n) = ceiling((sqrt(2)+1)^(2*n+1) - (sqrt(2)-1)^(2*n+1) - 2*(-1)^n)/4. - Lambert Klasen (Lambert.Klasen(AT)gmx.net), Nov 12 2004
a(n) is the k-th entry among the complete near-isosceles primitive Pythagorean triple A114336(n), where k = (3*(2n-1) - (-1)^n)/2, i.e., a(n) = A114336(A047235(n)), for positive n. - Lekraj Beedassy, Jun 04 2006
a(n) = A046727(n) - (-1)^n = 2*A114620(n). - Lekraj Beedassy, Aug 14 2006
From George F. Johnson, Aug 29 2012: (Start)
2*a(n)*(a(n) + (-1)^n) + 1 = (A000129(2*n+1))^2;
n > 0, 2*a(n)*(a(n) + (-1)^n) + 1 = ((a(n+1) - a(n-1))/4)^2, a perfect square.
a(n+1) = (3*a(n) + 2*(-1)^n) + 2*sqrt(2*a(n)*(a(n) + (-1)^n)+ 1).
a(n-1) = (3*a(n) + 2*(-1)^n) - 2*sqrt(2*a(n)*(a(n) + (-1)^n)+ 1).
a(n+1) = 6*a(n) - a(n-1) + 4*(-1)^n.
a(n+1) = 5*a(n) + 5*a(n-1) - a(n-2).
a(n+1) *a(n-1) = a(n)*(a(n) + 4*(-1)^n).
a(n) = (sqrt(1 + 8*A029549(n)) - (-1)^n)/2.
a(n) = A002315(n) - A084159(n) = A084159(n) - (-1)^n.
a(n) = A001652(n) + (1 - (-1)^n)/2 = A046090(n) - (1 + (-1)^n)/2.
Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2).
Limit_{n->oo} a(n)/a(n-2) = 17 + 12*sqrt(2).
Limit_{n->oo} a(n)/a(n-r) = (3 + 2*sqrt(2))^r.
Limit_{n->oo} a(n-r)/a(n) = (3 - 2*sqrt(2))^r. (End)
From G. C. Greubel, Feb 11 2023: (Start)
a(n) = (A001333(2*n+1) - 2*(-1)^n)/4.
a(n) = (1/2)*(A001109(n+1) + A001109(n) - (-1)^n). (End)
E.g.f.: exp(-x)*(exp(4*x)*(cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)) - 1)/2. - Stefano Spezia, Aug 03 2024

A047238 Numbers that are congruent to {0, 2} mod 6.

Original entry on oeis.org

0, 2, 6, 8, 12, 14, 18, 20, 24, 26, 30, 32, 36, 38, 42, 44, 48, 50, 54, 56, 60, 62, 66, 68, 72, 74, 78, 80, 84, 86, 90, 92, 96, 98, 102, 104, 108, 110, 114, 116, 120, 122, 126, 128, 132, 134, 138, 140, 144, 146, 150, 152, 156, 158, 162
Offset: 1

Views

Author

Keywords

Comments

Complement of A047251, or "Polyrhythmic Sequence" P(2,3); the present sequence represents where the "rests" occur in a "3 against 2" polyrhythm. (See A267027 for definition and description). - Bob Selcoe, Jan 12 2016

Crossrefs

Cf. A047270 [(6*n-(-1)^n-1)/2], A047235 [(6*n-(-1)^n-3)/2], A047241 [(6*n-(-1)^n-5)/2].

Programs

  • Magma
    [n: n in [0..200]|n mod 6 in {0,2}]; // Vincenzo Librandi, Jan 12 2016
  • Mathematica
    Select[Range[0,200],MemberQ[{0,2},Mod[#,6]]&] (* or *) LinearRecurrence[ {1,1,-1},{0,2,6},70] (* Harvey P. Dale, Jun 15 2011 *)
  • PARI
    forstep(n=0,200,[2,4],print1(n", ")) \\ Charles R Greathouse IV, Oct 17 2011
    

Formula

From Bruno Berselli, Jun 24 2010: (Start)
G.f.: 2*x*(1+2*x)/((1+x)*(1-x)^2).
a(n) = a(n-1) + a(n-2) - a(n-3), a(0)=0, a(1)=2, a(2)=6.
a(n) = (6*n - (-1)^n-7)/2.
a(n) = 2*A032766(n-1). (End)
a(n) = 6*n - a(n-1) - 10 (with a(1)=0). - Vincenzo Librandi, Aug 05 2010
a(n+1) = Sum_{k>=0} A030308(n,k)*A111286(k+2). - Philippe Deléham, Oct 17 2011
a(n) = 2*floor(3*n/2). - Enrique Pérez Herrero, Jul 04 2012
Sum_{n>=2} (-1)^n/a(n) = sqrt(3)*Pi/36 + log(3)/4. - Amiram Eldar, Dec 13 2021
E.g.f: 3*(x-1)*exp(x) - cosh(x) + 4. - David Lovler, Jul 11 2022
Showing 1-10 of 42 results. Next