A329613 Numbers k such that A071187(k) <> A329614(k).

324, 1296, 2500, 5625, 9604, 10000, 11664, 20736, 21609, 38416, 50625, 58564, 60025, 82944, 90000, 114244, 131769, 160000, 194481, 234256, 236196, 250000, 257049, 334084, 345744, 360000, 366025, 456976, 521284, 614656, 640000, 714025, 717409, 751689, 810000, 944784, 960400, 1119364, 1172889, 1185921, 1265625, 1327104, 1336336

Offset: 1

Author

Antti Karttunen, Nov 18 2019

Keywords

• nonn

Comments

All terms are squares; their square roots are in A329611.

Links

• David A. Corneth, Table of n, a(n) for n = 1..15072 (first 661 terms from Antti Karttunen, terms <= 10^14)
• David A. Corneth, PARI program

Crossrefs

Subsequence of A000290.

Cf. A055932, A071187, A329611, A329614.

Programs

PARI

A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A071187(n) = if(1==n, n, my(f = factor(numdiv(n))); vecmin(f[, 1]));
A329614(n) = A071187(A108951(n));
k=0; n=0; while(n<2^16, n++; u = n*n; if(A071187(u)!=A329614(u), k++; write("b329613.txt", k, " ", u); print1(u,", ")));

PARI

\\ See link