|
|
A346953
|
|
a(n) is the number of divisors of A346950(n) ending with 3.
|
|
2
|
|
|
1, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 4, 2, 2, 4, 2, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
EXAMPLE
|
a(17) = 4 since there are 4 divisors of A346950(17) = 429 ending with 3: 3, 13, 33 and 143.
|
|
MATHEMATICA
|
b={}; For[n=0, n<=450, n++, For[k=0, k<=n, k++, If[Mod[10*n+9, 10*k+3]==0 && Mod[(10*n+9)/(10*k+3), 10]==3 && 10*n+9>Max[b], AppendTo[b, 10*n+9]]]]; (* A346950 *) a={}; For[i =1, i<=Length[b], i++, AppendTo[a, Length[Drop[Select[Divisors[Part[b, i]], (Mod[#, 10]==3&)]]]]]; a
|
|
PROG
|
(Python)
from sympy import divisors
def f(n): return sum(d%10 == 3 for d in divisors(n)[1:-1])
def A346950upto(lim): return sorted(set(a*b for a in range(3, lim//3+1, 10) for b in range(a, lim//a+1, 10)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|