A378165 Differences between adjacent terms of A076467 that correspond to the locations of abc-quality records of A378164.
7, 44, 35, 10, 23440
Offset: 1
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.
Pairs b,c of consecutive terms of A076467 A378167 c-b b c = a(n) 8, 8, 16, 32, 32, 64, 343, 2401, 2744, 17576, 456976, 474552, 65610000, 157464000000, 157529610000, 11329982936, 407154266788096, 407165596771032, 26102469128, 1490999139060488, 1491025241529616, 315404039943, 173903379291252081, 173903694695292024, 152838610998696, 661905203228158707216, 661905356066769705912.
\\ Uses M. F. Hasler's A076467_vec from A076467 rad(x) = vecprod(factor(x)[,1]); a378166_7(upto) = {my(W=A076467_vec(upto)); for(k=2, #W, my(d=W[k]-W[k-1], q=rad(W[k])/rad(W[k]*W[k-1]*d)); if(q==1, print([d, W[k-1], W[k]])))}; \\ Alternative program not using rad, more efficient a378166_7(upto) = {my(W=A076467_vec(upto)); for(k=2, #W, my(C=Set(factor(W[k])[,1]), d=W[k]-W[k-1]); if(#setminus(Set(factor(d)[,1]), C)>0, , if(#setminus(Set(factor(W[k-1])[,1]), C)==0, print([d, W[k-1], W[k]]))))}; a378166_7(10^18)
Comments