UVA, 516 Prime Land

题意：给一个数的指数形式，例：5 1 2 1  num=5^1*2^1  求num-1的质因数分解结果

思路：素数筛法，质因数分解

代码如下：

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>

using namespace std;

const int N=1000000;
int cnt,tot;
long num;
long prime[N];
bool isprime[N];
string ss;
int a[N],b[N];

void getprime();
void showprime();
void getnum();
bool datecin();
void datecal(long);
void showres();


void getprime()
{
    memset(isprime,1,sizeof(isprime));
    isprime[0]=isprime[1]=cnt=0;
    for(int i=2;i<N;i++)
    {
        if(isprime[i])
            prime[cnt++]=i;
        for(int j=0;j<cnt&&i*prime[j]<N;j++)
        {
            isprime[i*prime[j]]=0;
            if(!(i%prime[j]))break;
        }
    }
}

void showprime()
{
    cout<<cnt<<endl;
    for(int i=0;i<cnt;i++)
        cout<<prime[i]<<' ';
    cout<<endl;
}

bool datecin()
{
    getline(cin,ss);
    if((ss[0]=='0'))
    {
        return false;
    }
    getnum();
    return true;
}

void getnum()//得到num
{
    int l=ss.length();//cout<<l<<endl;
    int k=1;
    double a,b,sum=1;
    for(int i=0;i<l;i++)
    {
        if(i==0)
        {
            a=atoi(ss.c_str());
        }
        else if(ss[i]==' ')
        {
            i++,k++;
            string s=ss.substr(i);
            int n=atoi(s.c_str());
            if(k==1)
            {
                a=n;
            }
            else if(k==2)
            {
                b=n;
                sum*=pow(a,b);
                k=0;
                num=(long)sum;
                //cout<<num<<endl;
            }
            //cout<<n<<endl;
        }
    }
    num--;
    //cout<<num<<endl;
}

void datecal(long num)
{
    double temp=sqrt(num)+1;
    tot=0;
    //cout<<'?'<<endl;
    for(int i=0;temp>prime[i];i++)
    {
        if(num%prime[i]==0)
        {
            a[++tot]=prime[i];
            b[tot]=0;
            while(num%prime[i]==0)
            {
                ++b[tot];
                num/=prime[i];
            }
        }
    }
    if(num!=1)
    {
        a[++tot]=num;
        b[tot]=1;
    }
}

void showres()
{
    for(int i=tot;i>=1;i--)
    {
        printf("%d %d",a[i],b[i]);
        if(i!=1)
            printf(" ");
    }
    printf("
");
}

int main()
{
    getprime();
    //showprime();
    while(datecin())
    {
        datecal(num);
        showres();
    }
    return 0;
}

关于输入，因为没有告诉你有多少个数字，用字符串输入转换比较好。