|
%I #12 Jul 06 2021 18:26:16
%S 1,1,1,1,1,2,1,2,3,2,1,3,1,2,3,8,1,3,1,4,3,2,1,6,5,2,3,4,1,6,1,4,3,2,
%T 5,6,1,2,3,8,1,6,1,4,15,2,1,8,7,5,3,4,1,6,5,8,3,2,1,12,1,2,9,8,5,6,1,
%U 4,3,10,1,12,1,2,15,4,7,6,1,40,9,2,1,12,5,2,3,8,1,15,7,4,3,2,5,12,1,7,9,10,1,6,1,8,15
%N a(n) = n / gcd(n, A002034(n)), where A002034(n) gives the smallest positive integer k such that n divides k!.
%H Antti Karttunen, <a href="/A345933/b345933.txt">Table of n, a(n) for n = 1..16384</a>
%H Antti Karttunen, <a href="/A345933/a345933.txt">Data supplement: n, a(n) computed for n = 1..65537</a>
%H <a href="/index/Fa#factorial">Index entries for sequences related to factorial numbers</a>
%F a(n) = n / A345931(n) = n / gcd(n, A002034(n)).
%t Table[n/GCD[n,m=1;While[Mod[m!,n]!=0,m++];m],{n,100}] (* _Giorgos Kalogeropoulos_, Jul 03 2021 *)
%o (PARI)
%o A002034(n) = if(1==n,n,my(s=factor(n)[, 1], k=s[#s], f=Mod(k!, n)); while(f, f*=k++); (k)); \\ After code in A002034.
%o A345933(n) = (n/gcd(n, A002034(n)));
%Y Cf. A002034, A345931, A345932.
%K nonn
%O 1,6
%A _Antti Karttunen_, Jul 01 2021
|
