OFFSET
1,1
COMMENTS
a(13) > 5*10^27.
EXAMPLE
Pairs b,c of consecutive terms of A076467
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.
PROG
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)
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Hugo Pfoertner, Nov 20 2024
STATUS
approved