|
|
A356158
|
|
a(n) = gcd(n, A347879(n)).
|
|
3
|
|
|
1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 3, 1, 1, 3, 2, 1, 2, 1, 10, 1, 2, 1, 6, 1, 1, 1, 28, 1, 2, 1, 2, 3, 2, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 3, 2, 1, 3, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 5, 2, 1, 1, 3, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 2, 1, 3, 1, 2, 1, 2, 1, 2, 1, 1, 3
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The fixed points of this sequence is given by the union of {2} and A336702.
|
|
LINKS
|
|
|
FORMULA
|
|
|
PROG
|
(PARI)
Abincompreflen(n, m) = { my(x=binary(n), y=binary(m), u=min(#x, #y)); for(i=1, u, if(x[i]!=y[i], return(i-1))); (u); };
Abinprefix(n, k) = { my(digs=binary(n)); fromdigits(vector(k, i, digs[i]), 2); };
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t); };
A156552(n) = {my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552
A347879(n) = A348041sq(n, sigma(n));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|