分解质因数

问题描述

　　求出区间[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

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
using namespace std;
int su[10010];
int book[10010];
int o;
int sum[10010];
void prim()
{
    memset(book,0,sizeof(book));
    book[0]=book[1]=1;
    for(int i=2; i<=10000; i++)
    {
        if(book[i]==0)
        {
            su[o++]=i;
            for(int j=i*2; j<=10000; j+=i)
                book[j]=1;
        }
    }

}
int main()
{
    int a,b;
    o=0;
    prim();
    while(~scanf("%d%d",&a,&b))
    {
        for(int i=a; i<=b; i++)
        {
            if(book[i]==0)//i为素数
            {
                printf("%d=%d\n",i,i);
                continue;
            }
            else
            {
                int num=i;
                int flag=0;
                int r=0;
                memset(sum,0,sizeof(sum));
                for(int j=0;j<o;j++)
                {
                    while(num%su[j]==0)
                    {
                        sum[r++]=su[j];
                        num=num/su[j];
                        if(book[num]==0)
                        {
                            sum[r++]=num;
                            num=0;
                            flag=1;
                            break;
                        }
                    }
                    if(num==0)
                    {
                        flag=1;
                        break;
                    }
                }
                if(flag==1)
                {
                    printf("%d=",i);
                    printf("%d",sum[0]);
                    for(int j=1;j<r;j++)
                        printf("*%d",sum[j]);
                    printf("\n");
                }
            }
        }
    }
    return 0;
}