cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A108298 Sum of the first 10^n terms in A097975. a(n) = sum_{m=1..10^n} t(m), where t(m) is the sum of the prime divisors of m that are greater than or equal to sqrt(m).

This page as a plain text file.
%I A108298 #8 Jul 07 2014 07:33:49
%S A108298 0,30,1797,132946,10034416,790688821,64867780292,5492352229154,
%T A108298 475943074590494,41984058676639733,3755707610763952011,
%U A108298 339758793864093720073,31019273006095379281810,2853680710328414627392965,264227600111858563511104972
%N A108298 Sum of the first 10^n terms in A097975. a(n) = sum_{m=1..10^n} t(m), where t(m) is the sum of the prime divisors of m that are greater than or equal to sqrt(m).
%C A108298 Does a(n+1)/a(n) converge?
%e A108298 The first 10^2 terms in A097975 sum to 1797, so a(2) = 1797.
%t A108298 s = 0; k = 1; Do[l = Select[Select[Divisors[n], PrimeQ], # >= Sqrt[n]&]; If[Length[l] > 0, s += l[[1]]]; If[n == k, Print[s]; s = 0; k *= 10], {n, 1, 10^7}]
%o A108298 (PARI) a(n) = sum(m=1, 10^n, sumdiv(m, d, d*isprime(d)*(d>=sqrt(m)))); \\ _Michel Marcus_, Jul 07 2014
%Y A108298 Cf. A097975.
%K A108298 more,nonn
%O A108298 0,2
%A A108298 _Ryan Propper_, Jul 24 2005
%E A108298 a(2)-a(7) and the example corrected and a(8)-a(14) from _Hiroaki Yamanouchi_, Jul 07 2014