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.
%I A181833 #8 Jun 28 2013 05:06:37 %S A181833 0,0,2,3,4,4,6,5,6,7,9,5,10,7,10,11,12,6,14,7,14,15,16,5,18,13,17,13, %T A181833 20,7,24,9,18,19,22,15,28,10,22,19,28,9,32,9,26,27,30,5,34,17,33,25, %U A181833 32,7,38,23,36,29,34,5,46 %N A181833 The number of positive integers <= n that are not strongly prime to n. %C A181833 k is strongly prime to n iff k is relatively prime to n and k does not divide n-1. %C A181833 a(n) = n - phi(n) + tau(n-1) if n > 0 and a(0) = 0. %C A181833 Here phi(n) = A000010(n) and tau(n) = A000005(n). %H A181833 Peter Luschny, <a href="http://www.oeis.org/wiki/User:Peter_Luschny/StrongCoprimality">Strong coprimality</a>. %e A181833 a(11) = 11 - card({3,4,6,7,8,9}) = 5. %p A181833 with(numtheory): %p A181833 A181833 := n -> `if`(n=0,0,n-phi(n)+tau(n-1)); %p A181833 A181833a := n -> n - A181830(n); %t A181833 a[n_] := Select[Range[n], Not[CoprimeQ[#, n] && !Divisible[n-1, #]] &] // Length; a[1] = 0; Table[a[n], {n, 0, 60}] (* _Jean-François Alcover_, Jun 28 2013 *) %Y A181833 Cf. A181830, A181831, A181832, A181834, A181835, A181836, A051953. %K A181833 nonn %O A181833 0,3 %A A181833 _Peter Luschny_, Nov 17 2010