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.

A196220 Integer quotients of k^2 by the sum of the prime distinct divisors of k^2+1, where k = A196219(n).

This page as a plain text file.
%I A196220 #14 Mar 09 2020 09:12:00
%S A196220 7,18,121,2268,13520,1377,8550,5157,7381,8496,76176,83521,161604,
%T A196220 284229,1028196,4092529,275804,274432,336985,1153476,962948,48841,
%U A196220 319225,276676,617796,3946827,684450,156349,632025,1256454,6368547,244917,2506180,2256004,5410947
%N A196220 Integer quotients of k^2 by the sum of the prime distinct divisors of k^2+1, where k = A196219(n).
%C A196220 Generated by  k = 7, 18, 187, 378, 1560, 1683, … (A196219).
%H A196220 Amiram Eldar, <a href="/A196220/b196220.txt">Table of n, a(n) for n = 1..500</a> (calculated from the b-file at A196219)
%F A196220 a(n) = A196219(n)^2/A008472(A196219(n)^2 + 1). - _Amiram Eldar_, Mar 09 2020
%e A196220 For k = 378, the prime distinct divisors of 378^2 + 1 are 5, 17, 41 and 378^2 /(5+17+41) = 2268. Hence 2268 is in the sequence.
%p A196220 with(numtheory):for k from 1 to 120000 do: y:=factorset(k^2+1): s:=sum(y[i],i=1..nops(y)):if irem(k^2,s)=0 then printf(`%d, `, k^2/s):else fi:od:
%t A196220 Select[Table[n^2/Total[Transpose[FactorInteger[n^2+1]][[1]]],{n,10^5}],IntegerQ] (* _Harvey P. Dale_, Apr 18 2015 *)
%Y A196220 Cf. A002522, A008472, A180278, A196219.
%K A196220 nonn
%O A196220 1,1
%A A196220 _Michel Lagneau_, Sep 29 2011