A264977 -id:A264977 - cp's OEIS frontend

A277826 a(n) = the least k for which A264977(k) = A264977(n).

0, 1, 2, 3, 4, 1, 6, 7, 8, 9, 2, 7, 12, 1, 14, 15, 16, 17, 18, 7, 4, 9, 14, 23, 24, 17, 2, 15, 28, 1, 30, 31, 32, 33, 34, 7, 36, 17, 14, 3, 8, 1, 18, 23, 28, 9, 46, 47, 48, 49, 34, 15, 4, 17, 30, 47, 56, 33, 2, 31, 60, 1, 62, 63, 64, 65, 66, 7, 68, 33, 14, 47, 72, 73, 34, 3, 28, 17, 6, 23, 16, 81, 2, 23, 36, 1, 46, 87, 56, 73, 18, 47, 92, 9

Offset: 0

Antti Karttunen, Nov 06 2016

Antti Karttunen, Table of n, a(n) for n = 0..65537

Cf. A264977, A277816.
Cf. A277701, A277712, A277713, A277715 (positions of 1, 2, 3 and 9 in this sequence).
Cf. A277824, A277884 and their scatter-plots.

(define (A277826 n) (let ((v (A264977 n))) (let loop ((k 0)) (if (= v (A264977 k)) k (loop (+ 1 k))))))

A283975 Odd bisection of A264977.

Original entry on oeis.org

1, 3, 1, 7, 5, 7, 1, 15, 13, 7, 5, 11, 13, 15, 1, 31, 29, 7, 13, 3, 1, 11, 5, 19, 21, 15, 13, 19, 29, 31, 1, 63, 61, 7, 29, 19, 25, 3, 13, 11, 9, 11, 1, 23, 25, 19, 5, 35, 37, 15, 21, 11, 9, 19, 13, 43, 37, 31, 29, 35, 61, 63, 1, 127, 125, 7, 61, 51, 41, 19, 29, 59, 49, 3, 25, 31, 17, 11, 13, 27, 25, 11, 9, 31, 21, 23, 1, 47, 33, 19, 25, 63, 41, 35, 5, 67, 69

Offset: 0

Antti Karttunen, Mar 25 2017

Antti Karttunen, Table of n, a(n) for n = 0..8190
Index entries for sequences related to binary expansion of n

Cf. A000120, A248663, A264977, A277324.

Cf. also A284266.

A264977[n_]:= If[n<2, n, If[EvenQ[n], 2 A264977[n/2], BitXor[A264977[(n - 1)/2], A264977[(n + 1)/2]]]]; Table[A264977[2n + 1], {n, 0, 150}] (* Indranil Ghosh, Mar 28 2017 *)

a(n) = if(n<2, n, if(n%2, bitxor(a((n - 1)/2), a((n + 1)/2)), 2*a(n/2))); for(n=0, 150, print1(a(2*n + 1),", ")) \\ Indranil Ghosh, Mar 28 2017

(define (A283975 n) (A264977 (+ n n 1)))

a(n) = A264977((2*n)+1).
a(n) = A248663(A277324(n)).
A000120(a(n)) = A284266(n).

A283979 a(n) = (n XOR A264977(n))/4, where XOR is bitwise-xor (A003987).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 3, 2, 3, 0, 3, 0, 0, 0, 7, 6, 5, 4, 4, 6, 7, 0, 5, 6, 5, 0, 7, 0, 0, 0, 15, 14, 9, 12, 10, 10, 9, 8, 10, 8, 8, 12, 10, 14, 15, 0, 9, 10, 15, 12, 14, 10, 9, 0, 9, 14, 9, 0, 15, 0, 0, 0, 31, 30, 17, 28, 22, 18, 21, 24, 20, 20, 18, 20, 16, 18, 17, 16, 22, 20, 22, 16, 21, 16, 16, 24, 16, 20, 18, 28, 22, 30, 31, 0, 17, 18, 27, 20, 28, 30

Offset: 0

Antti Karttunen, Mar 25 2017

Antti Karttunen, Table of n, a(n) for n = 0..8191
Index entries for sequences related to binary expansion of n

Cf. A003987, A264977, A265397, A284013.

A264977[n_]:= If[n<2, n, If[EvenQ[n], 2 A264977[n/2], BitXor[A264977[(n - 1)/2], A264977[(n + 1)/2]]]]; Table[BitXor[n, A264977[n]]/4, {n, 0, 100}] (* Indranil Ghosh, Mar 28 2017 *)

a(n) = if(n<2, n, if(n%2, bitxor(a((n - 1)/2), a((n + 1)/2)), 2*a(n/2)));
for(n=0, 150, print1(bitxor(n, a(n))/4,", ")) \\ Indranil Ghosh, Mar 28 2017

(define (A283979 n) (/ (A003987bi n (A264977 n)) 4)) ;; Where A003987bi implements bitwise-xor (A003987).

a(n) = (n XOR A264977(n))/4, where XOR is bitwise-xor (A003987).

A260443 Prime factorization representation of Stern polynomials: a(0) = 1, a(1) = 2, a(2n) = A003961(a(n)), a(2n+1) = a(n)*a(n+1).

Original entry on oeis.org

1, 2, 3, 6, 5, 18, 15, 30, 7, 90, 75, 270, 35, 450, 105, 210, 11, 630, 525, 6750, 245, 20250, 2625, 9450, 77, 15750, 3675, 47250, 385, 22050, 1155, 2310, 13, 6930, 5775, 330750, 2695, 3543750, 128625, 1653750, 847, 4961250, 643125, 53156250, 18865, 24806250, 202125, 727650, 143, 1212750, 282975, 57881250, 29645, 173643750, 1414875, 18191250, 1001

Offset: 0

Antti Karttunen, Jul 28 2015

The exponents in the prime factorization of term a(n) give the coefficients of the n-th Stern polynomial. See A125184 and the examples.
None of the terms have prime gaps in their factorization, i.e., all can be found in A073491.
Contains neither perfect squares nor prime powers with exponent > 1. A277701 gives the positions of the terms that are 2*square. - Antti Karttunen, Oct 27 2016
Many of the derived sequences (like A002487) have similar "Fir forest" or "Gaudian cathedrals" style scatter plot. - Antti Karttunen, Mar 21 2017

			n    a(n)   prime factorization    Stern polynomial
------------------------------------------------------------
0       1   (empty)                B_0(x) = 0
1       2   p_1                    B_1(x) = 1
2       3   p_2                    B_2(x) = x
3       6   p_2 * p_1              B_3(x) = x + 1
4       5   p_3                    B_4(x) = x^2
5      18   p_2^2 * p_1            B_5(x) = 2x + 1
6      15   p_3 * p_2              B_6(x) = x^2 + x
7      30   p_3 * p_2 * p_1        B_7(x) = x^2 + x + 1
8       7   p_4                    B_8(x) = x^3
9      90   p_3 * p_2^2 * p_1      B_9(x) = x^2 + 2x + 1

Alois P. Heinz, Table of n, a(n) for n = 0..4096 (first 1025 terms from Antti Karttunen)

Same sequence sorted into ascending order: A260442.
Cf. A000040, A000079, A000225, A001222, A002487, A003415, A003961, A005811, A007949, A046523, A056239, A073491, A090880, A097249, A101979, A125184, A178590, A186891, A206284, A277314, A277315, A277325, A277326, A277329, A277330, A277701, A277705, A277899, A278243, A278530, A278544, A284010, A284011.
Cf. also A048675, A277333 (left inverses).
Cf. A277323, A277324 (bisections), A277200 (even terms sorted), A277197 (first differences), A277198.
Cf. A277316 (values at primes), A277318.
Cf. A023758 (positions of squarefree terms), A101082 (of terms not squarefree), A277702 (positions of records), A277703 (their values).
Cf. A283992, A283993 (number of irreducible, reducible polynomials in range 1 .. n).
Cf. also A206296 (Fibonacci polynomials similarly represented).

b:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
a:= proc(n) option remember; `if`(n<2, n+1,
      `if`(irem(n, 2, 'h')=0, b(a(h)), a(h)*a(n-h)))
    end:
seq(a(n), n=0..56);  # Alois P. Heinz, Jul 04 2024

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[a@ n, {n, 0, 56}] (* Michael De Vlieger, Apr 05 2017 *)

A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From Michel Marcus
A260443(n) = if(n<2, n+1, if(n%2, A260443(n\2)*A260443(n\2+1), A003961(A260443(n\2)))); \\ After Charles R Greathouse IV's code for "ps" in A186891.
\\ Antti Karttunen, Oct 11 2016

from sympy import factorint, prime, primepi
from functools import reduce
from operator import mul
def a003961(n):
    F = factorint(n)
    return 1 if n==1 else reduce(mul, (prime(primepi(i) + 1)**F[i] for i in F))
def a(n): return n + 1 if n<2 else a003961(a(n//2)) if n%2==0 else a((n - 1)//2)*a((n + 1)//2)
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 21 2017

;; Uses memoization-macro definec:
(definec (A260443 n) (cond ((<= n 1) (+ 1 n)) ((even? n) (A003961 (A260443 (/ n 2)))) (else (* (A260443 (/ (- n 1) 2)) (A260443 (/ (+ n 1) 2))))))
;; A more standalone version added Oct 10 2016, requiring only an implementation of A000040 and the memoization-macro definec:
(define (A260443 n) (product_primes_to_kth_powers (A260443as_coeff_list n)))
(define (product_primes_to_kth_powers nums) (let loop ((p 1) (nums nums) (i 1)) (cond ((null? nums) p) (else (loop (* p (expt (A000040 i) (car nums))) (cdr nums) (+ 1 i))))))
(definec (A260443as_coeff_list n) (cond ((zero? n) (list)) ((= 1 n) (list 1)) ((even? n) (cons 0 (A260443as_coeff_list (/ n 2)))) (else (add_two_lists (A260443as_coeff_list (/ (- n 1) 2)) (A260443as_coeff_list (/ (+ n 1) 2))))))
(define (add_two_lists nums1 nums2) (let ((len1 (length nums1)) (len2 (length nums2))) (cond ((< len1 len2) (add_two_lists nums2 nums1)) (else (map + nums1 (append nums2 (make-list (- len1 len2) 0)))))))

a(0) = 1, a(1) = 2, a(2n) = A003961(a(n)), a(2n+1) = a(n)*a(n+1).
Other identities. For all n >= 0:
A001221(a(n)) = A277314(n). [#nonzero coefficients in each polynomial.]
A001222(a(n)) = A002487(n). [When each polynomial is evaluated at x=1.]
A048675(a(n)) = n.          [at x=2.]
A090880(a(n)) = A178590(n). [at x=3.]
A248663(a(n)) = A264977(n). [at x=2 over the field GF(2).]
A276075(a(n)) = A276081(n). ["at factorials".]
A156552(a(n)) = A277020(n). [Converted to "unary-binary" encoding.]
A051903(a(n)) = A277315(n). [Maximal coefficient.]
A277322(a(n)) = A277013(n). [Number of irreducible polynomial factors.]
A005361(a(n)) = A277325(n). [Product of nonzero coefficients.]
A072411(a(n)) = A277326(n). [And their LCM.]
A007913(a(n)) = A277330(n). [The squarefree part.]
A000005(a(n)) = A277705(n). [Number of divisors.]
A046523(a(n)) = A278243(n). [Filter-sequence.]
A284010(a(n)) = A284011(n). [True for n > 1. Another filter-sequence.]
A003415(a(n)) = A278544(n). [Arithmetic derivative.]
A056239(a(n)) = A278530(n). [Weighted sum of coefficients.]
A097249(a(n)) = A277899(n).
a(A000079(n)) = A000040(n+1).
a(A000225(n)) = A002110(n).
a(A000051(n)) = 3*A002110(n).
For n >= 1, a(A000918(n)) = A070826(n).
A007949(a(n)) is the interleaving of A000035 and A005811, probably A101979.
A061395(a(n)) = A277329(n).
Also, for all n >= 1:
A055396(a(n)) = A001511(n).
A252735(a(n)) = A061395(a(n)) - 1 = A057526(n).
a(A000040(n)) = A277316(n).
a(A186891(1+n)) = A277318(n). [Subsequence for irreducible polynomials].

More linking formulas added by Antti Karttunen, Mar 21 2017

A023758 Numbers of the form 2^i - 2^j with i >= j.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 12, 14, 15, 16, 24, 28, 30, 31, 32, 48, 56, 60, 62, 63, 64, 96, 112, 120, 124, 126, 127, 128, 192, 224, 240, 248, 252, 254, 255, 256, 384, 448, 480, 496, 504, 508, 510, 511, 512, 768, 896, 960, 992, 1008, 1016, 1020, 1022, 1023

Offset: 1

Olivier Gérard

Numbers whose digits in base 2 are in nonincreasing order.
Might be called "nialpdromes".
Subset of A077436. Proof: Since a(n) is of the form (2^i-1)*2^j, i,j >= 0, a(n)^2 = (2^(2i) - 2^(i+1))*2^(2j) + 2^(2j) where the first sum term has i-1 one bits and its 2j-th bit is zero, while the second sum term switches the 2j-th bit to one, giving i one bits, as in a(n). - Ralf Stephan, Mar 08 2004
Numbers whose binary representation contains no "01". - Benoit Cloitre, May 23 2004
Every polynomial with coefficients equal to 1 for the leading terms and 0 after that, evaluated at 2. For instance a(13) = x^4 + x^3 + x^2 at 2, a(14) = x^4 + x^3 + x^2 + x at 2. - Ben Paul Thurston, Jan 11 2008
From Gary W. Adamson, Jul 18 2008: (Start)
As a triangle by rows starting:
   1;
   2,  3;
   4,  6,  7;
   8, 12, 14, 15;
  16, 24, 28, 30, 31;
  ...,
equals A000012 * A130123 * A000012, where A130123 = (1, 0,2; 0,0,4; 0,0,0,8; ...). Row sums of this triangle = A000337 starting (1, 5, 17, 49, 129, ...). (End)
First differences are A057728 = 1; 1; 1; 1; 2,1; 1; 4,2,1; 1; 8,4,2,1; 1; ... i.e., decreasing powers of 2, separated by another "1". - M. F. Hasler, May 06 2009
Apart from first term, numbers that are powers of 2 or the sum of some consecutive powers of 2. - Omar E. Pol, Feb 14 2013
From Andres Cicuttin, Apr 29 2016: (Start)
Numbers that can be digitally generated with twisted ring (Johnson) counters. This is, the binary digits of a(n) correspond to those stored in a shift register where the input bit of the first bit storage element is the inverted output of the last storage element. After starting with all 0’s, each new state is obtained by rotating the stored bits but inverting at each state transition the last bit that goes to the first position (see link).
Examples: for a(n) represented by three bits
             Binary
a(5)= 4   ->  100   last bit = 0
a(6)= 6   ->  110   first bit = 1 (inverted last bit of previous number)
a(7)= 7   ->  111
and for a(n) represented by four bits
             Binary
a(8) = 8   -> 1000
a(9) = 12  -> 1100  last bit = 0
a(10)= 14  -> 1110  first bit = 1 (inverted last bit of previous number)
a(11)= 15  -> 1111
(End)
Powers of 2 represented in bases which are terms of this sequence must always contain at least one digit which is also a power of 2. This is because 2^i mod (2^i - 2^j) = 2^j, which means the last digit always cycles through powers of 2 (or if i=j+1 then the first digit is a power of 2 and the rest are trailing zeros). The only known non-member of this sequence with this property is 5. - Ely Golden, Sep 05 2017
Numbers k such that k = 2^(1 + A000523(k)) - 2^A007814(k). - Daniel Starodubtsev, Aug 05 2021
A002260(n) = v(a(n)/2^v(a(n))+1) and A002024(n) = A002260(n) + v(a(n)) where v is the dyadic valuation (i.e., A007814). - Lorenzo Sauras Altuzarra, Feb 01 2023

			a(22) = 64 = 32 + 32 = 2^5 + a(16) = 2^A003056(20) + a(22-5-1).
a(23) = 96 = 64 + 32 = 2^6 + a(16) = 2^A003056(21) + a(23-6-1).
a(24) = 112 = 64 + 48 = 2^6 + a(17) = 2^A003056(22) + a(24-6-1).

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000 (first 5051 terms from T. D. Noe)
Andreas M. Hinz and Paul K. Stockmeyer, Precious Metal Sequences and Sierpinski-Type Graphs, J. Integer Seq., Vol 25 (2022), Article 22.4.8.
S. M. Shabab Hossain, Md. Mahmudur Rahman and M. Sohel Rahman, Solving a Generalized Version of the Exact Cover Problem with a Light-Based Device, Optical Supercomputing, Lecture Notes in Computer Science, 2011, Volume 6748/2011, 23-31, DOI: 10.1007/978-3-642-22494-2_4.
Eric Weisstein's World of Mathematics, Digit.
Wikipedia, Ring counter.
Index entries for 2-automatic sequences.
Index entries for sequences related to binary expansion of n.

A000337(r) = sum of row T(r, c) with 0 <= c < r. See also A002024, A003056, A140129, A140130, A221975.
Cf. A007088, A130123, A101082 (complement), A340375 (characteristic function).
This is the base-2 version of A064222. First differences are A057728.
Subsequence of A077436, of A129523, of A277704, and of A333762.
Subsequences: A043569 (nonzero even terms, or equally, nonzero terms doubled), A175332, A272615, A335431, A000396 (its even terms only), A324200.
Positions of zeros in A049502, A265397, A277899, A284264.
Positions of ones in A283983, A283989.
Positions of nonzero terms in A341509 (apart from the initial zero).
Positions of squarefree terms in A260443.
Fixed points of A264977, A277711, A283165, A334666.
Distinct terms in A340632.
Cf. also A002024, A002260, A007814, A133797, A152449, A181666, A211705, A277699, A306441.
Cf. also A309758, A309759, A309761 (for analogous sequences).

import Data.Set (singleton, deleteFindMin, insert)
a023758 n = a023758_list !! (n-1)
a023758_list = 0 : f (singleton 1) where
f s = x : f (if even x then insert z s' else insert z $ insert (z+1) s')
where z = 2*x; (x, s') = deleteFindMin s
-- Reinhard Zumkeller, Sep 24 2014, Dec 19 2012

a:=proc(n) local n2,d: n2:=convert(n,base,2): d:={seq(n2[j]-n2[j-1],j=2..nops(n2))}: if n=0 then 0 elif n=1 then 1 elif d={0,1} or d={0} or d={1} then n else fi end: seq(a(n),n=0..2100); # Emeric Deutsch, Apr 22 2006

Union[Flatten[Table[2^i - 2^j, {i, 0, 100}, {j, 0, i}]]] (* T. D. Noe, Mar 15 2011 *)
Select[Range[0, 2^10], NoneTrue[Differences@ IntegerDigits[#, 2], # > 0 &] &] (* Michael De Vlieger, Sep 05 2017 *)

for(n=0,2500,if(prod(k=1,length(binary(n))-1,component(binary(n),k)+1-component(binary(n),k+1))>0,print1(n,",")))

A023758(n)= my(r=round(sqrt(2*n--))); (1<<(n-r*(r-1)/2)-1)<<(r*(r+1)/2-n)
/* or, to illustrate the "decreasing digit" property and analogy to A064222: */
A023758(n,show=0)={ my(a=0); while(n--, show & print1(a","); a=vecsort(binary(a+1)); a*=vector(#a,j,2^(j-1))~); a} \\ M. F. Hasler, May 06 2009

is(n)=if(n<5,1,n>>=valuation(n,2);n++;n>>valuation(n,2)==1) \\ Charles R Greathouse IV, Jan 04 2016

list(lim)=my(v=List([0]),t); for(i=1,logint(lim\1+1,2), t=2^i-1; while(t<=lim, listput(v,t); t*=2)); Set(v) \\ Charles R Greathouse IV, May 03 2016

def a_next(a_n): return (a_n | (a_n >> 1)) + (a_n & 1)
a_n = 1; a = [0]
for i in range(55): a.append(a_n); a_n = a_next(a_n) # Falk Hüffner, Feb 19 2022

from math import isqrt
def A023758(n): return (1<<(m:=isqrt(n-1<<3)+1>>1))-(1<<(m*(m+1)-(n-1<<1)>>1)) # Chai Wah Wu, Feb 23 2025

a(n) = 2^s(n) - 2^((s(n)^2 + s(n) - 2n)/2) where s(n) = ceiling((-1 + sqrt(1+8n))/2). - Sam Alexander, Jan 08 2005
a(n) = 2^k + a(n-k-1) for 1 < n and k = A003056(n-2). The rows of T(r, c) = 2^r-2^c for 0 <= c < r read from right to left produce this sequence: 1; 2, 3; 4, 6, 7; 8, 12, 14, 15; ... - Frank Ellermann, Dec 06 2001
For n > 0, a(n) mod 2 = A010054(n). - Benoit Cloitre, May 23 2004
A140130(a(n)) = 1 and for n > 1: A140129(a(n)) = A002262(n-2). - Reinhard Zumkeller, May 14 2008
a(n+1) = (2^(n - r(r-1)/2) - 1) 2^(r(r+1)/2 - n), where r=round(sqrt(2n)). - M. F. Hasler, May 06 2009
Start with A000225. If k is in the sequence, then so is 2k. - Ralf Stephan, Aug 16 2013
G.f.: (x^2/((2-x)*(1-x)))*(1 + Sum_{k>=0} x^((k^2+k)/2)*(1 + x*(2^k-1))). The sum is related to Jacobi theta functions. - Robert Israel, Feb 24 2015
A049502(a(n)) = 0. - Reinhard Zumkeller, Jun 17 2015
a(n) = a(n-1) + a(n-d)/a(d*(d+1)/2 + 2) if n > 1, d > 0, where d = A002262(n-2). - Yuchun Ji, May 11 2020
A277699(a(n)) = a(n)^2, A306441(a(n)) = a(n+1). - Antti Karttunen, Feb 15 2021 (the latter identity from A306441)
Sum_{n>=2} 1/a(n) = A211705. - Amiram Eldar, Feb 20 2022

Definition changed by N. J. A. Sloane, Jan 05 2008

A248663 Binary encoding of the prime factors of the squarefree part of n.

Original entry on oeis.org

0, 1, 2, 0, 4, 3, 8, 1, 0, 5, 16, 2, 32, 9, 6, 0, 64, 1, 128, 4, 10, 17, 256, 3, 0, 33, 2, 8, 512, 7, 1024, 1, 18, 65, 12, 0, 2048, 129, 34, 5, 4096, 11, 8192, 16, 4, 257, 16384, 2, 0, 1, 66, 32, 32768, 3, 20, 9, 130, 513, 65536, 6, 131072, 1025, 8, 0, 36, 19

Offset: 1

Peter Kagey, Jan 11 2015

The binary digits of a(n) encode the prime factorization of A007913(n), where the i-th digit from the right is 1 if and only if prime(i) divides A007913(n), otherwise 0. - Robert Israel, Jan 12 2015
Old name: a(1) = 0; a(A000040(n)) = 2^(n-1), and a(n*m) = a(n) XOR a(m).
XOR is the bitwise exclusive or operation (A003987).
a(k^2) = 0 for a natural number k.
Equivalently, the i-th binary digit from the right is 1 iff prime(i) divides n an odd number of times, otherwise zero. - Ethan Beihl, Oct 15 2016
When a polynomial with nonnegative integer coefficients is encoded with the prime factorization of n (e.g., as in A206296, A260443, with scheme explained in A206284), then A048675(n) gives the evaluation of that polynomial at x=2. This sequence is otherwise similar, except the polynomial is evaluated over the field GF(2), which implies also that all its coefficients are essentially reduced modulo 2. - Antti Karttunen, Dec 11 2015
Squarefree numbers (A005117) give the positions k where a(k) = A048675(k). - Antti Karttunen, Oct 29 2016
From Peter Munn, Jun 07 2021: (Start)
When we encode polynomials with nonnegative integer coefficients as described by Antti Karttunen above, polynomial addition is represented by integer multiplication, multiplication is represented by A297845(.,.), and this sequence represents a surjective semiring homomorphism to polynomials in GF(2)[x] (encoded as described in A048720). The mapping of addition operations by this homomorphism is part of the sequence definition: "a(n*m) = a(n) XOR a(m)". The mapping of multiplication is given by a(A297845(n, k)) = A048720(a(n), a(k)).
In a related way, A329329 defines a representation of a different set of polynomials as positive integers, namely polynomials in GF(2)[x,y].
Let P_n(x,y) denote the polynomial represented, as in A329329, by n >= 1. If 0 is substituted for y in P_n(x,y), we get a polynomial P'_n(x,y) (in which y does not appear, of course) that is equivalent to a polynomial P'_n(x) in GF(2)[x]. a(n) is the integer encoding of P'_n(x) (described in A048720).
Viewed as above, this sequence represents another surjective homomorphism, a homomorphism between polynomial rings, with A329329(.,.)/A059897(.,.) and A048720(.,.)/A003987(.,.) as the respective ring operations.
a(n) can be composed as a(n) = A048675(A007913(n)) and the effect of the A007913(.) component corresponds to different operations on the respective polynomial domains of the two homomorphisms described above. In the first homomorphism, coefficients are reduced modulo 2; in the second, 0 is substituted for y. This is illustrated in the examples.
(End)

			a(3500) = a(2^2 * 5^3 * 7) = a(2) XOR a(2) XOR a(5) XOR a(5) XOR a(5) XOR a(7) = 1 XOR 1 XOR 4 XOR 4 XOR 4 XOR 8 = 0b0100 XOR 0b1000 = 0b1100 = 12.
From _Peter Munn_, Jun 07 2021: (Start)
The examples in the table below illustrate the homomorphisms (between polynomial structures) represented by this sequence.
The staggering of the rows is to show how the mapping n -> A007913(n) -> A048675(A007913(n)) = a(n) relates to the encoded polynomials, as not all encodings are relevant at each stage.
For an explanation of each polynomial encoding, see the sequence referenced in the relevant column heading. (Note also that A007913 generates squarefree numbers, and with these encodings, all squarefree numbers represent equivalent polynomials in N[x] and GF(2)[x,y].)
                     |<-----    encoded polynomials    ----->|
  n  A007913(n) a(n) |         N[x]    GF(2)[x,y]    GF(2)[x]|
                     |Cf.:  A206284       A329329     A048720|
--------------------------------------------------------------
  24                            x+3         x+y+1
          6                     x+1           x+1
                  3                                       x+1
--------------------------------------------------------------
  36                           2x+2          xy+y
          1                       0             0
                  0                                         0
--------------------------------------------------------------
  60                        x^2+x+2       x^2+x+y
         15                   x^2+x         x^2+x
                  6                                     x^2+x
--------------------------------------------------------------
  90                       x^2+2x+1      x^2+xy+1
         10                   x^2+1         x^2+1
                  5                                     x^2+1
--------------------------------------------------------------
This sequence is a left inverse of A019565. A019565(.) maps a(n) to A007913(n) for all n, effectively reversing the second stage of the mapping from n to a(n) shown above. So, with the encodings used here, A019565(.) represents each of two injective homomorphisms that map polynomials in GF(2)[x] to equivalent polynomials in N[x] and GF(2)[x,y] respectively.
(End)

Peter Kagey, Table of n, a(n) for n = 1..5000
Eric Weisstein's World of Mathematics, Squarefree Part
Wikipedia, Polynomial ring
Index entries for sequences operating on GF(2)[X]-polynomials

Cf. A000040, A000079, A000720, A005117, A020639, A032742, A048672, A055396, A100112, A206284.
A048675 composed with A007913. A007814 composed with A225546.
A homomorphism from A297845/A003991 and from A329329/A059897 to A048720/A003987, mapping A206296 to A168081, A260443 to A264977, A265408 to A265407, A275734 to A275808, A276076 to A276074, A283477 to A006068.
A left inverse of A019565.
Other sequences used to express relationship between terms of this sequence: A003961, A007913, A331590, A334747.
Cf. also A099884, A277330.
A087207 is the analogous sequence with OR.
A277417 gives the positions where coincides with A277333.
A000290 gives the positions of zeros.

import Data.Bits (xor)
a248663 = foldr (xor) 0 . map (\i -> 2^(i - 1)) . a112798_row
-- Peter Kagey, Sep 16 2016

f:= proc(n)
local F,f;
F:= select(t -> t[2]::odd, ifactors(n)[2]);
add(2^(numtheory:-pi(f[1])-1), f = F)
end proc:
seq(f(i),i=1..100); # Robert Israel, Jan 12 2015

a[1] = 0; a[n_] := a[n] = If[PrimeQ@ n, 2^(PrimePi@ n - 1), BitXor[a[#], a[n/#]] &@ FactorInteger[n][[1, 1]]]; Array[a, 66] (* Michael De Vlieger, Sep 16 2016 *)

A248663(n) = vecsum(apply(p -> 2^(primepi(p)-1),factor(core(n))[,1])); \\ Antti Karttunen, Feb 15 2021

from sympy import factorint, primepi
from sympy.ntheory.factor_ import core
def a048675(n):
    f=factorint(n)
    return 0 if n==1 else sum([f[i]*2**(primepi(i) - 1) for i in f])
def a(n): return a048675(core(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 21 2017

require 'prime'
def f(n)
  a = 0
  reverse_primes = Prime.each(n).to_a.reverse
  reverse_primes.each do |prime|
    a <<= 1
    while n % prime == 0
      n /= prime
      a ^= 1
    end
  end
  a
end
(Scheme, with memoizing-macro definec)
(definec (A248663 n) (cond ((= 1 n) 0) ((= 1 (A010051 n)) (A000079 (- (A000720 n) 1))) (else (A003987bi (A248663 (A020639 n)) (A248663 (A032742 n)))))) ;; Where A003987bi computes bitwise-XOR as in A003987.
;; Alternatively:
(definec (A248663 n) (cond ((= 1 n) 0) (else (A003987bi (A000079 (- (A055396 n) 1)) (A248663 (A032742 n))))))
;; Antti Karttunen, Dec 11 2015

a(1) = 0; for n > 1, if n is a prime, a(n) = 2^(A000720(n)-1), otherwise a(A020639(n)) XOR a(A032742(n)). [After the definition.] - Antti Karttunen, Dec 11 2015
For n > 1, this simplifies to: a(n) = 2^(A055396(n)-1) XOR a(A032742(n)). [Where A055396(n) gives the index of the smallest prime dividing n and A032742(n) gives the largest proper divisor of n. Cf. a similar formula for A048675.]
Other identities and observations. For all n >= 0:
a(n) = A048672(A100112(A007913(n))). - Peter Kagey, Dec 10 2015
From Antti Karttunen, Dec 11 2015, Sep 19 & Oct 27 2016, Feb 15 2021: (Start)
a(n) = a(A007913(n)). [The result depends only on the squarefree part of n.]
a(n) = A048675(A007913(n)).
a(A206296(n)) = A168081(n).
a(A260443(n)) = A264977(n).
a(A265408(n)) = A265407(n).
a(A275734(n)) = A275808(n).
a(A276076(n)) = A276074(n).
a(A283477(n)) = A006068(n).
(End)
From Peter Munn, Jan 09 2021 and Apr 20 2021: (Start)
a(n) = A007814(A225546(n)).
a(A019565(n)) = n; A019565(a(n)) = A007913(n).
a(A003961(n)) = 2 * a(n).
a(A297845(n, k)) = A048720(a(n), a(k)).
a(A329329(n, k)) = A048720(a(n), a(k)).
a(A059897(n, k)) = A003987(a(n), a(k)).
a(A331590(n, k)) = a(n) + a(k).
a(A334747(n)) = a(n) + 1.
(End)

New name from Peter Munn, Nov 01 2023

A277330 a(0)=1, a(1)=2, a(2n) = A003961(a(n)), a(2n+1) = lcm(a(n),a(n+1))/gcd(a(n),a(n+1)).

Original entry on oeis.org

1, 2, 3, 6, 5, 2, 15, 30, 7, 10, 3, 30, 35, 2, 105, 210, 11, 70, 21, 30, 5, 10, 105, 42, 77, 70, 3, 210, 385, 2, 1155, 2310, 13, 770, 231, 30, 55, 70, 105, 6, 7, 2, 21, 42, 385, 10, 165, 66, 143, 110, 231, 210, 5, 70, 1155, 66, 1001, 770, 3, 2310, 5005, 2, 15015, 30030, 17, 10010, 3003, 30, 715, 770, 105, 66, 91, 154, 231, 6, 385, 70, 15, 42, 11, 14, 3, 42, 55, 2

Offset: 0

Antti Karttunen, Oct 27 2016

nonn

Each term is a squarefree number, A005117.

Antti Karttunen, Table of n, a(n) for n = 0..8191

Cf. A001511, A003961, A005117, A007913, A019565, A048675, A055396, A260443, A264977, A277707.
Cf. A023758 (positions where coincides with A260443).
Cf. A277701, A277712, A277713 for the positions of 2's, 3's and 6's in this sequence, which are also the first three rows of array A277710.
Cf. also A255483.

a(0) = 1, a(1) = 2, a(2n) = A003961(a(n)), a(2n+1) = lcm(a(n),a(n+1))/gcd(a(n),a(n+1)).
Other identities. For all n >= 0:
a(n) = A007913(A260443(n)).
a(n) = A019565(A264977(n)), A048675(a(n)) = A264977(n).
A055396(a(n)) = A277707(A260443(n)) = A001511(n).

A277817 Positions of zeros in A277815; sequence A277711 sorted into ascending order.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 9, 12, 14, 15, 16, 17, 18, 23, 24, 28, 30, 31, 32, 33, 34, 36, 46, 47, 48, 49, 56, 60, 62, 63, 64, 65, 66, 68, 72, 73, 81, 87, 92, 94, 95, 96, 97, 98, 111, 112, 120, 124, 126, 127, 128, 129, 130, 132, 135, 136, 137, 143, 144, 145, 146, 153, 159, 162, 174, 175, 177, 184, 188, 190, 191, 192, 193, 194, 196, 222

Offset: 0

Antti Karttunen, Nov 06 2016

nonn

Numbers n for which there is no such k < n that A264977(k) = A264977(n).

After zero, numbers not in the range of A277816.

Antti Karttunen, Table of n, a(n) for n = 0..10000

Cf. A264977, A277696, A277815, A277816.
Sequence A277711 sorted into ascending order.
Positions of ones in A277824.
Cf. A277884 (a left inverse).

Other identities. For all n >= 0:

A277884(a(n)) = n.

A277695 Permutation of natural numbers: a(1) = 1; for n > 1, if A277815(n) = 0, a(n) = 2a(A277814(n)-1), otherwise a(n) = 1 + 2a(A277815(n)).

Original entry on oeis.org

1, 2, 4, 8, 3, 16, 6, 32, 12, 5, 13, 64, 7, 24, 10, 26, 128, 14, 27, 17, 25, 49, 48, 20, 257, 11, 21, 52, 15, 256, 28, 54, 34, 50, 55, 98, 515, 99, 9, 65, 31, 29, 97, 105, 51, 96, 40, 514, 22, 101, 43, 35, 1031, 513, 81, 42, 69, 23, 57, 104, 63, 30, 512, 56, 108, 68, 111, 100, 139, 199, 163, 110, 196, 203, 19, 211, 2063, 33, 195, 53, 1030, 47

Offset: 1

Antti Karttunen, Nov 06 2016

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences that are permutations of the natural numbers

Inverse: A277696.

Cf. A264977, A277710, A277814, A277815, A277816, A277817.

a(1) = 1; for n > 1, if A277815(n) = 0 [when n is in A277817], a(n) = 2*a(A277814(n)-1), otherwise a(n) = 1 + 2*a(A277815(n)).

a(1) = 1, a(A277817(1+n)) = 2*a(n), a(A277816(n)) = 1 + 2*a(n). [Implicit form.]

A265407 Spironacci-style recurrence: a(0)=0, a(1)=1, a(n) = 2*a(n) XOR a(A265409(n)).

Original entry on oeis.org

0, 1, 2, 4, 8, 16, 32, 64, 129, 259, 519, 1036, 2074, 4150, 8296, 16600, 33208, 66424, 132832, 265696, 531424, 1062880, 2125696, 4251521, 8502785, 17005825, 34011905, 68023301, 136047622, 272093206, 544188470, 1088378998, 2176753882, 4353515996, 8707015520, 17414063992, 34828160840, 69656354600, 139312643368

Offset: 0

Antti Karttunen, Dec 13 2015

nonn

Spironacci-polynomials evaluated at X=2 over the field GF(2).

This is otherwise computed like A078510, which starts with a(0)=0 placed in the center of spiral (in square grid), followed by a(1) = 1, after which each term is a sum of two previous terms that are nearest when terms are arranged in a spiral, that is terms a(n-1) and a(A265409(n)), except here we first multiply the term a(n-1) by 2, and use carryless XOR (A003987) instead of normal addition.

Antti Karttunen, Table of n, a(n) for n = 0..256

Cf. A003987, A248663, A265408, A265409.

Cf. also A078510, A264977.

a(0)=0, a(1)=1; after which, a(n) = 2*a(n) XOR a(A265409(n)).

a(n) = A248663(A265408(n)).

cp's OEIS Frontend

A277826 a(n) = the least k for which A264977(k) = A264977(n).

Views

Author

Keywords

Links

Crossrefs

Programs

Scheme

A283975 Odd bisection of A264977.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

A283979 a(n) = (n XOR A264977(n))/4, where XOR is bitwise-xor (A003987).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

A260443 Prime factorization representation of Stern polynomials: a(0) = 1, a(1) = 2, a(2n) = A003961(a(n)), a(2n+1) = a(n)*a(n+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Scheme

Formula

Extensions

A023758 Numbers of the form 2^i - 2^j with i >= j.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Python

Formula

Extensions

A248663 Binary encoding of the prime factors of the squarefree part of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

A277695 Permutation of natural numbers: a(1) = 1; for n > 1, if A277815(n) = 0, a(n) = 2a(A277814(n)-1), otherwise a(n) = 1 + 2a(A277815(n)).