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

A378226 XOR-Moebius transform of A318457, where A318457(n) = n XOR (sigma(n)-n).

Original entry on oeis.org

1, 2, 3, 4, 5, 0, 7, 8, 15, 4, 11, 24, 13, 0, 1, 16, 17, 8, 19, 4, 27, 16, 23, 40, 27, 4, 27, 0, 29, 52, 31, 32, 39, 36, 45, 8, 37, 32, 57, 16, 41, 0, 43, 24, 5, 32, 47, 80, 63, 0, 53, 20, 53, 104, 41, 112, 63, 4, 59, 124, 61, 0, 7, 64, 91, 48, 67, 76, 75, 36, 71, 0, 73, 68, 103, 56, 83, 36, 79, 48, 111, 84, 83
Offset: 1

Views

Author

Antti Karttunen, Nov 26 2024

Keywords

Comments

Unique sequence satisfying SumXOR_{d divides n} a(d) = A318457(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

Cf. A000203, A001065, A003987, A318457, A378227, A378230 (positions of 0's), A378441 (fixed points).

Programs

  • PARI
    A318457(n) = bitxor(n,sigma(n)-n);
    A378226(n) = { my(v=0); fordiv(n, d, if(issquarefree(n/d), v=bitxor(v, A318457(d)))); (v); }

A378230 Positions of 0's in A378226, where A278226 is XOR-Moebius transform of A318457, and A318457(n) = n XOR (sigma(n)-n).

Original entry on oeis.org

6, 14, 28, 42, 50, 62, 72, 114, 124, 150, 186, 248, 254, 376, 402, 412, 426, 434, 474, 496, 498, 508, 762, 784, 786, 796, 868, 938, 978, 994, 1002, 1016, 1302, 1528, 1568, 1578, 1626, 1778, 1834, 1852, 1888, 1948, 1988, 2032, 3056, 3064, 3094, 3282, 3350, 3556, 3568, 3644, 3682, 3794, 3858, 3868, 3882, 3954, 4064
Offset: 1

Views

Author

Antti Karttunen, Nov 26 2024

Keywords

Crossrefs

Positions of 0's in A378226.
Subsequences: A000396 (at least the even terms).

Programs

A378441 Fixed points of A378226, where A278226 is XOR-Moebius transform of A318457, and A318457(n) = n XOR (sigma(n)-n).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 11, 13, 16, 17, 19, 23, 27, 29, 31, 32, 37, 41, 43, 47, 53, 59, 61, 64, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 125, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 243, 251, 256, 257, 263, 269, 271, 277, 281
Offset: 1

Views

Author

Antti Karttunen, Nov 26 2024

Keywords

Crossrefs

Cf. A000040 (subsequence), A001065, A003987, A318457, A378226.

Programs

Formula

{k such that A378226(k) is equal to k}.

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

A318458 a(n) = n AND A001065(n), where AND is bitwise-and (A004198) & A001065 = sum of proper divisors.

Original entry on oeis.org

0, 0, 1, 0, 1, 6, 1, 0, 0, 8, 1, 0, 1, 10, 9, 0, 1, 16, 1, 20, 1, 6, 1, 0, 0, 16, 9, 28, 1, 10, 1, 0, 1, 0, 1, 36, 1, 6, 1, 32, 1, 34, 1, 40, 33, 10, 1, 0, 0, 34, 17, 36, 1, 2, 17, 0, 17, 32, 1, 44, 1, 34, 41, 0, 1, 66, 1, 0, 1, 66, 1, 72, 1, 8, 1, 64, 1, 74, 1, 64, 0, 0, 1, 4, 21, 6, 1, 88, 1, 16, 17, 76, 1, 18, 25, 0, 1, 64, 33, 100, 1, 98, 1, 104, 65
Offset: 1

Views

Author

Antti Karttunen, Aug 26 2018

Keywords

Comments

The peculiar look of the scatterplot is partly an artifact of the logarithmic scale. Compare also to the scatterplot of A318468.

Crossrefs

Programs

  • Magma
    [SumOfDivisors(n)-BitwiseOr(n, SumOfDivisors(n)-n): n in [1..100]]; // Vincenzo Librandi, Aug 29 2018
  • Mathematica
    Table[BitAnd[n, DivisorSigma[1, n] - n], {n, 100}] (* Vincenzo Librandi, Aug 29 2018 *)
  • PARI
    A318458(n) = bitand(n,sigma(n)-n);
    

Formula

a(n) = A004198(n, A001065(n)).
a(n) = A000203(n) - A318456(n) = (A000203(n)-A318457(n))/2.

A318456 a(n) = n OR A001065(n), where OR is bitwise-or (A003986) and A001065 = sum of proper divisors.

Original entry on oeis.org

1, 3, 3, 7, 5, 6, 7, 15, 13, 10, 11, 28, 13, 14, 15, 31, 17, 23, 19, 22, 31, 30, 23, 60, 31, 26, 31, 28, 29, 62, 31, 63, 47, 54, 47, 55, 37, 54, 55, 58, 41, 62, 43, 44, 45, 62, 47, 124, 57, 59, 55, 62, 53, 118, 55, 120, 63, 58, 59, 124, 61, 62, 63, 127, 83, 78, 67, 126, 95, 78, 71, 123, 73, 106, 123, 76, 95, 94, 79, 122, 121, 126, 83
Offset: 1

Views

Author

Antti Karttunen, Aug 26 2018

Keywords

Crossrefs

Programs

  • Mathematica
    Array[BitOr[#, DivisorSigma[1, #] - #] &, 100] (* Paolo Xausa, Mar 11 2024 *)
  • PARI
    A318456(n) = bitor(n,sigma(n)-n);
    
  • Python
    from sympy import divisor_sigma
    def A318456(n): return n|divisor_sigma(n)-n # Chai Wah Wu, Jul 01 2022

Formula

a(n) = A003986(n, A001065(n)).
a(n) = A000203(n) - A318458(n).

A318467 a(n) = 2*n XOR A000203(n), where XOR is bitwise-xor (A003987) and A000203 = sum of divisors.

Original entry on oeis.org

3, 7, 2, 15, 12, 0, 6, 31, 31, 6, 26, 4, 20, 4, 6, 63, 48, 3, 50, 2, 10, 8, 54, 12, 45, 30, 30, 0, 36, 116, 30, 127, 114, 114, 118, 19, 108, 112, 118, 10, 120, 52, 122, 12, 20, 20, 110, 28, 91, 57, 46, 10, 92, 20, 38, 8, 34, 46, 74, 208, 68, 28, 22, 255, 214, 20, 194, 246, 234, 28, 198, 83, 216, 230, 234, 20, 250, 52, 206, 26
Offset: 1

Views

Author

Antti Karttunen, Aug 26 2018

Keywords

Crossrefs

Cf. A000396 (positions of zeros), A378227 (XOR-Moebius transform), A379234 (fixed points), A379236.
Cf. also A294899, A318457, A378988.

Programs

  • Mathematica
    Table[BitXor[2n,DivisorSigma[1,n]],{n,80}] (* Harvey P. Dale, Oct 30 2022 *)
  • PARI
    A318467(n) = bitxor(2*n,sigma(n));

Formula

a(n) = A003987(2*n, A000203(n)).
a(n) = A224880(n) - 2*A318468(n).
a(n) = 2*n XOR (A318457(n)+2*A318458(n)). - Antti Karttunen, Jan 08 2025

A318508 a(n) = A032742(n) AND A001065(n)-A032742(n), where AND is bitwise-and (A004198) and A001065 = sum of proper divisors and A032742 = the largest proper divisor of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 3, 0, 0, 1, 1, 0, 2, 0, 3, 4, 0, 0, 8, 0, 8, 4, 3, 0, 8, 1, 1, 0, 14, 0, 11, 0, 0, 0, 1, 6, 0, 0, 3, 4, 20, 0, 1, 0, 18, 2, 3, 0, 16, 1, 16, 0, 16, 0, 3, 2, 4, 0, 1, 0, 14, 0, 3, 20, 0, 4, 33, 0, 0, 4, 35, 0, 4, 0, 1, 24, 2, 8, 35, 0, 0, 9, 1, 0, 34, 0, 3, 4, 32, 0, 33, 8, 14, 4, 3, 2, 32, 0, 16, 0, 2, 0, 51, 0, 52, 32
Offset: 1

Views

Author

Antti Karttunen, Aug 28 2018

Keywords

Crossrefs

Programs

Formula

a(n) = A004198(A032742(n), A318505(n)).
For n > 1, a(n) = A001065(n) - A318506(n) = (A001065(n) - A318507(n))/2.

A318507 a(n) = A032742(n) XOR A001065(n)-A032742(n), where XOR is bitwise-or (A003987) and A001065 = sum of proper divisors and A032742 = the largest proper divisor of n.

Original entry on oeis.org

1, 1, 1, 3, 1, 0, 1, 7, 2, 6, 1, 12, 1, 4, 1, 15, 1, 5, 1, 6, 3, 8, 1, 20, 4, 14, 13, 0, 1, 20, 1, 31, 15, 18, 1, 55, 1, 16, 9, 10, 1, 52, 1, 4, 29, 20, 1, 44, 6, 11, 21, 14, 1, 60, 13, 56, 23, 30, 1, 80, 1, 28, 1, 63, 11, 12, 1, 58, 19, 4, 1, 115, 1, 38, 1, 60, 3, 20, 1, 106, 22, 42, 1, 72, 23, 40, 25, 28, 1, 78, 5, 48, 27, 44, 21, 92, 1
Offset: 1

Views

Author

Antti Karttunen, Aug 28 2018

Keywords

Comments

Note that here zeros occur only on even perfect numbers (even terms of A000396), in contrast to A318457, which would be zero also for any hypothetical odd perfect number. - Antti Karttunen, Aug 29 2018

Crossrefs

Programs

Formula

a(n) = A003987(A032742(n), A318505(n)).
For n > 1, a(n) = A001065(n) - 2*A318508(n).

A324867 a(n) = A156552(n) XOR A324865(n), where XOR is bitwise-xor, A003987.

Original entry on oeis.org

0, 1, 3, 2, 7, 4, 15, 6, 0, 13, 31, 10, 63, 16, 2, 6, 127, 12, 255, 18, 7, 46, 511, 22, 28, 82, 4, 46, 1023, 30, 2047, 30, 54, 174, 2, 22, 4095, 256, 12, 54, 8191, 36, 16383, 66, 24, 798, 32767, 46, 60, 31, 248, 130, 65535, 28, 19, 70, 12, 1300, 131071, 42, 262143, 2734, 48, 22, 126, 94, 524287, 302, 774, 40, 1048575, 38, 2097151, 4354, 10
Offset: 1

Views

Author

Antti Karttunen, Mar 18 2019

Keywords

Crossrefs

Programs

  • PARI
    A156552(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}; \\ From A156552
    A318457(n) = bitxor(n,sigma(n)-n);
    A324867(n) = if(1==n,0,A318457(A156552(n)));
    
  • PARI
    A324867(n) = { my(k=A156552(n)); bitxor(k,(A323243(n)-k)); }; \\ Needs also code from A323243.

Formula

a(1) = 0; for n > 1, a(n) = A318457(A156552(n)).
a(n) = A156552(n) XOR (A323243(n) - A156552(n)).
a(n) = A323243(n) - 2*A324398(n).
Showing 1-10 of 12 results. Next