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.

Previous Showing 11-20 of 28 results. Next

A091205 Factorization and index-recursion preserving isomorphism from binary codes of GF(2) polynomials to integers.

Original entry on oeis.org

0, 1, 2, 3, 4, 9, 6, 5, 8, 15, 18, 7, 12, 23, 10, 27, 16, 81, 30, 13, 36, 25, 14, 69, 24, 11, 46, 45, 20, 21, 54, 19, 32, 57, 162, 115, 60, 47, 26, 63, 72, 61, 50, 33, 28, 135, 138, 17, 48, 35, 22, 243, 92, 39, 90, 37, 40, 207, 42, 83, 108, 29, 38, 75, 64, 225, 114, 103
Offset: 0

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

This "deeply multiplicative" bijection is one of the deep variants of A091203 which satisfy most of the same identities as the latter, but it additionally preserves also the structures where we recurse on irreducible polynomial's A014580-index. E.g., we have: A091238(n) = A061775(a(n)). The reason this holds is that when the permutation is restricted to the binary codes for irreducible polynomials over GF(2) (A014580), it induces itself: a(n) = A049084(a(A014580(n))).
On the other hand, when this permutation is restricted to the union of {1} and reducible polynomials over GF(2) (A091242), permutation A245813 is induced.

Crossrefs

Programs

  • PARI
    allocatemem(123456789);
    v091226 = vector(2^22);
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
    n=2; while((n < 2^22), if(isA014580(n), v091226[n] = v091226[n-1]+1, v091226[n] = v091226[n-1]); n++)
    A091226(n) = v091226[n];
    A091205(n) = if(n<=1,n,if(isA014580(n),prime(A091205(A091226(n))),{my(irfs,t); irfs=subst(lift(factor(Mod(1,2)*Pol(binary(n)))),x,2); irfs[,1]=apply(t->A091205(t),irfs[,1]); factorback(irfs)}));
    for(n=0, 8192, write("b091205.txt", n, " ", A091205(n)));
    \\ Antti Karttunen, Aug 16 2014

Formula

a(0)=0, a(1)=1. For n that is coding an irreducible polynomial, that is if n = A014580(i), we have a(n) = A000040(a(i)) and for reducible polynomials a(ir_i X ir_j X ...) = a(ir_i) * a(ir_j) * ..., where ir_i = A014580(i), X stands for carryless multiplication of polynomials over GF(2) (A048720) and * for the ordinary multiplication of integers (A004247).
As a composition of related permutations:
a(n) = A245821(A245704(n)).
Other identities.
For all n >= 0, the following holds:
a(A091230(n)) = A007097(n). [Maps iterates of A014580 to the iterates of primes. Permutation A245704 has the same property.]
For all n >= 1, the following holds:
A010051(a(n)) = A091225(n). [After a(1)=1, maps binary representations of irreducible GF(2) polynomials, A014580, bijectively to primes and the binary representations of corresponding reducible polynomials, A091242, to composite numbers, in some order. The permutations A091203, A106443, A106445, A106447, A235042 and A245704 have the same property.]

Extensions

Name changed by Antti Karttunen, Aug 16 2014

A091203 Factorization-preserving isomorphism from binary codes of GF(2) polynomials to integers.

Original entry on oeis.org

0, 1, 2, 3, 4, 9, 6, 5, 8, 15, 18, 7, 12, 11, 10, 27, 16, 81, 30, 13, 36, 25, 14, 33, 24, 17, 22, 45, 20, 21, 54, 19, 32, 57, 162, 55, 60, 23, 26, 63, 72, 29, 50, 51, 28, 135, 66, 31, 48, 35, 34, 243, 44, 39, 90, 37, 40, 99, 42, 41, 108, 43, 38, 75, 64, 225, 114, 47, 324
Offset: 0

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Comments

E.g. we have the following identities: A000040(n) = a(A014580(n)), A091219(n) = A008683(a(n)), A091220(n) = A000005(a(n)), A091221(n) = A001221(a(n)), A091222(n) = A001222(a(n)), A091225(n) = A010051(a(n)), A091227(n) = A049084(a(n)), A091247(n) = A066247(a(n)).

Crossrefs

Programs

  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
    A091225(n) = polisirreducible(Pol(binary(n))*Mod(1, 2));
    A305419(n) = if(n<3,1, my(k=n-1); while(k>1 && !A091225(k),k--); (k));
    A305422(n) = { my(f = subst(lift(factor(Pol(binary(n))*Mod(1, 2))),x,2)); for(i=1,#f~,f[i,1] = Pol(binary(A305419(f[i,1])))); fromdigits(Vec(factorback(f))%2,2); };
    A091203(n) = if(n<=1,n,if(!(n%2),2*A091203(n/2),A003961(A091203(A305422(n))))); \\ Antti Karttunen, Jun 10 2018

Formula

a(0)=0, a(1)=1. For n's coding an irreducible polynomial ir_i, that is if n=A014580(i), we have a(n) = A000040(i) and for composite polynomials a(ir_i X ir_j X ...) = p_i * p_j * ..., where p_i = A000040(i) and X stands for carryless multiplication of GF(2)[X] polynomials (A048720) and * for the ordinary multiplication of integers (A004247).
Other identities. For all n >= 1, the following holds:
A010051(a(n)) = A091225(n). [After a(1)=1, maps binary representations of irreducible GF(2) polynomials, A014580, to primes and the binary representations of corresponding reducible polynomials, A091242, to composite numbers. The permutations A091205, A106443, A106445, A106447, A235042 and A245704 have the same property.]
From Antti Karttunen, Jun 10 2018: (Start)
For n <= 1, a(n) = n, for n > 1, a(n) = 2*a(n/2) if n is even, and if n is odd, then a(n) = A003961(a(A305422(n))).
a(n) = A005940(1+A305418(n)) = A163511(A305428(n)).
A046523(a(n)) = A278233(n).
(End)

A235042 Factorization-preserving bijection from GF(2)[X]-polynomials to nonnegative integers, version which fixes the elements that are irreducible in both semirings.

Original entry on oeis.org

0, 1, 2, 3, 4, 9, 6, 7, 8, 21, 18, 11, 12, 13, 14, 27, 16, 81, 42, 19, 36, 49, 22, 39, 24, 5, 26, 63, 28, 33, 54, 31, 32, 93, 162, 91, 84, 37, 38, 99, 72, 41, 98, 15, 44, 189, 78, 47, 48, 77, 10, 243, 52, 57, 126, 17, 56, 117, 66, 59, 108, 61, 62, 147, 64, 441
Offset: 0

Views

Author

Antti Karttunen, Jan 02 2014

Keywords

Comments

Like A091203 this is a factorization-preserving isomorphism from GF(2)[X]-polynomials to integers. The former are encoded in the binary representation of n like this: n=11, '1011' in binary, stands for polynomial x^3+x+1, n=25, '11001' in binary, stands for polynomial x^4+x^3+1. However, this version does not map the irreducible GF(2)[X] polynomials (A014580) straight to the primes (A000040), but instead fixes the intersection of those two sets (A091206), and maps the elements in their set-wise difference A014580 \ A000040 (= A091214) in numerical order to the set-wise difference A000040 \ A014580 (= A091209).
The composite values are defined by the multiplicativity. E.g., we have a(A048724(n)) = 3*a(n) and a(A001317(n)) = A000244(n) = 3^n for all n.
This map satisfies many of the same identities as A091203, e.g., we have A091219(n) = A008683(a(n)), A091220(n) = A000005(a(n)), A091221(n) = A001221(a(n)), A091222(n) = A001222(a(n)), A091225(n) = A010051(a(n)) and A091247(n) = A066247(a(n)) for all n >= 1.

Examples

			Here (t X u) = A048720(t,u):
a(2)=2, a(3)=3 and a(7)=7, as 2, 3 and 7 are all in A091206.
a(4) = a(2 X 2) = a(2)*a(2) = 2*2 = 4.
a(5) = a(3 X 3) = a(3)*a(3) = 3*3 = 9.
a(9) = a(3 X 7) = a(3)*a(7) = 3*7 = 21.
a(10) = a(2 X 3 X 3) = a(2)*a(3)*a(3) = 2*3*3 = 18.
a(15) = a(3 X 3 X 3) = a(3)^3 = 3^3 = 27.
a(17) = a(3 X 3 X 3 X 3) = a(3)^4 = 3^4 = 81.
a(21) = a(7 X 7) = a(7)*a(7) = 7*7 = 49.
a(25) = 5, as 25 is the first term of A091214 and 5 is the first term of A091209.
a(50) = a(2 X 25) = a(2)*a(25) = 2*5 = 10.
		

Crossrefs

Inverse: A235041. Fixed points: A235045.
Similar cross-multiplicative permutations: A091203, A091205, A106443, A106445, A106447.

Formula

a(0)=0, a(1)=1, a(p) = p for those irreducible GF(2)[X]-polynomials whose binary encoding is a prime (i.e., p is in A091206), and for the rest of irreducible GF(2)[X]-polynomials (those which are encoded by a composite natural number, i.e., q is in A091214), a(q) = A091209(A235044(q)), and for reducible polynomials, a(i X j X k X ...) = a(i) * a(j) * a(k) * ..., where each i, j, k, ... is in A014580, X stands for carryless multiplication of GF(2)[X] polynomials (A048720) and * for the ordinary multiplication of integers (A004247).

A176642 Triangle T(n, k) = 8^(k*(n-k)), read by rows.

Original entry on oeis.org

1, 1, 1, 1, 8, 1, 1, 64, 64, 1, 1, 512, 4096, 512, 1, 1, 4096, 262144, 262144, 4096, 1, 1, 32768, 16777216, 134217728, 16777216, 32768, 1, 1, 262144, 1073741824, 68719476736, 68719476736, 1073741824, 262144, 1, 1, 2097152, 68719476736, 35184372088832, 281474976710656, 35184372088832, 68719476736, 2097152, 1
Offset: 0

Views

Author

Roger L. Bagula, Apr 22 2010

Keywords

Examples

			Triangle begins as:
  1;
  1,      1;
  1,      8,          1;
  1,     64,         64,           1;
  1,    512,       4096,         512,           1;
  1,   4096,     262144,      262144,        4096,          1;
  1,  32768,   16777216,   134217728,    16777216,      32768,      1;
  1, 262144, 1073741824, 68719476736, 68719476736, 1073741824, 262144, 1;
		

Crossrefs

Cf. this sequence (q=2), A176643 (q=3), A176644 (q=4).
Cf. A117401 (m=0), A118180 (m=1), A118185 (m=2), A118190 (m=3), A158116 (m=4), this sequence (m=6), A158117 (m=8), A176627 (m=10), A176639 (m=13), A156581 (m=15), A176643 (m=19), A176631 (m=20), A176641 (m=26).

Programs

  • Magma
    [8^(k*(n-k)): k in [0..n], n in [0..12]]; // G. C. Greubel, Jun 30 2021
    
  • Mathematica
    T[n_, k_, q_]:= (q*(3*q-2))^(k*(n-k)); Table[T[n, k, 2], {n, 0, 12}, {k, 0, n}]//Flatten
    With[{m=6}, Table[(m+2)^(k*(n-k)), {n,0,12}, {k,0,n}]//Flatten] (* G. C. Greubel, Jun 30 2021 *)
  • Sage
    flatten([[8^(k*(n-k)) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jun 30 2021

Formula

T(n, k, q) = c(n,q)/(c(k, q)*c(n-k, q)) where c(n, q) = (q*(3*q - 2))^binomial(n+1,2) and q = 2.
T(n, k, q) = (q*(3*q-2))^(k*(n-k)) with q = 2.
T(n, k) = 8^A004247(n,k), where A004247 is interpreted as a triangle. [relation detected by sequencedb.net]. - R. J. Mathar, Jun 30 2021
T(n, k, m) = (m+2)^(k*(n-k)) with m = 6. - G. C. Greubel, Jun 30 2021

Extensions

Edited by R. J. Mathar and G. C. Greubel, Jun 30 2021

A067138 OR-numbral multiplication table, read by antidiagonals.

Original entry on oeis.org

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

Views

Author

Jens Voß, Jan 02 2002

Keywords

Comments

See A048888 for the definition of OR-numbral arithmetic

Examples

			The top left 0..16 x 0..16 corner of the array:
  0,  0,  0,  0,  0,  0,  0,   0,   0,   0,   0,   0,   0,   0,   0,   0,
  0,  1,  2,  3,  4,  5,  6,   7,   8,   9,  10,  11,  12,  13,  14,  15,
  0,  2,  4,  6,  8, 10, 12,  14,  16,  18,  20,  22,  24,  26,  28,  30,
  0,  3,  6,  7, 12, 15, 14,  15,  24,  27,  30,  31,  28,  31,  30,  31,
  0,  4,  8, 12, 16, 20, 24,  28,  32,  36,  40,  44,  48,  52,  56,  60,
  0,  5, 10, 15, 20, 21, 30,  31,  40,  45,  42,  47,  60,  61,  62,  63,
  0,  6, 12, 14, 24, 30, 28,  30,  48,  54,  60,  62,  56,  62,  60,  62,
  0,  7, 14, 15, 28, 31, 30,  31,  56,  63,  62,  63,  60,  63,  62,  63,
  0,  8, 16, 24, 32, 40, 48,  56,  64,  72,  80,  88,  96, 104, 112, 120,
  0,  9, 18, 27, 36, 45, 54,  63,  72,  73,  90,  91, 108, 109, 126, 127,
  0, 10, 20, 30, 40, 42, 60,  62,  80,  90,  84,  94, 120, 122, 124, 126,
  0, 11, 22, 31, 44, 47, 62,  63,  88,  91,  94,  95, 124, 127, 126, 127,
  0, 12, 24, 28, 48, 60, 56,  60,  96, 108, 120, 124, 112, 124, 120, 124,
  0, 13, 26, 31, 52, 61, 62,  63, 104, 109, 122, 127, 124, 125, 126, 127,
  0, 14, 28, 30, 56, 62, 60,  62, 112, 126, 124, 126, 120, 126, 124, 126,
  0, 15, 30, 31, 60, 63, 62,  63, 120, 127, 126, 127, 124, 127, 126, 127,
  0, 16, 32, 48, 64, 80, 96, 112, 128, 144, 160, 176, 192, 208, 224, 240,
.
Multiplying 3 ("11" in binary) with itself in this system means taking bitwise-or of "11" with itself, when shifted one bit-position left:
       11
      110
  -------
OR:   111 = 7 in decimal = A(3,3).
.
Multiplying 10 (= "1010" in binary) and 11 (= "1011" in binary) in this system means taking bitwise-or of binary number 1011 when shifted once left with the same binary number when shifted three bit-positions left:
      10110
    1011000
    -------
OR: 1011110 = 94 in decimal = A(10,11) = A(11,10).
		

Crossrefs

Cf. A003986, A067139, A048888, A007059, A067398 (main diagonal).
Cf. also A004247, A048720 for analogous multiplication tables.

Programs

  • PARI
    t(n, k) = {res = 0; for (i=0, length(binary(n))-1, if (bittest(n, i), res = bitor(res, shift(k, i)));); return (res);} \\ Michel Marcus, Apr 14 2013

Formula

From Rémy Sigrist, Mar 17 2021: (Start)
T(n, 0) = 0.
T(n, 1) = n.
T(n, 2^k) = n*2^k for any k >= 0.
T(n, n) = A067398(n).
(End)
For all n, k: A048720(n,k) <= A(n,k) <= A004247(n,k). - Antti Karttunen, Mar 17 2021

Extensions

Example-section rewritten by Antti Karttunen, Mar 17 2021

A325820 Multiplication table for carryless product i X j in base 3 for i >= 0 and j >= 0, read by antidiagonals.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 1, 3, 0, 0, 4, 6, 6, 4, 0, 0, 5, 8, 9, 8, 5, 0, 0, 6, 7, 12, 12, 7, 6, 0, 0, 7, 3, 15, 16, 15, 3, 7, 0, 0, 8, 5, 18, 11, 11, 18, 5, 8, 0, 0, 9, 4, 21, 24, 13, 24, 21, 4, 9, 0, 0, 10, 18, 24, 19, 21, 21, 19, 24, 18, 10, 0, 0, 11, 20, 27, 23, 26, 9, 26, 23, 27, 20, 11, 0, 0, 12, 19, 30, 36, 19, 15, 15, 19, 36, 30, 19, 12, 0
Offset: 0

Views

Author

Antti Karttunen, May 22 2019

Keywords

Examples

			The array begins as:
  0,  0,  0,  0,  0,  0,  0,  0,  0,   0,   0,   0,   0, ...
  0,  1,  2,  3,  4,  5,  6,  7,  8,   9,  10,  11,  12, ...
  0,  2,  1,  6,  8,  7,  3,  5,  4,  18,  20,  19,  24, ...
  0,  3,  6,  9, 12, 15, 18, 21, 24,  27,  30,  33,  36, ...
  0,  4,  8, 12, 16, 11, 24, 19, 23,  36,  40,  44,  48, ...
  0,  5,  7, 15, 11, 13, 21, 26, 19,  45,  50,  52,  33, ...
  0,  6,  3, 18, 24, 21,  9, 15, 12,  54,  60,  57,  72, ...
  0,  7,  5, 21, 19, 26, 15, 13, 11,  63,  70,  68,  57, ...
  0,  8,  4, 24, 23, 19, 12, 11, 16,  72,  80,  76,  69, ...
  0,  9, 18, 27, 36, 45, 54, 63, 72,  81,  90,  99, 108, ...
  0, 10, 20, 30, 40, 50, 60, 70, 80,  90, 100,  83, 120, ...
  0, 11, 19, 33, 44, 52, 57, 68, 76,  99,  83,  91, 132, ...
  0, 12, 24, 36, 48, 33, 72, 57, 69, 108, 120, 132, 144, ...
  etc.
A(2,2) = 2*2 mod 3 = 1.
		

Crossrefs

Cf. A169999 (the main diagonal).
Row/Column 0: A000004, Row/Column 1: A001477, Row/Column 2: A004488, Row/Column 3: A008585, Row/Column 4: A242399, Row/Column 9: A008591.
Cf. A325821 (same table without the zero row and column).
Cf. A048720 (binary), A059692 (decimal), A004247 (full multiply).

Programs

  • PARI
    up_to = 105;
    A325820sq(b, c) = fromdigits(Vec(Pol(digits(b,3))*Pol(digits(c,3)))%3, 3);
    A325820list(up_to) = { my(v = vector(up_to), i=0); for(a=0,oo, for(col=0,a, if(i++ > up_to, return(v)); v[i] = A325820sq(a-col,col))); (v); };
    v325820 = A325820list(up_to);
    A325820(n) = v325820[1+n];

A279777 Numbers k such that the sum of digits of 9k is 27.

Original entry on oeis.org

111, 211, 221, 222, 311, 321, 322, 331, 332, 333, 411, 421, 422, 431, 432, 433, 441, 442, 443, 444, 511, 521, 522, 531, 532, 533, 541, 542, 543, 544, 551, 552, 553, 554, 555, 611, 621, 622, 631, 632, 633, 641, 642, 643, 644, 651, 652, 653, 654, 655, 661
Offset: 1

Views

Author

M. F. Hasler, Dec 23 2016

Keywords

Comments

The digital sum of 9k is always a multiple of 9. For most numbers below 100 it is actually equal to 9. Numbers such that the digital sum of 9k is 18 are listed in A279769. Only every third term of the present sequence is divisible by 3.
The sequence of record gaps [and upper end of the gap] is: 100 [a(2) = 211], 101 [a(221) = 1211], 111 [a(4841) = 11211], 111 [a(10121) = 22311], 111 [a(15752) = 33411], ..., 111 [a(45133) = 88911], 111 [a(50413) = 100011], 211 [a(55253) = 111211], 311 [a(110000) = 222311], ..., 911 [a(380557) = 888911], 1011 [a(411049) = 1000011], 1211 [a(436976) = 1111211], 2311 [a(840281) = 2222311], ..., 8911 [a(2451241) = 8888911], ...

Crossrefs

Cf. A008591, A084854, A003991, A004247, A279769 (sumdigits(9n) = 18).
Digital sum of m*n equals m: A088404 = A069537/2, A088405 = A052217/3, A088406 = A063997/4, A088407 = A069540/5, A088408 = A062768/6, A088409 = A063416/7, A088410 = A069543/8.
Numbers with given digital sum: A011557 (1), A052216 (2), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225 (14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A007953 (digital sum), A005349 (Niven or Harshad numbers), A245062 (arranged in rows by digit sums).
Cf. A082259.

Programs

  • Mathematica
    Select[Range@ 661, Total@ IntegerDigits[9 #] == 27 &] (* Michael De Vlieger, Dec 23 2016 *)
  • PARI
    is(n)=sumdigits(9*n)==27

A003988 Triangle with subscripts (1,1),(2,1),(1,2),(3,1),(2,2),(1,3), etc. in which entry (i,j) is [ i/j ].

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Another version of A010766.

Crossrefs

Row sums are in A006218. Antidiagonal sums are in A002541.

Programs

  • Haskell
    a003988 n k = (n + 1 - k) `div` k
    a003988_row n = zipWith div [n,n-1..1] [1..n]
    a003988_tabl = map a003988_row [1..]
    -- Reinhard Zumkeller, Apr 13 2012
  • Mathematica
    t[n_, k_] := Quotient[n, k]; Table[t[n-k+1, k], {n, 1, 14}, {k, 1, n}] // Flatten (* Jean-François Alcover, Nov 21 2013 *)

Formula

From Franklin T. Adams-Watters, Jan 28 2006: (Start)
T(n,k) = Sum_{i=1..k} A077049(n,i).
G.f.: (1/(1-x))*Sum_{k>0} x^k*y^k/(1-x^k) = (1/(1-x))*Sum_{k>0} x^k * y / (1 - x^k y) = (1/(1-x)) * Sum_{k>0} x^k * Sum_{d|k} y^d = A(x,y)/(1-x) where A(x,y) is the g.f. of A077049. (End)
T(n,k) = floor((n + 1 - k) / k). - Reinhard Zumkeller, Apr 13 2012

Extensions

More terms from James Sellers

A353109 Array read by antidiagonals: A(n, k) is the digital root of n*k with n >= 0 and k >= 0.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 4, 3, 0, 0, 4, 6, 6, 4, 0, 0, 5, 8, 9, 8, 5, 0, 0, 6, 1, 3, 3, 1, 6, 0, 0, 7, 3, 6, 7, 6, 3, 7, 0, 0, 8, 5, 9, 2, 2, 9, 5, 8, 0, 0, 9, 7, 3, 6, 7, 6, 3, 7, 9, 0, 0, 1, 9, 6, 1, 3, 3, 1, 6, 9, 1, 0, 0, 2, 2, 9, 5, 8, 9, 8, 5, 9, 2, 2, 0
Offset: 0

Views

Author

Stefano Spezia, Apr 24 2022

Keywords

Examples

			The array begins:
    0, 0, 0, 0, 0, 0, 0, 0, ...
    0, 1, 2, 3, 4, 5, 6, 7, ...
    0, 2, 4, 6, 8, 1, 3, 5, ...
    0, 3, 6, 9, 3, 6, 9, 3, ...
    0, 4, 8, 3, 7, 2, 6, 1, ...
    0, 5, 1, 6, 2, 7, 3, 8, ...
    0, 6, 3, 9, 6, 3, 9, 6, ...
    0, 7, 5, 3, 1, 8, 6, 4, ...
    ...
		

Crossrefs

Cf. A003991, A004247, A010888, A056992 (diagonal), A073636, A139413, A180592, A180593, A180594, A180595, A180596, A180597, A180598, A180599, A303296, A336225, A353128 (antidiagonal sums), A353933, A353974 (partial sum of the main diagonal).

Programs

  • Mathematica
    A[i_,j_]:=If[i*j==0,0,1+Mod[i*j-1,9]];Flatten[Table[A[n-k,k],{n,0,12},{k,0,n}]]
  • PARI
    T(n,k) = if (n && k, (n*k-1)%9+1, 0); \\ Michel Marcus, May 12 2022

Formula

A(n, k) = A010888(A004247(n, k)).
A(n, k) = A010888(A003991(n, k)) for n*k > 0.

A336225 Table read by antidiagonals: T(n, k) = digitsum(n*k) with n >= 0 and k >= 0.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 4, 3, 0, 0, 4, 6, 6, 4, 0, 0, 5, 8, 9, 8, 5, 0, 0, 6, 1, 3, 3, 1, 6, 0, 0, 7, 3, 6, 7, 6, 3, 7, 0, 0, 8, 5, 9, 2, 2, 9, 5, 8, 0, 0, 9, 7, 3, 6, 7, 6, 3, 7, 9, 0, 0, 1, 9, 6, 10, 3, 3, 10, 6, 9, 1, 0, 0, 2, 2, 9, 5, 8, 9, 8, 5, 9, 2, 2, 0
Offset: 0

Views

Author

Stefano Spezia, Jul 12 2020

Keywords

Examples

			The table T(n, k) begins
0   0   0   0   0   0   0   0 ...
0   1   2   3   4   5   6   7 ...
0   2   4   6   8   1   3   5 ...
0   3   6   9   3   6   9   3 ...
0   4   8   3   7   2   6  10 ...
0   5   1   6   2   7   3   8 ...
0   6   3   9   6   3   9   6 ...
0   7   5   3  10   8   6  13 ...
...
		

Crossrefs

Cf. A003991, A004092, A004159 (diagonal), A004164 (digitsum of n^3), A004247, A007953, A055565 (digitsum of n^4), A055566 (digitsum of n^5), A055567 (digitsum of n^6).

Programs

  • Mathematica
    T[n_,k_]:=Total[IntegerDigits[n*k]]; Table[T[n-k,k],{n,0,12},{k,0,n}]//Flatten
  • PARI
    T(n, k) = sumdigits(n*k);

Formula

T(n, k) = A007953(A004247(n, k)).
T(n, 1) = T(1, n) = A007953(n).
T(n, 2) = T(2, n) = A004092(n).
T(n, k) = A007953(A003991(n, k)) for n*k > 0. - Michel Marcus, Jul 13 2020.
Previous Showing 11-20 of 28 results. Next