OFFSET
1,1
COMMENTS
PROG
(PARI)
allocatemem(2^30);
A061395(n) = if(1==n, 0, primepi(vecmax(factor(n)[, 1]))); \\ After M. F. Hasler's code for A006530.
A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ Michel Marcus, Oct 10 2016
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From Michel Marcus
i=1; n=0; while(i < 100, n++; if(isA277200(n), write("b277200.txt", i, " ", n); i++));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 14 2016
STATUS
approved