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

A349444 Dirichlet convolution of A003602 (Kimberling's paraphrases) with A092673 (Dirichlet inverse of A001511).

Original entry on oeis.org

1, -1, 1, 0, 2, -1, 3, 0, 3, -2, 5, 0, 6, -3, 4, 0, 8, -3, 9, 0, 6, -5, 11, 0, 10, -6, 9, 0, 14, -4, 15, 0, 10, -8, 12, 0, 18, -9, 12, 0, 20, -6, 21, 0, 12, -11, 23, 0, 21, -10, 16, 0, 26, -9, 20, 0, 18, -14, 29, 0, 30, -15, 18, 0, 24, -10, 33, 0, 22, -12, 35, 0, 36, -18, 20, 0, 30, -12, 39, 0, 27, -20, 41, 0, 32
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2021

Keywords

Crossrefs

Cf. A001511, A003602, A008683, A092673, A349445 (Dirichlet inverse), A349446 (sum with it).
Cf. also A349431, A349447.

Programs

  • Mathematica
    s[n_] := MoebiusMu[n] - If[OddQ[n], 0, MoebiusMu[n/2]]; k[n_] := (n/2^IntegerExponent[n, 2] + 1)/2; a[n_] := DivisorSum[n, s[#]*k[n/#] &]; Array[a, 100] (* Amiram Eldar, Nov 19 2021 *)
  • PARI
    A003602(n) = (1+(n>>valuation(n,2)))/2;
    A092673(n) = if(n<1, 0, moebius(n) - if( n%2, 0, moebius(n/2))); \\ From A092673
    A349444(n) = sumdiv(n,d,A003602(n/d)*A092673(d));

Formula

a(n) = Sum_{d|n} A003602(n/d) * A092673(d).

A092149 Partial sums of A092673.

Original entry on oeis.org

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

Views

Author

Jon Perry, Mar 31 2004

Keywords

Crossrefs

Programs

  • Mathematica
    Accumulate@ Array[MoebiusMu[#] - If[OddQ@ #, 0, MoebiusMu[#/2]] &, 106] (* Michael De Vlieger, Mar 31 2021 *)
  • PARI
    a(n)=my(s);forstep(k=bitor(n\4+1,1),n\2,2,s-=moebius(k));forstep(k=bitor(n\2+1,1),n,2,s+=moebius(k)); s \\ Charles R Greathouse IV, Feb 07 2013
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A092149(n):
        if n == 1:
            return 1
        c, j = n+1, 2
        k1 = n//j
        while k1 > 1:
            j2 = n//k1 + 1
            c += (j2-j)*A092149(k1)
            j, k1 = j2, n//j2
        return j-c # Chai Wah Wu, Mar 31 2021

Formula

G.f. Sum_{n >= 1} a(n)*(x^n)/((1-x^n)*(x^(n+1)-1))*x = -(x^2) and -1/x. [Mats Granvik, Oct 11 2010]
On the Riemann hypothesis, |a(n)| = O(n^(1/2+e)) for any e > 0. - Charles R Greathouse IV, Feb 07 2013
a(1)=1, then for n>=2, Sum_{k=1..n} a(floor(n/k)) = 0. - Benoit Cloitre, Feb 21 2013
G.f. A(x) satisfies x * (1 - x) = Sum_{k>=1} (1 - x^k) * A(x^k). - Seiichi Manyama, Mar 31 2023

A091428 Numbers k such that abs(A092673(k)) = 1.

Original entry on oeis.org

1, 3, 4, 5, 7, 11, 12, 13, 15, 17, 19, 20, 21, 23, 28, 29, 31, 33, 35, 37, 39, 41, 43, 44, 47, 51, 52, 53, 55, 57, 59, 60, 61, 65, 67, 68, 69, 71, 73, 76, 77, 79, 83, 84, 85, 87, 89, 91, 92, 93, 95, 97, 101, 103, 105, 107, 109, 111, 113, 115, 116, 119, 123, 124, 127, 129
Offset: 1

Views

Author

Jon Perry, Mar 02 2004

Keywords

Comments

The asymptotic density of this sequence is 5/Pi^2 = 0.506605... . - Amiram Eldar, Jan 11 2023

Crossrefs

Disjoint union of A056911 and A081770.
Cf. A092673, A359592 (characteristic function).
Distinct from A020486.

Programs

  • Mathematica
    Select[Range[130], Abs[MoebiusMu[#] - If[EvenQ[#], MoebiusMu[#/2], 0]] == 1 &] (* Amiram Eldar, Jan 11 2023 *)

A092152 Sign of A092673(n).

Original entry on oeis.org

1, -1, -1, 1, -1, 1, -1, 0, 0, 1, -1, -1, -1, 1, 1, 0, -1, 0, -1, -1, 1, 1, -1, 0, 0, 1, 0, -1, -1, -1, -1, 0, 1, 1, 1, 0, -1, 1, 1, 0, -1, -1, -1, -1, 0, 1, -1, 0, 0, 0, 1, -1, -1, 0, 1, 0, 1, 1, -1, 1, -1, 1, 0, 0, 1, -1, -1, -1, 1, -1, -1, 0, -1, 1, 0, -1, 1, -1, -1, 0, 0, 1, -1, 1, 1, 1, 1, 0, -1, 0, 1, -1, 1, 1, 1, 0, -1, 0, 0, 0, -1, -1, -1, 0, -1, 1
Offset: 1

Views

Author

Jon Perry, Mar 31 2004

Keywords

Comments

Multiplicative because A092673 is. - Andrew Howroyd, Jul 27 2018

Crossrefs

Cf. A092673.

Programs

  • PARI
    a(n)=sign(moebius(n) - if(n%2, 0, moebius(n/2))) \\ Andrew Howroyd, Jul 27 2018

Formula

a(2*n) = sign(mu(2*n) - mu(n)); a(2*n + 1) = mu(2*n + 1). - Andrew Howroyd, Jul 27 2018

A092151 A092673(n)*A092674(n).

Original entry on oeis.org

0, -6, -3, 1, -5, 0, -7, 0, 0, 4, -11, -5, -13, 8, 7, 0, -17, 0, -19, -9, 11, 16, -23, 0, 0, 20, 0, -13, -29, -20, -31, 0, 19, 28, 23, 0, -37, 32, 23, 0, -41, -28, -43, -21, 0, 40, -47, 0, 0, 0
Offset: 1

Views

Author

Jon Perry, Mar 31 2004

Keywords

A092155 First differences of A092673.

Original entry on oeis.org

-3, 1, 2, -2, 3, -3, 1, 0, 2, -3, 0, 0, 3, -1, -1, -1, 1, -1, 0, 2, 1, -3, 1, 0, 2, -2, -1, 0, -1, 1, 1, 1, 1, -1, -1, -1, 3, -1, -1, -1, -1, 1, 0, 1, 2, -3, 1, 0, 0, 1, -2, 0, 1, 1, -1, 1, 1, -3, 2, -2, 3, -2, 0, 1, -3, 1, 0, 2, -3, 1, 1, -1, 3, -2, -1, 2, -3, 1, 1, 0, 2, -3, 2, 0, 1, -1, -1, -1, 1, 1, -2, 2, 1, -1, -1, -1, 1, 0, 0, -1, -1, 1, 1, -1, 3
Offset: 1

Views

Author

Jon Perry, Mar 31 2004

Keywords

Comments

Empirical: For n > 1, letting M(n) denote the n X n matrix whereby the (i,j)-entry of M(n) is -Sum_{k=1..j} floor(i/k), we have that the first n-1 entries in the first column of the inverse of M(n) are a(1) through a(n-1). - John M. Campbell, Aug 30 2017
John Campbell's empirical data is proved at the Mathematics Stack Exchange link. - Sungjin Kim, Jul 17 2019

Crossrefs

Cf. A092673.

A001511 The ruler function: exponent of the highest power of 2 dividing 2n. Equivalently, the 2-adic valuation of 2n.

Original entry on oeis.org

1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 5, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 6, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 5, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 5, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1, 6, 1, 2, 1, 3, 1, 2, 1, 4, 1
Offset: 1

Views

Author

Keywords

Comments

Number of 2's dividing 2*n.
a(n) is equivalently the exponent of the smallest power of 2 which does not divide n. - David James Sycamore, Oct 02 2023
a(n) - 1 is the number of trailing zeros in the binary expansion of n.
If you are counting in binary and the least significant bit is numbered 1, the next bit is 2, etc., a(n) is the bit that is incremented when increasing from n-1 to n. - Jud McCranie, Apr 26 2004
Number of steps to reach an integer starting with (n+1)/2 and using the map x -> x*ceiling(x) (cf. A073524).
a(n) is the number of the disk to be moved at the n-th step of the optimal solution to Towers of Hanoi problem (comment from Andreas M. Hinz).
Shows which bit to flip when creating the binary reflected Gray code (bits are numbered from the right, offset is 1). This is essentially equivalent to Hinz's comment. - Adam Kertesz, Jul 28 2001
a(n) is the Hamming distance between n and n-1 (in binary). This is equivalent to Kertesz's comments above. - Tak-Shing Chan (chan12(AT)alumni.usc.edu), Feb 25 2003
Let S(0) = {1}, S(n) = {S(n-1), S(n-1)-{x}, x+1} where x = last term of S(n-1); sequence gives S(infinity). - Benoit Cloitre, Jun 14 2003
The sum of all terms up to and including the first occurrence of m is 2^m-1. - Donald Sampson (marsquo(AT)hotmail.com), Dec 01 2003
m appears every 2^m terms starting with the 2^(m-1)th term. - Donald Sampson (marsquo(AT)hotmail.com), Dec 08 2003
Sequence read mod 4 gives A092412. - Philippe Deléham, Mar 28 2004
If q = 2n/2^A001511(n) and if b(m) is defined by b(0)=q-1 and b(m)=2*b(m-1)+1, then 2n = b(A001511(n)) + 1. - Gerald McGarvey, Dec 18 2004
Repeating pattern ABACABADABACABAE ... - Jeremy Gardiner, Jan 16 2005
Relation to C(n) = Collatz function iteration using only odd steps: a(n) is the number of right bits set in binary representation of A004767(n) (numbers of the form 4*m+3). So for m=A004767(n) it follows that there are exactly a(n) recursive steps where m
Between every two instances of any positive integer m there are exactly m distinct values (1 through m-1 and one value greater than m). - Franklin T. Adams-Watters, Sep 18 2006
Number of divisors of n of the form 2^k. - Giovanni Teofilatto, Jul 25 2007
Every prefix up to (but not including) the first occurrence of some k >= 2 is a palindrome. - Gary W. Adamson, Sep 24 2008
1 interleaved with (2 interleaved with (3 interleaved with ( ... ))). - Eric D. Burgess (ericdb(AT)gmail.com), Oct 17 2009
A054525 (Möbius transform) * A001511 = A036987 = A047999^(-1) * A001511. - Gary W. Adamson, Oct 26 2009
Equals A051731 * A036987, (inverse Möbius transform of the Fredholm-Rueppel sequence) = A047999 * A036987. - Gary W. Adamson, Oct 26 2009
Cf. A173238, showing links between generalized ruler functions and A000041. - Gary W. Adamson, Feb 14 2010
Given A000041, P(x) = A(x)/A(x^2) with P(x) = (1 + x + 2x^2 + 3x^3 + 5x^4 + 7x^5 + ...), A(x) = (1 + x + 3x^2 + 4x^3 + 10x^4 + 13x^5 + ...), A(x^2) = (1 + x^2 + 3x^4 + 4x^6 + 10x^8 + ...), where A092119 = (1, 1, 3, 4, 10, ...) = Euler transform of the ruler sequence, A001511. - Gary W. Adamson, Feb 11 2010
Subtracting 1 from every term and deleting any 0's yields the same sequence, A001511. - Ben Branman, Dec 28 2011
In the listing of the compositions of n as lists in lexicographic order, a(k) is the last part of composition(k) for all k <= 2^(n-1) and all n, see example. - Joerg Arndt, Nov 12 2012
According to Hinz, et al. (see links), this sequence was studied by Louis Gros in his 1872 pamphlet "Théorie du Baguenodier" and has therefore been called the Gros sequence.
First n terms comprise least squarefree word of length n using positive integers, where "squarefree" means that the word contains no consecutive identical subwords; e.g., 1 contains no square; 11 contains a square but 12 does not; 121 contains no square; both 1211 and 1212 have squares but 1213 does not; etc. - Clark Kimberling, Sep 05 2013
Length of 0-run starting from 2 (10, 100, 110, 1000, 1010, ...), or length of 1-run starting from 1 (1, 11, 101, 111, 1001, 1011, ...) of every second number, from right to left in binary representation. - Armands Strazds, Apr 13 2017
a(n) is also the frequency of the largest part in the integer partition having viabin number n. The viabin number of an integer partition is defined in the following way. Consider the southeast border of the Ferrers board of the integer partition and consider the binary number obtained by replacing each east step with 1 and each north step, except the last one, with 0. The corresponding decimal form is, by definition, the viabin number of the given integer partition. "Viabin" is coined from "via binary". For example, consider the integer partition [2,2,2,1]. The southeast border of its Ferrers board yields 10100, leading to the viabin number 20. - Emeric Deutsch, Jul 24 2017
As A000005(n) equals the number of even divisors of 2n and A001227(n) = A001227(2n), the formula A001511(n) = A000005(n)/A001227(n) might be read as "The number of even divisors of 2n is always divisible by the number of odd divisors of 2n" (where number of divisors means sum of zeroth powers of divisors). Conjecture: For any nonnegative integer k, the sum of the k-th powers of even divisors of n is always divisible by the sum of the k-th powers of odd divisors of n. - Ivan N. Ianakiev, Jul 06 2019
From Benoit Cloitre, Jul 14 2022: (Start)
To construct the sequence, start from 1's separated by a place 1,,1,,1,,1,,1,,1,,1,,1,,1,,1,,1,,1,,1,,1,...
Then put the 2's in every other remaining place
1,2,1,,1,2,1,,1,2,1,,1,2,1,,1,2,1,,1,2,1,,1,2,1,...
Then the 3's in every other remaining place
1,2,1,3,1,2,1,,1,2,1,3,1,2,1,,1,2,1,3,1,2,1,,1,2,1,...
Then the 4's in every other remaining place
1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,,1,2,1,3,1,2,1,4,1,2,1,...
By iterating this process, we get the ruler function 1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,5,1,2,1,3,1,2,1,4,1,2,1,... (End)
a(n) is the least positive integer k for which there does not exist i+j=n and a(i)=a(j)=k (cf. A322523). - Rémy Sigrist and Jianing Song, Aug 23 2022
a(n) is the smallest positive integer that does not occur in the coincidences of the sequence so far a(1..n-1) and its reverse. - Neal Gersh Tolunsky, Jan 18 2023
The geometric mean of this sequence approaches the Somos constant (A112302). - Jwalin Bhatt, Jan 31 2025

Examples

			For example, 2^1|2, 2^2|4, 2^1|6, 2^3|8, 2^1|10, 2^2|12, ... giving the initial terms 1, 2, 1, 3, 1, 2, ...
From _Omar E. Pol_, Jun 12 2009: (Start)
Triangle begins:
1;
2,1;
3,1,2,1;
4,1,2,1,3,1,2,1;
5,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1;
6,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,5,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1;
7,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,5,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,6,1,2,1,3,...
(End)
S(0) = {} S(1) = 1 S(2) = 1, 2, 1 S(3) = 1, 2, 1, 3, 1, 2, 1 S(4) = 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 2, 1. - Yann David (yann_david(AT)hotmail.com), Mar 21 2010
From _Joerg Arndt_, Nov 12 2012: (Start)
The 16 compositions of 5 as lists in lexicographic order:
[ n]  a(n)  composition
[ 1]  [ 1]  [ 1 1 1 1 1 ]
[ 2]  [ 2]  [ 1 1 1 2 ]
[ 3]  [ 1]  [ 1 1 2 1 ]
[ 4]  [ 3]  [ 1 1 3 ]
[ 5]  [ 1]  [ 1 2 1 1 ]
[ 6]  [ 2]  [ 1 2 2 ]
[ 7]  [ 1]  [ 1 3 1 ]
[ 8]  [ 4]  [ 1 4 ]
[ 9]  [ 1]  [ 2 1 1 1 ]
[10]  [ 2]  [ 2 1 2 ]
[11]  [ 1]  [ 2 2 1 ]
[12]  [ 3]  [ 2 3 ]
[13]  [ 1]  [ 3 1 1 ]
[14]  [ 2]  [ 3 2 ]
[15]  [ 1]  [ 4 1 ]
[16]  [ 5]  [ 5 ]
a(n) is the last part in each list.
(End)
From _Omar E. Pol_, Aug 20 2013: (Start)
Also written as a triangle in which the right border gives A000027 and row lengths give A011782 and row sums give A000079 the sequence begins:
1;
2;
1,3;
1,2,1,4;
1,2,1,3,1,2,1,5;
1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,6;
1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,5,1,2,1,3,1,2,1,4,1,2,1,3,1,2,1,7;
(End)
G.f. = x + 2*x^2 + x^3 + 3*x^4 + x^5 + 2*x^6 + x^7 + 4*x^8 + x^9 + 2*x^10 + ...
		

References

  • J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003.
  • E. R. Berlekamp, J. H. Conway and R. K. Guy, Winning Ways, Academic Press, NY, 2 vols., 2nd ed., 2001-2003; see Dim- and Dim+ on p. 98; Dividing Rulers, on pp. 436-437; The Ruler Game, pp. 469-470; Ruler Fours, Fives, ... Fifteens on p. 470.
  • L. Gros, Théorie du Baguenodier, Aimé Vingtrinier, Lyon, 1872.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section E22.
  • A. M. Hinz, The Tower of Hanoi, in Algebras and combinatorics (Hong Kong, 1997), 277-289, Springer, Singapore, 1999.
  • D. E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, Section 7.1.3, Problem 41, p. 589.
  • Andrew Schloss, "Towers of Hanoi" composition, in The Digital Domain. Elektra/Asylum Records 9 60303-2, 1983. Works by Jaffe (Finale to "Silicon Valley Breakdown"), McNabb ("Love in the Asylum"), Schloss ("Towers of Hanoi"), Mattox ("Shaman"), Rush, Moorer ("Lions are Growing") and others.
  • 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

Column 1 of table A050600.
Sequence read mod 2 gives A035263.
Sequence is bisection of A007814, A050603, A050604, A067029, A089309.
This is Guy Steele's sequence GS(4, 2) (see A135416).
Cf. A005187 (partial sums), A085058 (bisection), A112302 (geometric mean), A171977 (2^a(n)).
Cf. A287896, A002487, A209229 (Mobius trans.), A092673 (Dirichlet inv.).
Cf. generalized ruler functions for k=3,4,5: A051064, A115362, A055457.

Programs

  • Haskell
    a001511 n = length $ takeWhile ((== 0) . (mod n)) a000079_list
    -- Reinhard Zumkeller, Sep 27 2011
    
  • Haskell
    a001511 n | odd n = 1 | otherwise = 1 + a001511 (n `div` 2)
    -- Walt Rorie-Baety, Mar 22 2013
    
  • MATLAB
    nmax=5;r=1;for n=2:nmax;r=[r n r];end % Adriano Caroli, Feb 26 2016
    
  • Magma
    [Valuation(2*n,2): n in [1..105]]; // Bruno Berselli, Nov 23 2015
    
  • Maple
    A001511 := n->2-wt(n)+wt(n-1); # where wt is defined in A000120
    # This is the binary logarithm of the denominator of (256^n-1)B_{8n}/n, in Maple parlance a := n -> log[2](denom((256^n-1)*bernoulli(8*n)/n)). - Peter Luschny, May 31 2009
    A001511 := n -> padic[ordp](2*n,2): seq(A001511(n), n=1..105);  # Peter Luschny, Nov 26 2010
    a:= n-> ilog2((Bits[Xor](2*n, 2*n-1)+1)/2): seq(a(n), n=1..50);  # Gary Detlefs, Dec 13 2018
  • Mathematica
    Array[ If[ Mod[ #, 2] == 0, FactorInteger[ # ][[1, 2]], 0] &, 105] + 1 (* or *)
    Nest[ Flatten[ # /. a_Integer -> {1, a + 1}] &, {1}, 7] (* Robert G. Wilson v, Mar 04 2005 *)
    IntegerExponent[2*n, 2] (* Alexander R. Povolotsky, Aug 19 2011 *)
    myHammingDistance[n_, m_] := Module[{g = Max[m, n], h = Min[m, n]}, b1 = IntegerDigits[g, 2]; b2 = IntegerDigits[h, 2, Length[b1]]; HammingDistance[b1, b2]] (* Vladimir Shevelev A206853 *) Table[ myHammingDistance[n, n - 1], {n, 111}] (* Robert G. Wilson v, Apr 05 2012 *)
    Table[Position[Reverse[IntegerDigits[n,2]],1,1,1],{n,110}]//Flatten (* Harvey P. Dale, Aug 18 2017 *)
  • PARI
    a(n) = sum(k=0,floor(log(n)/log(2)),floor(n/2^k)-floor((n-1)/2^k)) /* Ralf Stephan */
    
  • PARI
    a(n)=if(n%2,1,factor(n)[1,2]+1) /* Jon Perry, Jun 06 2004 */
    
  • PARI
    {a(n) = if( n, valuation(n, 2) + 1, 0)}; /* Michael Somos, Sep 30 2006 */
    
  • PARI
    {a(n)=if(n==1,1,polcoeff(x-sum(k=1, n-1, a(k)*x^k*(1-x^k)*(1-x+x*O(x^n))), n))} /* Paul D. Hanna, Jun 22 2007 */
    
  • Python
    def a(n): return bin(n)[2:][::-1].index("1") + 1 # Indranil Ghosh, May 11 2017
    
  • Python
    A001511 = lambda n: (n&-n).bit_length() # M. F. Hasler, Apr 09 2020
    
  • Python
    def A001511(n): return (~n & n-1).bit_length()+1 # Chai Wah Wu, Jul 01 2022
    
  • Sage
    [valuation(2*n,2) for n in (1..105)]  # Bruno Berselli, Nov 23 2015
    
  • Scheme
    (define (A001511 n) (let loop ((n n) (e 1)) (if (odd? n) e (loop (/ n 2) (+ 1 e))))) ;; Antti Karttunen, Oct 06 2017

Formula

a(n) = A007814(n) + 1 = A007814(2*n).
a(2*n+1) = 1; a(2*n) = 1 + a(n). - Philippe Deléham, Dec 08 2003
a(n) = 2 - A000120(n) + A000120(n-1), n >= 1. - Daniele Parisse
a(n) = 1 + log_2(abs(A003188(n) - A003188(n-1))).
Multiplicative with a(p^e) = e+1 if p = 2; 1 if p > 2. - David W. Wilson, Aug 01 2001
For any real x > 1/2: lim_{N->infinity} (1/N)*Sum_{n=1..N} x^(-a(n)) = 1/(2*x-1); also lim_{N->infinity} (1/N)*Sum_{n=1..N} 1/a(n) = log(2). - Benoit Cloitre, Nov 16 2001
s(n) = Sum_{k=1..n} a(k) is asymptotic to 2*n since s(n) = 2*n - A000120(n). - Benoit Cloitre, Aug 31 2002
For any n >= 0, for any m >= 1, a(2^m*n + 2^(m-1)) = m. - Benoit Cloitre, Nov 24 2002
a(n) = Sum_{d divides n and d is odd} mu(d)*tau(n/d). - Vladeta Jovovic, Dec 04 2002
G.f.: A(x) = Sum_{k>=0} x^(2^k)/(1-x^(2^k)). - Ralf Stephan, Dec 24 2002
a(1) = 1; for n > 1, a(n) = a(n-1) + (-1)^n*a(floor(n/2)). - Vladeta Jovovic, Apr 25 2003
A fixed point of the mapping 1->12; 2->13; 3->14; 4->15; 5->16; ... . - Philippe Deléham, Dec 13 2003
Product_{k>0} (1+x^k)^a(k) is g.f. for A000041(). - Vladeta Jovovic, Mar 26 2004
G.f. A(x) satisfies A(x) = A(x^2) + x/(1-x). - Franklin T. Adams-Watters, Feb 09 2006
a(A118413(n,k)) = A002260(n,k); = a(A118416(n,k)) = A002024(n,k); a(A014480(n)) = A003602(A014480(n)). - Reinhard Zumkeller, Apr 27 2006
Ordinal transform of A003602. - Franklin T. Adams-Watters, Aug 28 2006 (The ordinal transform of a sequence b_0, b_1, b_2, ... is the sequence a_0, a_1, a_2, ... where a_n is the number of times b_n has occurred in {b_0 ... b_n}.)
Could be extended to n <= 0 using a(-n) = a(n), a(0) = 0, a(2*n) = a(n)+1 unless n=0. - Michael Somos, Sep 30 2006
A094267(2*n) = A050603(2*n) = A050603(2*n + 1) = a(n). - Michael Somos, Sep 30 2006
Sequence = A129360 * A000005 = M*V, where M = an infinite lower triangular matrix and V = d(n) as a vector: [1, 2, 2, 3, 2, 4, ...]. - Gary W. Adamson, Apr 15 2007
Row sums of triangle A130093. - Gary W. Adamson, May 13 2007
Dirichlet g.f.: zeta(s)*2^s/(2^s-1). - Ralf Stephan, Jun 17 2007
a(n) = -Sum_{d divides n} mu(2*d)*tau(n/d). - Benoit Cloitre, Jun 21 2007
G.f.: x/(1-x) = Sum_{n>=1} a(n)*x^n*( 1 - x^n ). - Paul D. Hanna, Jun 22 2007
2*n = 2^a(n)* A000265(n). - Eric Desbiaux, May 14 2009 [corrected by Alejandro Erickson, Apr 17 2012]
Multiplicative with a(2^k) = k + 1, a(p^k) = 1 for any odd prime p. - Franklin T. Adams-Watters, Jun 09 2009
With S(n): 2^n - 1 first elements of the sequence then S(0) = {} (empty list) and if n > 0, S(n) = S(n-1), n, S(n-1). - Yann David (yann_david(AT)hotmail.com), Mar 21 2010
a(n) = log_2(A046161(n)/A046161(n-1)). - Johannes W. Meijer, Nov 04 2012
a((2*n-1)*2^p) = p+1, p >= 0 and n >= 1. - Johannes W. Meijer, Feb 05 2013
a(n+1) = 1 + Sum_{j=0..ceiling(log_2(n+1))} (j * (1 - abs(sign((n mod 2^(j + 1)) - 2^j + 1)))). - Enrico Borba, Oct 01 2015
Conjecture: a(n) = A181988(n)/A003602(n). - L. Edson Jeffery, Nov 21 2015
a(n) = log_2(A006519(n)) + 1. - Doug Bell, Jun 02 2017
Inverse Moebius transform of A209229. - Andrew Howroyd, Aug 04 2018
a(n) = 1 + (A183063(n)/A001227(n)). - Omar E. Pol, Nov 06 2018 (after Franklin T. Adams-Watters)
a(n) = log_2((Xor(2*n,2*n-1)+1)/2). - Gary Detlefs, Dec 13 2018
(2^(a(n)-1)-1)*(n mod 4) = 2*floor(((n+1) mod 4)/3). - Gary Detlefs, Dec 14 2018
a(n) = A000005(n)/A001227(n). - Ivan N. Ianakiev, Jul 05 2019
a(n) = Sum_{j=1..r} (j/2^j)*(Product_{k=1..j} (1 - (-1)^floor( (n+2^(j-1))/2^(k-1) ))), for n < a predefined 2^r. - Adriano Caroli, Sep 30 2019

Extensions

Name edited following suggestion by David James Sycamore, Oct 05 2023

A039956 Even squarefree numbers.

Original entry on oeis.org

2, 6, 10, 14, 22, 26, 30, 34, 38, 42, 46, 58, 62, 66, 70, 74, 78, 82, 86, 94, 102, 106, 110, 114, 118, 122, 130, 134, 138, 142, 146, 154, 158, 166, 170, 174, 178, 182, 186, 190, 194, 202, 206, 210, 214, 218, 222, 226, 230, 238, 246, 254, 258, 262
Offset: 1

Author

Keywords

Comments

Sum of even divisors = 2 * the sum of odd divisors. - Amarnath Murthy, Sep 07 2002
From Daniel Forgues, May 27 2009: (Start)
a(n) = n * (3/1) * zeta(2) + O(n^(1/2)) = n * (3/1) * (Pi^2 / 6) + O(n^(1/2)).
For any prime p_i, the n-th squarefree number even to p_i (divisible by p_i) is:
n * ((p_i + 1)/1) * zeta(2) + O(n^(1/2)) = n * ((p_i + 1)/1) * (Pi^2 / 6) + O(n^(1/2)).
For any prime p_i, there are as many squarefree numbers having p_i as a factor as squarefree numbers not having p_i as a factor amongst all the squarefree numbers (one-to-one correspondence, both cardinality aleph_0).
E.g., there are as many even squarefree numbers as there are odd squarefree numbers.
For any prime p_i, the density of squarefree numbers having p_i as a factor is 1/p_i of the density of squarefree numbers not having p_i as a factor.
E.g., the density of even squarefree numbers is 1/p_i = 1/2 of the density of odd squarefree numbers (which means that 1/(p_i + 1) = 1/3 of the squarefree numbers are even and p_i/(p_i + 1) = 2/3 are odd) and as a consequence the n-th even squarefree number is very nearly p_i = 2 times the n-th odd squarefree number (which means that the n-th even squarefree number is very nearly (p_i + 1) = 3 times the n-th squarefree number while the n-th odd squarefree number is very nearly (p_i + 1)/ p_i = 3/2 the n-th squarefree number).
(End)
Apart from first term, these are the tau2-atoms as defined in [Anderson, Frazier] and [Lanterman]. - Michel Marcus, May 15 2019

References

  • Richard A. Mollin, Quadratics, CRC Press, 1996, Tables B1-B3.

Crossrefs

Programs

  • Haskell
    a039956 n = a039956_list !! (n-1)
    a039956_list = filter even a005117_list  -- Reinhard Zumkeller, Aug 15 2011
    
  • Magma
    [n: n in [2..262 by 2] | IsSquarefree(n)];  // Bruno Berselli, Mar 03 2011
    
  • Maple
    select(numtheory:-issqrfree,[seq(i,i=2..1000,4)]); # Robert Israel, Dec 23 2015
  • Mathematica
    Select[Range[2,270,2],SquareFreeQ] (* Harvey P. Dale, Jul 23 2011 *)
  • PARI
    is(n)=n%4==2 && issquarefree(n) \\ Charles R Greathouse IV, Sep 13 2013

Formula

Numbers k such that A092673(k) = +- 2. - Jon Perry, Mar 02 2004
Sum_{n>=1} 1/a(n)^s = zeta(s)/((1+2^s)*zeta(2*s)). - Enrique Pérez Herrero, Sep 15 2012 [corrected by Amiram Eldar, Sep 26 2023]
a(n) = 2*A056911(n). - Robert Israel, Dec 23 2015
a(n) = 2*(1+2*A264387(n)), n >= 1. - Wolfdieter Lang, Dec 24 2015

A001747 2 together with primes multiplied by 2.

Original entry on oeis.org

2, 4, 6, 10, 14, 22, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, 146, 158, 166, 178, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 386, 394, 398, 422, 446, 454, 458, 466, 478, 482, 502
Offset: 1

Keywords

Comments

When supplemented with 8, may be considered the "even primes", since these are the even numbers n = 2k which are divisible just by 1, 2, k and 2k. - Louis Zuckerman (louis(AT)trapezoid.com), Sep 12 2000
Sequence gives solutions of sigma(n) - phi(n) = n + tau(n) where tau(n) is the number of divisors of n.
Numbers n such that sigma(n) = 3*(n - phi(n)).
Except for 2, orders of non-cyclic groups k (in A060679(n)) such that x^k==1 (mod k) has only 1 solution 2<=x<=k. - Benoit Cloitre, May 10 2002
Numbers n such that A092673(n) = 2. - Jon Perry, Mar 02 2004
Except for initial terms, this sequence = A073582 = A074845 = A077017. Starting with the term 10, they are identical. - Robert G. Wilson v, Jun 15 2004
Together with 8 and 16, even numbers n such that n^2 does not divide (n/2)!. - Arkadiusz Wesolowski, Jul 16 2011
Twice noncomposite numbers. - Omar E. Pol, Jan 30 2012

Crossrefs

Equals {2} UNION {A100484}.

Programs

  • GAP
    Concatenation([2], List([1..60], n-> 2*Primes[n])); # G. C. Greubel, May 18 2019
  • Magma
    [2] cat [2*NthPrime(n): n in [1..60]]; // G. C. Greubel, May 18 2019
    
  • Mathematica
    Join[{2},2*Prime[Range[60]]] (* Harvey P. Dale, Jul 23 2013 *)
  • PARI
    print1(2);forprime(p=2,97,print1(", "2*p)) \\ Charles R Greathouse IV, Jan 31 2012
    
  • Sage
    [2]+[2*nth_prime(n) for n in (1..60)] # G. C. Greubel, May 18 2019
    

Formula

a(n) = A001043(n) - A001223(n+1), except for initial term.
a(n) = A116366(n-2,n-2) for n>2. - Reinhard Zumkeller, Feb 06 2006
A006093(n) = A143201(a(n+1)) for n>1. - Reinhard Zumkeller, Aug 12 2008
a(n) = 2*A008578(n). - Omar E. Pol, Jan 30 2012, and Reinhard Zumkeller, Feb 16 2012

A154955 a(1) = 1, a(2) = -1, followed by 0, 0, 0, ... .

Original entry on oeis.org

1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Author

Mats Granvik, Jan 18 2009

Keywords

Comments

Matrix inverse of A000012.
Moebius transform of the sequence A000035. Dirichlet inverse of A209229. Partial sums of a(n) is characteristic function of 1 (A063524). a(n)=(-1)^(n+1)*A019590(n). a(n) for n >= 1 is Dirichlet convolution of following functions b(n), c(n), a(n) = Sum_{d|n} b(d)*c(n/d): a(n) = A000012(n) * A092673(n). Examples of Dirichlet convolutions with function a(n), i.e. b(n) = Sum_{d|n} a(d)*c(n/d): a(n) * A000012(n) = A000035(n), a(n) * A000027(n) = A026741(n), a(n) * A008683(n) = A092673(n), a(n) * A036987(n-1) = A063524(n), a(n) * A000005(n) = A001227(n). - Jaroslav Krizek, Mar 21 2009
The Kn21 sums, see A180662, of triangle A108299 equal the terms of this sequence. - Johannes W. Meijer, Aug 14 2011
{a(n-1)}A132393.%20-%20_Wolfdieter%20Lang">{n>=1}, gives the alternating row sums of A132393. - _Wolfdieter Lang, May 09 2017
With offset 0 the alternating row sums of A097805. - Peter Luschny, Sep 07 2017

Programs

Formula

G.f.: A(x) = x - x^2 = x / (1 + x / (1 - x)). - Michael Somos, Jan 03 2013
a(n) = (2/sqrt(3))*sin((2*Pi/3)*n!). - Lorenzo Pinlac, Jan 16 2022
a(n) = [n = 1] - [n = 2], where [] is the Iverson bracket. - Wesley Ivan Hurt, Jun 22 2024
Multiplicative with a(2) = -1, a(2^e) = 0 if e > 1, a(p^e) = 0 if p > 2. - Antti Karttunen, Dec 17 2024

Extensions

Keyword:mult added by Antti Karttunen, Dec 17 2024
Showing 1-10 of 19 results. Next