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.

A343249 a(n) is the least k0 <= n such that v_2(n), the 2-adic order of n, can be obtained by the formula: v_2(n) = log_2(n / L_2(k0, n)), where L_2(k0, n) is the lowest common denominator of the elements of the set S_2(k0, n) = {(1/n)*binomial(n, k), with 0 < k <= k0 such that k is not divisible by 2} or 0 if no such k0 exists.

This page as a plain text file.
%I A343249 #80 Jul 12 2022 12:22:45
%S A343249 1,1,3,1,5,3,7,1,9,5,11,3,13,7,5,1,17,9,19,5,7,11,23,3,25,13,27,7,29,
%T A343249 5,31,1,11,17,7,9,37,19,13,5,41,7,43,11,9,23,47,3,49,25,17,13,53,27,
%U A343249 11,7,19,29,59,5,61,31,9,1,13,11,67,17,23,7,71,9,73,37,25,19,11,13,79,5
%N A343249 a(n) is the least k0 <= n such that v_2(n), the 2-adic order of n, can be obtained by the formula: v_2(n) = log_2(n / L_2(k0, n)), where L_2(k0, n) is the lowest common denominator of the elements of the set S_2(k0, n) = {(1/n)*binomial(n, k), with 0 < k <= k0 such that k is not divisible by 2} or 0 if no such k0 exists.
%C A343249 Conjecture: a(n) is the greatest power of a prime different from 2 that divides n.
%H A343249 Dario T. de Castro, <a href="/A343249/b343249.txt">Table of n, a(n) for n = 1..1000</a>
%H A343249 Dario T. de Castro, <a href="http://math.colgate.edu/~integers/w61/w61.pdf">P-adic Order of Positive Integers via Binomial Coefficients</a>, INTEGERS, Electronic J. of Combinatorial Number Theory, Vol. 22, Paper A61, 2022.
%e A343249 For n = 15, a(15) = 5. To understand this result, consider the largest set S_2, which is the S_2(k0=15, 15). According to the definition, S_2(n, n) is the set of elements of the form (1/n)*binomial(n, k), where k goes from 1 to n, skipping the multiples of 2. The elements of S_2(15, 15) are: {1, 0, 91/3, 0, 1001/5, 0, 429, 0, 1001/3, 0, 91, 0, 7, 0, 1/15}, where the zeros were put pedagogically to identify the skipped terms, i.e., when k is divisible by 2. At this point we verify which of the nested subsets {1}, {1, 0}, {1, 0, 91/3}, {1, 0, 91/3, 0}, {1, 0, 91/3, 0, 1001/5},... will match for the first time the p-adic order’s formula. If k vary from 1 to 5 (instead of 15) we see that the lowest common denominator of the set S_2(5, 15) will be 15. So, L_2(5, 15) = 15 and the equation v_2(15) = log_2(15/15) yields a True result. Then we may say that a(15) = 5 specifically because 5 was the least k0.
%t A343249 j = 1;
%t A343249 Nmax = 250;
%t A343249 Array[val, Nmax];
%t A343249 Do[val[i] = 0, {i, 1, Nmax}];
%t A343249 Do[flag = 0;
%t A343249   Do[If[(flag == 0 &&
%t A343249       Prime[j]^IntegerExponent[n, Prime[j]] ==
%t A343249        n/LCM[Table[
%t A343249            If[Divisible[k, Prime[j]], 1,
%t A343249             Denominator[(1/n) Binomial[n, k]]], {k, 1, k}] /.
%t A343249           List -> Sequence]), val[n] = k; flag = 1;, Continue], {k, 1,
%t A343249      n, 1}], {n, 1, Nmax}];
%t A343249 tabseq = Table[val[i], {i, 1, Nmax}];
%o A343249 (PARI) Lp(k, n, p) = {my(list = List()); for (i=1, k, if (i%p, listput(list, binomial(n,i)/n));); lcm(apply(denominator, Vec(list)));}
%o A343249 isok(k, n, v, p) = p^v == n/Lp(k, n, p);
%o A343249 a(n, p=2) = {my(k=1, v=valuation(n, p)); for (k=1, n, if (isok(k, n, v, p), return(k));); n;} \\ _Michel Marcus_, Apr 22 2021
%Y A343249 Cf. A007814, A343250, A343251, A343252, A343253, A345531.
%K A343249 nonn
%O A343249 1,3
%A A343249 _Dario T. de Castro_, Apr 09 2021