OFFSET
1,2
COMMENTS
Starts to differ from A322791 in row n=16, where 4 is an exponential divisor but not an exponential unitary divisor.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..15331 (first 10000 rows, flattened)
Nicusor Minculete and László Tóth, Exponential unitary divisors, Annales Univ. Sci. Budapest., Sect. Comp., Vol. 35 (2011), pp. 205-216.
MAPLE
A361255 := proc(n)
local expundivs , d, isue, p, ai, bi;
expudvs := {} ;
for d in numtheory[divisors](n) do
isue := true ;
for p in numtheory[factorset](n) do
ai := padic[ordp](n, p) ;
bi := padic[ordp](d, p) ;
if bi > 0 then
if modp(ai, bi) <>0 or igcd(bi, ai/bi) <> 1 then
isue := false;
end if;
else
isue := false ;
end if;
end do;
if isue then
expudvs := expudvs union {d} ;
end if;
end do:
sort(expudvs) ;
end proc:
seq(op(A361255(n)), n=1..60) ;
MATHEMATICA
udivQ[n_, m_] := (n > 0 && m > 0 && Divisible[n, m] && CoprimeQ[m, n/m]);
expuDivQ[n_, d_] := Module[{f = FactorInteger[n]}, And @@ MapThread[udivQ, {f[[;; , 2]], IntegerExponent[d, f[[;; , 1]]]}]]; expuDivs[1] = {1};
expuDivs[n_] := Module[{d = Rest[Divisors[n]]}, Select[d, expuDivQ[n, #] &]];
Table[expuDivs[n], {n, 1, 70}] // Flatten (* Amiram Eldar, Mar 11 2023 *)
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
R. J. Mathar, Mar 06 2023
STATUS
approved