function?uns?=?UDall（ns）
h?=?1:n;
ind?=?find（gcd（hn）==1）;
hm?=?h（ind）;?
m?=?length（hm）;?
udt?=?mod（h‘*hmn）;?
ind0?=?find（udt==0）;
udt（ind0）=n;
if?s>m
disp（‘s?必須小于或等于?m‘）;
return;
else
mcs?=nchoosek（ms）;
if?mcs<1e5
tind?=?nchoosek（1:ms）;
[pq]?=?size（tind）;cd2?=?zeros（p1）;
for?k=1:p
UT?=?udt（1:ntind（k:））;
cd2（k1）?=?UDCD2（UT）;
end
tc=tind（find（abs（cd2?-?min（cd2））?<1e1?-?5）:）;
for?r=1:size（tc1）;
uns?（::r）=?udt（:tc（r:））;
end
else
for?k?=?1:n
a?=?k;
UT?=?mod（h‘*a.^（0:s-?1）n）;
cd2（k1）?=?UDCD2（UT）;
end
tc?=?find（abs（cd2?-?min（cd2））for?r=1:size（tc1）;
uns?（::r）=?mod（h‘*tc（r）.^（0:s-?1）n）;
end
ind0?=?find（uns==0）;?uns（ind0）=n;

end

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????472??2018-09-10?21:18??代碼\UDCD2.m
?????文件?????????707??2018-09-10?21:40??代碼\UDall.m
?????目錄???????????0??2018-09-10?21:32??代碼\