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

View Code

 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 }

posted @ 2012-06-18 09:02 Hug_Sea 阅读(124) 评论(0) 收藏举报

刷新页面返回顶部