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 12 results. Next

A001651 Numbers not divisible by 3.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 49, 50, 52, 53, 55, 56, 58, 59, 61, 62, 64, 65, 67, 68, 70, 71, 73, 74, 76, 77, 79, 80, 82, 83, 85, 86, 88, 89, 91, 92, 94, 95, 97, 98, 100, 101, 103, 104
Offset: 1

Views

Author

Keywords

Comments

Inverse binomial transform of A084858. - Benoit Cloitre, Jun 12 2003
Earliest monotonic sequence starting with (1,2) and satisfying the condition: "a(n)+a(n-1) is not in the sequence." - Benoit Cloitre, Mar 25 2004. [The numbers of the form a(n)+a(n-1) form precisely the complement with respect to the positive integers. - David W. Wilson, Feb 18 2012]
a(1) = 1; a(n) is least number which is relatively prime to the sum of all the previous terms. - Amarnath Murthy, Jun 18 2001
For n > 3, numbers having 3 as an anti-divisor. - Alexandre Wajnberg, Oct 02 2005
Also numbers n such that (n+1)*(n+2)/6 = A000292(n)/n is an integer. - Ctibor O. Zizka, Oct 15 2010
Notice the property described by Gary Detlefs in A113801: more generally, these numbers are of the form (2*h*n + (h-4)*(-1)^n-h)/4 (h, n natural numbers), therefore ((2*h*n + (h-4)*(-1)^n - h)/4)^2 - 1 == 0 (mod h); in this case, a(n)^2 - 1 == 0 (mod 3). - Bruno Berselli, Nov 17 2010
A001651 mod 9 gives A141425. - Paul Curtz, Dec 31 2010. (Correct for the modified offset 1. - M. F. Hasler, Apr 07 2015)
The set of natural numbers (1, 2, 3, ...), sequence A000027; represents the numbers of ordered compositions of n using terms in the signed set: (1, 2, -4, -5, 7, 8, -10, -11, 13, 14, ...). This follows from (1, 2, 3, ...) being the INVERT transform of A011655, signed and beginning: (1, 1, 0, -1, -1, 0, 1, 1, 0, ...). - Gary W. Adamson, Apr 28 2013
Union of A047239 and A047257. - Wesley Ivan Hurt, Dec 19 2013
Numbers whose sum of digits (and digital root) is != 0 (mod 3). - Joerg Arndt, Aug 29 2014
The number of partitions of 3*(n-1) into at most 2 parts. - Colin Barker, Apr 22 2015
a(n) is the number of partitions of 3*n into two distinct parts. - L. Edson Jeffery, Jan 14 2017
Conjectured (and like even easily proved) to be the graph bandwidth of the complete bipartite graph K_{n,n}. - Eric W. Weisstein, Apr 24 2017
Numbers k such that Fibonacci(k) mod 4 = 1 or 3. Equivalently, sequence lists the indices of the odd Fibonacci numbers (see A014437). - Bruno Berselli, Oct 17 2017
Minimum value of n_3 such that the "rectangular spiral pattern" is the optimal solution for Ripà's n_1 X n_2 x n_3 Dots Problem, for any n_1 = n_2. For example, if n_1 = n_2 = 5, n_3 = floor((3/2)*(n_1 - 1)) + 1 = a(5). - Marco Ripà, Jul 23 2018
For n >= 54, a(n) = sat(n, P_n), the minimum number of edges in a P_n-saturated graph on n vertices, where P_n is the n-vertex path (see Dudek, Katona, and Wojda, 2003; Frick and Singleton, 2005). - Danny Rorabaugh, Nov 07 2017
From Roger Ford, May 09 2021: (Start)
a(n) is the smallest sum of arch lengths for the top arches of a semi-meander with n arches. An arch length is the number of arches covered + 1.
/\ The top arch has a length of 3. /\ The top arch has a length of 3.
/ \ Both bottom arches have a //\\ The middle arch has a length of 2.
//\/\\ length of 1. ///\\\ The bottom arch has a length of 1.
Example: a(6) = 8 /\ /\
//\\ /\ //\\ /\ 2 + 1 + 1 + 2 + 1 + 1 = 8. (End)
This is the lexicographically earliest increasing sequence of positive integers such that no polynomial of degree d can be fitted to d+2 consecutive terms (equivalently, such that no iterated difference is zero). - Pontus von Brömssen, Dec 26 2021

Examples

			G.f.: x + 2*x^2 + 4*x^3 + 5*x^4 + 7*x^5 + 8*x^6 + 10*x^7 + 11*x^8 + 13*x^9 + ...
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • GAP
    Filtered([0..110],n->n mod 3<>0); # Muniru A Asiru, Jul 24 2018
    
  • Haskell
    a001651 = (`div` 2) . (subtract 1) . (* 3)
    a001651_list = filter ((/= 0) . (`mod` 3)) [1..]
    -- Reinhard Zumkeller, Jul 07 2012, Aug 23 2011
    
  • Magma
    [3*(2*n-1)/4-(-1)^n/4: n in [1..80]]; // Vincenzo Librandi, Jun 07 2011
    
  • Maple
    A001651 := n -> 3*floor(n/2) - (-1)^n; # Corrected by M. F. Hasler, Apr 07 2015
    A001651:=(1+z+z**2)/(z+1)/(z-1)**2; # Simon Plouffe in his 1992 dissertation
    a[1]:=1:a[2]:=2:for n from 3 to 100 do a[n]:=a[n-2]+3 od: seq(a[n], n=1..69); # Zerinvary Lajos, Mar 16 2008, offset corrected by M. F. Hasler, Apr 07 2015
  • Mathematica
    Select[Table[n,{n,200}],Mod[#,3]!=0&] (* Vladimir Joseph Stephan Orlovsky, Feb 18 2011 *)
    Drop[Range[200 + 1], {1, -1, 3}] - 1 (* József Konczer, May 24 2016 *)
    Floor[(3 Range[70] - 1)/2] (* Eric W. Weisstein, Apr 24 2017 *)
    CoefficientList[Series[(x^2 + x + 1)/((x - 1)^2 (x + 1)), {x, 0, 70}],
      x] (* or *)
    LinearRecurrence[{1, 1, -1}, {1, 2, 4}, 70] (* Robert G. Wilson v, Jul 25 2018 *)
  • PARI
    {a(n) = n + (n-1)\2}; /* Michael Somos, Jan 15 2011 */
    
  • PARI
    x='x+O('x^100); Vec(x*(1+x+x^2)/((1-x)*(1-x^2))) \\ Altug Alkan, Oct 22 2015
    
  • Python
    print([k for k in range(1, 105) if k%3]) # Michael S. Branicky, Sep 06 2021
    
  • Python
    def A001651(n): return (n<<1)-(n>>1)-1 # Chai Wah Wu, Mar 05 2024

Formula

a(n) = 3 + a(n-2) for n > 2.
a(n) = a(n-1) + a(n-2) - a(n-3) for n > 3.
a(2*n+1) = 3*n+1, a(2*n) = 3*n-1.
G.f.: x * (1 + x + x^2) / ((1 - x) * (1 - x^2)). - Michael Somos, Jun 08 2000
a(n) = (4-n)*a(n-1) + 2*a(n-2) + (n-3)*a(n-3) (from the Carlitz et al. article).
a(n) = floor((3*n-1)/2). [Corrected by Gary Detlefs]
a(1) = 1, a(n) = 2*a(n-1) - 3*floor(a(n-1)/3). - Benoit Cloitre, Aug 17 2002
a(n+1) = 1 + n - n mod 2 + (n + n mod 2)/2. - Reinhard Zumkeller, Dec 17 2002
a(1) = 1, a(n+1) = a(n) + (a(n) mod 3). - Reinhard Zumkeller, Mar 23 2003
a(1) = 1, a(n) = 3*(n-1) - a(n-1). - Benoit Cloitre, Apr 12 2003
a(n) = 3*(2*n-1)/4 - (-1)^n/4. - Benoit Cloitre, Jun 12 2003
Nearest integer to (Sum_{k>=n} 1/k^3)/(Sum_{k>=n} 1/k^4). - Benoit Cloitre, Jun 12 2003
Partial sums of A040001. a(n) = A032766(n-1)+1. - Paul Barry, Sep 02 2003
a(n) = T(n, 1) = T(n, n-1), where T is the array in A026386. - Emeric Deutsch, Feb 18 2004
a(n) = sqrt(3*A001082(n)+1). - Zak Seidov, Dec 12 2007
a(n) = A077043(n) - A077043(n-1). - Reinhard Zumkeller, Dec 28 2007
a(n) = A001477(n-1) + A008619(n-1). - Yosu Yurramendi, Aug 10 2008
Euler transform of length 3 sequence [2, 1, -1]. - Michael Somos, Sep 06 2008
A011655(a(n)) = 1. - Reinhard Zumkeller, Nov 30 2009
a(n) = n - 1 + ceiling(n/2). - Michael Somos, Jan 15 2011
a(n) = 3*A000217(n)+1 - 2*Sum_{i=1..n-1} a(i), for n>1. - Bruno Berselli, Nov 17 2010
a(n) = 3*floor(n/2) + (-1)^(n+1). - Gary Detlefs, Dec 29 2011
A215879(a(n)) > 0. - Reinhard Zumkeller, Dec 28 2012 [More precisely, A215879 is the characteristic function of A001651. - M. F. Hasler, Apr 07 2015]
a(n) = 2n - 1 - floor(n/2). - Wesley Ivan Hurt, Oct 25 2013
a(n) = (3n - 2 + (n mod 2)) / 2. - Wesley Ivan Hurt, Mar 31 2014
a(n) = A000217(n) - A000982(n-1). - Bui Quang Tuan, Mar 28 2015
1/1^3 - 1/2^3 + 1/4^3 - 1/5^3 + 1/7^3 - 1/8^3 + ... = 4 Pi^3/(3 sqrt(3)). - M. F. Hasler, Mar 29 2015
E.g.f.: (4 + sinh(x) - cosh(x) + 3*(2*x - 1)*exp(x))/4. - Ilya Gutkovskiy, May 24 2016
a(n) = a(n+k-1) + a(n-k) - a(n-1) for n > k >= 0. - Bob Selcoe, Feb 03 2017
a(n) = -a(1-n) for all n in Z. - Michael Somos, Jul 31 2018
a(n) = n + A004526(n-1). - David James Sycamore, Sep 06 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(3)) (A073010). - Amiram Eldar, Dec 04 2021
From Amiram Eldar, Nov 22 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 1.
Product_{n>=2} (1 + (-1)^n/a(n)) = 2*Pi/(3*sqrt(3)) (A248897). (End)

Extensions

This is a list, so the offset should be 1. I corrected this and adjusted some of the comments and formulas. Other lines probably also need to be adjusted. - N. J. A. Sloane, Jan 01 2011
Offset of pre-2011 formulas verified or corrected by M. F. Hasler, Apr 07-18 2015 and by Danny Rorabaugh, Oct 23 2015

A000325 a(n) = 2^n - n.

Original entry on oeis.org

1, 1, 2, 5, 12, 27, 58, 121, 248, 503, 1014, 2037, 4084, 8179, 16370, 32753, 65520, 131055, 262126, 524269, 1048556, 2097131, 4194282, 8388585, 16777192, 33554407, 67108838, 134217701, 268435428, 536870883, 1073741794, 2147483617
Offset: 0

Views

Author

Rosario Salamone (Rosario.Salamone(AT)risc.uni-linz.ac.at)

Keywords

Comments

Number of permutations of degree n with at most one fall; called "Grassmannian permutations" by Lascoux and Schützenberger. - Axel Kohnert (Axel.Kohnert(AT)uni-bayreuth.de)
Number of different permutations of a deck of n cards that can be produced by a single shuffle. [DeSario]
Number of Dyck paths of semilength n having at most one long ascent (i.e., ascent of length at least two). Example: a(4)=12 because among the 14 Dyck paths of semilength 4, the only paths that have more than one long ascent are UUDDUUDD and UUDUUDDD (each with two long ascents). Here U = (1, 1) and D = (1, -1). Also number of ordered trees with n edges having at most one branch node (i.e., vertex of outdegree at least two). - Emeric Deutsch, Feb 22 2004
Number of {12,1*2*,21*}-avoiding signed permutations in the hyperoctahedral group.
Number of 1342-avoiding circular permutations on [n+1].
2^n - n is the number of ways to partition {1, 2, ..., n} into arithmetic progressions, where in each partition all the progressions have the same common difference and have lengths at least 1. - Marty Getz (ffmpg1(AT)uaf.edu) and Dixon Jones (fndjj(AT)uaf.edu), May 21 2005
if b(0) = x and b(n) = b(n-1) + b(n-1)^2*x^(n-2) for n > 0, then b(n) is a polynomial of degree a(n). - Michael Somos, Nov 04 2006
The chromatic invariant of the Mobius ladder graph M_n for n >= 2. - Jonathan Vos Post, Aug 29 2008
Dimension sequence of the dual alternative operad (i.e., associative and satisfying the identity xyz + yxz + zxy + xzy + yzx + zyx = 0) over the field of characteristic 3. - Pasha Zusmanovich, Jun 09 2009
An elephant sequence, see A175654. For the corner squares six A[5] vectors, with decimal values between 26 and 176, lead to this sequence (without the first leading 1). For the central square these vectors lead to the companion sequence A168604. - Johannes W. Meijer, Aug 15 2010
a(n+1) is also the number of order-preserving and order-decreasing partial isometries (of an n-chain). - Abdullahi Umar, Jan 13 2011
A040001(n) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, May 12 2012
A130103(n+1) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, May 12 2012
The number of labeled graphs with n vertices whose vertex set can be partitioned into a clique and a set of isolated points. - Alex J. Best, Nov 20 2012
For n > 0, a(n) is a B_2 sequence. - Thomas Ordowski, Sep 23 2014
See coefficients of the linear terms of the polynomials of the table on p. 10 of the Getzler link. - Tom Copeland, Mar 24 2016
Consider n points lying on a circle, then for n>=2 a(n-1) is the maximum number of ways to connect two points with non-intersecting chords. - Anton Zakharov, Dec 31 2016
Also the number of cliques in the (n-1)-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017
From Eric M. Schmidt, Jul 17 2017: (Start)
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) != e(j) < e(k). [Martinez and Savage, 2.7]
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i), e(j), e(k) pairwise distinct. [Martinez and Savage, 2.7]
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(j) >= e(k) and e(i) != e(k) pairwise distinct. [Martinez and Savage, 2.7]
(End)
Number of F-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are F-equivalent iff the positions of pattern F are identical in these paths. - Sergey Kirgizov, Apr 08 2018
From Gus Wiseman, Feb 10 2019: (Start)
Also the number of connected partitions of an n-cycle. For example, the a(1) = 1 through a(4) = 12 connected partitions are:
{{1}} {{12}} {{123}} {{1234}}
{{1}{2}} {{1}{23}} {{1}{234}}
{{12}{3}} {{12}{34}}
{{13}{2}} {{123}{4}}
{{1}{2}{3}} {{124}{3}}
{{134}{2}}
{{14}{23}}
{{1}{2}{34}}
{{1}{23}{4}}
{{12}{3}{4}}
{{14}{2}{3}}
{{1}{2}{3}{4}}
(End)
Number of subsets of n-set without the single-element subsets. - Yuchun Ji, Jul 16 2019
For every prime p, there are infinitely many terms of this sequence that are divisible by p (see IMO Compendium link and Doob reference). Corresponding indices n are: for p = 2, even numbers A299174; for p = 3, A047257; for p = 5, A349767. - Bernard Schott, Dec 10 2021
Primes are in A081296 and corresponding indices in A048744. - Bernard Schott, Dec 12 2021

Examples

			G.f. = 1 + x + 2*x^2 + 5*x^3 + 12*x^4 + 27*x^5 + 58*x^6 + 121*x^7 + ...
		

References

  • Michael Doob, The Canadian Mathematical Olympiad & L'Olympiade Mathématique du Canada 1969-1993, Canadian Mathematical Society & Société Mathématique du Canada, Problem 4, 1983, page 158, 1993.

Crossrefs

Column 1 of triangle A008518.
Row sum of triangles A184049 and A184050.

Programs

  • Haskell
    a000325 n = 2 ^ n - n
    a000325_list = zipWith (-) a000079_list [0..]
    -- Reinhard Zumkeller, Jul 17 2012
    
  • Magma
    [2^n - n: n in [0..35]]; // Vincenzo Librandi, May 13 2011
    
  • Maple
    A000325 := proc(n) option remember; if n <=1 then n+1 else 2*A000325(n-1)+n-1; fi; end;
    g:=1/(1-2*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)-n, n=0..31); # Zerinvary Lajos, Jan 09 2009
  • Mathematica
    Table[2^n - n, {n, 0, 39}] (* Alonso del Arte, Sep 15 2014 *)
    LinearRecurrence[{4, -5, 2}, {1, 2, 5}, {0, 20}] (* Eric W. Weisstein, Jul 14 2017 *)
  • PARI
    {a(n) = 2^n - n}; /* Michael Somos, Nov 04 2006 */
    
  • Python
    def A000325(n): return (1<Chai Wah Wu, Jan 11 2023

Formula

a(n+1) = 2*a(n) + n - 1, a(0) = 1. - Reinhard Zumkeller, Apr 12 2003
Binomial transform of 1, 0, 1, 1, 1, .... The sequence starting 1, 2, 5, ... has a(n) = 1 + n + 2*Sum_{k=2..n} binomial(n, k) = 2^(n+1) - n - 1. This is the binomial transform of 1, 1, 2, 2, 2, 2, .... a(n) = 1 + Sum_{k=2..n} C(n, k). - Paul Barry, Jun 06 2003
G.f.: (1-3x+3x^2)/((1-2x)*(1-x)^2). - Emeric Deutsch, Feb 22 2004
A107907(a(n+2)) = A000051(n+2) for n > 0. - Reinhard Zumkeller, May 28 2005
a(n+1) = sum of n-th row of the triangle in A109128. - Reinhard Zumkeller, Jun 20 2005
Row sums of triangle A133116. - Gary W. Adamson, Sep 14 2007
G.f.: 1 / (1 - x / (1 - x / ( 1 - x / (1 + x / (1 - 2*x))))). - Michael Somos, May 12 2012
First difference is A000225. PSUM transform is A084634. - Michael Somos, May 12 2012
a(n) = [x^n](B(x)^n-B(x)^(n-1)), n>0, a(0)=1, where B(x) = (1+2*x+sqrt(1+4*x^2))/2. - Vladimir Kruchinin, Mar 07 2014
E.g.f.: (exp(x) - x)*exp(x). - Ilya Gutkovskiy, Aug 07 2016
a(n) = A125128(n) - A000225(n) + 1. - Miquel Cerda, Aug 12 2016
a(n) = 2*A125128(n) - A095151(n) + 1. - Miquel Cerda, Aug 12 2016
a(n) = A079583(n-1) - A000225(n-1). - Miquel Cerda, Aug 15 2016
a(n)^2 - 4*a(n-1)^2 = (n-2)*(a(n)+2*a(n-1)). - Yuchun Ji, Jul 13 2018
a(n) = 2^(-n) * A186947(n) = 2^n * A002064(-n) for all n in Z. - Michael Somos, Jul 18 2018
a(2^n) = (2^a(n) - 1)*2^n. - Lorenzo Sauras Altuzarra, Feb 01 2022

A327859 a(n) = A276086(A003415(n)), where A003415 is the arithmetic derivative, and A276086 is the primorial base exp-function.

Original entry on oeis.org

1, 1, 2, 2, 9, 2, 18, 2, 25, 5, 10, 2, 225, 2, 30, 15, 21, 2, 750, 2, 625, 45, 50, 2, 525, 45, 150, 3750, 21, 2, 14, 2, 18375, 75, 250, 25, 49, 2, 750, 225, 735, 2, 630, 2, 875, 210, 1250, 2, 385875, 75, 1050, 375, 13125, 2, 36750, 225, 1029, 1125, 14, 2, 1029, 2, 42, 5250, 2941225, 125, 98, 2, 1225, 1875, 78750
Offset: 0

Views

Author

Antti Karttunen, Sep 30 2019

Keywords

Comments

Sequence contains only terms of A048103.
Are there fixed points other than 1, 2, 10, 15, 5005? (There are none in the range 5006 .. 402653184.) See A369650.
Records occur at n = 0, 2, 4, 6, 8, 12, 18, 27, 32, 48, 64, 80, 144, 224, 256, 336, 448, 480, 512, 1728, ... (see also A131117).
a(n) and n are never multiples of 9 at the same time, thus the fixed points certainly exclude any terms of A008591. For a proof, consider my comment in A047257 and that A003415(9*n) is always a multiple of 3. - Antti Karttunen, Feb 08 2024

Crossrefs

Cf. A003415, A008591, A048103, A131117, A276086, A327858, A327860, A341517 [= mu(a(n))], A341518 (k where a(k) is squarefree), A369641 (composite k where a(k) is squarefree), A369642.
Cf. A370114 (where a(k) is a multiple of k), A370115 (where k is a multiple of a(k)), A369650.

Programs

  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
    A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
    A327859(n) = A276086(A003415(n));

Formula

a(n) = A276086(A003415(n)).
a(p) = 2 for all primes p.

A299174 The positive even integers.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, 124, 126, 128, 130, 132, 134, 136, 138, 140, 142, 144
Offset: 1

Views

Author

Joss Langford, Feb 04 2018

Keywords

Comments

Possible periods of Post's {00, 1101} tag system. - Charles R Greathouse IV, Dec 13 2021
Numbers m such that 2^m - m is divisible by 2. - Bernard Schott, Dec 15 2021

Crossrefs

Equals A005843 without the leading zero.
Bisection of A000027. Complement of A004273. - Omar E. Pol, Feb 25 2018
First row of A083140.
Cf. A005408.
Essentially the same as A163300, A103517, A051755, A005843 and A004277.

Programs

Formula

a(n) = 2*n, n >= 1.
G.f.: 2*x/(1 - x)^2; corrected by Ilya Gutkovskiy, Mar 29 2018
a(n) = 2*a(n-1) - a(n-2). - Wesley Ivan Hurt, Jul 17 2025

A328110 Fixed points of A327860: numbers k for which A003415(A276086(k)) = k, where A003415 is the arithmetic derivative, and A276086 is the primorial base exp-function.

Original entry on oeis.org

0, 1, 7, 8, 2556
Offset: 1

Views

Author

Antti Karttunen, Oct 08 2019

Keywords

Comments

Applying A276086 to these terms gives the fixed points of A327859: 1, 2, 10, 15, 5005, ..., i.e., A369650 without any of the terms of A100716.
No more terms below <= 2550136832.
From Antti Karttunen, Feb 09 2024: (Start)
The known five terms are all members of A276156, which is equal to the claim that the intersection of A048103 and A369650 is squarefree. See the example, and also comments in A351088 and in A380527.
Even terms here must be multiples of 4, see comment in A327860.
No terms of A047257 may occur in this sequence, which is equal to the claim that A276086(a(n)) is never a multiple of 9. See comment in A327859.
(End)

Examples

			Computing A327860(2556) is easy, because it is a member of A276156, as 2556 = 6 + 30 + 210 + 2310. Therefore A327860(2556) = A003415(5*7*11*13) = (5*7*11) + (5*7*13) + (5*11*13) + (7*11*13) = 2556, and 2556 is included in this sequence. - _Antti Karttunen_, Feb 04 2024
		

Crossrefs

After 0, the intersection of A351087 and A380527, thus like the latter, also this is conjectured to be a subsequence of A276156.
After two initial terms (0 & 1), a subsequence of A328118. Subsequence of A351088.

Programs

  • PARI
    A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
    isA328110(n) = (A327860(n) == n);

A047246 Numbers that are congruent to {0, 1, 2, 3} mod 6.

Original entry on oeis.org

0, 1, 2, 3, 6, 7, 8, 9, 12, 13, 14, 15, 18, 19, 20, 21, 24, 25, 26, 27, 30, 31, 32, 33, 36, 37, 38, 39, 42, 43, 44, 45, 48, 49, 50, 51, 54, 55, 56, 57, 60, 61, 62, 63, 66, 67, 68, 69, 72, 73, 74, 75, 78, 79, 80, 81, 84, 85, 86, 87, 90, 91, 92, 93, 96, 97, 98
Offset: 1

Views

Author

Keywords

Comments

The sequence is the interleaving of A047238 with A047241. - Guenther Schrack, Feb 12 2019

Crossrefs

Cf. A045331 (primes congruent to {1,2,3} mod 6), A047238, A047241.
Complement: A047257.

Programs

  • GAP
    Filtered([0..100],n->n mod 6 = 0 or n mod 6 = 1 or n mod 6 = 2 or n mod 6 = 3); # Muniru A Asiru, Feb 20 2019
  • Haskell
    a047246 n = a047246_list !! (n-1)
    a047246_list = [0..3] ++ map (+ 6) a047246_list
    -- Reinhard Zumkeller, Jan 15 2013
    
  • Magma
    [Floor((6/5)*Floor(5*(n-1)/4)) : n in [1..100]]; // Wesley Ivan Hurt, May 21 2016
    
  • Maple
    A047246:=n->(6*n-9-I^(2*n)-(1-I)*I^(-n)-(1+I)*I^n)/4: seq(A047246(n), n=1..100); # Wesley Ivan Hurt, May 21 2016
  • Mathematica
    Table[(6n-9-I^(2n)-(1-I)*I^(-n)-(1+I)*I^n)/4, {n, 80}] (* Wesley Ivan Hurt, May 21 2016 *)
  • PARI
    my(x='x+O('x^70)); concat([0], Vec(x^2*(1+x+x^2+3*x^3)/((1-x)*(1-x^4)))) \\ G. C. Greubel, Feb 16 2019
    
  • Sage
    a=(x^2*(1+x+x^2+3*x^3)/((1-x)*(1-x^4))).series(x, 72).coefficients(x, sparse=False); a[1:] # G. C. Greubel, Feb 16 2019
    

Formula

G.f.: x^2*(1+x+x^2+3*x^3) / ((1+x)*(1-x)^2*(1+x^2)). - R. J. Mathar, Oct 08 2011
a(n) = floor((6/5)*floor(5*(n-1)/4)). - Bruno Berselli, May 03 2016
From Wesley Ivan Hurt, May 21 2016: (Start)
a(n) = a(n-1) + a(n-4) - a(n-5) for n>5.
a(n) = (6*n - 9 - i^(2*n) - (1-i)*i^(-n) - (1+i)*i^n)/4 where i=sqrt(-1).
a(2*n) = A047241(n), a(2*n-1) = A047238(n). (End)
E.g.f.: (6 + sin(x) - cos(x) + (3*x - 4)*sinh(x) + (3*x - 5)*cosh(x))/2. - Ilya Gutkovskiy, May 21 2016
From Guenther Schrack, Feb 12 2019: (Start)
a(n) = (6*n - 9 - (-1)^n - 2*(-1)^(n*(n+1)/2))/4.
a(n) = a(n-4) + 6, a(1)=0, a(2)=1, a(3)=2, a(4)=3, for n > 4. (End)
Sum_{n>=2} (-1)^n/a(n) = Pi/(6*sqrt(3)) + 2*log(2)/3. - Amiram Eldar, Dec 16 2021
a(n)-a(n-1) = A093148(n-2). - R. J. Mathar, May 01 2024

Extensions

More terms from Wesley Ivan Hurt, May 21 2016

A047247 Numbers that are congruent to {2, 3, 4, 5} (mod 6).

Original entry on oeis.org

2, 3, 4, 5, 8, 9, 10, 11, 14, 15, 16, 17, 20, 21, 22, 23, 26, 27, 28, 29, 32, 33, 34, 35, 38, 39, 40, 41, 44, 45, 46, 47, 50, 51, 52, 53, 56, 57, 58, 59, 62, 63, 64, 65, 68, 69, 70, 71, 74, 75, 76, 77, 80, 81, 82, 83, 86, 87, 88, 89, 92, 93, 94, 95, 98, 99
Offset: 1

Views

Author

Keywords

Comments

The sequence is the interleaving of A047235 with A047270. - Guenther Schrack, Feb 10 2019
Numbers k for which A276076(k) and A276086(k) are multiples of three. For a simple proof, consider the penultimate digit in the factorial and primorial base expansions of n, A007623 and A049345. - Antti Karttunen, Feb 08 2024

Crossrefs

Programs

  • Magma
    [n : n in [0..100] | n mod 6 in [2, 3, 4, 5]]; // Wesley Ivan Hurt, May 21 2016
    
  • Maple
    A047247:=n->(6*n-1-I^(2*n)-(1-I)*I^(-n)-(1+I)*I^n)/4: seq(A047247(n), n=1..100); # Wesley Ivan Hurt, May 21 2016
  • Mathematica
    Table[(6n-1-I^(2n)-(1-I)*I^(-n)-(1+I)*I^n)/4, {n, 80}] (* Wesley Ivan Hurt, May 21 2016 *)
    LinearRecurrence[{1,0,0,1,-1},{2,3,4,5,8},70] (* Harvey P. Dale, May 25 2024 *)
  • PARI
    my(x='x+O('x^70)); Vec(x*(2+x+x^2+x^3+x^4)/((1-x)*(1-x^4))) \\ G. C. Greubel, Feb 16 2019
    
  • Sage
    a=(x*(2+x+x^2+x^3+x^4)/((1-x)*(1-x^4))).series(x, 72).coefficients(x, sparse=False); a[1:] # G. C. Greubel, Feb 16 2019

Formula

G.f.: x*(2+x+x^2+x^3+x^4) / ( (1+x)*(1+x^2)*(1-x)^2 ). - R. J. Mathar, Oct 08 2011
From Wesley Ivan Hurt, May 21 2016: (Start)
a(n) = a(n-1) + a(n-4) - a(n-5) for n>5.
a(n) = (6*n - 1 - i^(2*n) - (1-i)*i^(-n) - (1+i)*i^n)/4 where i = sqrt(-1).
a(2*n) = A047270(n), a(2*n-1) = A047235(n).
a(n) = A047227(n) + 1, a(1-n) = - A047227(n). (End)
From Guenther Schrack, Feb 10 2019: (Start)
a(n) = (6*n - 1 - (-1)^n -2*(-1)^(n*(n+1)/2))/4.
a(n) = a(n-4) + 6, a(1)=2, a(2)=3, a(3)=4, a(4)=5, for n > 4.
a(n) = A047227(n) + 1. a(n) = A047246(n) + 2. (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = sqrt(3)*Pi/12 - 2*log(2)/3 + log(3)/4. - Amiram Eldar, Dec 17 2021

Extensions

More terms from Wesley Ivan Hurt, May 21 2016

A380474 Numbers k such that A380459(k) has no divisors of the form p^p, while A003415(k) has such a divisor or is 0.

Original entry on oeis.org

1, 158, 482, 698, 914, 1238, 1346, 1454, 1994, 2102, 2426, 2642, 2858, 2966, 3398, 3506, 3722, 4262, 4478, 4586, 4694, 5234, 5342, 5666, 5774, 6098, 6638, 6746, 7286, 7394, 7934, 8042, 8258, 9014, 9122, 9446, 9662, 9986, 10202, 10418, 10958, 11282, 11498, 11714, 12146, 12686, 12794, 12902
Offset: 1

Views

Author

Antti Karttunen, Feb 02 2025

Keywords

Comments

Because all terms k of A380468 are squarefree, they are also in A048103, so A003415(k) is outside of A048103 only if k is in A327934 or k = 1.
Like A380478, also this (after the initial 1) is a subsequence of A039956. The first three terms k with A001222(k) > 2 are: 11362082, 16782482, 20965982.
Contains 2*A141964 as a subsequence, because all primes p congruent to 25 mod 27 are also congruent to 1 mod 6, therefore A276086(p) is a nonmultiple of 3 in those cases and thus coprime with A276086(2) = 3.
In contrast, neither all 2*(primes congruent to 3123 mod 3125) nor all 2*(primes congruent to (7^7)-2 mod 7^7, like 1647082) are present. The missing ones are those for which A276086(p) is a multiple of 9, i.e., when p is in A047257.

Crossrefs

Intersection of A380468 and ({1} U A327929), or equally of A380478 and ({1} U A327934).
After initial 1, a subsequence of A039956.

Programs

A380475 a(n) is the least term in A380468 that has exactly n prime factors.

Original entry on oeis.org

1, 2, 6, 186, 4686
Offset: 0

Views

Author

Antti Karttunen, Feb 03 2025

Keywords

Comments

If it exists, a(5) > 2^35. - Antti Karttunen, Feb 19 2025
Conjecture: Sequence is finite and a(4) is the last term. This is equivalent to the claim that no arithmetic derivative (A003415) of a product of five or more distinct primes, (i.e., the value of the (n-1)-st elementary symmetric polynomial formed from those n distinct primes) can be formed as a carry-free sum of those n summands in primorial base (A049345). See also A380476 and A380528, A380530.

Examples

			186 = 2*3*31 and A276086(186/2) = 2058 = 2 * 3 * 7^3, A276086(186/3) = 3 * 7^2, A276086(186/31) = 5, whose product =  2^1 * 3^2 * 5^1 * 7^5 = 1512630 = A380459(186), and as all the exponents are less than the corresponding primes, the product is in A048103, and because there are no any smaller number with three prime factors satisfying the same condition (of A380468), 186 is the term a(3) of this sequence. Note that A049345(A003415(186)) = 5121, where the digits are the exponents in the product read from the largest to the smallest prime factor.
See also the example in A380476 about 4686.
		

Crossrefs

Cf. A001221, A001222, A003415, A048103, A049345, A276086, A380459, A380468, A380476 (terms of A380468 with more than three prime factors).

Formula

a(n) = Min_{k in A380468} for which A001221(k) = n.

A380476 Numbers k with at least 4 prime factors such that A380459(k) is in A048103, i.e., has no divisors of the form p^p.

Original entry on oeis.org

4686, 32406, 184866, 209166, 388086, 1099626, 1714866, 2111406, 2166846, 2356206, 3081606, 3303366, 6445806, 11366106, 21621606, 23022366, 39824466, 39826986, 42882846, 43197846, 46043826, 58216686, 61265886, 63603546, 66496506, 66611166, 87941706, 88968246, 92086746, 97117026, 101108706, 103367886, 118743306, 119658066
Offset: 1

Views

Author

Antti Karttunen, Feb 04 2025

Keywords

Comments

Numbers m with four or more distinct prime factors such that their arithmetic derivative (A003415) can be formed as a carryless (or "carry-free") sum (in the primorial base, A049345) of the respective summands. See the example.
The terms are all squarefree and even (see A380468 and A380478 to find out why). Moreover, they are all multiples of six, because A380459(n) = Product_{d|n} A276086(n/d)^A349394(d) applied to a product of 2*p*q*r, with p, q, r three odd primes > 3 would yield three subproducts which would be multiples of 3 (consider A047247), so the 3-adic valuation of the whole product would be >= 3; hence the second smallest prime factor must be 3. For a similar reason, with terms that are product of four primes, the two remaining prime factors are either both of the form 6m+1 (A002476), or they are both of the form 6m-1 (A007528).
It is conjectured that there are no terms with more than four prime factors. See A380475 and A380528, A380530, also A380526.

Examples

			4686 = 2*3*11*71 and taking subproducts of three primes at time, we obtain 2*3*11 = 66, 2*3*71 = 426, 2*11*71 = 1562, 3*11*71 = 2343. Then A380459(4686) = A276086(66) * A276086(426) * A276086(1562) * A276086(2343) = 1622849599205985150 = 2^1 * 3^2 * 5^2 * 7^6 * 11^9 * 13^1, and because all the exponents are less than the corresponding primes, the product is in A048103.
Considering the primorial base expansions of the same summands (subproducts), we obtain
    2100  = A049345(66)
   20100  = A049345(426)
   73010  = A049345(1562)
  101011  = A049345(2343)
  ------
  196221  = A049345(A003415(4686)), with the summands adding together cleanly without any carries.
Note how the primorial base digits at the bottom are the exponents in the product A380459(4686) given above, read from the largest to the smallest prime factor
		

Crossrefs

Intersection of A033987 and A380468.
Subsequence of A005117, A358673, A380478.
Conjectured to be a subsequence of A046386.

Programs

  • PARI
    is_A380476(n) = (issquarefree(n) && (omega(n)>=4) && A380467(n)); \\ Note that issquarefree here is just an optimization as A380467(n) = 1 implies squarefreeness of n.
Showing 1-10 of 12 results. Next