A273102
Difference table of the divisors of the positive integers.
Original entry on oeis.org
1, 1, 2, 1, 1, 3, 2, 1, 2, 4, 1, 2, 1, 1, 5, 4, 1, 2, 3, 6, 1, 1, 3, 0, 2, 2, 1, 7, 6, 1, 2, 4, 8, 1, 2, 4, 1, 2, 1, 1, 3, 9, 2, 6, 4, 1, 2, 5, 10, 1, 3, 5, 2, 2, 0, 1, 11, 10, 1, 2, 3, 4, 6, 12, 1, 1, 1, 2, 6, 0, 0, 1, 4, 0, 1, 3, 1, 2, 1, 1, 13, 12, 1, 2, 7, 14, 1, 5, 7, 4, 2, -2, 1, 3, 5, 15, 2, 2, 10, 0, 8, 8
Offset: 1
For n = 18 the divisors of 18 are 1, 2, 3, 6, 9, 18, so the difference triangle of the divisors of 18 is
1 . 2 . 3 . 6 . 9 . 18
1 . 1 . 3 . 3 . 9
0 . 2 . 0 . 6
2 .-2 . 6
-4 . 8
12
and the 18th slice is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2,-2, 6;
-4, 8;
12;
The tetrahedron begins:
1;
1, 2;
1;
1, 3;
2;
1, 2, 4;
1, 2;
1;
...
This is also an irregular triangle T(n,r) read by rows in which row n lists the difference triangle of the divisors of n flattened. Row lengths are the terms of A184389. Row sums give A273103.
Triangle begins:
1;
1, 2, 1;
1, 3, 2;
1, 2, 4, 1, 2, 1;
...
-
Table[Drop[FixedPointList[Differences, Divisors@ n], -2], {n, 15}] // Flatten (* Michael De Vlieger, May 16 2016 *)
-
def A273102_DTD(n): # DTD = Difference Table of Divisors
D = divisors(n)
T = matrix(ZZ, len(D))
for (m, d) in enumerate(D):
T[0, m] = d
for k in range(m-1, -1, -1) :
T[m-k, k] = T[m-k-1, k+1] - T[m-k-1, k]
return [T.row(k)[:len(D)-k] for k in range(len(D))]
# Keeps the rows of the DTD, for instance
# A273102_DTD(18)[1] = 1,1,3,3,9 (see the example above).
for n in range(1,19): print(A273102_DTD(n)) # Peter Luschny, May 18 2016
A187215
Sum of the elements of the absolute difference table of the divisors of n.
Original entry on oeis.org
1, 4, 6, 11, 10, 21, 14, 26, 25, 31, 22, 52, 26, 45, 54, 57, 34, 82, 38, 82, 72, 73, 46, 119, 71, 87, 90, 108, 58, 161, 62, 120, 108, 115, 134, 181, 74, 129, 126, 193, 82, 221, 86, 172, 218, 157, 94, 252, 141, 190, 162, 204, 106, 285, 202, 233
Offset: 1
For n = 14 the divisors of 14 are 1, 2, 7, 14, and the absolute difference triangle of the divisors is
1 . 2 . 7 . 14
. 1 . 5 . 7
. . 4 . 2
. . . 2
The sum of all elements of the triangle is 1 + 2 + 7 + 14 + 1 + 5 + 7 + 4 + 2 + 2 = 45, so a(14) = 45.
-
with(numtheory):
DD:= l-> [seq(abs(l[i]-l[i-1]), i=2..nops(l))]:
a:= proc(n) local l;
l:= sort([divisors(n)[]], `>`);
add(j, j=[seq((DD@@i)(l)[], i=0..nops(l)-1)]);
end:
seq(a(n), n=1..100); # Alois P. Heinz, Aug 02 2011
-
Table[Total@ Flatten@ NestWhileList[Abs@ Differences@ # &, Divisors@ n, Length@ # > 1 &], {n, 56}] (* Michael De Vlieger, May 18 2016 *)
A273109
Numbers n such that in the difference triangle of the divisors of n (including the divisors of n) the diagonal from the bottom entry to n gives the divisors of n.
Original entry on oeis.org
1, 2, 4, 8, 12, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824, 2147483648, 4294967296, 8589934592
Offset: 1
For n = 12 the difference triangle of the divisors of 12 is
1 . 2 . 3 . 4 . 6 . 12
. 1 . 1 . 1 . 2 . 6
. . 0 . 0 . 1 . 4
. . . 0 . 1 . 3
. . . . 1 . 2
. . . . . 1
The bottom entry is 1, and the diagonal from the bottom entry to 12 is [1, 2, 3, 4, 6, 12] hence the diagonal gives the divisors of 12, so 12 is in the sequence.
Note that for n = 12 and the powers of 2 the descending diagonals, from left to right, are symmetrics, for example: the first diagonal is 1, 1, 0, 0, 1, 1.
-
aQ[n_] := Module[{d=Divisors[n]}, nd = Length[d]; vd = d; ans = True; Do[ vd = Differences[vd]; If[Max[vd] != d[[nd-k]], ans=False; Break[]], {k,1,nd-1}]; ans]; Select[Range[100000], aQ] (* Amiram Eldar, Feb 23 2019 *)
-
isok(n) = {my(d = divisors(n)); my(nd = #d); my(vd = d); for (k=1, nd-1, vd = vector(#vd-1, j, vd[j+1] - vd[j]); if (vecmax(vd) != d[nd-k], return (0));); return (1);} \\ Michel Marcus, May 16 2016
A273135
Difference table of the divisors of the positive integers (with every table read by antidiagonals downwards).
Original entry on oeis.org
1, 1, 2, 1, 1, 3, 2, 1, 2, 1, 4, 2, 1, 1, 5, 4, 1, 2, 1, 3, 1, 0, 6, 3, 2, 2, 1, 7, 6, 1, 2, 1, 4, 2, 1, 8, 4, 2, 1, 1, 3, 2, 9, 6, 4, 1, 2, 1, 5, 3, 2, 10, 5, 2, 0, 1, 11, 10, 1, 2, 1, 3, 1, 0, 4, 1, 0, 0, 6, 2, 1, 1, 1, 12, 6, 4, 3, 2, 1, 1, 13, 12, 1, 2, 1, 7, 5, 4, 14, 7, 2, -2, 1, 3, 2, 5, 2, 0, 15, 10, 8, 8
Offset: 1
The tables of the first nine positive integers are
1; 1, 2; 1, 3; 1, 2, 4; 1, 5; 1, 2, 3, 6; 1, 7; 1, 2, 4, 8; 1, 3, 9;
1; 2; 1, 2; 4; 1, 1, 3; 6; 1, 2, 4; 2, 6;
1; 0, 2; 1, 2; 4;
2; 1;
For n = 18 the difference table of the divisors of 18 is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2, -2, 6;
-4, 8;
12;
This table read by antidiagonals downwards gives the finite subsequence [1], [2, 1], [3, 1, 0], [6, 3, 2, 2], [9, 3, 0, -2, -4], [18, 9, 6, 6, 8, 12].
Cf.
A000005,
A000217,
A027750,
A161700,
A184389,
A187202,
A272210,
A272121,
A273102,
A273103,
A273262.
-
Table[Table[#[[m - k + 1, k]], {m, Length@ #}, {k, m, 1, -1}] &@ NestWhileList[Differences, Divisors@ n, Length@ # > 1 &], {n, 15}] // Flatten (* Michael De Vlieger, Jun 26 2016 *)
A272210
Difference table of the divisors of the positive integers (with every table read by antidiagonals upwards).
Original entry on oeis.org
1, 1, 1, 2, 1, 2, 3, 1, 1, 2, 1, 2, 4, 1, 4, 5, 1, 1, 2, 0, 1, 3, 2, 2, 3, 6, 1, 6, 7, 1, 1, 2, 1, 2, 4, 1, 2, 4, 8, 1, 2, 3, 4, 6, 9, 1, 1, 2, 2, 3, 5, 0, 2, 5, 10, 1, 10, 11, 1, 1, 2, 0, 1, 3, 0, 0, 1, 4, 1, 1, 1, 2, 6, 1, 2, 3, 4, 6, 12, 1, 12, 13, 1, 1, 2, 4, 5, 7, -2, 2, 7, 14, 1, 2, 3, 0, 2, 5, 8, 8, 10, 15
Offset: 1
The tables of the first nine positive integers are
1; 1, 2; 1, 3; 1, 2, 4; 1, 5; 1, 2, 3, 6; 1, 7; 1, 2, 4, 8; 1, 3, 9;
. 1; 2; 1, 2; 4; 1, 1, 3; 6; 1, 2, 4; 2, 6;
. 1; 0, 2; 1, 2; 4;
. 2; 1;
.
For n = 18 the difference table of the divisors of 18 is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2, -2, 6;
-4, 8;
12;
This table read by antidiagonals upwards gives the finite subsequence [1], [1, 2], [0, 1, 3], [2, 2, 3, 6], [-4, -2, 0, 3, 9], [12, 8, 6, 6, 9, 18].
Cf.
A000005,
A000217,
A027750,
A161700,
A184389,
A187202,
A273102,
A273103,
A273109,
A273135,
A273132,
A273136,
A273261,
A273262,
A273263.
-
Table[Table[#[[m - k + 1, k]], {m, Length@ #}, {k, m}] &@ NestWhileList[Differences, Divisors@ n, Length@ # > 1 &], {n, 15}] // Flatten (* Michael De Vlieger, Jun 29 2016 *)
A273130
Numbers which have only positive entries in the difference table of their divisors.
Original entry on oeis.org
1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 21, 23, 25, 27, 29, 31, 32, 33, 37, 39, 41, 43, 47, 49, 51, 53, 55, 57, 59, 61, 64, 65, 67, 69, 71, 73, 79, 81, 83, 85, 87, 89, 93, 95, 97, 101, 103, 107, 109, 111, 113, 115, 119, 121, 123, 125, 127, 128, 129, 131, 133
Offset: 1
85 is in the sequence because the difference table of the divisors of 85 has only entries greater than 0:
[1, 5, 17, 85]
[4, 12, 68]
[8, 56]
[48]
-
Select[Range@ 1000, {} == NestWhile[ Differences, Divisors @ #, # != {} && Min[#] > 0 &] &] (* Giovanni Resta, May 16 2016 *)
-
has(v)=if(#v<2, v[1]>0, if(vecmin(v)<1, 0, has(vector(#v-1,i,v[i+1]-v[i]))))
is(n)=has(divisors(n)) \\ Charles R Greathouse IV, May 16 2016
-
def sf(z):
D = divisors(z)
T = matrix(ZZ, len(D))
for m, d in enumerate(D):
T[0, m] = d
for k in range(m-1, -1, -1) :
T[m-k, k] = T[m-k-1, k+1] - T[m-k-1, k]
if T[m-k, k] <= 0: return False
return True
print([z for z in range(1,100) if sf(z)])
A273262
Irregular triangle read by rows: T(n,k) = sum of the elements of the k-th antidiagonal of the difference table of the divisors of n.
Original entry on oeis.org
1, 1, 3, 1, 5, 1, 3, 7, 1, 9, 1, 3, 4, 13, 1, 13, 1, 3, 7, 15, 1, 5, 19, 1, 3, 10, 17, 1, 21, 1, 3, 4, 5, 11, 28, 1, 25, 1, 3, 16, 21, 1, 5, 7, 41, 1, 3, 7, 15, 31, 1, 33, 1, 3, 4, 13, 6, 59, 1, 37, 1, 3, 7, 3, 31, 21, 1, 5, 13, 53, 1, 3, 28, 29, 1, 45, 1, 3, 4, 5, 11, 4, 36, 39, 1, 9, 61, 1, 3, 34, 33, 1, 5, 19, 65
Offset: 1
Triangle begins:
1;
1, 3;
1, 5;
1, 3, 7;
1, 9;
1, 3, 4, 13;
1, 13;
1, 3, 7, 15;
1, 5, 19;
1, 3, 10, 17;
1, 21;
1, 3, 4, 5, 11, 28;
1, 25;
1, 3, 16, 21;
1, 5, 7, 41;
1, 3, 7, 15, 31;
1, 33;
1, 3, 4, 13, 6, 59;
1, 37;
1, 3, 7, 3, 31, 21;
1, 5, 13, 53;
1, 3, 28, 29;
1, 45;
1, 3, 4, 5, 11, 4, 36, 39;
1, 9, 61;
1, 3, 34, 33;
1, 5, 19, 65;
...
For n = 18 the divisors of 18 are 1, 2, 3, 6, 9, 18, and the difference triangle of the divisors is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2, -2, 6;
-4, 8;
12;
The antidiagonal sums give [1, 3, 4, 13, 6, 59] which is also the 18th row of the irregular triangle.
-
Table[Map[Total, Table[#[[m - k + 1, k]], {m, Length@ #}, {k, m}], {1}] &@ NestWhileList[Differences, Divisors@ n, Length@ # > 1 &], {n, 27}] (* Michael De Vlieger, Jun 26 2016 *)
-
row(n) = {my(d = divisors(n)); my(nd = #d); my(m = matrix(#d, #d)); for (j=1, nd, m[1,j] = d[j];); for (i=2, nd, for (j=1, nd - i +1, m[i,j] = m[i-1,j+1] - m[i-1,j];);); vector(nd, i, sum(k=0, i-1, m[i-k, k+1]));}
tabf(nn) = for (n=1, nn, print(row(n)););
lista(nn) = for (n=1, nn, v = row(n); for (j=1, #v, print1(v[j], ", "));); \\ Michel Marcus, Jun 25 2016
A273263
Irregular triangle read by rows: T(n,k) is the sum of the elements of the k-th column of the difference table of the divisors of n.
Original entry on oeis.org
1, 2, 2, 3, 3, 3, 4, 4, 5, 5, 4, 5, 6, 6, 7, 7, 4, 6, 8, 8, 7, 9, 9, 4, 7, 10, 10, 11, 11, 4, 6, 8, 10, 12, 12, 13, 13, 4, 9, 14, 14, 11, 13, 15, 15, 5, 8, 12, 16, 16, 17, 17, 12, 11, 12, 15, 18, 18, 19, 19, -3, 4, 10, 15, 20, 20, 13, 17, 21, 21, 4, 13, 22, 22, 23, 23, -4, 3, 8, 12, 16, 20, 24, 24, 21, 25, 25
Offset: 1
Triangle begins:
1;
2, 2;
3, 3;
3, 4, 4;
5, 5;
4, 5, 6, 6;
7, 7;
4, 6, 8, 8;
7, 9, 9;
4, 7, 10, 10;
11, 11;
4, 6, 8, 10, 12, 12;
13, 13;
4, 9, 14, 14;
11, 13, 15, 15;
5, 8, 12, 16, 16;
17, 17;
12, 11, 12, 15, 18, 18;
19, 19;
-3, 4, 10, 15, 20, 20;
13, 17, 21, 21;
4, 13, 22, 22;
23, 23;
-4, 3, 8, 12, 16, 20, 24, 24;
21, 25, 25;
4, 15, 26, 26;
...
For n = 18 the divisors of 18 are 1, 2, 3, 6, 9, 18, and the difference triangle of the divisors is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2, -2, 6;
-4, 8;
12;
The column sums give [12, 11, 12, 15, 18, 18] which is also the 18th row of the irregular triangle.
-
Table[Total /@ Transpose@ Map[Function[w, PadRight[w, Length@ #]], NestWhileList[Differences, #, Length@ # > 1 &]] &@ Divisors@ n, {n, 25}] // Flatten (* Michael De Vlieger, Jun 26 2016 *)
-
row(n) = {my(d = divisors(n)); my(nd = #d); my(m = matrix(#d, #d)); for (j=1, nd, m[1,j] = d[j];); for (i=2, nd, for (j=1, nd - i +1, m[i,j] = m[i-1,j+1] - m[i-1,j];);); vector(nd, j, sum(i=1, nd, m[i, j]));}
tabf(nn) = for (n=1, nn, print(row(n)););
lista(nn) = for (n=1, nn, v = row(n); for (j=1, #v, print1(v[j], ", "));); \\ Michel Marcus, Jun 25 2016
A273261
Irregular triangle read by rows: T(n,k) = sum of the elements of the k-th row of the difference table of the divisors of n.
Original entry on oeis.org
1, 3, 1, 4, 2, 7, 3, 1, 6, 4, 12, 5, 2, 2, 8, 6, 15, 7, 3, 1, 13, 8, 4, 18, 9, 4, 0, 12, 10, 28, 11, 5, 4, 3, 1, 14, 12, 24, 13, 6, -2, 24, 14, 8, 8, 31, 15, 7, 3, 1, 18, 16, 39, 17, 8, 6, 4, 12, 20, 18, 42, 19, 9, 4, 3, -11, 32, 20, 12, 8, 36, 21, 10, -6, 24, 22, 60, 23, 11, 8, 6, 3, 4, -12, 31, 24, 16, 42, 25, 12, -8
Offset: 1
Triangle begins:
1;
3, 1;
4, 2;
7, 3, 1;
6, 4;
12, 5, 2, 2;
8, 6;
15, 7, 3, 1;
13, 8, 4;
18, 9, 4, 0;
12, 10;
28, 11, 5, 4, 3, 1;
14, 12;
24, 13, 6, -2;
24, 14, 8, 8;
31, 15, 7, 3, 1;
18, 16;
39, 17, 8, 6, 4, 12;
20, 18;
42, 19, 9, 4, 3, -11;
32, 20, 12, 8;
36, 21, 10, -6;
24, 22;
60, 23, 11, 8, 6, 3, 4, -12;
31, 24, 16;
42, 25, 12, -8;
...
For n = 14 the divisors of 14 are 1, 2, 7, 14, and the difference triangle of the divisors is
1, 2, 7, 14;
1, 5, 7;
4, 2;
-2;
The row sums give [24, 13, 6, -2] which is also the 14th row of the irregular triangle.
In the first row, the last element is 14, the first is 1. So the sum of the second row is 14 - 1 is 13. Similarly, the sum of the third row is 7 - 1 = 6, and of the last row, 2 - 4 = -2. - _David A. Corneth_, Jun 25 2016
-
Map[Total, Table[NestWhileList[Differences, Divisors@ n, Length@ # > 1 &], {n, 26}], {2}] // Flatten (* Michael De Vlieger, Jun 26 2016 *)
-
row(n) = {my(d = divisors(n));my(nd = #d); my(m = matrix(#d, #d)); for (j=1, nd, m[1,j] = d[j];); for (i=2, nd, for (j=1, nd - i +1, m[i,j] = m[i-1,j+1] - m[i-1,j];);); vector(nd, i, sum(j=1, nd, m[i, j]));}
tabf(nn) = for (n=1, nn, print(row(n)););
lista(nn) = for (n=1, nn, v = row(n); for (j=1, #v, print1(v[j], ", "));); \\ Michel Marcus, Jun 25 2016
A273136
Difference table of the divisors of the positive integers (with every table read by columns).
Original entry on oeis.org
1, 1, 1, 2, 1, 2, 3, 1, 1, 1, 2, 2, 4, 1, 4, 5, 1, 1, 0, 2, 2, 1, 2, 3, 3, 6, 1, 6, 7, 1, 1, 1, 1, 2, 2, 2, 4, 4, 8, 1, 2, 4, 3, 6, 9, 1, 1, 2, 0, 2, 3, 2, 5, 5, 10, 1, 10, 11, 1, 1, 0, 0, 1, 1, 2, 1, 0, 1, 2, 3, 1, 1, 3, 4, 2, 4, 6, 6, 12, 1, 12, 13, 1, 1, 4, -2, 2, 5, 2, 7, 7, 14, 1, 2, 0, 8, 3, 2, 8, 5, 10, 15
Offset: 1
The tables of the first nine positive integers are
1; 1, 2; 1, 3; 1, 2, 4; 1, 5; 1, 2, 3, 6; 1, 7; 1, 2, 4, 8; 1, 3, 9;
. 1; 2; 1, 2; 4; 1, 1, 3; 6; 1, 2, 4; 2, 6;
. 1; 0, 2; 1, 2; 4;
. 2; 1;
.
For n = 18 the difference table of the divisors of 18 is
1, 2, 3, 6, 9, 18;
1, 1, 3, 3, 9;
0, 2, 0, 6;
2, -2, 6;
-4, 8;
12;
This table read by columns gives the finite subsequence [1, 1, 0, 2, -4, 12], [2, 1, 2, -2, 8], [3, 3, 0, 6], [6, 3, 6], [9, 9], [18].
Cf.
A000005,
A000217,
A027750,
A184389,
A187202,
A272210,
A273102,
A273103,
A273135,
A273137,
A273263.
-
Table[Transpose@ Map[Function[w, PadRight[w, Length@ #]], NestWhileList[Differences, #, Length@ # > 1 &]] &@ Divisors@ n, {n, 15}] /. 0 -> {} // Flatten (* Michael De Vlieger, Jun 26 2016 *)
Showing 1-10 of 13 results.
Comments