分解质因数

问题描述
　　求出区间[a,b]中所有整数的质因数分解。
输入格式
　　输入两个整数a，b。
输出格式
　　每行输出一个数的分解，形如k=a1*a2*a3…(a1<=a2<=a3…，k也是从小到大的)(具体可看样例)
样例输入
3 10
样例输出
3=3
4=2*2
5=5
6=2*3
7=7
8=2*2*2
9=3*3
10=2*5

数据规模和约定
　　2<=a<=b<=10000

提示
　　先打表所有素数，然后再分解。

#include <iostream>
using namespace std;
int main()
{
    int i,j,a,b,k,m,flag,num[10000];//储存10000以内所有质数
    cin>>a>>b;
    num[0]=2;
    j=1;
    for(i=1;i<10000;i++)//打表10000内所有质数
    {
        for(k=2;k<i;k++)
        {
            if(i%k==0)
                break;
        }
        if(k==i)
            num[j++]=i;
    }
    for(i=a;i<=b;i++)
    {
        flag=0;
        for(k=2;k<i;k++)
        {
            if(i%k==0)
            {
                flag=1;
                break;
            }
        }
        if(flag==0)//i是质数，分解为其本身
            cout<<i<<"="<<i<<endl;
        else//i不是质数，将其分解
        {
            cout<<i<<"=";
            m=i;
            while(m!=1)//m==1说明已分解完
            {
            for(k=0;k<10000;k++)//循环10000内的质数
            {
                if(m%num[k]==0)
                {
                    cout<<num[k];
                    m=m/num[k];
                    if(m!=1)
                        cout<<"*";
                    break;
                }
            }
            }
            cout<<endl;
        }
    }
    return 0;
}