|
|
A077662
|
|
a(n) = n-th positive integer not relatively prime to n, with a(1)=1.
|
|
3
|
|
|
1, 4, 9, 8, 25, 9, 49, 16, 27, 16, 121, 18, 169, 24, 33, 32, 289, 27, 361, 34, 49, 40, 529, 36, 125, 48, 81, 49, 841, 40, 961, 64, 84, 64, 112, 54, 1369, 72, 102, 66, 1681, 58, 1849, 80, 96, 88, 2209, 72, 343, 84, 136, 96, 2809, 81, 200, 98, 153, 112, 3481, 82
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MATHEMATICA
|
Join[{1}, Table[Select[Range[5000], !CoprimeQ[#, n]&][[n]], {n, 2, 60}]] (* Harvey P. Dale, May 26 2023 *)
|
|
PROG
|
(PARI) a(n) = local(i, k); if(n<=1, n, k=n-eulerphi(n); i=n\k; k*=i; i*=n; while(k<n, k+=(gcd(i++, n)!=1)); i)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|