题意:
给一个整数x的素数分解形式,求出x-1的素数分解形式。
先把x算出来,再把x-1素数分解就可以了,注意一下pow的精度损失
代码:
#include <cstdlib>
#include <cctype>
#include <cstring>
#include <cstdio>
#include <cmath>
#include<climits>
#include <algorithm>
#include <vector>
#include <string>
#include <iostream>
#include <sstream>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <fstream>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <list>
#include <stdexcept>
#include <functional>
#include <utility>
#include <ctime>
using namespace std;
#define PB push_back
#define MP make_pair
#define REP(i,x,n) for(int i=x;i<(n);++i)
#define FOR(i,l,h) for(int i=(l);i<=(h);++i)
#define FORD(i,h,l) for(int i=(h);i>=(l);--i)
#define SZ(X) ((int)(X).size())
#define ALL(X) (X).begin(), (X).end()
#define RI(X) scanf("%d", &(X))
#define RII(X, Y) scanf("%d%d", &(X), &(Y))
#define RIII(X, Y, Z) scanf("%d%d%d", &(X), &(Y), &(Z))
#define DRI(X) int (X); scanf("%d", &X)
#define DRII(X, Y) int X, Y; scanf("%d%d", &X, &Y)
#define DRIII(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z)
#define OI(X) printf("%d",X);
#define RS(X) scanf("%s", (X))
#define MS0(X) memset((X), 0, sizeof((X)))
#define MS1(X) memset((X), -1, sizeof((X)))
#define LEN(X) strlen(X)
#define F first
#define S second
#define Swap(a, b) (a ^= b, b ^= a, a ^= b)
#define Dpoint strcut node{int x,y}
#define cmpd int cmp(const int &a,const int &b){return a>b;}
/*#ifdef HOME
freopen("in.txt","r",stdin);
#endif*/
const int MOD = 1e9+7;
typedef vector<int> VI;
typedef vector<string> VS;
typedef vector<double> VD;
typedef long long LL;
typedef pair<int,int> PII;
//#define HOME
int Scan()
{
int res = 0, ch, flag = 0;
if((ch = getchar()) == '-') //判断正负
flag = 1;
else if(ch >= '0' && ch <= '9') //得到完整的数
res = ch - '0';
while((ch = getchar()) >= '0' && ch <= '9' )
res = res * 10 + ch - '0';
return flag ? -res : res;
}
/*----------------PLEASE-----DO-----NOT-----HACK-----ME--------------------*/
int in[40000];
char s[100000];
int prime[10000];
int vis[40005];
int ans[40005];
void getprime()
{
int cnt=0;
for(int i=2;i<=40000;i++)
{
if(!vis[i])
{
prime[cnt++]=i;
for(int j=0;j<cnt&&prime[j]<=40000/i;j++)
{
prime[prime[j]*i]=1;
if(i%prime[j]==0)
break;
}
}
}
}
int main()
{getprime();
while(gets(s))
{
if(s[0]=='0')
break;
int len=strlen(s);
int t=0;
int cnt=0;
REP(i,0,len)
{if(s[i]==' ')
{
in[cnt++]=t;
t=0;
continue;
}
t=t*10+s[i]-'0';
}
in[cnt++]=t;
int x=1;
for(int i=0;i<cnt;i+=2)
{
x=x*(int)(pow((double)in[i],(double)in[i+1])+0.05);
}
x-=1;
int p=(int)sqrt((double)x);
int q=0;
for(int i=2;i<=p;i++)
{
int tmp=0;
if(x%i==0)
{while(x%i==0)
{
x=x/i;
tmp++;
}
ans[q++]=tmp;
ans[q++]=i;
}
}
if(x!=1)
{
ans[q++]=1;
ans[q++]=x;
}
for(int i=q-1;i>0;i--)
printf("%d ",ans[i]);
printf("%d\n",ans[0]);
}
return 0;
}