[算法笔记]1059 Prime Factors (25 分)

1059 Prime Factors (25 分)

Input Specification:

Each input file contains one test case which gives a positive integer NNN in the range of long int.

Output Specification:

Sample Input:

97532468

Sample Output:

97532468=2^2*11*17*101*1291

思路： 找一个数的质因子

1. 首先建立素数表
2. 用结构体表示质因子（质因子和个数）
3. 如果是质因子就加入数组，并且在循环中除以质因子，找到质因子个数
4. 如果n无法被小于sqrt(n)的质因子除尽，那么一定有一个大于sqrt(n)的质因子

 

#include<bits/stdc++.h>
using namespace std;

#define maxn 101000

struct factor
{
    int x;  //质因子
    int cnt;    //个数
}fac[10];

//判断素数
bool isprime(int n)
{
//    if(n==1)return 0;
    for(int i=2;i<(int)sqrt(n*1.0);i++)
    {
        if(n%i==0)
            return 0;
    }
    return 1;
}

//建立素数表
int prime[maxn],j=0;
void findPrime()
{
    for(int i=2;i<maxn;i++)         //1本来就不是素数，为什么还需要判断呢？
    {
        if(isprime(i))
            prime[j++]=i;
    }
}

int findfactor(int n)
{
    int num=0;
    int sqr=(int)sqrt(n*1.0);
    for(int i=0;i<j&&prime[i]<=sqr;i++)     //如果大于n只会浪费时间
    {
        if(n%prime[i]==0)           //如果是质因子就加入数组并初始化数量为0
        {
            fac[num].x=prime[i];
            fac[num].cnt=0;
            while(n%prime[i]==0)
            {
                fac[num].cnt++;
                n/=prime[i];
            }
            num++;
        }
        if(n==1) break;         //节省时间
    }
    if(n!=1)    /*如果无法被小于sqrt(n)的质因子除尽，那么一定有一个大于sqrt(n)的质因子*/
    {
        fac[num].x=n;
        fac[num++].cnt=1;
    }
    return num;
}

int print(int n)
{
    if(n==1)            //1的情况需要特判
    {
        printf("1=1");
        return 0;
    }
    int d=findfactor(n);
    printf("%d=",n);
    for(int i=0;i<d;i++)
    {
        if(i>0)
            cout<<"*";
        printf("%d",fac[i].x);
        if(fac[i].cnt>1)
             printf("^%d",fac[i].cnt);
    }
}
int main()
{
    int n;

    cin>>n;
    findPrime();
    print(n);

    return 0;
}