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

A115872 Square array where row n gives all solutions k > 0 to the cross-domain congruence n*k = A048720(A065621(n),k), zero sequence (A000004) if no such solutions exist.

Original entry on oeis.org

1, 2, 1, 3, 2, 3, 4, 3, 6, 1, 5, 4, 7, 2, 7, 6, 5, 12, 3, 14, 3, 7, 6, 14, 4, 15, 6, 7, 8, 7, 15, 5, 28, 7, 14, 1, 9, 8, 24, 6, 30, 12, 15, 2, 15, 10, 9, 28, 7, 31, 14, 28, 3, 30, 7, 11, 10, 30, 8, 56, 15, 30, 4, 31, 14, 3, 12, 11, 31, 9, 60, 24, 31, 5, 60, 15, 6, 3, 13, 12, 48, 10, 62, 28, 56, 6, 62, 28, 12, 6, 5, 14, 13, 51, 11, 63, 30, 60, 7, 63, 30, 15, 7, 10, 7
Offset: 1

Views

Author

Antti Karttunen, Feb 07 2006

Keywords

Comments

Here * stands for ordinary multiplication and X means carryless (GF(2)[X]) multiplication (A048720).
Square array is read by descending antidiagonals, as A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), etc.
Rows at positions 2^k are 1, 2, 3, ..., (A000027). Row 2n is equal to row n.
Numbers on each row give a subset of positions of zeros at the corresponding row of A284270. - Antti Karttunen, May 08 2019

Examples

			Fifteen initial terms of rows 1 - 19 are listed below:
   1:  1,  2,  3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15, ...
   2:  1,  2,  3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15, ...
   3:  3,  6,  7,  12,  14,  15,  24,  28,  30,  31,  48,  51,  56,  60,  62, ...
   4:  1,  2,  3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15, ...
   5:  7, 14, 15,  28,  30,  31,  56,  60,  62,  63, 112, 120, 124, 126, 127, ...
   6:  3,  6,  7,  12,  14,  15,  24,  28,  30,  31,  48,  51,  56,  60,  62, ...
   7:  7, 14, 15,  28,  30,  31,  56,  60,  62,  63, 112, 120, 124, 126, 127, ...
   8:  1,  2,  3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15, ...
   9: 15, 30, 31,  60,  62,  63, 120, 124, 126, 127, 240, 248, 252, 254, 255, ...
  10:  7, 14, 15,  28,  30,  31,  56,  60,  62,  63, 112, 120, 124, 126, 127, ...
  11:  3,  6, 12,  15,  24,  27,  30,  31,  48,  51,  54,  60,  62,  63,  96, ...
  12:  3,  6,  7,  12,  14,  15,  24,  28,  30,  31,  48,  51,  56,  60,  62, ...
  13:  5, 10, 15,  20,  21,  30,  31,  40,  42,  45,  47,  60,  61,  62,  63, ...
  14:  7, 14, 15,  28,  30,  31,  56,  60,  62,  63, 112, 120, 124, 126, 127, ...
  15: 15, 30, 31,  60,  62,  63, 120, 124, 126, 127, 240, 248, 252, 254, 255, ...
  16:  1,  2,  3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15, ...
  17: 31, 62, 63, 124, 126, 127, 248, 252, 254, 255, 496, 504, 508, 510, 511, ...
  18: 15, 30, 31,  60,  62,  63, 120, 124, 126, 127, 240, 248, 252, 254, 255, ...
  19:  7, 14, 28,  31,  56,  62,  63, 112, 119, 124, 126, 127, 224, 238, 248, ...

Links

Crossrefs

Transpose: A114388. First column: A115873.
Cf. A048720, A065621, A115857, A115871, A325565, A325566, A325567, A325568, A325570, A325571.
Cf. also arrays A277320, A277810, A277820, A284270.
A few odd-positioned rows: row 1: A000027, Row 3: A048717, Row 5: A115770 (? Checked for all values less than 2^20), Row 7: A115770, Row 9: A115801, Row 11: A115803, Row 13: A115772, Row 15: A115801 (? Checked for all values less than 2^20), Row 17: A115809, Row 19: A115874, Row 49: A114384, Row 57: A114386.

Programs

  • Mathematica
    X[a_, b_] := Module[{A, B, C, x},
     A = Reverse@IntegerDigits[a, 2];
     B = Reverse@IntegerDigits[b, 2];
     C = Expand[
        Sum[A[[i]]*x^(i-1), {i, 1, Length[A]}]*
        Sum[B[[i]]*x^(i-1), {i, 1, Length[B]}]];
     PolynomialMod[C, 2] /. x -> 2];
T[n_, k_] := Module[{x = BitXor[n-1, 2n-1], k0 = k},
     For[i = 1, True, i++, If[n*i == X[x, i],
     If[k0 == 1, Return[i], k0--]]]];
Table[T[n-k+1, k], {n, 1, 14}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jan 04 2022 *)
  • PARI
    up_to = 120;
A048720(b,c) = fromdigits(Vec(Pol(binary(b))*Pol(binary(c)))%2, 2);
A065621(n) = bitxor(n-1,n+n-1);
A115872sq(n, k) = { my(x = A065621(n)); for(i=1,oo,if((n*i)==A048720(x,i),if(1==k,return(i),k--))); };
A115872list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A115872sq(col,(a-(col-1))))); (v); };
v115872 = A115872list(up_to);
A115872(n) = v115872[n]; \\ (Slow) - Antti Karttunen, May 08 2019

Extensions

Example section added and the data section extended up to n=105 by Antti Karttunen, May 08 2019

A325570 Numbers n that have no divisor d > 1 such that A048720(A065621(d),n/d) = n.

Original entry on oeis.org

1, 3, 5, 7, 11, 13, 15, 17, 19, 23, 25, 27, 29, 31, 37, 39, 41, 43, 47, 51, 53, 55, 57, 59, 61, 63, 67, 69, 71, 73, 77, 79, 81, 83, 85, 87, 89, 91, 95, 97, 99, 101, 103, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 131, 137, 139, 141, 143, 145, 147, 149, 151, 157, 159, 163, 167, 169, 171, 173, 175, 177, 179, 181
Offset: 1

Views

Author

Antti Karttunen, May 10 2019

Keywords

Links

Crossrefs

Cf. A048720, A065621, A115872, A115873.
Positions of ones in A325565 and A325566.
Cf. A065091 (a subsequence), A325571 (the composite terms), A325572 (complement).
Subsequence of A005408 (odd numbers).

Programs

  • PARI
    A048720(b,c) = fromdigits(Vec(Pol(binary(b))*Pol(binary(c)))%2, 2);
A065621(n) = bitxor(n-1,n+n-1);
isA325570(n) = fordiv(n,d,if(A048720(A065621(n/d),d)==n,return(d==n)));

A325573 Odd numbers n that have divisor d > 1 such that A048720(A065621(d),n/d) = n.

Original entry on oeis.org

9, 21, 33, 35, 45, 49, 65, 75, 93, 105, 129, 133, 135, 153, 155, 161, 165, 189, 195, 217, 225, 259, 273, 279, 297, 309, 315, 341, 345, 381, 385, 403, 441, 465, 513, 525, 527, 561, 567, 585, 589, 597, 611, 621, 635, 645, 651, 681, 693, 705, 713, 729, 765, 775, 793, 819, 837, 889, 899, 945, 961, 1025, 1029, 1035, 1057, 1065
Offset: 1

Views

Author

Antti Karttunen, May 10 2019

Keywords

Links

Crossrefs

Cf. A048720, A065621, A115872, A325566, A325567, A325571.
Subsequence of A071904 and of A325572.

Programs

  • PARI
    A048720(b,c) = fromdigits(Vec(Pol(binary(b))*Pol(binary(c)))%2, 2);
A065621(n) = bitxor(n-1,n+n-1);
isA325573(n) = ((n%2)&&fordiv(n,d,if(A048720(A065621(n/d),d)==n,return(d

A383363 Composite numbers k all of whose proper divisors have binary weights that are not equal to the binary weight of k.

Original entry on oeis.org

15, 25, 27, 39, 51, 55, 57, 63, 69, 77, 81, 85, 87, 91, 95, 99, 111, 115, 117, 119, 121, 123, 125, 141, 143, 145, 147, 159, 169, 171, 175, 177, 183, 185, 187, 201, 203, 205, 207, 209, 213, 215, 219, 221, 231, 235, 237, 243, 245, 247, 249, 253, 255, 261, 265, 275
Offset: 1

Views

Author

Amiram Eldar, Apr 24 2025

Keywords

Comments

First differs from A325571 at n = 56: A325571(56) = 267 is not a term of this sequence. Differs from A325571 by having the terms 16849, 35235, 101265, 268357, 295717, ..., and not having the terms 267, 295, 327, 387, 395, ... .
Composite numbers k such that A380844(k) = 1.
All the odd primes p have A380844(p) = 1.
All the terms are odd numbers since for an even number k, A000120(k) = A000120(k/2).

Examples

			15 = 3 * 5 is a term since it is composite, and its binary weight, A000120(15) = 4 is different from the binary weights of its proper divisors: A000120(1) = 1, A000120(3) = 2, and A000120(5) = 2.

Links

Crossrefs

Cf. A000120, A325571, A380844, A383364, A383365.

Programs

  • Mathematica
    q[k_] := CompositeQ[k] && DivisorSum[k, 1 &, DigitCount[#, 2, 1] == DigitCount[k, 2, 1] &] == 1; Select[Range[1, 300, 2], q]
  • PARI
    isok(k) = if(k == 1 || isprime(k), 0, my(h = hammingweight(k)); sumdiv(k, d, hammingweight(d) == h) == 1);
Showing 1-4 of 4 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.