PTA 7-9 素因子分解

给定某个正整数 N，求其素因子分解结果，即给出其因式分解表达式 N=p1的k1次方*p2的k2次方…pm的km次方；

输入格式：

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

输出格式：

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

输入样例：

123

输出样例：

1323=3^3*72；

这里次方没弄好，就是从头往后读。

#include 
 int main()
 {
     long long n;
     int i;
     scanf("%lld",&n);
     printf("%lld=",n);
     if (n>1) {
         int flag = 0;
         for (i=2;i<=n;i++) {
             int k = 0;
             if (n%i==0) {
                 while (n%i==0) {
                     k++;
                     n /= i;
                 }
                 if (flag)
                     printf("*");
                 printf("%d",i);
                 flag = 1;
                 if (k>1)
                     printf("^%d", k);
             }
             
         }
     }
   	 else {
        printf("%d", 1);
    }
     printf("\n");
     return 0; }