login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A344760 "Inverted variant" of A342002: a(n) = A342002(A289234(n)). 6
0, 1, 1, 5, 2, 7, 1, 7, 8, 31, 13, 41, 3, 11, 14, 43, 19, 53, 2, 9, 11, 37, 16, 47, 4, 13, 17, 49, 22, 59, 1, 9, 10, 41, 17, 55, 12, 59, 71, 247, 106, 317, 26, 87, 113, 331, 148, 401, 19, 73, 92, 289, 127, 359, 33, 101, 134, 373, 169, 443, 4, 15, 19, 59, 26, 73, 27, 89, 116, 337, 151, 407, 41, 117, 158, 421, 193, 491, 34 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
LINKS
FORMULA
a(n) = A342001(A346101(n)) = A342002(A289234(n)).
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A342001(n) = (A003415(n) / A003557(n));
A346101(n) = { my(pr=1, p=2, m=1); while(n>0, my(d=n%p); if(d>0, m *= p^lift(1/Mod(d, p))); pr *= p; n \= p; p = nextprime(1+p)); return(m); };
(PARI)
A289234(n) = { my(pr=1, p=2, v=0); while(n>0, my (d=n%p); if(d>0, v += pr * lift(1/Mod(d, p))); pr *= p; n \= p; p = nextprime(1+p)); return(v); }; \\ From A289234.
A342002(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= p^(e>0); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
CROSSREFS
Cf. also A344761, A344762.
Sequence in context: A263832 A351952 A342002 * A343422 A214969 A093591
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Jul 10 2021
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 19 21:09 EDT 2024. Contains 371798 sequences. (Running on oeis4.)