[go: up one dir, main page]

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!)
A353574 a(n) = A342002(n) / gcd(A342002(n), A342002(A276154(n))). 5
1, 1, 5, 1, 7, 1, 7, 2, 31, 13, 41, 1, 9, 11, 37, 2, 47, 1, 11, 7, 43, 19, 53, 1, 13, 17, 49, 11, 59, 1, 9, 5, 41, 17, 55, 2, 59, 71, 247, 53, 317, 19, 73, 46, 289, 127, 359, 13, 87, 113, 331, 74, 401, 11, 101, 67, 373, 169, 443, 1, 11, 13, 47, 5, 61, 17, 69, 43, 277, 121, 347, 2, 83, 107, 319, 71, 389, 31, 97, 16 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
Compare the scatter plot to that of A342002.
LINKS
FORMULA
a(n) = A342002(n) / A353573(n) = A342002(n) / gcd(A342002(n), A353572(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]));
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A342001(n) = (A003415(n) / A003557(n));
A353574(n) = { my(u=A342001(A276086(n))); u/gcd(u, A353571(A276086(n))); };
CROSSREFS
Sequence in context: A318676 A346475 A346474 * A265293 A348502 A340089
KEYWORD
nonn,base,easy,look
AUTHOR
Antti Karttunen, Apr 29 2022
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 August 29 19:56 EDT 2024. Contains 375518 sequences. (Running on oeis4.)