OFFSET
1,2
COMMENTS
a(n) = 2 iff n is even and is a term of A226872. - Daniel Suteu, Jul 28 2019
From Bernard Schott, Jul 30 2019: (Start)
a(n) = n if n = 1, 2, 6, 42, 1806.
a(n) = 6 if n is of the form 2^i*3^j, i and j >= 1, so if n is a term of A033845.
a(n) = 10 if n is of the form 2^i*5^j, i >= 2 and j >= 1.
a(n) = 30 if n is of the form 2^i*3^j*5^k, i >=2, j >= 1 and k >= 1. (End)
LINKS
Eric M. Schmidt, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
fa=FactorInteger; d[m_]:= Product[If[IntegerQ[m/(fa[m][[i, 1]]-1)], fa[m][[i, 1]], 1], {i, Length@fa@m}]; Table[d[n], {n, 1, 333}]
PROG
(Sage) def A225821(n) : return prod(p for (p, m) in factor(n) if n%(p-1)==0) # Eric M. Schmidt, Jul 31 2013
(PARI) a(n)=my(f=factor(n)[, 1]); prod(i=1, #f, if(n%(f[i]-1)==0, f[i], 1)) \\ Charles R Greathouse IV, Nov 13 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Jul 30 2013
STATUS
approved