function y=jie(n)
y=1;
disp('n!=');
i=2;
l=0;
k=1;
while i<=n;
    if pan(i)==1;
        j=1;
        a=i^j;
        x=0;
        
        while a<=n;
            x=fix(n/a)+x;
            j=j+1;
            a=i^j;
        end
        fprintf('%d^%d*',i,x);
        if rem(k,20)==0
            %disp('\n');
            fprintf('\n');
        end
        k=k+1;
        
        y=y*i^x;
    end
        if l==0;
            i=i+1;
        else
            i=i+2;
        end
      l=1;  
end