function [flag,p,q]=factors_resolve(n,B) 
%p-1因子分解法,n必須為奇整數(shù) 
a=2;
for j=2:B
    a=big_mod(a,j,n);
end
d=gcd(a-1,n);
if 1<d&d<n
    flag=1; %分解成功,n可以分解成兩個素數(shù)的乘積 
    p=d;q=n/d;
else
    flag=0;display('分解失敗,n不能分解成兩個素數(shù)的乘積');    
    p=0;q=0;
end