POJ 1365 因子分解

题意：给一个整数n的素数因子指数序列：p0^e0,p1^e1,……pi^ei.它表示n=p0^e0*p1^e1*……*pi^ei（其中pi为素因子）。现在让你求n-1这个数的素因子指数序列，而且是按素因子降序输出。

#include<iostream>
#include<string>
using namespace std;
string s;
int a[1000];
int num;
void init(){
    int sum=0;
    num=1;
    for(int i=0;i<s.size();i++){
        if(s[i]==' '){
            a[num++]=sum;
            sum=0;
            continue;
        }
        sum = sum*10 + ( s[i]-'0' );
    }
    if(sum) a[num++]=sum; 
}
void fan(long long sum){
    if(sum==2||sum==3) { printf("%d 1\n",sum); return ;}
    int b[100],k=0;
    memset(b,0,sizeof(b));
    for(int i=2;i*i<=sum;i++){
        if(sum%i==0){
            b[k++]=i;
            while(sum%i==0){
                b[k]++;
                sum /= i;
            }
            k++;
        }
    }
    if(sum>1) { b[k++]=sum; b[k++]=1; } /// 最后的数可能是素数，要注意了，其他的没什么 
    for(int i=k-2;i>=0;i-=2){
        if(i==k-2) printf("%d %d",b[i],b[i+1]);
        else printf(" %d %d",b[i],b[i+1]);
    }
    printf("\n");
}
int main(){
    while(getline(cin,s)){
        if(s[0]=='0') break;
        init();
        long long sum=1;
        for(int i=2;i<num;i+=2){
            while(a[i]--){
                sum *= a[i-1];
            }
        }
        sum--;
        fan(sum);
    }
}