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

A275809 Indices of zeros in A275808.

Original entry on oeis.org

0, 5, 14, 19, 22, 54, 59, 74, 84, 89, 93, 97, 100, 111, 114, 119, 264, 269, 278, 283, 286, 366, 371, 408, 413, 422, 427, 430, 440, 463, 466, 482, 492, 497, 501, 536, 552, 557, 566, 571, 574, 579, 589, 592, 596, 601, 604, 615, 618, 623, 655, 658, 675, 685, 688, 692, 696, 701, 710, 715, 718, 1560, 1565, 1574, 1579, 1582, 1614, 1619, 1634, 1644
Offset: 0

Views

Author

Antti Karttunen, Aug 09 2016

Keywords

Comments

Indexing begins from zero, because a(0) = 0 is a special case in this sequence.
Terms A009445(n)-1, [0, 5, 119, 5039, 362879, 39916799, ...] form a subsequence, and also the terms of A010050(n)-2, [0, 22, 718, 40318, 3628798, ...] form a subsequence.

Links

Crossrefs

Cf. A275736, A275808.
Cf. A275810 (first differences).
Subsequence of A275813 and of A275805 (after the initial 0).
Cf. also A009445, A010050.

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

Views

Author

Peter Kagey, Jan 11 2015

Keywords

Comments

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)

Examples

			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)

Links

Crossrefs

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.

Programs

  • Haskell
    import Data.Bits (xor)
a248663 = foldr (xor) 0 . map (\i -> 2^(i - 1)) . a112798_row
-- Peter Kagey, Sep 16 2016
  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    A248663(n) = vecsum(apply(p -> 2^(primepi(p)-1),factor(core(n))[,1])); \\ Antti Karttunen, Feb 15 2021
  • Python
    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
  • Ruby
    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

Formula

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)

Extensions

New name from Peter Munn, Nov 01 2023

A275734 Prime-factorization representations of "factorial base slope polynomials": a(0) = 1; for n >= 1, a(n) = A275732(n) * a(A257684(n)).

Original entry on oeis.org

1, 2, 3, 6, 2, 4, 5, 10, 15, 30, 10, 20, 3, 6, 9, 18, 6, 12, 2, 4, 6, 12, 4, 8, 7, 14, 21, 42, 14, 28, 35, 70, 105, 210, 70, 140, 21, 42, 63, 126, 42, 84, 14, 28, 42, 84, 28, 56, 5, 10, 15, 30, 10, 20, 25, 50, 75, 150, 50, 100, 15, 30, 45, 90, 30, 60, 10, 20, 30, 60, 20, 40, 3, 6, 9, 18, 6, 12, 15, 30, 45, 90, 30, 60, 9, 18, 27
Offset: 0

Views

Author

Antti Karttunen, Aug 08 2016

Keywords

Comments

These are prime-factorization representations of single-variable polynomials where the coefficient of term x^(k-1) (encoded as the exponent of prime(k) in the factorization of n) is equal to the number of nonzero digits that occur on the slope (k-1) levels below the "maximal slope" in the factorial base representation of n. See A275811 for the definition of the "digit slopes" in this context.

Examples

			For n=23 ("321" in factorial base representation, A007623), all three nonzero digits are maximal for their positions (they all occur on "maximal slope"), thus a(23) = prime(1)^3 = 2^3 = 8.
For n=29 ("1021"), there are three nonzero digits, where both 2 and the rightmost 1 are on the "maximal slope", while the most significant 1 is on the "sub-sub-sub-maximal", thus a(29) = prime(1)^2 * prime(4)^1 = 2*7 = 28.
For n=37 ("1201"), there are three nonzero digits, where the rightmost 1 is on the maximal slope, 2 is on the sub-maximal, and the most significant 1 is on the "sub-sub-sub-maximal", thus a(37) = prime(1) * prime(2) * prime(4) = 2*3*7 = 42.
For n=55 ("2101"), the least significant 1 is on the maximal slope, and the digits "21" at the beginning are together on the sub-sub-maximal slope (as they are both two less than the maximal digit values 4 and 3 allowed in those positions), thus a(55) = prime(1)^1 * prime(3)^2 = 2*25 = 50.

Links

Crossrefs

Cf. A001221, A001222, A002110, A007489, A007814, A048675, A051903, A056169, A056170, A060130, A060502, A225901.
Cf. A257684, A275732.
Cf. A275811.
Cf. A260736, A275728, A275808, A275812, A275946, A275947, A275962.
Cf. A275804 (indices of squarefree terms), A275805 (of terms not squarefree).
Cf. also A275725, A275733, A275735, A276076 for other such prime factorization encodings of A060117/A060118-related polynomials.

Programs

  • Python
    from operator import mul
from sympy import prime, factorial as f
def a007623(n, p=2): return n if n
    0 else '0' for i in x)[::-1]
    return 0 if n==1 else sum(int(y[i])*f(i + 1) for i in range(len(y)))
def a(n): return 1 if n==0 else a275732(n)*a(a257684(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 19 2017

Formula

a(0) = 1; for n >= 1, a(n) = A275732(n) * a(A257684(n)).
Other identities and observations. For all n >= 0:
a(n) = A275735(A225901(n)).
a(A007489(n)) = A002110(n).
A001221(a(n)) = A060502(n).
A001222(a(n)) = A060130(n).
A007814(a(n)) = A260736(n).
A051903(a(n)) = A275811(n).
A048675(a(n)) = A275728(n).
A248663(a(n)) = A275808(n).
A056169(a(n)) = A275946(n).
A056170(a(n)) = A275947(n).
A275812(a(n)) = A275962(n).

A275728 a(0) = 0, for n >= 1, a(n) = A275736(n) + a(A257684(n)); a(n) = A048675(A275734(n)).

Original entry on oeis.org

0, 1, 2, 3, 1, 2, 4, 5, 6, 7, 5, 6, 2, 3, 4, 5, 3, 4, 1, 2, 3, 4, 2, 3, 8, 9, 10, 11, 9, 10, 12, 13, 14, 15, 13, 14, 10, 11, 12, 13, 11, 12, 9, 10, 11, 12, 10, 11, 4, 5, 6, 7, 5, 6, 8, 9, 10, 11, 9, 10, 6, 7, 8, 9, 7, 8, 5, 6, 7, 8, 6, 7, 2, 3, 4, 5, 3, 4, 6, 7, 8, 9, 7, 8, 4, 5, 6, 7, 5, 6, 3, 4, 5, 6, 4, 5, 1, 2, 3, 4
Offset: 0

Views

Author

Antti Karttunen, Aug 09 2016

Keywords

Comments

See graph: pine trees on a snowy mountain. - N. J. A. Sloane, Aug 12 2016

Links

Crossrefs

Cf. A048675, A257684, A275734, A275736.
Cf. also A275727, A275808.

Formula

a(0) = 0, for n >= 1, a(n) = A275736(n) + a(A257684(n)).
a(n) = A048675(A275734(n)).

A276074 A276076-polynomials evaluated at X=2 over the field GF(2): a(n) = A248663(A276076(n)).

Original entry on oeis.org

0, 1, 2, 3, 0, 1, 4, 5, 6, 7, 4, 5, 0, 1, 2, 3, 0, 1, 4, 5, 6, 7, 4, 5, 8, 9, 10, 11, 8, 9, 12, 13, 14, 15, 12, 13, 8, 9, 10, 11, 8, 9, 12, 13, 14, 15, 12, 13, 0, 1, 2, 3, 0, 1, 4, 5, 6, 7, 4, 5, 0, 1, 2, 3, 0, 1, 4, 5, 6, 7, 4, 5, 8, 9, 10, 11, 8, 9, 12, 13, 14, 15, 12, 13, 8, 9, 10, 11, 8, 9, 12, 13, 14, 15, 12, 13, 0, 1, 2, 3, 0, 1, 4, 5, 6, 7, 4, 5, 0, 1, 2, 3, 0, 1
Offset: 0

Views

Author

Antti Karttunen, Aug 18 2016

Keywords

Links

Crossrefs

Cf. A248663, A276076.
Cf. also A276073, A275808

Programs

Formula

a(n) = A248663(A276076(n)).

A276010 a(0) = 0, for n >= 1, a(n) = A275736(n) OR a(A257684(n)), where OR is given by A003986.

Original entry on oeis.org

0, 1, 2, 3, 1, 1, 4, 5, 6, 7, 5, 5, 2, 3, 2, 3, 3, 3, 1, 1, 3, 3, 1, 1, 8, 9, 10, 11, 9, 9, 12, 13, 14, 15, 13, 13, 10, 11, 10, 11, 11, 11, 9, 9, 11, 11, 9, 9, 4, 5, 6, 7, 5, 5, 4, 5, 6, 7, 5, 5, 6, 7, 6, 7, 7, 7, 5, 5, 7, 7, 5, 5, 2, 3, 2, 3, 3, 3, 6, 7, 6, 7, 7, 7, 2, 3, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 1, 1, 3, 3, 1, 1, 5, 5, 7, 7, 5, 5, 3
Offset: 0

Views

Author

Antti Karttunen, Aug 17 2016

Keywords

Links

Crossrefs

Cf. A003986, A087207, A257684, A275734.
Cf. also A000120, A060502, A275736, A275728, A275808.

Formula

a(0) = 0, for n >= 1, a(n) = A275736(n) OR a(A257684(n)), where OR is given by A003986.
a(n) = A087207(A275734(n)).
Other identities. For all n >= 1:
A000120(a(n)) = A060502(n).
Showing 1-6 of 6 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.