|
|
A098192
|
|
Terms in a specific cycle of length 29 of the map x->A098189(x).
|
|
4
|
|
|
2206, 2210, 2932, 2950, 3394, 3398, 3400, 3402, 3520, 3768, 3844, 4240, 4440, 4516, 4778, 4782, 4828, 4870, 4884, 6428, 6484, 6840, 7056, 7624, 7680, 7984, 8560, 9072, 10264
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Iteration of the map x -> A098189(x) enters cycles of various lengths.
The 29 terms in the first cycle where A098190(x)=29 are listed here in ascending order; the example section of A098190 shows them in the mapping order.
|
|
LINKS
|
|
|
EXAMPLE
|
The cycle follows the mapping 3768 -> 4440 -> ... -> 2210-> 3768 (returning to the first element).
This cycle is entered, for example, if the mapping is started at x=246 (Cf. A098190).
|
|
MATHEMATICA
|
Function[s, Union@ Drop[s, #] &@ Count[DeleteDuplicates@ Map[Position[s, #] &, s], k_ /; Length@ k == 1]]@ NestList[Function[n, DivisorSum[n, # &, CoprimeQ[#, n/#] &] - EulerPhi@n], 246, 10^3] (* Michael De Vlieger, Mar 01 2017 *)
|
|
PROG
|
(PARI) A098189(n)=my(f=factor(n)); prod(k=1, #f~, f[k, 1]^f[k, 2]+1) - eulerphi(f)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,fini,full
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|