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-8 of 8 results.

A374354 Irregular table T(n, k), n >= 0, 0 <= k < A277561(n), read by rows; the n-th row lists the fibbinary numbers f <= n such that n - f is also a fibbinary number whose binary expansion has no common 1's with that of f (where fibbinary numbers correspond to A003714).

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 2, 0, 4, 0, 1, 4, 5, 2, 4, 2, 5, 0, 8, 0, 1, 8, 9, 0, 2, 8, 10, 1, 2, 9, 10, 4, 8, 4, 5, 8, 9, 4, 10, 5, 10, 0, 16, 0, 1, 16, 17, 0, 2, 16, 18, 1, 2, 17, 18, 0, 4, 16, 20, 0, 1, 4, 5, 16, 17, 20, 21, 2, 4, 18, 20, 2, 5, 18, 21, 8, 16, 8, 9, 16, 17
Offset: 0

Views

Author

Rémy Sigrist, Jul 06 2024

Keywords

Comments

In other words, we partition n into pairs of fibbinary numbers whose binary expansions have no common 1's and list the corresponding fibbinary numbers to get the n-th row.

Examples

			Triangle T(n, k) begins:
  n   n-th row
  --  -----------
   0  0
   1  0, 1
   2  0, 2
   3  1, 2
   4  0, 4
   5  0, 1, 4, 5
   6  2, 4
   7  2, 5
   8  0, 8
   9  0, 1, 8, 9
  10  0, 2, 8, 10
  11  1, 2, 9, 10
  12  4, 8
  13  4, 5, 8, 9
  14  4, 10
  15  5, 10
  16  0, 16
		

Crossrefs

See A295989 and A374361 for similar sequences.

Programs

  • PARI
    row(n) = { my (r = [0], e, x, y, b); while (n, x = y = 0; e = valuation(n, 2); for (k = 0, oo, if (bittest(n, e+k), n -= b = 2^(e+k); [x, y] = [y + b, x], r = concat([v + y | v <- r], [v + x | v <- r]); break;););); return (r); }

Formula

T(n, 0) = 0 iff n is a fibbinary number.
T(n, k) + T(n, A277561(n)-1-k) = n.
T(n, 0) = A374355(n).
T(n, A277561(n)-1) = A374356(n).
Sum_{k = 0..A277561(n)-1} T(n, k) = n * 2^A037800(n).

A374394 Irregular table T(n, k), n >= 0, 0 <= k < A277561(1+A003754(n)), read by rows; the n-th row lists the numbers z <= n such that the Zeckendorf representations of z and n-z have no common Fibonacci numbers and when combined together correspond to the lazy Fibonacci representation of n.

Original entry on oeis.org

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

Views

Author

Rémy Sigrist, Jul 07 2024

Keywords

Examples

			Triangle T(n, k) begins:
  n   n-th row
  --  ------------------------
   0  0
   1  0, 1
   2  0, 2
   3  1, 2
   4  0, 1, 3, 4
   5  2, 3
   6  2, 4
   7  0, 2, 5, 7
   8  1, 2, 6, 7
   9  3, 4, 5, 6
  10  3, 7
  11  4, 7
  12  0, 1, 3, 4, 8, 9, 11, 12
  13  2, 3, 10, 11
  14  2, 4, 10, 12
		

Crossrefs

Programs

  • PARI
    \\ See Links section.

Formula

T(n, k) = A022290(A374354(1+A003754(n)), k).

A034444 a(n) is the number of unitary divisors of n (d such that d divides n, gcd(d, n/d) = 1).

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 2, 2, 4, 2, 4, 2, 4, 4, 2, 2, 4, 2, 4, 4, 4, 2, 4, 2, 4, 2, 4, 2, 8, 2, 2, 4, 4, 4, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 4, 4, 4, 4, 4, 2, 8, 2, 4, 4, 2, 4, 8, 2, 4, 4, 8, 2, 4, 2, 4, 4, 4, 4, 8, 2, 4, 2, 4, 2, 8, 4, 4, 4, 4, 2, 8, 4, 4, 4, 4, 4, 4, 2, 4, 4, 4, 2, 8, 2, 4, 8
Offset: 1

Views

Author

Keywords

Comments

If n = Product p_i^a_i, d = Product p_i^c_i is a unitary divisor of n if each c_i is 0 or a_i.
Also the number of squarefree divisors of n. - Labos Elemer
Also number of divisors of the squarefree kernel of n: a(n) = A000005(A007947(n)). - Reinhard Zumkeller, Jul 19 2002
Also shadow transform of pronic numbers A002378.
For n >= 1 define an n X n (0,1) matrix A by A[i,j] = 1 if lcm(i,j) = n, A[i,j] = 0 if lcm(i,j) <> n for 1 <= i,j <= n. a(n) is the rank of A. - Yuval Dekel (dekelyuval(AT)hotmail.com), Aug 11 2003
a(n) is also the number of solutions to x^2 - x == 0 (mod n). - Yuval Dekel (dekelyuval(AT)hotmail.com), Sep 21 2003
a(n) is the number of squarefree divisors of n, but in general the set of unitary divisors of n is not the set of squarefree divisors (compare the rows of A077610 and A206778). - Jaroslav Krizek, May 04 2009
Row lengths of the triangles in A077610 and in A206778. - Reinhard Zumkeller, Feb 12 2012
a(n) is also the number of distinct residues of k^phi(n) (mod n), k=0..n-1. - Michel Lagneau, Nov 15 2012
a(n) is the number of irreducible fractions y/x that satisfy x*y=n (and gcd(x,y)=1), x and y positive integers. - Luc Rousseau, Jul 09 2017
a(n) is the number of (x,y) lattice points satisfying both x*y=n and (x,y) is visible from (0,0), x and y positive integers. - Luc Rousseau, Jul 10 2017
Conjecture: For any nonnegative integer k and positive integer n, the sum of the k-th powers of the unitary divisors of n is divisible by the sum of the k-th powers of the odd unitary divisors of n (note that this sequence lists the sum of the 0th powers of the unitary divisors of n). - Ivan N. Ianakiev, Feb 18 2018
a(n) is the number of one-digit numbers, k, when written in base n such that k and k^2 end in the same digit. - Matthew Scroggs, Jun 01 2018
Dirichlet convolution of A271102 and A000005. - Vaclav Kotesovec, Apr 08 2019
Conjecture: Let b(i; n), n > 0, be multiplicative sequences for some fixed integer i >= 0 with b(i; p^e) = (Sum_{k=1..i+1} A164652(i, k) * e^(k-1)) * (i+2) / (i!) for prime p and e > 0. Then we have Dirichlet generating functions: Sum_{n > 0} b(i; n) / n^s = (zeta(s))^(i+2) / zeta((i+2) * s). Examples for i=0 this sequence, for i=1 A226602, and for i=2 A286779. - Werner Schulte, Feb 17 2022
The smallest integer with 2^m unitary divisors, or equivalently, the smallest integer with 2^m squarefree divisors, is A002110(m). - Bernard Schott, Oct 04 2022

Examples

			a(12) = 4 because the four unitary divisors of 12 are 1, 3, 4, 12, and also because the four squarefree divisors of 12 are 1, 2, 3, 6.
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, Sect. B3.

Crossrefs

Sum of the k-th powers of the squarefree divisors of n for k=0..10: this sequence (k=0), A048250 (k=1), A351265 (k=2), A351266 (k=3), A351267 (k=4), A351268 (k=5), A351269 (k=6), A351270 (k=7), A351271 (k=8), A351272 (k=9), A351273 (k=10).
Sequences of the form n^k * Product_ {p|n, p prime} (1 + 1/p^k) for k=0..10: this sequence (k=0), A001615 (k=1), A065958 (k=2), A065959 (k=3), A065960 (k=4), A351300 (k=5), A351301 (k=6), A351302 (k=7), A351303 (k=8), A351304 (k=9), this sequence (k=10).
Cf. A020821 (Dgf at s=2), A177057 (Dgf at s=4).

Programs

  • Haskell
    a034444 = length . a077610_row  -- Reinhard Zumkeller, Feb 12 2012
    
  • Magma
    [#[d:d in Divisors(n)|Gcd(d,n div d) eq 1]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020
    
  • Magma
    [&+[Abs(MoebiusMu(d)):d in Divisors(n)]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020
  • Maple
    with(numtheory): for n from 1 to 200 do printf(`%d,`,2^nops(ifactors(n)[2])) od:
    with(numtheory);
    # returns the number of unitary divisors of n and a list of them
    f:=proc(n)
    local ct,i,t1,ans;
    ct:=0; ans:=[];
    t1:=divisors(n);
    for i from 1 to nops(t1) do
    d:=t1[i];
    if igcd(d,n/d)=1 then ct:=ct+1; ans:=[op(ans),d]; fi;
    od:
    RETURN([ct,ans]);
    end;
    # N. J. A. Sloane, May 01 2013
    # alternative Maple program:
    a:= n-> 2^nops(ifactors(n)[2]):
    seq(a(n), n=1..105);  # Alois P. Heinz, Jan 23 2024
    a := n -> 2^NumberTheory:-NumberOfPrimeFactors(n, distinct):  # Peter Luschny, May 13 2025
  • Mathematica
    a[n_] := Count[Divisors[n], d_ /; GCD[d, n/d] == 1]; a /@ Range[105] (* Jean-François Alcover, Apr 05 2011 *)
    Table[2^PrimeNu[n],{n,110}] (* Harvey P. Dale, Jul 14 2011 *)
  • PARI
    a(n)=1<Charles R Greathouse IV, Feb 11 2011
    
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (1+X)/(1-X))[n], ", ")) \\ Vaclav Kotesovec, Sep 26 2020
    
  • Python
    from sympy import divisors, gcd
    def a(n):
        return sum(1 for d in divisors(n) if gcd(d, n//d)==1)
    # Indranil Ghosh, Apr 16 2017
    
  • Python
    from sympy import primefactors
    def a(n): return 2**len(primefactors(n))
    print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Apr 16 2017
    
  • Scheme
    (define (A034444 n) (if (= 1 n) n (* 2 (A034444 (A028234 n))))) ;; Antti Karttunen, May 29 2017
    

Formula

a(n) = Sum_{d|n} abs(mu(n)) = 2^(number of different primes dividing n) = 2^A001221(n), with mu(n) = A008683(n). [Added Möbius formula. - Wolfdieter Lang, Jan 11 2020]
a(n) = Product_{ primes p|n } (1 + Legendre(1, p)).
Multiplicative with a(p^k)=2 for p prime and k>0. - Henry Bottomley, Oct 25 2001
a(n) = Sum_{d|n} tau(d^2)*mu(n/d), Dirichlet convolution of A048691 and A008683. - Benoit Cloitre, Oct 03 2002
Dirichlet generating function: zeta(s)^2/zeta(2s). - Franklin T. Adams-Watters, Sep 11 2005
Inverse Mobius transform of A008966. - Franklin T. Adams-Watters, Sep 11 2005
Asymptotically [Finch] the cumulative sum of a(n) = Sum_{n=1..N} a(n) ~ (6/(Pi^2))*N*log(N) + (6/(Pi^2))*(2*gamma - 1 - (12/(Pi^2))*zeta'(2))*N + O(sqrt(N)). - Jonathan Vos Post, May 08 2005 [typo corrected by Vaclav Kotesovec, Sep 13 2018]
a(n) = Sum_{d|n} floor(rad(d)/d), where rad is A007947 and floor(rad(n)/n) = A008966(n). - Enrique Pérez Herrero, Nov 13 2009
a(n) = A000005(n) - A048105(n); number of nonzero terms in row n of table A225817. - Reinhard Zumkeller, Jul 30 2013
G.f.: Sum_{n>0} A008966(n)*x^n/(1-x^n). - Mircea Merca, Feb 25 2014
a(n) = Sum_{d|n} lambda(d)*mu(d), where lambda is A008836. - Enrique Pérez Herrero, Apr 27 2014
a(n) = A277561(A156552(n)). - Antti Karttunen, May 29 2017
a(n) = A005361(n^2)/A005361(n). - Velin Yanev, Jul 26 2017
L.g.f.: -log(Product_{k>=1} (1 - mu(k)^2*x^k)^(1/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Jul 30 2018
a(n) = Sum_{d|n} A001615(d) * A023900(n/d). - Torlach Rush, Jan 20 2020
Sum_{d|n, gcd(d, n/d) = 1} a(d) * (-1)^omega(n/d) = 1. - Amiram Eldar, May 29 2020
a(n) = lim_{k->oo} A000005(n^(2*k))/A000005(n^k). - Velin Yanev and Amiram Eldar, Jan 10 2025

Extensions

More terms from James Sellers, Jun 20 2000

A156552 Unary-encoded compressed factorization of natural numbers.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 8, 7, 6, 9, 16, 11, 32, 17, 10, 15, 64, 13, 128, 19, 18, 33, 256, 23, 12, 65, 14, 35, 512, 21, 1024, 31, 34, 129, 20, 27, 2048, 257, 66, 39, 4096, 37, 8192, 67, 22, 513, 16384, 47, 24, 25, 130, 131, 32768, 29, 36, 71, 258, 1025, 65536, 43, 131072, 2049, 38, 63, 68, 69, 262144
Offset: 1

Views

Author

Leonid Broukhis, Feb 09 2009

Keywords

Comments

The primes become the powers of 2 (2 -> 1, 3 -> 2, 5 -> 4, 7 -> 8); the composite numbers are formed by taking the values for the factors in the increasing order, multiplying them by the consecutive powers of 2, and summing. See the Example section.
From Antti Karttunen, Jun 27 2014: (Start)
The odd bisection (containing even terms) halved gives A244153.
The even bisection (containing odd terms), when one is subtracted from each and halved, gives this sequence back.
(End)
Question: Are there any other solutions that would satisfy the recurrence r(1) = 0; and for n > 1, r(n) = Sum_{d|n, d>1} 2^A033265(r(d)), apart from simple variants 2^k * A156552(n)? See also A297112, A297113. - Antti Karttunen, Dec 30 2017

Examples

			For 84 = 2*2*3*7 -> 1*1 + 1*2 + 2*4 + 8*8 =  75.
For 105 = 3*5*7 -> 2*1 + 4*2 + 8*4 = 42.
For 137 = p_33 -> 2^32 = 4294967296.
For 420 = 2*2*3*5*7 -> 1*1 + 1*2 + 2*4 + 4*8 + 8*16 = 171.
For 147 = 3*7*7 = p_2 * p_4 * p_4 -> 2*1 + 8*2 + 8*4 = 50.
		

Crossrefs

One less than A005941.
Inverse permutation: A005940 with starting offset 0 instead of 1.
Cf. also A297106, A297112 (Möbius transform), A297113, A153013, A290308, A300827, A323243, A323244, A323247, A324201, A324812 (n for which a(n) is a square), A324813, A324822, A324823, A324398, A324713, A324815, A324819, A324865, A324866, A324867.

Programs

  • Mathematica
    Table[Floor@ Total@ Flatten@ MapIndexed[#1 2^(#2 - 1) &, Flatten[ Table[2^(PrimePi@ #1 - 1), {#2}] & @@@ FactorInteger@ n]], {n, 67}] (* Michael De Vlieger, Sep 08 2016 *)
  • PARI
    a(n) = {my(f = factor(n), p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ David A. Corneth, Mar 08 2019
    
  • PARI
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n)))); \\ (based on the given recurrence) - Antti Karttunen, Mar 08 2019
    
  • Perl
    # Program corrected per instructions from Leonid Broukhis. - Antti Karttunen, Jun 26 2014
    # However, it gives correct answers only up to n=136, before corruption by a wrap-around effect.
    # Note that the correct answer for n=137 is A156552(137) = 4294967296.
    $max = $ARGV[0];
    $pow = 0;
    foreach $i (2..$max) {
    @a = split(/ /, `factor $i`);
    shift @a;
    $shift = 0;
    $cur = 0;
    while ($n = int shift @a) {
    $prime{$n} = 1 << $pow++ if !defined($prime{$n});
    $cur |= $prime{$n} << $shift++;
    }
    print "$cur, ";
    }
    print "\n";
    (Scheme, with memoization-macro definec from Antti Karttunen's IntSeq-library, two different implementations)
    (definec (A156552 n) (cond ((= n 1) 0) (else (+ (A000079 (+ -2 (A001222 n) (A061395 n))) (A156552 (A052126 n))))))
    (definec (A156552 n) (cond ((= 1 n) (- n 1)) ((even? n) (+ 1 (* 2 (A156552 (/ n 2))))) (else (* 2 (A156552 (A064989 n))))))
    ;; Antti Karttunen, Jun 26 2014
    
  • Python
    from sympy import primepi, factorint
    def A156552(n): return sum((1<Chai Wah Wu, Mar 10 2023

Formula

From Antti Karttunen, Jun 26 2014: (Start)
a(1) = 0, a(n) = A000079(A001222(n)+A061395(n)-2) + a(A052126(n)).
a(1) = 0, a(2n) = 1+2*a(n), a(2n+1) = 2*a(A064989(2n+1)). [Compare to the entanglement recurrence A243071].
For n >= 0, a(2n+1) = 2*A244153(n+1). [Follows from the latter clause of the above formula.]
a(n) = A005941(n) - 1.
As a composition of related permutations:
a(n) = A003188(A243354(n)).
a(n) = A054429(A243071(n)).
For all n >= 1, A005940(1+a(n)) = n and for all n >= 0, a(A005940(n+1)) = n. [The offset-0 version of A005940 works as an inverse for this permutation.]
This permutations also maps between the partition-lists A112798 and A125106:
A056239(n) = A161511(a(n)). [The sums of parts of each partition (the total sizes).]
A003963(n) = A243499(a(n)). [And also the products of those parts.]
(End)
From Antti Karttunen, Oct 09 2016: (Start)
A161511(a(n)) = A056239(n).
A029837(1+a(n)) = A252464(n). [Binary width of terms.]
A080791(a(n)) = A252735(n). [Number of nonleading 0-bits.]
A000120(a(n)) = A001222(n). [Binary weight.]
For all n >= 2, A001511(a(n)) = A055396(n).
For all n >= 2, A000120(a(n))-1 = A252736(n). [Binary weight minus one.]
A252750(a(n)) = A252748(n).
a(A250246(n)) = A252754(n).
a(A005117(n)) = A277010(n). [Maps squarefree numbers to a permutation of A003714, fibbinary numbers.]
A085357(a(n)) = A008966(n). [Ditto for their characteristic functions.]
For all n >= 0:
a(A276076(n)) = A277012(n).
a(A276086(n)) = A277022(n).
a(A260443(n)) = A277020(n).
(End)
From Antti Karttunen, Dec 30 2017: (Start)
For n > 1, a(n) = Sum_{d|n, d>1} 2^A033265(a(d)). [See comments.]
More linking formulas:
A106737(a(n)) = A000005(n).
A290077(a(n)) = A000010(n).
A069010(a(n)) = A001221(n).
A136277(a(n)) = A181591(n).
A132971(a(n)) = A008683(n).
A106400(a(n)) = A008836(n).
A268411(a(n)) = A092248(n).
A037011(a(n)) = A010052(n) [conjectured, depends on the exact definition of A037011].
A278161(a(n)) = A046951(n).
A001316(a(n)) = A061142(n).
A277561(a(n)) = A034444(n).
A286575(a(n)) = A037445(n).
A246029(a(n)) = A181819(n).
A278159(a(n)) = A124859(n).
A246660(a(n)) = A112624(n).
A246596(a(n)) = A069739(n).
A295896(a(n)) = A053866(n).
A295875(a(n)) = A295297(n).
A284569(a(n)) = A072411(n).
A286574(a(n)) = A064547(n).
A048735(a(n)) = A292380(n).
A292272(a(n)) = A292382(n).
A244154(a(n)) = A048673(n), a(A064216(n)) = A244153(n).
A279344(a(n)) = A279339(n), a(A279338(n)) = A279343(n).
a(A277324(n)) = A277189(n).
A037800(a(n)) = A297155(n).
For n > 1, A033265(a(n)) = 1+A297113(n).
(End)
From Antti Karttunen, Mar 08 2019: (Start)
a(n) = A048675(n) + A323905(n).
a(A324201(n)) = A000396(n), provided there are no odd perfect numbers.
The following sequences are derived from or related to the base-2 expansion of a(n):
A000265(a(n)) = A322993(n).
A002487(a(n)) = A323902(n).
A005187(a(n)) = A323247(n).
A324288(a(n)) = A324116(n).
A323505(a(n)) = A323508(n).
A079559(a(n)) = A323512(n).
A085405(a(n)) = A323239(n).
The following sequences are obtained by applying to a(n) a function that depends on the prime factorization of its argument, which goes "against the grain" because a(n) is the binary code of the factorization of n, which in these cases is then factored again:
A000203(a(n)) = A323243(n).
A033879(a(n)) = A323244(n) = 2*a(n) - A323243(n),
A294898(a(n)) = A323248(n).
A000005(a(n)) = A324105(n).
A000010(a(n)) = A324104(n).
A083254(a(n)) = A324103(n).
A001227(a(n)) = A324117(n).
A000593(a(n)) = A324118(n).
A001221(a(n)) = A324119(n).
A009194(a(n)) = A324396(n).
A318458(a(n)) = A324398(n).
A192895(a(n)) = A324100(n).
A106315(a(n)) = A324051(n).
A010052(a(n)) = A324822(n).
A053866(a(n)) = A324823(n).
A001065(a(n)) = A324865(n) = A323243(n) - a(n),
A318456(a(n)) = A324866(n) = A324865(n) OR a(n),
A318457(a(n)) = A324867(n) = A324865(n) XOR a(n),
A318458(a(n)) = A324398(n) = A324865(n) AND a(n),
A318466(a(n)) = A324819(n) = A323243(n) OR 2*a(n),
A318467(a(n)) = A324713(n) = A323243(n) XOR 2*a(n),
A318468(a(n)) = A324815(n) = A323243(n) AND 2*a(n).
(End)

Extensions

More terms from Antti Karttunen, Jun 28 2014

A278219 Filter-sequence related to base-2 run-length encoding: a(n) = A046523(A243353(n)).

Original entry on oeis.org

1, 2, 4, 2, 4, 8, 6, 2, 4, 12, 16, 8, 6, 12, 6, 2, 4, 12, 36, 12, 16, 32, 24, 8, 6, 30, 24, 12, 6, 12, 6, 2, 4, 12, 36, 12, 36, 72, 60, 12, 16, 48, 64, 32, 24, 72, 24, 8, 6, 30, 60, 30, 24, 48, 60, 12, 6, 30, 24, 12, 6, 12, 6, 2, 4, 12, 36, 12, 36, 72, 60, 12, 36, 180, 144, 72, 60, 180, 60, 12, 16, 48, 144, 48, 64, 128, 96, 32, 24, 120, 216, 72, 24, 72
Offset: 0

Views

Author

Antti Karttunen, Nov 16 2016

Keywords

Crossrefs

Other base-2 related filter sequences: A278217, A278222.
Sequences that (seem to) partition N into same or coarser equivalence classes are at least these: A005811, A136004, A033264, A037800, A069010, A087116, A090079 and many others like A105500, A106826, A166242, A246960, A277561, A037834, A225081 although these have not been fully checked yet.

Programs

  • Mathematica
    f[n_, i_, x_] := Which[n == 0, x, EvenQ@ n, f[n/2, i + 1, x], True, f[(n - 1)/2, i, x Prime@ i]]; g[n_] := If[n == 1, 1, Times @@ MapIndexed[ Prime[First@ #2]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]]];
    Table[g@ f[BitXor[n, Floor[n/2]], 1, 1], {n, 0, 93}] (* Michael De Vlieger, May 09 2017 *)
  • Python
    from sympy import prime, factorint
    import math
    def A(n): return n - 2**int(math.floor(math.log(n, 2)))
    def b(n): return n + 1 if n<2 else prime(1 + (len(bin(n)[2:]) - bin(n)[2:].count("1"))) * b(A(n))
    def a005940(n): return b(n - 1)
    def P(n):
        f = factorint(n)
        return sorted([f[i] for i in f])
    def a046523(n):
        x=1
        while True:
            if P(n) == P(x): return x
            else: x+=1
    def a003188(n): return n^int(n/2)
    def a243353(n): return a005940(1 + a003188(n))
    def a(n): return a046523(a243353(n)) # Indranil Ghosh, May 07 2017
  • Scheme
    (define (A278219 n) (A046523 (A243353 n)))
    

Formula

a(n) = A046523(A243353(n)).
a(n) = A278222(A003188(n)).
a(n) = A278220(1+A075157(n)).

A278243 Filter-sequence for Stern polynomials: Least number with the same prime signature as A260443(n).

Original entry on oeis.org

1, 2, 2, 6, 2, 12, 6, 30, 2, 60, 12, 120, 6, 180, 30, 210, 2, 420, 60, 1080, 12, 2160, 120, 2520, 6, 2520, 180, 7560, 30, 6300, 210, 2310, 2, 4620, 420, 37800, 60, 90720, 1080, 75600, 12, 226800, 2160, 544320, 120, 453600, 2520, 138600, 6, 138600, 2520, 756000, 180, 2268000, 7560, 831600, 30, 415800, 6300, 2079000, 210, 485100, 2310, 30030, 2, 60060, 4620
Offset: 0

Views

Author

Antti Karttunen, Nov 16 2016

Keywords

Comments

This sequence can be used for filtering certain Stern polynomial (see A125184, A260443) related sequences, because it matches only with any such sequence b that can be computed as b(n) = f(A260443(n)), where f(n) is any function that depends only on the prime signature of n (some of these are listed under the index entry for "sequences computed from exponents in ...").
Matching in this context means that the sequence a matches with the sequence b iff for all i, j: a(i) = a(j) => b(i) = b(j). In other words, iff the sequence b partitions the natural numbers to the same or coarser equivalence classes (as/than the sequence a) by the distinct values it obtains.
Some of these are listed on the last line ("Sequences that partition N into ...") of Crossrefs section.

Crossrefs

Sequences that partition or seem to partition N into same or coarser equivalence classes: A002487, A126606, A277314, A277315, A277325, A277326, A277700, A277705.
The following are less certain: A007302 (not proved, but the first 10000 terms match), A072453, A110955 (uncertain, but related to A007302), A218799, A218800.
Note that the base-2 related sequences A069010 and A277561 (= 2^A069010(n)) do not match, although at first it seems so, up to for at least 139 initial terms. Also A028928 belongs to a different family.

Programs

  • Mathematica
    a[n_] := a[n] = Which[n < 2, n + 1, EvenQ@ n, Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1] &@ a[n/2], True, a[#] a[# + 1] &[(n - 1)/2]]; Table[Times @@ MapIndexed[Prime[First@ #2]^#1 &, Sort[FactorInteger[#][[All, -1]], Greater]] - Boole[# == 1] &@ a@ n, {n, 0, 66}] (* Michael De Vlieger, May 12 2017 *)
  • Scheme
    (define (A278243 n) (A046523 (A260443 n)))

Formula

a(n) = A046523(A260443(n)).

A374356 a(n) is the greatest fibbinary number f <= n such that n - f is also a fibbinary number whose binary expansion has no common 1's with that of f (where fibbinary numbers correspond to A003714).

Original entry on oeis.org

0, 1, 2, 2, 4, 5, 4, 5, 8, 9, 10, 10, 8, 9, 10, 10, 16, 17, 18, 18, 20, 21, 20, 21, 16, 17, 18, 18, 20, 21, 20, 21, 32, 33, 34, 34, 36, 37, 36, 37, 40, 41, 42, 42, 40, 41, 42, 42, 32, 33, 34, 34, 36, 37, 36, 37, 40, 41, 42, 42, 40, 41, 42, 42, 64, 65, 66, 66
Offset: 0

Views

Author

Rémy Sigrist, Jul 06 2024

Keywords

Comments

To compute a(n): replace every other bit with zero (starting with the second bit) in each run of consecutive 1's in the binary expansion of n.
From Gus Wiseman, Jul 11 2025: (Start)
This is the greatest binary rank of a sparse subset of the binary indices of n, where:
1. The binary indices of a nonnegative integer are the positions of 1 in its reversed binary expansion.
2. A set is sparse iff 1 is not a first difference.
3. The binary rank of a set {S_1,S_2,...} is Sum_i 2^(S_i-1).
(End)

Examples

			The first terms, in decimal and in binary, are:
  n   a(n)  bin(n)  bin(a(n))
  --  ----  ------  ---------
   0     0       0          0
   1     1       1          1
   2     2      10         10
   3     2      11         10
   4     4     100        100
   5     5     101        101
   6     4     110        100
   7     5     111        101
   8     8    1000       1000
   9     9    1001       1001
  10    10    1010       1010
  11    10    1011       1010
  12     8    1100       1000
  13     9    1101       1001
  14    10    1110       1010
  15    10    1111       1010
  16    16   10000      10000
		

Crossrefs

The union is A003714 (Fibbinary numbers).
For prime instead of binary indices we have A385216.
A034839 counts subsets by number of maximal runs, for strict partitions A116674.
A166469 counts sparse submultisets of prime indices, maximal A385215.
A245564 counts sparse subsets of binary indices, maximal case A384883.
A319630 ranks sparse submultisets of prime indices, complement A104210.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    fbi[q_]:=If[q=={},0,Total[2^q]/2];
    Table[Max@@fbi/@Select[Subsets[bpe[n]],FreeQ[Differences[#],1]&],{n,0,100}] (* Gus Wiseman, Jul 11 2025 *)
  • PARI
    a(n) = { my (v = 0, e, x, y, b); while (n, x = y = 0; e = valuation(n, 2); for (k = 0, oo, if (bittest(n, e+k), n -= b = 2^(e+k); [x, y] = [y + b, x], v += x; break;););); return (v); }

Formula

a(n) = A374354(n, A277561(n)-1).
a(n) = n - A374355(n).
a(n) <= n with equality iff n is a fibbinary number.

A356397 a(n) is the product of the terms in the n-th row of triangle A343835; a(0) = 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 6, 7, 8, 8, 16, 24, 12, 12, 14, 15, 16, 16, 32, 48, 64, 64, 96, 112, 24, 24, 48, 72, 28, 28, 30, 31, 32, 32, 64, 96, 128, 128, 192, 224, 256, 256, 512, 768, 384, 384, 448, 480, 48, 48, 96, 144, 192, 192, 288, 336, 56, 56, 112, 168, 60, 60, 62
Offset: 0

Views

Author

Rémy Sigrist, Aug 05 2022

Keywords

Examples

			For n = 11:
- row 11 of A343835 is (8, 3),
- so a(11) = 8 * 3 = 24.
		

Crossrefs

Programs

  • PARI
    a(n) = { my (v=1); while (n, my (z=valuation(n, 2), o=valuation(n/2^z+1, 2), r=(2^o-1)*2^z); n-=r; v*=r); v }

Formula

a(2*n) = a(n) * A277561(n).
a(n) = n iff n belongs to A023758 \ {0}.
A246674(n) divides a(n).
Showing 1-8 of 8 results.