|
%I #40 Nov 09 2022 00:05:21
%S 1,2,3,4,5,6,7,8,9,10,12,14,15,16,18,20,21,24,25,27,28,30,32,35,36,40,
%T 42,45,48,49,50,54,56,60,63,64,70,72,75,80,81,84,88,90,91,92,96,98,
%U 100,105,108,112,115,120,125,126,128,135,140,144,147,150,160,161,162,168
%N Numbers k that divide tau(k), where tau(k)=A000594(k) is Ramanujan's tau function.
%C Although most small numbers are in the sequence, it becomes sparser for larger values; e.g., only 504 numbers up to 10000 and only 184 numbers from 10001 to 20000 are in the sequence.
%H Seiichi Manyama, <a href="/A063938/b063938.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/TauFunction.html">Tau Function.</a>
%t (* First do <<NumberTheory`Ramanujan` *) test[n_] := Mod[RamanujanTau[n], n]==0; Select[Range[200], test]
%t (* Second program: *)
%t Select[Range@ 168, Divisible[RamanujanTau@ #, #] &] (* _Michael De Vlieger_, Dec 23 2017 *)
%o (PARI 2.8) for (n=1,1000,if(Mod(ramanujantau(n),n)==0,print1(n", "))) \\ _Dana Jacobsen_, Sep 06 2015
%o (Perl) use ntheory ":all"; my @p = grep { !(ramanujan_tau($_) % $_) } 1..1000; say "@p"; # _Dana Jacobsen_, Sep 06 2015
%o (Python)
%o from itertools import count, islice
%o from sympy import divisor_sigma
%o def A063938_gen(startvalue=1): # generator of terms >= startvalue
%o return filter(lambda n: not -840*(pow(m:=n+1>>1,2,n)*(0 if n&1 else pow(m*divisor_sigma(m),2,n))+(sum(pow(i,4,n)*divisor_sigma(i)*divisor_sigma(n-i) for i in range(1,m))<<1)) % n, count(max(startvalue,1)))
%o A063938_list = list(islice(A063938_gen(),25)) # _Chai Wah Wu_, Nov 08 2022
%Y For the sequence when n is prime see A007659.
%Y Cf. A063940, A000594, A079334, A296991, A296993.
%K nonn,easy
%O 1,2
%A _Robert G. Wilson v_, Aug 31 2001
%E More terms from _Dean Hickerson_, Jan 03 2003
|
