POJ-1365 Prime Land
题意:已知任意一个大于1的数可以表示成一些素数的乘积,即x=p1^e1*p2^e2……pn^en (pi 为素数,ei 为对应素数的个数),现给你x的表示,要你求x-1的表示。
例:输入:5 1 2 1 则x=5^1*2^1=10,所以x-1=9,9可以表示成:9=3^2
输出:3 2
思路:打出素数表、然后从小到大除即可。
题目链接:http://poj.org/problem?id=1365
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include <cstdio> 2 #include <cstring> 3 #include <cstdlib> 4 #include <cmath> 5 #include <string> 6 #include <algorithm> 7 #include <iostream> 8 using namespace std; 9 const int N=40000; 10 11 bool flag[N]; 12 int prime[N]; 13 typedef struct In{ 14 int val; 15 int cnt; 16 }In; 17 In b[N]; 18 19 int main(){ 20 21 // freopen("data.in","r",stdin); 22 // freopen("data.out","w",stdout); 23 24 memset(flag,false,sizeof(flag)); 25 for(int i=2;i<=N;i++){ 26 if(!flag[i]) 27 for(int j=2;i*j<=N;j++) 28 flag[i*j]=true; 29 } 30 int cnt=0; 31 for(int i=2;i<=N;i++) 32 if(!flag[i]) prime[++cnt]=i; 33 long long n,m,sum; 34 char c; 35 while(scanf("%lld",&n),n){ 36 sum=1; 37 scanf("%c%lld",&c,&m); 38 while(m--) sum*=n; 39 while(scanf("%c",&c)!=EOF){ 40 if(c=='\n') break; 41 scanf("%lld%c%lld",&n,&c,&m); 42 while(m--) sum*=n; 43 } 44 sum--; 45 cnt=1; 46 int ct=1; 47 memset(b,0,sizeof(b)); 48 while(sum!=1){ 49 while(sum%prime[cnt]==0){ 50 sum/=prime[cnt]; 51 b[ct].val=prime[cnt]; 52 b[ct].cnt++; 53 } 54 cnt++; 55 if(b[ct].cnt!=0) ct++; 56 } 57 for(int i=ct-1;i>=1;i--) printf(i==1?"%d %d\n":"%d %d ",b[i].val,b[i].cnt); 58 } 59 return 0; 60 }