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.

A071189 Smallest prime factor of sum of divisors of n.

This page as a plain text file.
%I A071189 #23 Mar 24 2024 03:55:31
%S A071189 1,3,2,7,2,2,2,3,13,2,2,2,2,2,2,31,2,3,2,2,2,2,2,2,31,2,2,2,2,2,2,3,2,
%T A071189 2,2,7,2,2,2,2,2,2,2,2,2,2,2,2,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,127,2,2,
%U A071189 2,2,2,2,2,3,2,2,2,2,2,2,2,2,11,2,2,2,2,2,2,2
%N A071189 Smallest prime factor of sum of divisors of n.
%H A071189 Antti Karttunen, <a href="/A071189/b071189.txt">Table of n, a(n) for n = 1..10000</a>
%F A071189 a(n) = A020639(A000203(n)).
%F A071189 a(n) = 2 if and only if n is in A028983. - _Amiram Eldar_, Mar 24 2024
%t A071189 Table[FactorInteger[DivisorSigma[1,n]][[1,1]],{n,90}] (* _Harvey P. Dale_, May 15 2011 *)
%o A071189 (PARI) A071189(n) = if(1==n, n, my(f = factor(sigma(n))); vecmin(f[, 1])); \\ _Antti Karttunen_, Jul 24 2017
%o A071189 (PARI) first(n) = {my(v = vector(n, i, 2), sq = List()); for(i=1, sqrtint(n), listput(sq, i^2); listput(sq, 2*i^2)); listsort(sq); v[1]=1; for(i=2, #sq, if(sq[i]>n,break); v[sq[i]] = factor(sigma(sq[i]))[, 1]~[1]);v} \\ _David A. Corneth_, Jul 24 2017
%Y A071189 Cf. A000203, A020639, A028982, A028983, A071187, A071190.
%K A071189 nonn
%O A071189 1,2
%A A071189 _Reinhard Zumkeller_, May 15 2002