素因子分解

给定某个正整数 N，求其素因子分解结果，即给出其因式分解表达式 N=p​1​^k1*p​2​k​2​​⋯p​m^​k​m，其中pi为素因子并要求由小到大输出，指数ki为pi的个数；当ki为1即因子pi只有一个时不输出ki。

输入格式：
输入long int范围内的正整数 N。

输出格式：
按给定格式输出N的素因式分解表达式，即 N=p1^k1 *p2^k2*…*pmkm，其中pi为素因子并要求由小到大输出，指数ki为pi的个数；当ki为1即因子pi只有一个时不输出ki。

输入样例：
1323
输出样例：
1323=3^3*7^2

#include
using namespace std;

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

int main()
{
	long int n,cnt=0;  int flag=0,flag2=0,flag3=0;
	cin>>n;
	cout<<n<<"=";
	if(n==1)
	{
		cout<<1;
		return 0;
	}
	if(prime(n))
	{
		cout<<n;
		return 0;
	}
	for(int i=2;i<=n;i++)
	{
		while(n%i==0&&prime(i))
		{
			cnt++;
			n=n/i;
			flag2=1;
			flag3=1;
		}
		if(flag3==0)
		continue;
		if(flag==0&&flag2)
		{
			if(cnt==1)
			cout<<i;
			else
			cout<<i<<"^"<<cnt;
		}
		else if(flag&&flag2)
		{
			if(cnt==1)
			cout<<"*"<<i;
			else
			cout<<"*"<<i<<"^"<<cnt;
		}
		flag=1;
		cnt=0;
		flag2=0;
			
	} 
	return 0;
}