OFFSET
1,2
COMMENTS
The sequence can also be obtained by taking a(1) = 1 and then multiplying the previous term by n if n is coprime to the previous term a(n-1), dividing the previous term by n if n divides the previous term a(n-1), taking a(n) = a(n-1) if n is unrelated to a(n-1). - Amarnath Murthy, Oct 30 2002 (corrected by Franklin T. Adams-Watters, Dec 13 2006)
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..1000
MATHEMATICA
a = {1}; Do[AppendTo[a, If[GCD[a[[-1]], n] == 1, a[[-1]]*n, If[Divisible[a[[-1]], n], a[[-1]]/n, a[[-1]]]]], {n, 2, 32}]; a (* Ivan Neretin, May 21 2015 *)
PROG
(PARI) print1(k=1); for(n=2, 99, if(gcd(k, n)==1, k*=n, if(k%n==0, k/=n)); print1(", "k)) \\ Charles R Greathouse IV, May 21 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Amarnath Murthy, Feb 26 2002
EXTENSIONS
a(26)-a(32) corrected by Ivan Neretin, May 21 2015
STATUS
approved