OFFSET
1,1
COMMENTS
This equivalence criterion splits a set of composite numbers into two classes and can be used to count certain combinatorial objects.
PROG
(MATLAB)
n=100; % gives all terms of the sequence not exceeding n
A=[];
for i=1:n
dn=divisors(i);
if size(dn, 2)>2 && mod(totient(i)/totient(dn(2)), totient(i)/totient(dn(end-1)))==0
A=[A i];
end
end
function [res] = totient(n)
res=0;
for i=1:n
if gcd(i, n)==1
res=res+1;
end
end
end
(PARI) isok(c) = if ((c>1) && !isprime(c), my(t=eulerphi(c), d=divisors(c)); ((t/eulerphi(d[2])) % (t/eulerphi(d[#d-1]))) == 0); \\ Michel Marcus, Dec 28 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Maxim Karimov, Dec 27 2020
STATUS
approved