【C/C++】任意大于1的整数分解成素数因子乘积的形式

//

#include<stdio.h>
#include<math.h>
#include<malloc.h>
int isprime(long n);
void decompose_to_primes(int n);
int main()
{
    decompose_to_primes(3);
    return 0;
}

void decompose_to_primes(int n)
{
    int num;
    int *prime[2];
    int pi=0;
    int temp,i,j,exp;
    temp=(int)sqrt(n);
    num=temp;
    prime[0]=(int*)malloc(num*sizeof(int));
    prime[1]=(int*)malloc(num*sizeof(int));
    if(n<4)// n==2 or n==3
    {prime[0][0]=n;prime[1][0]=1;pi++;}
    for(i=2;i<=temp;i++)
        if(isprime(i))
        {
            exp=0;
            j=i;
            while(n%j==0 && j<=n)
            {exp++;j*=i;}
            j/=i;
            //if(exp)
            {prime[0][pi]=i;prime[1][pi]=exp;
            pi++;
            n/=j;
            }          
        }
        //output primefactors whose exp!=0
        for(i=0;i<pi;i++)
            if(prime[1][i]!=0)
                printf("<%d,%d> ",prime[0][i],prime[1][i]);
            printf("
");
            free(prime[0]);
            free(prime[1]);
}

int isprime(long n)
{
    int bound;
    int i;
    bound=(int)sqrt(n);
    for(i=2;i<=bound;i++)
        if(n%i==0)
            return 0;
        return 1;
}