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

A295901 Unique sequence satisfying SumXOR_{d divides n} a(d) = n^2 for any n > 0, where SumXOR is the analog of summation under the binary XOR operation.

Original entry on oeis.org

1, 5, 8, 20, 24, 40, 48, 80, 88, 120, 120, 160, 168, 240, 240, 320, 288, 312, 360, 480, 384, 408, 528, 640, 616, 520, 648, 960, 840, 816, 960, 1280, 1072, 1440, 1248, 1248, 1368, 1224, 1360, 1920, 1680, 1920, 1848, 1632, 1872, 2640, 2208, 2560, 2384, 3016
Offset: 1

Views

Author

Rémy Sigrist, Nov 29 2017

Keywords

Comments

This sequence is a variant of A256739; both sequences have nice graphical features.
Replacing "SumXOR" by "Sum" in the name leads to the Jordan function J_2 (A007434).
For any sequence f of nonnegative integers with positive indices:
- let x_f be the unique sequence satisfying SumXOR_{d divides n} x_f(d) = f(n) for any n > 0,
- in particular, x_A000027 = A256739 and x_A000290 = a (this sequence),
- also, x_A178910 = A000027 and x_A055895 = A000079,
- see the links section for a gallery of x_f plots for some classic f functions,
- x_f(1) = f(1),
- x_f(p) = f(1) XOR f(p) for any prime p,
- x_f(n) = SumXOR_{d divides n and n/d is squarefree} f(d) for any n > 0,
- the function x: f -> x_f is a bijection,
- A000004 is the only fixed point of x (i.e. x_f = f if and only if f = A000004),
- for any sequence f, x_{2*f} = 2 * x_f,
- for any sequences g and f, x_{g XOR f} = x_g XOR x_f.
From Antti Karttunen, Dec 29 2017: (Start)
The transform x_f described above could be called "Xor-Moebius transform of f" because of its analogous construction to Möbius transform with A008683 replaced by A008966 and the summation replaced by cumulative XOR.
(End)

Crossrefs

Programs

  • PARI
    a(n{, f=k->k^2}) = my (v=0); fordiv(n,d,if (issquarefree(n/d), v=bitxor(v,f(d)))); return (v)

Formula

a(n) = SumXOR_{d divides n and n/d is squarefree} d^2.

A296203 Xor-Moebius transform of A000203, the sum of divisors.

Original entry on oeis.org

1, 2, 5, 4, 7, 10, 9, 8, 9, 22, 13, 20, 15, 18, 27, 16, 19, 34, 21, 60, 45, 42, 25, 40, 25, 38, 37, 36, 31, 78, 33, 32, 57, 38, 63, 108, 39, 42, 51, 120, 43, 90, 45, 116, 95, 82, 49, 80, 49, 86, 95, 76, 55, 122, 67, 72, 65, 70, 61, 204, 63, 66, 65, 64, 93, 130, 69, 76, 125, 166, 73, 184, 75, 86, 125, 180, 101, 190, 81, 240, 81
Offset: 1

Views

Author

Antti Karttunen, Dec 25 2017, from Rémy Sigrist's "gallery of scatter plots" in A295901

Keywords

Comments

Unique sequence satisfying SumXOR_{d divides n} a(d) = sigma(n) for all n > 0, where SumXOR is the analog of summation under the binary XOR operation. See A295901 for a list of some of the properties of Xor-Moebius transform.

Crossrefs

Programs

  • PARI
    A296203(n) = { my(v=0); fordiv(n, d, if(issquarefree(n/d), v=bitxor(v, sigma(d)))); (v); } \\ after code in A295901.

A297110 Xor-Moebius transform of A006068, inverse of the binary Gray code.

Original entry on oeis.org

1, 2, 3, 4, 7, 4, 4, 8, 12, 8, 12, 8, 8, 12, 15, 16, 31, 20, 28, 16, 31, 20, 27, 16, 23, 24, 28, 24, 23, 16, 20, 32, 48, 32, 63, 40, 56, 36, 48, 32, 48, 32, 51, 40, 48, 44, 52, 32, 36, 56, 63, 48, 39, 36, 47, 48, 48, 56, 44, 32, 40, 60, 63, 64, 112, 80, 124, 64, 96, 64, 123, 80, 112, 72, 111, 72, 127, 80, 116, 64
Offset: 1

Views

Author

Antti Karttunen, Dec 25 2017

Keywords

Comments

Unique sequence satisfying SumXOR_{d divides n} a(d) = A006068(n) for all n > 0, where SumXOR is the analog of summation under the binary XOR operation. See A295901 for a list of some of the properties of the Xor-Moebius transform.

Crossrefs

Programs

  • PARI
    A006068(n)= { my(s=1, ns); while(1, ns = n >> s; if(0==ns, break()); n = bitxor(n, ns); s <<= 1; ); return (n); } \\ Essentially Joerg Arndt's Jul 19 2012 code.
    A297110(n) = { my(v=0); fordiv(n, d, if(issquarefree(n/d), v=bitxor(v, A006068(d)))); (v); };
Showing 1-3 of 3 results.