UVA-10375 唯一分解定理

#include<iostream>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<stdio.h>
#define rep(i,j,k) for(int i=j;i<=k;i++) using namespace std;
const int N = ;
bool is_prime[N];
int prime[N];
int cnt=;
void get_prime(){
int m=sqrt(N);
memset(is_prime,,sizeof(is_prime));
rep(i,,m){
if (!is_prime[i]){
for (int j=i*i;j<=N;j+=i){
is_prime[j]=;
}
}
}
rep(i,,N){
if (is_prime[i]==){
prime[cnt++]=i;
}
}
}
int e[N];
void add_integer(int num,int d){
for (int i=;i<cnt;i++){
while(num%prime[i]==){
num/=prime[i];
e[i]+=d;
}
if (num==)break;
}
}
void add_factorial(int n,int d){
for(int i=;i<=n;i++)
add_integer(i,d);
}
int main(){
int p,q,r,s;
get_prime();
while(~scanf("%d%d%d%d",&p,&q,&r,&s)){
memset(e,,sizeof(e));
add_factorial(p,);
add_factorial(q,-);
add_factorial(p-q,-);
add_factorial(r,-);
add_factorial(s,);
add_factorial(r-s,);
int maxx = max(r,p);
double ans=;
for (int i=;i<=maxx;i++){
ans*=pow(prime[i],e[i]);
}
printf("%.5f\n",ans);
}
return ;
}
上一篇:oracle的一知半解


下一篇:HTML5入门总结 HTML5API