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 32 results. Next

A272619 Irregular array read by rows: n-th row contains (in ascending order) the numbers 1 <= k < n such that at least one prime divisor p of k also divides n and at least one prime divisor q of k is coprime to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 6, 6, 6, 0, 10, 0, 6, 10, 12, 6, 10, 12, 6, 10, 12, 14, 0, 10, 14, 15, 0, 6, 12, 14, 15, 18, 6, 12, 14, 15, 18, 6, 10, 12, 14, 18, 20, 0, 10, 14, 15, 20, 21, 22, 10, 15, 20, 6, 10, 12, 14, 18, 20, 22, 24, 6, 12, 15, 18, 21, 24, 6, 10, 12, 18, 20, 21, 22, 24, 26, 0, 14, 21, 22
Offset: 1

Views

Author

Michael De Vlieger, May 03 2016

Keywords

Comments

The k are the "semitotatives" of n as counted by A243823(n).
All nonzero terms k are composite and pertain to composite rows n. This is because prime k must either divide or be coprime to n, and k = 1 is both a divisor of and coprime to n. Further, the terms k must have at least two distinct prime divisors p and q.
Row n for prime p contains zero, since numbers 1 <= k < p must either divide or be coprime to prime p.
Row n for prime powers p^e contains all the numbers k in the corresponding row of A133995. There is only one prime divisor p of p^e and every power 1 <= m <= e of p divides p^e, thus none of the terms of the corresponding row of A133995 are in A272618(n).
Rows n = 4 and 6 are special cases of composite n that contains zero. 4 is the smallest composite number; there are no composites k < n. 6 has the prime divisors 2 and 3, thus 5 is the smallest prime coprime to 6; the product of the minimum prime divisor and minimum prime coprime to 6 is 10, which exceeds 6 and falls outside the considered range. The situation is not so for composite n > 6. Thus rows n for composite n > 6 contain at least 1 nonzero value.
The smallest k of row n = A096014(n) < n, i.e., those values of A096014(n) pertaining to composite n > 6, a product of the smallest prime divisor p of n and the smallest prime q coprime to n. The smallest k of n are even squarefree semiprimes since 2 either divides n or is coprime to n and k is by definition a number with at least two distinct primes. The smallest k = 2p for p^2 sets record values for A096014(n) when we ignore values pertaining to prime n, n = 4, and n = 6.
In base n, 1/a(n) has a mixed recurrent expansion.

Examples

			For n = 12, the numbers 1 <= k < n such that the prime divisors p of k also divide n are {2, 3, 4, 6, 8, 9}; {2, 3, 4, 6} divide n = 12, thus row n = 12 is {8, 9}.
n:   k
1:   0
2:   0
3:   0
4:   0
5:   0
6:   0
7:   0
8:   6
9:   6
10:  6
11:  0
12: 10
13:  0
14:  6 10 12
15:  6 10 12
16:  6 10 12 14
17:  0
18: 10 14 15
19:  0
20:  6 12 14 15 18
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, pp. 144-5, Theorem 136.

Crossrefs

The union of nonzero terms of a(n) and A272618 = A133995, thus A243822(n) + A243823(n) = A045763(n).

Programs

  • Mathematica
    Table[With[{r = First /@ FactorInteger@ n}, Select[Range@ n, Function[m, And[! SubsetQ[r, First /@ FactorInteger@ m], 1 < GCD[m, n] < n]]]], {n, 30}] /. {} -> {0} // Flatten (* Michael De Vlieger, May 03 2016 *)

A243823 Quantity of "semitotatives," numbers m < n that are products of at least one prime divisor p of n and one prime q coprime to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 3, 3, 4, 0, 3, 0, 5, 5, 6, 0, 6, 3, 8, 6, 9, 0, 5, 0, 11, 8, 11, 7, 11, 0, 13, 10, 14, 0, 12, 0, 16, 14, 17, 0, 18, 5, 19, 14, 20, 0, 21, 11, 22, 16, 23, 0, 19, 0, 25, 20, 26, 13, 25, 0, 27, 20, 27, 0, 31, 0, 30, 27, 31, 13, 32, 0, 35, 23, 34, 0, 33, 17, 36, 25, 38, 0, 35, 15, 39, 27, 40, 19, 45, 0, 44, 32, 46
Offset: 1

Views

Author

Michael De Vlieger, Jun 11 2014

Keywords

Comments

Semitotatives m < n have a regular factor that is the product of prime divisors of n, and a coprime factor that is the product of primes q that are coprime to n.
The unit fractions of semitotatives have a mixed recurrent expansion in base n (See Hardy & Wright).

Examples

			For n = 10 with prime divisors {2, 5} and prime totatives {3, 7}, the only semitotative is 6. For n = 16, with the prime divisor 2 and the prime totatives {3, 5, 7, 11, 13}, there are four semitotatives {6, 10, 12, 14}.
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Sixth Edition, Oxford University Press, 2008, pages 144-145 (last part of Theorem 136).

Crossrefs

Programs

  • Maple
    f:= n -> n + 1 - numtheory:-phi(n) - add(numtheory:-mobius(k)*floor(n/k), k=select(t -> igcd(n,t)=1, [$1..n])):
    map(f, [$1..100]); # Robert Israel, May 10 2016
  • Mathematica
    Table[n + 1 - EulerPhi@ n - Total[MoebiusMu[#] Floor[n/#] &@ Select[Range@ n, CoprimeQ[#, n] &]], {n, 120}] (* Michael De Vlieger, May 10 2016 *)

Formula

a(n) = A045763(n) - A243822(n).
a(n) = n + 1 - phi(n) - Sum_{1 <= k <= n, gcd(n, k) = 1} mu(k)*floor(n/k). - Michael De Vlieger, May 10 2016, after Benoit Cloitre at A010846.

A355432 a(n) = number of k < n such that rad(k) = rad(n) and k does not divide n, where rad(k) = A007947(k).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 2, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 4, 0, 2, 0, 2, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Michael De Vlieger, Feb 22 2023

Keywords

Comments

a(n) = 0 for prime powers and squarefree numbers.

Examples

			a(1) = 18, since 18/6 >= 3. We note that rad(12) = rad(18) = 6, yet 12 does not divide 18.
a(2) = 24, since 24/6 >= 3. rad(18) = rad(24) = 6 and 24 mod 18 = 6.
a(3) = 36, since 36/6 >= 3. rad(24) = rad(36) = 6 and 36 mod 24 = 12.
a(6) = 54, since 54/6 >= 3. m in {12, 24, 36, 48} are such that rad(m) = rad(54) = 6, but none divides 54, etc.
		

Crossrefs

Programs

  • Mathematica
    rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]]; Table[Which[PrimePowerQ[n], 0, SquareFreeQ[n], 0, True, r = rad[n]; Count[Select[Range[n], Nor[PrimePowerQ[#], SquareFreeQ[#]] &], _?(And[rad[#] == r, Mod[n, #] != 0] &)]], {n, 120}]
  • PARI
    rad(n) = factorback(factorint(n)[, 1]); \\ A007947
    a(n) = my(rn=rad(n)); sum(k=1, n-1, if (n % k, rad(k)==rn)); \\ Michel Marcus, Feb 23 2023

Formula

a(n) > 0 for n in A360768.
a(n) < A243822(n) < A010846(n).
a(n) = A008479(n) - A005361(n). - Amiram Eldar, Oct 25 2024

A361098 Intersection of A360765 and A360768.

Original entry on oeis.org

36, 48, 50, 54, 72, 75, 80, 96, 98, 100, 108, 112, 135, 144, 147, 160, 162, 189, 192, 196, 200, 216, 224, 225, 240, 242, 245, 250, 252, 270, 288, 294, 300, 320, 324, 336, 338, 350, 352, 360, 363, 375, 378, 384, 392, 396, 400, 405, 416, 432, 441, 448, 450, 468, 480, 484, 486, 490, 500, 504, 507, 525
Offset: 1

Views

Author

Michael De Vlieger, Mar 15 2023

Keywords

Comments

Numbers k that are neither prime powers nor squarefree, such that rad(k) * A053669(k) < k and k/rad(k) >= A119288(k), where rad(k) = A007947(k).
Numbers k such that A360480(k), A360543(k), A361235(k), and A355432(k) are positive.
Subset of A126706. All terms are neither prime powers nor squarefree.
From Michael De Vlieger, Aug 03 2023: (Start)
Superset of A286708 = A001694 \ {{1} U A246547}, which in turn is a superset of A303606. We may write k in A286708 as m*rad(k)^2, m >= 1. Since omega(k) > 1, it is clear both k/rad(k) > A053669(k) and k/rad(k) >= A119288(k). Also superset of A359280 = A286708 \ A303606.
This sequence contains {A002182 \ A168263}. (End)

Examples

			For prime p, A360480(p) = A360543(p) = A361235(p) = A355432(p) = 0, since k < p is coprime to p.
For prime power n = p^e > 4, e > 0, A360543(n) = p^(e-1) - e, but A360480(n) = A361235(n) = A355432(n) = 0, since the other sequences require omega(n) > 1.
For squarefree composite n, A360480(n) >= 1 and A361235(n) >= 1 (the latter for n > 6), but A360543(n) = A355432(n) = 0, since the other sequences require at least 1 prime power factor p^e | n with e > 0.
For n = 18, A360480(n) = | {10, 14, 15} | = 3,
            A360543(n) = | {} | = 0,
            A361235(n) = | {4, 8, 16} | = 3,
            A355432(n) = | {12} | = 1.
Therefore 18 is not in the sequence.
For n = 36, A360480(n) = | {10, 14, 15, 20, 21, 22, 26, 28, 33, 34} | = 10,
            A360543(n) = | {30} | = 1,
            A361235(n) = | {8, 16, 27, 32} | = 4,
            A355432(n) = | {24} | = 1.
Therefore 36 is the smallest term in the sequence.
Table pertaining to the first 12 terms:
Key: a = A360480, b = A360543, c = A243823; d = A361235, e = A355432, f = A243822;
g = A046753 = f + c, tau = A000005, phi = A000010.
    n |  a + b =  c | d + e = f | g + tau + phi - 1 =  n
  ------------------------------------------------------
   36 | 10 + 1 = 11 | 4 + 1 = 5 | 16 +  9 + 12 - 1 =  36
   48 | 16 + 2 = 18 | 3 + 2 = 5 | 23 + 10 + 16 - 1 =  48
   50 | 18 + 1 = 19 | 4 + 2 = 6 | 25 +  6 + 20 - 1 =  50
   54 | 19 + 2 = 21 | 4 + 4 = 8 | 29 +  8 + 18 - 1 =  54
   72 | 27 + 4 = 31 | 4 + 2 = 6 | 37 + 12 + 24 - 1 =  72
   75 | 25 + 2 = 27 | 2 + 1 = 3 | 30 +  6 + 40 - 1 =  75
   80 | 32 + 3 = 35 | 3 + 1 = 4 | 39 + 10 + 32 - 1 =  80
   96 | 38 + 7 = 45 | 4 + 4 = 8 | 53 + 12 + 32 - 1 =  96
   98 | 41 + 3 = 44 | 5 + 2 = 7 | 51 +  6 + 42 - 1 =  98
  100 | 42 + 4 = 46 | 4 + 2 = 6 | 52 +  9 + 40 - 1 = 100
  108 | 44 + 8 = 52 | 5 + 4 = 9 | 61 + 12 + 36 - 1 = 108
  112 | 48 + 3 = 51 | 3 + 1 = 4 | 55 + 10 + 48 - 1 = 112
		

Crossrefs

Programs

  • Mathematica
    nn = 2^16;
    a053669[n_] := If[OddQ[n], 2, p = 2; While[Divisible[n, p], p = NextPrime[p]]; p];
    s = Select[Range[nn], Nor[PrimePowerQ[#], SquareFreeQ[#]] &];
    Reap[ Do[n = s[[j]];
        If[And[#1*a053669[n] < n, n/#1 >= #2] & @@ {Times @@ #, #[[2]]} &@
          FactorInteger[n][[All, 1]], Sow[n]], {j, Length[s]}]][[-1, -1]]

A244974 Sum of numbers m <= n whose set of prime divisors is a subset of the set of prime divisors of n.

Original entry on oeis.org

1, 3, 4, 7, 6, 16, 8, 15, 13, 30, 12, 45, 14, 36, 33, 31, 18, 79, 20, 66, 41, 64, 24, 103, 31, 70, 40, 80, 30, 235, 32, 63, 84, 114, 73, 198, 38, 120, 92, 163, 42, 310, 44, 140, 130, 132, 48, 246, 57, 213, 108, 154, 54, 300, 97, 217, 116, 150, 60, 600, 62, 156, 180, 127, 109, 540, 68, 246
Offset: 1

Views

Author

Michael De Vlieger, Jul 08 2014

Keywords

Comments

a(n) = A000203(n) when n is prime or a perfect prime power (A000961). This is because all products of the prime divisor p in such numbers produce divisors.
a(n) > A000203(n) when n is composite and not a perfect prime power.

Examples

			For n = 4, A162306(4) = {1, 2, 4} and a(4) = 7.
For n = 5, A162306(5) = {1, 5} and a(5) = 6.
For n = 6, A162306(6) = {1, 2, 3, 4, 6} and a(6) = 16.
		

Crossrefs

a(n) = sum of terms of n-th row of triangle A162306(n,k).

Programs

  • Mathematica
    Table[Total@ Union[{1}, Function[d, Select[Range@ n, Union[d, First /@ FactorInteger@ #] == d &]][First /@ FactorInteger@ n]], {n, 68}] (* or *)
    Table[Sum[k (Floor[n^k/k] - Floor[(n^k - 1)/k]), {k, n}], {n, 68}] (* Michael De Vlieger, May 26 2016 *)
  • PARI
    a(n) = {summ = 0; spn = factor(n)[,1]~; for (m=1, n, spm = factor(m)[,1]~; if (setintersect(spm, spn) == spm, summ += m);); summ;} \\ Michel Marcus, Jul 17 2014

Formula

a(n) = Sum_{k=1..n} k*( floor(n^k/k)-floor((n^k - 1)/k) ). - Anthony Browne, May 25 2016
a(n) = Sum_{j=1..n} Sum_{i=j..gcd(n^j,j)} i. - Wesley Ivan Hurt, Apr 05 2021

A360543 a(n) = number of numbers k < n, gcd(k, n) > 1, such that omega(k) > omega(n) and rad(n) | rad(k), where omega(n) = A001221(n) and rad(n) = A007947(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 6, 0, 0, 0, 0, 11, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 5, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 0, 0, 0, 1, 26, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 2, 0, 0, 0, 0, 3, 23, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 7, 0, 3, 1, 4
Offset: 1

Views

Author

Michael De Vlieger, Mar 06 2023

Keywords

Examples

			a(4) = 0 since k = 1..3 are prime powers.
a(8) = 1 since only k = 6 is such that p = 3, q = 5, but gcd(6, 10) = 2.
a(9) = 1 since the following satisfies definition: {6},
a(16) = 4, i.e., {6, 10, 12, 14},
a(25) = 3, i.e., {10, 15, 20},
a(27) = 6, i.e., {6, 12, 15, 18, 21, 24},
a(32) = 11, i.e., {6, 10, 12, 14, 18, 20, 22, 24, 26, 28, 30},
a(36) = 1, i.e., {30},
a(40) = 1, i.e., {30},
a(45) = 1, i.e., {30}, etc.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]]; c = Select[Range[4, nn], CompositeQ]; Table[Function[{q, r}, Count[TakeWhile[c, # <= n &], _?(And[PrimeNu[#] > q, Divisible[rad[#], r]] &)]] @@ {PrimeNu[n], rad[n]}, {n, nn}]

Formula

a(n) = A243823(n) - A360480(n).
a(n) = A045763(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A000005(n) - A243822(n) - A360480(n).
a(n) = A051953(n) - A010846(n) - A360480(n).
a(n) = A243823(n) = A045763(n) for n in A246547.
For prime power n = p^e, n > 1, a(n) = p^(e-1) - e.
For n in A360765, a(n) > 0.

A275280 Irregular triangle listing numbers m of n that have prime divisors p that also divide n, in order of appearance in a matrix of products that arranges the powers of prime divisors p of n along independent axes.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 2, 4, 1, 5, 1, 2, 4, 3, 6, 1, 7, 1, 2, 4, 8, 1, 3, 9, 1, 2, 4, 8, 5, 10, 1, 11, 1, 2, 4, 8, 3, 6, 12, 9, 1, 13, 1, 2, 4, 8, 7, 14, 1, 3, 9, 5, 15, 1, 2, 4, 8, 16, 1, 17, 1, 2, 4, 8, 16, 3, 6, 12, 9, 18, 1, 19, 1, 2, 4, 8, 16, 5, 10, 20, 1, 3, 9, 7, 21, 1, 2, 4, 8, 16, 11, 22, 1, 23
Offset: 1

Views

Author

Michael De Vlieger, Jul 28 2016

Keywords

Comments

Product matrix of tensors T = 1,p,p^2,...,p^e that include the powers 1 <= e of prime divisors p such that p^e <= n.
This sequence is analogous to A275055 but differs from it in that the tensors T include not only powers p^e that divide n but all powers p^e <= n.
The matrix a(n) is bounded by n, thus all products m <= n.
Let omega(n) = A001221(n). The matrix a(n) has omega(n) dimensions and is an omega(n)-dimensional simplex with (omega(n)-1) "right-angle: sides and 1 irregular surface that is bounded by n.
A027750(n) is a subset of A162306(n) and in a(n), the terms of A275055(n) appear in an contiguous omega(n)-dimensional parallelepiped (parallelotope) with 1 at the origin and n at the opposite corner. Thus the omega(n)-dimensional array described by A275055(n) is fully contained in the simplex-like matrix described by a(n). Divisors appear within the parallelepiped while nondivisors appear in the field outside the parallelepiped (see examples). Terms within the parallelepiped appear in A027750(n) while those outside appear in A272618(n).
For a(2^x + 2) there is a term m = (n-2); m != (n-1) except for n=2, since GCD(n, n-1)=1.
a(p^e) = A027750(p^e) = A162306(p^e) = A275055(p^e) for e >= 1.

Examples

			Triangle begins:
1;
1, 2;
1, 3;
1, 2, 4;
1, 5;
1, 2, 4, 3, 6;
1, 7;
1, 2, 4, 8;
1, 3, 9;
1, 2, 4, 8, 5, 10;
1, 11;
1, 2, 4, 8, 3, 6, 12, 9;
1, 13;
1, 2, 4, 8, 7, 14;
1, 3, 9, 5, 15;
1  2, 4, 8, 16;
1, 17;
1, 2, 4, 8, 16, 3, 6, 12, 9, 18;
...
2 prime divisors: n = 96
   1  2  4  8 16 32 64
   3  6 12 24 48 96
   9 18 36 72
  27 54
  81
thus a(96) = {1,2,4,8,16,32,64,3,6,12,24,48,96,9,18,36,72,27,54,81}.
The divisors of 72 (thus the terms of A275055(72)) appear in a rectangle delimited by 1 at top left and 72 at bottom right.
3 prime divisors: n = 60
(the 3 dimensional levels correspond with powers of 5)
   level 5^0:            level 5^1:         level 5^2:
   1  2  4  8 16 32  |    5 10 20 40    |   25 50
   3  6 12 24 48     |   15 30 60       |
   9 18 36           |   45             |
  27 54              |                  |
thus a(60) = {1,2,4,8,16,32,3,6,12,24,48,9,18,36,27,54,5,10,20,40,15,30,60,45,25,50}.
The divisors of 60 (thus the terms of A275055(60)) appear in a parallelepiped delimited by 1 at top left of level 5^0 and 60 at bottom right of level 5^1.
		

Crossrefs

Cf. A162306, A010846 (row length), A243103 (row product), A027750 (divisors of n), A000005 (number of divisors of n), A272618 (nondivisors m <= n that have prime divisors p that also divide n), A243822 (number of such nondivisors of n), A275055 (Product of tensor of prime divisor powers that are also divisors).

Programs

  • Mathematica
    f[n_] := If[n == 1, 1, Function[w, ToExpression@ StringJoin["With[{n=", ToString@ n, "}, Table[", ToString@ InputForm[Times @@ Map[Power @@ # &, w]], ", ", Most@ Flatten@ Map[{#, ", "} &, #], "]]"] &@ MapIndexed[Function[p, StringJoin["{", ToString@ Last@ p, ", 0, Log[", ToString@ First@ p, ", n/(", ToString@ InputForm[Times @@ Map[Power @@ # &, Take[w, First@ #2 - 1]]], ")]}"]]@ w[[First@ #2]] &, w]]@ Map[{#, ToExpression["p" <> ToString@ PrimePi@ #]} &, Reverse[FactorInteger[n][[All, 1]]]] ]; Array[f, 24] // Flatten (* Michael De Vlieger, Mar 08 2017 *)

A279907 Triangle read by rows: T(n,k) is the smallest power of n that is divisible by k, or -1 if no such power exists.

Original entry on oeis.org

0, 0, 1, 0, -1, 1, 0, 1, -1, 1, 0, -1, -1, -1, 1, 0, 1, 1, 2, -1, 1, 0, -1, -1, -1, -1, -1, 1, 0, 1, -1, 1, -1, -1, -1, 1, 0, -1, 1, -1, -1, -1, -1, -1, 1, 0, 1, -1, 2, 1, -1, -1, 3, -1, 1, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 0, 1, 1, 1, -1, 1, -1, 2, 2, -1, -1, 1, 0, -1, -1, -1, -1, -1, -1
Offset: 1

Views

Author

Michael De Vlieger, Dec 26 2016

Keywords

Comments

T(n,1) = 0 since 1 | n^0.
T(n,p) = 1 for prime divisors p of n since p | n^1.
T(n,d) = 1 for divisors d > 1 of n since d | n^1.
Row n for prime p have maximum value 1, since all k < p are coprime to p, and k | p^1 only when k = p.
Values greater than 1 pertain only to composite k of composite n > 4, but not in all cases. T(n,k) = 1 for squarefree kernels k of composite n.
T(n,k) = -1 for numbers k > 1 coprime to n and for numbers that are products of at least one prime q coprime to n and one prime p | n.
T(n,k) is nonnegative for all numbers k for which n^k (mod k) = 0, i.e., all the prime divisors p of k also divide n.
The largest possible value s in row n of T = floor(log_2(n)), since the largest possible multiplicity of any number m <= n pertains to perfect powers of 2, as 2 is the smallest prime. This number s first appears at T(2^s + 2, 2^s) for s > 1.
If T(n,k) is positive, 1/k terminates T(n,k) digits after the radix point in base n. If T(n,k) is negative, 1/k is recurrent in base n.
From Robert Israel, Dec 28 2016: (Start)
T(a*b,c*d) = max(T(a,c),T(b,d)) if GCD(a,b)=1, GCD(b,d)=1,T(a,c)>=0 and T(b,d)>=0.
T(n,a*b) = max(T(n,a),T(n,b)) if GCD(a,b)=1 and T(n,a)>=0 and T(n,b)>=0. (End)

Examples

			The triangle T(n,k) begins (with -1 shown as "." for clarity):
n\k 1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 ...
1:  0
2:  0  1
3:  0  .  1
4:  0  1  .  1
5:  0  .  .  .  1
6:  0  1  1  2  .  1
7:  0  .  .  .  .  .  1
8:  0  1  .  1  .  .  .  1
9:  0  .  1  .  .  .  .  .  1
10: 0  1  .  2  1  .  .  3  .  1
11: 0  .  .  .  .  .  .  .  .  .  1
12: 0  1  1  1  .  1  .  2  2  .  .  1
13: 0  .  .  .  .  .  .  .  .  .  .  .  1
14: 0  1  .  2  .  .  1  3  .  .  .  .  .  1
15: 0  .  1  .  1  .  .  .  2  .  .  .  .  .  1
...
		

Crossrefs

Cf.: A010846 (number of nonnegative k in row n), A162306 (k with nonnegative values in a(n)), A051731 (k with values 0 or 1), A000005 (number of k in row n with values 0 or 1), A272618 (k with values > 1), A243822 (number of k in row n with values > 1), A007947.

Programs

  • Maple
    f:= proc(n,k) local Fk,Fn,i;
       if k = 1 then return 0 fi;
       Fk:= ifactors(k)[2];
       Fn:= map(t -> padic:-ordp(n,t[1]),Fk);
       if min(Fn) = 0 then -1 else max(seq(ceil(Fk[i,2]/Fn[i]),i=1..nops(Fk))) fi
    end proc:
    seq(seq(f(n,k),k=1..n),n=1..20); # Robert Israel, Dec 28 2016
  • Mathematica
    Table[Boole[k == 1] + (Boole[#[[-1, 1]] == 1] (-1 + Length@ #) /. 0 -> -1) &@ NestWhileList[Function[s, {#1/s, s}]@ GCD[#1, #2] & @@ # &, {k, n}, And[First@ # != 1, ! CoprimeQ @@ #] &], {n, 16}, {k, n}] // Flatten (* or *)
    Table[SelectFirst[Range[0, Floor@ Log2@ n], PowerMod[n, #, k] == 0 &] /. k_ /; MissingQ@ k -> -1, {n, 12}, {k, n}] // TableForm (* Version 10.2 *)

A280269 Irregular triangle T(n,m) read by rows: smallest power e of n that is divisible by m = term k in row n of A162306.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 2, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 2, 1, 3, 1, 0, 1, 0, 1, 1, 1, 1, 2, 2, 1, 0, 1, 0, 1, 2, 1, 3, 1, 0, 1, 1, 2, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 2, 1, 3, 1, 2, 4, 1, 0, 1, 0, 1, 1, 1, 2, 1, 2, 1, 0, 1, 1, 2, 1, 0, 1, 2, 3, 1, 4, 1, 0, 1, 0, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Michael De Vlieger, Dec 30 2016

Keywords

Comments

This table eliminates the negative values in row n of A279907.
Let k = A162306(n,m), i.e., the value in column m of row n.
T(n,1) = 0 since 1 | n^0.
T(n,p) = 1 for prime divisors p of n since p | n^1.
T(n,d) = 1 for divisors d > 1 of n since d | n^1.
Row n for prime p have two terms, {0,1}, the maximum value 1, since all k < p are coprime to p, and k | p^1 only when k = p.
Row n for prime power p^i have (i+1) terms, one zero and i ones, since all k that appear in corresponding row n of A162306 are divisors d of p^i.
Values greater than 1 pertain only to composite k of composite n > 4, but not in all cases. T(n,k) = 1 for squarefree kernels k of composite n.
Numbers k > 1 coprime to n and numbers that are products of at least one prime q coprime to n and one prime p | n do not appear in A162306; these do not divide n^e evenly.
T(n,k) is nonnegative for all numbers k for which n^k (mod k) = 0, i.e., all the prime divisors p of k also divide n.
The largest possible value s in row n of T = floor(log_2(n)), since the largest possible multiplicity of any number m <= n pertains to perfect powers of 2, as 2 is the smallest prime. This number s first appears at T(2^s + 2, 2^s) for s > 1.
1/k terminates T(n,k) digits after the radix point in base n for values of k that appear in row n of A162306.
Originally from Robert Israel at A279907: (Start)
T(a*b,c*d) = max(T(a,c),T(b,d)) if GCD(a,b)=1, GCD(b,d)=1,T(a,c)>=0 and T(b,d)>=0.
T(n,a*b) = max(T(n,a),T(n,b)) if GCD(a,b)=1 and T(n,a)>=0 and T(n,b)>=0.
(End)

Examples

			Triangle T(n,m) begins:  Triangle A162306(n,k):
1:  0                    1
2:  0  1                 1  2
3:  0  1                 1  3
4:  0  1  1              1  2  4
5:  0  1                 1  5
6:  0  1  1  2  1        1  2  3  4  6
7:  0  1                 1  7
8:  0  1  1  1           1  2  4  8
9:  0  1  1              1  3  9
10: 0  1  2  1  3  1     1  2  4  5  8  10
...
		

Crossrefs

Cf. A162306, A279907 (T(n,k) with values for all 1 <= k <= n), A280274 (maximum values in row n), A010846 (number of nonnegative k in row n), A051731 (k with e <= 1), A000005 (number of k in row n with e <= 1), A272618 (k with e > 1), A243822 (number of k in row n with e > 1), A007947.

Programs

  • Mathematica
    Table[SelectFirst[Range[0, #], PowerMod[n, #, k] == 0 &] /. m_ /; MissingQ@ m -> Nothing &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* Version 10.2, or *)
    DeleteCases[#, -1] & /@ Table[If[# == {}, -1, First@ #] &@ Select[Range[0, #], PowerMod[n, #, k] == 0 &] &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* or *)
    DeleteCases[#, -1] & /@ Table[Boole[k == 1] + (Boole[#[[-1, 1]] == 1] (-1 + Length@ #) /. 0 -> -1) &@ NestWhileList[Function[s, {#1/s, s}]@ GCD[#1, #2] & @@ # &, {k, n}, And[First@# != 1, ! CoprimeQ @@ #] &], {n, 24}, {k, n}] // Flatten

A280274 a(n) is the maximum value in row n of A279907.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 3, 2, 1, 1, 4, 1, 2, 2, 4, 1, 2, 1, 4, 1, 2, 1, 4, 1, 1, 3, 5, 2, 3, 1, 5, 3, 2, 1, 5, 1, 3, 2, 5, 1, 3, 1, 5, 3, 3, 1, 5, 2, 2, 3, 5, 1, 3, 1, 5, 2, 1, 2, 6, 1, 3, 3, 6, 1, 2, 1, 6, 3, 3, 2, 6, 1, 2, 1, 6, 1, 4, 2, 6, 4, 2, 1, 6, 2, 3, 4, 6, 2, 4, 1, 6, 2, 3
Offset: 1

Views

Author

Michael De Vlieger, Dec 30 2016

Keywords

Comments

Consider integers 1 <= r <= n where all the prime divisors p of r also divide n. Call such r of n "regulars" of n, or "r regular to n".
Consider rho = smallest power e of n such that r | n^e. a(n) is the largest value of rho found among regular 1 <= r <= n of n.
For n=1, r=1 | n^0; since it is the only integer in the range 1 <= k <= n and since 1 is the empty product, a(1) = 0.
a(p) = 1 for prime p, since all 1 <= k <= n must either divide or be coprime to p; if k is coprime to p it is nonregular and does not qualify. Thus only k=1 and k=p divide some power of n; 1 | p^0 and p | p^1 by definition. Thus the largest value of rho among r={1,p} is 1.
a(p^x) = 1 for prime powers p^x with x >= 1, since the only regular 1 <= r <= p^x are divisors that are powers {1, p^1, p^2, ... p^(x-1), p^x}; since these r are all divisors d, d | n^1 by definition, thus the largest rho among these r is 1. Thus, a(n) = 1 for n with omega(n) = 1.
a(4) = 1 since all regular 1 <= r <= 4 divide 4; all divisors d divide n^1 by definition, thus the largest rho among these r is 1.
a(n) > 1 for composite n > 4, since there is at least one nondivisor regular ("semidivisor") 1 <= r <= n. (See A243822.) If r does not divide n, then it divides some power n^e with e > 1, since all the prime divisors p of r divide n. Another way to look at this is that the set of prime divisors p of semidivisor r is a subset of the set of prime divisors p of n, yet r does not itself divide n. The multiplicity of at least one prime divisor p of r exceeds that of the corresponding prime divisor p of n. The maximum possible multiplicity of any number m < n pertains to the largest power of 2 < n, thus the greatest possible value of rho = floor(log_2(n)).
a(n) for n = 2 (mod 4) = floor(log_2(n)), since such n is an odd multiple of 2. Thus the largest power of 2^x less than n has multiplicity x for 2 while that in n is 1. Any other prime divisor q of n has multiplicity less than that of 2. Thus 2^x | n^x, and the largest rho among 1 <= r <= n = x = floor(log_2(n)).
a(n) for squarefree n = floor(log_p(n)) = A280363(n), where p is the smallest distinct prime divisor of n.
Consider the standard form prime decomposition of n = p_1^e_1 * p_2^3_2 * ... p_i^e_i. a(n) for all other n is the maximum value of ceiling(floor(log_p_i(n))/e_i), i.e., ceiling(A280363(n)/e_i).
a(n) < n.
a(n) is the longest terminating base-n expansion of 1/r for 1 <= r <= n. Example: in base 10, the longest terminating expansion of 1/r is 3 for r = 8. 1/8 = .125. a(10) = 3.
Also maximum value in row n of A280269.

Examples

			Row n of A280269               a(n)
1:  0                          0
2:  0  1                       1
3:  0  1                       1
4:  0  1  1                    1
5:  0  1                       1
6:  0  1  1  2  1              2
7:  0  1                       1
8:  0  1  1  1                 1
9:  0  1  1                    1
10: 0  1  2  1  3  1           3
11: 0  1                       1
12: 0  1  1  1  1  2  2  1     2
13: 0  1                       1
14: 0  1  2  1  3  1           3
15: 0  1  1  2  1              2
16: 0  1  1  1  1              1
...
		

Crossrefs

Cf.: A162306, A279907 (T(n,k) with values for all 1 <= k <= n), A280269, A007947, A280363.

Programs

  • Mathematica
    Table[Function[f, Which[n == 1, 0, Length@ f == 1, 1, Max[f[[All, -1]]] == 1, Floor[Log[f[[1, 1]], n]], True, Max@ Map[With[{p = #1, e = #2}, Ceiling[Floor[Log[p, n]]/e]] & @@ # &, f]]][FactorInteger[n]], {n, 120}] (* most efficient, or *)
    Table[If[PrimeNu@ n == 1, 1, Max@Map[Function[k, SelectFirst[Range[0, #], PowerMod[n, #, k] == 0 &] /. m_ /; MissingQ@ m -> Nothing], Range@ n] &@ Floor@ Log2@ n], {n, 120}] (* Version 10.2, or *)
    Max@ DeleteCases[#, -1] & /@ Table[If[# == {}, -1, First@ #] &@ Select[Range[0, #], PowerMod[n, #, k] == 0 &] &@ Floor@ Log2@ n, {n, 120}, {k, n}] (* or *)
    Max@ DeleteCases[#, -1] & /@ Table[Boole[k == 1] + (Boole[#[[-1, 1]] == 1] (-1 + Length@#) /. 0 -> -1) &@ NestWhileList[Function[s, {#1/s, s}]@ GCD[#1, #2] & @@ # &, {k, n}, And[First@ # != 1, ! CoprimeQ @@ #] &], {n, 120}, {k, n}]
Previous Showing 11-20 of 32 results. Next