A249714 -id:A249714 - cp's OEIS frontend

A249441 a(n) is the smallest prime whose square divides at least one entry in the n-th row of Pascal's triangle, or 0 if there is no such prime.

0, 0, 0, 0, 2, 0, 2, 0, 2, 2, 2, 0, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 0, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2

Offset: 0

Vladimir Shevelev, Oct 28 2014

a(n) = 3 for 15, 31, 47, 63, 95, 127, 191, 255, 383, 511, 767, 1023, 1535, 2047, 3071, etc.
The above values all occur in A249723 and from 31 onward seem to be given by A052955(n>=8). (Cf. also A249714 & A249715). - Antti Karttunen, Nov 04 2014
Using the Kummer theorem on carries, one can prove that, if a(n)>3 or 0, then n>23 takes the form of either 1...1 or 101...1 in base 2 and simultaneously 212...2 in base 3. However, it is easy to see that this leads to a contradiction. Thus there are no terms greater than 3 and only 8 zeros, i.e., there are only 8 rows in Pascal's triangle that contain all squarefree numbers. It turns out that the latter result has been known for a long time (see A048278).

Charles R Greathouse IV, Table of n, a(n) for n = 0..10000
E. E. Kummer, Über die Ergänzungssätze zu den allgemeinen Reciprocitätsgesetzen, J. Reine Angew. Math. 44 (1852), 93-146.
Mihai Prunescu, Sign-reductions, p-adic valuations, binomial coefficients modulo p^k and triangular symmetries. Preprint 2013.

Cf. A005117, A007913, A048278, A052955, A249695, A249714, A249715, A249723.

a_list := proc(len) local s; s := proc(L,p) local n; seq(max(op(map(b-> padic[ordp](b,p),{seq(binomial(n,k),k=0..n)}))),n=0..L); map(k-> `if`(k<2,0,p),[%]) end: zip((x,y)-> `if`(x=0,y,x),s(len,2),s(len,3)) end: a_list(86); # Peter Luschny, Nov 01 2014
# alternative
A249441 := proc(n)
    local p,wrks,bi,k;
    if n in [0,1,2,3,5,7,11,23] then
        return 0 ;
    end if;
    p :=2 ;
    while true do
        wrks := false;
        bi := 1 ;
        for k from 0 to n do
            if modp(bi,p^2) = 0 then
                wrks := true;
                break;
            end if;
            bi := bi*(n-k)/(1+k) ;
        end do:
        if wrks then
            return p;
        end if;
        p := nextprime(p) ;
    end do:
end proc: # R. J. Mathar, Nov 04 2014

row[n_] := Table[Binomial[n, k], {k, 1, (n-Mod[n, 2])/2}];
a[n_] := If[MemberQ[{0, 1, 2, 3, 5, 7, 11, 23}, n], 0, For[p = 2, True, p = NextPrime[p], If[AnyTrue[row[n], Divisible[#, p^2]&], Return[p]]]];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Jul 30 2018 *)

a(n) = my(o=0); for(k=1,n\2, o+=valuation((n-k+1)/k, 2); if(o>1, return(2))); if(n<24 && n!=15, 0, 3) \\ Charles R Greathouse IV, Nov 03 2014

A249441(n) = { forprime(p=2,3,for(k=0,n\2,if((0==(binomial(n,k)%(p*p))),return(p)))); return(0); } \\ This is more straightforward, but a slower implementation - Antti Karttunen, Nov 03 2014

a(n)=if((n+1)>>valuation(n+1,2)<5, if(n<24 && setsearch([1,2,3,5,7,11,23],n), 0, 3), 2) \\ Charles R Greathouse IV, Nov 06 2014

More terms from Peter J. C. Moses, Oct 28 2014

A249695 a(n)=0, if A249441(n)=0; otherwise, a(n) is the smallest i such that A249441(n)^2 divides binomial(n,i).

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 3, 0, 1, 2, 3, 0, 1, 6, 3, 7, 1, 2, 3, 4, 1, 6, 3, 0, 1, 2, 3, 12, 1, 6, 3, 5, 1, 2, 3, 4, 1, 6, 3, 8, 1, 2, 3, 12, 1, 6, 3, 21, 1, 2, 3, 4, 1, 6, 3, 24, 1, 2, 3, 12, 1, 6, 3, 1, 1, 2, 3, 4, 1, 6, 3, 8, 1, 2, 3, 12, 1, 6, 3, 16, 1, 2, 3, 4, 1, 6, 3

Offset: 0

Vladimir Shevelev, Nov 04 2014

nonn

After a(0) = 0, A048278 gives the positions of seven other zeros in the sequence. - Antti Karttunen, Nov 04 2014

Antti Karttunen, Table of n, a(n) for n = 0..10000

A249714 and A249715 give the record values and their positions.
Cf. A048278, A249441, A249722, A249723, A249726.
Differs from A249442 for the first time at n=9.

A249695 := proc(n)
    a41n := A249441(n) ;
    if a41n = 0 then
        return 0;
    end if;
    bi := 1;
    for i from 0 do
        if modp(bi,a41n^2)= 0 then
            return i;
        end if;
        bi := bi*(n-i)/(1+i) ;
    end do:
end proc: # R. J. Mathar, Nov 04 2014

bb[n_] := Table[Binomial[n, k], {k, 1, (n - Mod[n, 2])/2}];
a41[n_] := If[MemberQ[{0, 1, 2, 3, 5, 7, 11, 23}, n], 0, For[p = 2, True, p = NextPrime[p], If[AnyTrue[bb[n], Divisible[#, p^2]&], Return[p]]]];
a[n_] := If[(a41n = a41[n]) == 0, 0, For[i = 1, True, i++, If[Divisible[ Binomial[n, i], a41n^2], Return[i]]]];
a /@ Range[0, 100] (* Jean-François Alcover, Mar 27 2020 *)

A249695(n) = { forprime(p=2,3,for(k=0,floor(n/2),if((0==(binomial(n,k)%(p*p))),return(k)))); return(0); } \\ Straightforward and unoptimized version. But fast enough for 10000 terms.
A249695(n) = { for(p=2,3, my(o=0); for(k=1, n\2, o+=valuation((n-k+1)/k, p); if(o>1, return(k)))); return(0); } \\ This version is based on Charles R Greathouse IV's code for A249441.
for(n=0, 10000, write("b249695.txt", n, " ", A249695(n)));
\\ Antti Karttunen, Nov 04 2014

cp's OEIS Frontend

A249441 a(n) is the smallest prime whose square divides at least one entry in the n-th row of Pascal's triangle, or 0 if there is no such prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Extensions

A249695 a(n)=0, if A249441(n)=0; otherwise, a(n) is the smallest i such that A249441(n)^2 divides binomial(n,i).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A249715 Positions of records in A249695.

Views

Author

Keywords

Comments

Crossrefs