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.

A243822 Number of k < n such that rad(k) | n but k does not divide n, where rad = A007947.

This page as a plain text file.
%I A243822 #63 Jul 13 2025 11:06:09
%S A243822 0,0,0,0,0,1,0,0,0,2,0,2,0,2,1,0,0,4,0,2,1,3,0,3,0,3,0,2,0,10,0,0,2,4,
%T A243822 1,5,0,4,2,3,0,11,0,3,2,4,0,5,0,6,2,3,0,8,1,3,2,4,0,14,0,4,2,0,1,14,0,
%U A243822 4,2,12,0,6,0,5,3,4,1,15,0,4,0,5,0,16,1,5,3,3,0,20,1,4,3,5,1,8,0,7,2,6
%N A243822 Number of k < n such that rad(k) | n but k does not divide n, where rad = A007947.
%C A243822 Former name: number of "semidivisors" of n, numbers m < n that do not divide n but divide n^e for some integer e > 1. See ACM Inroads paper.
%H A243822 Michael De Vlieger, <a href="/A243822/b243822.txt">Table of n, a(n) for n = 1..10000</a>
%H A243822 Michael De Vlieger, <a href="https://doi.org/10.1145/2077808.2077809">Exploring Number Bases as Tools</a>, ACM Inroads, March 2012, Vol. 3, No. 1, pp. 4-12.
%H A243822 Michael De Vlieger, <a href="https://doi.org/10.13140/RG.2.2.28311.18084">Regular and coregular numbers</a>, ResearchGate, 2024.
%H A243822 Michael De Vlieger, <a href="/A243822/a243822.png">Log log scatterplot of a(n)</a>, n = 1..2^20
%F A243822 a(n) = A010846(n) - A000005(n) = card({row n of A162306} \ {row n of A027750}).
%F A243822 a(n) = A045763(n) - A243823(n).
%F A243822 a(n) = (Sum_{1<=k<=n, gcd(n,k)=1} mu(k)*floor(n/k)) - tau(n). - _Michael De Vlieger_, May 10 2016, after _Benoit Cloitre_ at A010846.
%F A243822 From _Michael De Vlieger_, Aug 11 2024" (Start)
%F A243822 a(n) = 0 for n in A000961, a(n) > 0 for n in A024619.
%F A243822 a(n) = A051953(n) - A000005(n) + 1 = n - A000010(n) - A000005(n) - A243823(n) + 1.
%F A243822 a(n) = A355432(n) + A361235(n).
%F A243822 a(n) = A355432(n) for n in A360768.
%F A243822 a(n) = A361235(n) for n not in A360768.
%F A243822 a(n) = number of terms in row n of A272618.
%F A243822 a(n) = sum of row n of A304570. (End)
%e A243822 From _Michael De Vlieger_, Aug 11 2024: (Start)
%e A243822 Let S(n) = row n of A162306 and let D(n) = row n of A027750.a(2) = 0 since S(2) \ D(2) = {1, 2} \ {1, 2} is null.
%e A243822 a(10) = 2 since S(10) \ D(10) = {1, 2, 4, 5, 8, 10} \ {1, 2, 5, 10} = {4, 8}.a(16) = 0 since S(16) \ D(16) = {1, 2, 4, 8, 16} \ {1, 2, 4, 8, 16} is null, etc.Table of a(n) and S(n) \ D(n):
%e A243822    n  a(n)  row n of A272618.
%e A243822   ---------------------------
%e A243822    6    1   {4}
%e A243822   10    2   {4, 8}
%e A243822   12    2   {8, 9}
%e A243822   14    2   {4, 8}
%e A243822   15    1   {9}
%e A243822   18    4   {4, 8, 12*, 16}
%e A243822   20    2   {8, 16}
%e A243822   21    1   {9}
%e A243822   22    3   {4, 8, 16}
%e A243822   24    3   {9, 16, 18*}
%e A243822   26    3   {4, 8, 16}
%e A243822   28    2   {8, 16}
%e A243822   30   10   {4, 8, 9, 12, 16, 18, 20, 24, 25, 27}
%e A243822 Terms in A272618 marked with an asterisk are counted by A355432. All other terms are counted by A361235. (End)
%t A243822 Table[Count[Range[n], _?(And[Divisible[n, Times @@ FactorInteger[#][[All, 1]]], ! Divisible[n, #]] &)], {n, 120}] (* _Michael De Vlieger_, Aug 11 2024 *)
%Y A243822 Cf. A000005, A000961, A024619, A027750, A010846, A045763, A162306, A243823, A272618, A304570, A355432, A361235.
%K A243822 nonn
%O A243822 1,10
%A A243822 _Michael De Vlieger_, Jun 11 2014
%E A243822 New name from _David James Sycamore_, Aug 11 2024