A121920 a(n) = n-th integer from among those positive integers k where GCD(k,n) = prime power >1. a(1) = 1.

1, 4, 9, 8, 25, 10, 49, 16, 27, 18, 121, 22, 169, 26, 36, 32, 289, 34, 361, 38, 54, 42, 529, 46, 125, 50, 81, 54, 841, 63, 961, 64, 88, 66, 120, 70, 1369, 74, 105, 78, 1681, 87, 1849, 86, 111, 90, 2209, 94, 343, 98, 141, 102, 2809, 106, 210, 110, 159, 114, 3481, 125

Offset

1,2

Links

Table of n, a(n) for n=1..60.

Example

2,4,5,6,8,12,14,15,16,18,22,... is the sequence of positive integers k where GCD(k, 10) = p^j, j = positive integer, p = prime. 18 is the 10th of these integers, so a(10) = 18.

Mathematica

f[n_] := Block[{c = n, k = 0}, If[n == 1, 1, While[c > 0, k++; While[Length[FactorInteger[GCD[k, n]]] != 1, k++ ]; c--; ]; k]]; Table[f[n], {n, 65}] (* Ray Chandler, Sep 06 2006 *)

Crossrefs

Sequence in context: A077662 A063718 A063748 * A318279 A065642 A285109

Adjacent sequences: A121917 A121918 A121919 * A121921 A121922 A121923

Keyword

nonn

Author

Leroy Quet, Sep 02 2006

Extensions

Extended by Ray Chandler, Sep 06 2006

Status

approved