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.

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.

This page as a plain text file.
%I A280269 #16 Jul 15 2024 03:21:13
%S A280269 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,
%T A280269 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,
%U A280269 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
%N 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.
%C A280269 This table eliminates the negative values in row n of A279907.
%C A280269 Let k = A162306(n,m), i.e., the value in column m of row n.
%C A280269 T(n,1) = 0 since 1 | n^0.
%C A280269 T(n,p) = 1 for prime divisors p of n since p | n^1.
%C A280269 T(n,d) = 1 for divisors d > 1 of n since d | n^1.
%C A280269 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.
%C A280269 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.
%C A280269 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.
%C A280269 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.
%C A280269 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.
%C A280269 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.
%C A280269 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.
%C A280269 Originally from _Robert Israel_ at A279907: (Start)
%C A280269 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.
%C A280269 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.
%C A280269 (End)
%H A280269 Michael De Vlieger, <a href="/A280269/b280269.txt">Table of n, a(n) for n = 1..10202</a> (rows 1 <= n <= 660)
%e A280269 Triangle T(n,m) begins:  Triangle A162306(n,k):
%e A280269 1:  0                    1
%e A280269 2:  0  1                 1  2
%e A280269 3:  0  1                 1  3
%e A280269 4:  0  1  1              1  2  4
%e A280269 5:  0  1                 1  5
%e A280269 6:  0  1  1  2  1        1  2  3  4  6
%e A280269 7:  0  1                 1  7
%e A280269 8:  0  1  1  1           1  2  4  8
%e A280269 9:  0  1  1              1  3  9
%e A280269 10: 0  1  2  1  3  1     1  2  4  5  8  10
%e A280269 ...
%t A280269 Table[SelectFirst[Range[0, #], PowerMod[n, #, k] == 0 &] /. m_ /; MissingQ@ m -> Nothing &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* Version 10.2, or *)
%t A280269 DeleteCases[#, -1] & /@ Table[If[# == {}, -1, First@ #] &@ Select[Range[0, #], PowerMod[n, #, k] == 0 &] &@ Floor@ Log2@ n, {n, 24}, {k, n}] // Flatten (* or *)
%t A280269 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
%Y A280269 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.
%K A280269 nonn,tabf,easy
%O A280269 1,14
%A A280269 _Michael De Vlieger_, Dec 30 2016