poj 1365 Prime Land 【簡單數論】

#include<iostream>
#include<sstream>
#include<cstdio>
#include<vector>
using namespace std;
bool not_prime[40000];
vector<int>prime,ans;
int main()
{
    for(int i=2;i<40000;i++)
        if(!not_prime[i])
        {
            prime.push_back(i);
            for(int j=i+i;j<40000;j+=i) not_prime[j]=1;
        }
    string s;
    while(getline(cin,s),s[0]!=48)
    {
        int sum=1,x,y;
        ans.clear();
        istringstream ss(s);
        while(ss>>x>>y)
             while(y--) sum*=x;
        sum--;
        for(int i=prime.size()-1;i>=0;i--)
        {
            if(sum%prime[i]==0)
            {
                ans.push_back(prime[i]);
                int cnt=0;
                while(sum%prime[i]==0) cnt++,sum/=prime[i];
                ans.push_back(cnt);
            }
        }
        for(int i=0;i<ans.size();i++)
            cout<<ans[i]<<(i==ans.size()-1?'\n':' ');
    }
    return 0;
}