poj 2034 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; }