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.

A381096 Number of k <= n such that k is neither coprime to n and rad(k) != rad(n), where rad = A007947.

This page as a plain text file.
%I A381096 #4 Feb 16 2025 23:02:47
%S A381096 0,0,0,0,0,3,0,1,1,5,0,6,0,7,6,4,0,10,0,10,8,11,0,13,3,13,6,14,0,21,0,
%T A381096 11,12,17,10,20,0,19,14,21,0,29,0,22,19,23,0,28,5,28,18,26,0,33,14,29,
%U A381096 20,29,0,42,0,31,25,26,16,45,0,34,24,45,0,42,0,37
%N A381096 Number of k <= n such that k is neither coprime to n and rad(k) != rad(n), where rad = A007947.
%C A381096 Number of k <= n in the cototient of n that do not share the same squarefree kernel as n.
%C A381096 Define a number k "neutral" to n to be such that 1 < gcd(k,n) < k, that is, k neither divides n nor is coprime to n. A045763(n) is the number of k < n such that k is neutral to n.
%C A381096 Define quality Q(k) to be true if k is such that 1 < gcd(k,n) and rad(k) != rad(n).
%C A381096 Then for k <= n and n > 1, a(n) = A045763(n), but admitting divisors k | n such that rad(k) != rad(n), and eliminating occasional nondivisors k such that rad(k) = rad(n), i.e., k listed in row n of A359929 for n = A360768(i).
%H A381096 Michael De Vlieger, <a href="/A381096/b381096.txt">Table of n, a(n) for n = 1..10000</a>
%H A381096 Michael De Vlieger, <a href="/A381096/a381096.png">Log log scatterplot of a(n)</a>, n = 1..2^20.
%F A381096 a(1) = 0, a(p) = a(4) = 0.
%F A381096 a(n) = A045763(n) - A005361(n).
%F A381096 For n > 1, a(n) = n - phi(n) - tau(n/rad(n)) = A000010(n) - A005361(n).
%F A381096 For n > 1, a(n) = n - A000010(n) - A008479(n) + A355432(n).
%e A381096 a(6) = 3 since {2, 3, 4} are neither coprime to 6 and do not have the squarefree kernel 6.
%e A381096 a(8) = 1 since only 6 is neither coprime to 8 and does not share the squarefree kernel 2 with 8.
%e A381096 a(10) = 5 since {2, 4, 5, 6, 8} are neither coprime to 10 nor have the squarefree kernel 10.
%e A381096 a(12) = 6 since {2, 3, 4, 8, 9, 10} are neither coprime to 12 nor have the squarefree kernel 6.
%e A381096 a(14) = 7 since {2, 4, 6, 7, 8, 10, 12} are neither coprime to 14 nor have the squarefree kernel 14, etc.
%t A381096 {0}~Join~Table[n - EulerPhi[n] - DivisorSigma[0, n/rad[n]], {n, 2, 120}]
%Y A381096 Cf. A000010, A005361, A008479, A355432, A359929, A381094.
%K A381096 nonn,easy
%O A381096 1,6
%A A381096 _Michael De Vlieger_, Feb 14 2025