Choose and divide UVA - 10375 唯一分解定理+Eratosthenes筛法

Eratosthenes筛法筛一遍素数 直接朴素判也能过 然后 写一个add_factorial算阶乘   add_integer用来唯一分解定理

#include<bits/stdc++.h>
using namespace std;
bool is_prime[100000+5];
vector<int>prime;
int e[100000+5];
void init(){
	memset(is_prime,0,sizeof(is_prime));
	int m=sqrt(100000+0.5);
	for(int i=2;i<=m;i++)if(!is_prime[i])
		for(int j=i*i;j<=100000;j+=i)is_prime[j]=1;
}
/*bool is_prime(int n) {
  int m = floor(sqrt(n) + 0.5);
  for(int a = 2; a <= m; a++)
    if(n % a == 0) return false;
  return true;
}*/

void add_integer(int n,int d){
	for(int i=0;i<prime.size();i++){
		while(n%prime[i]==0){
			e[i]+=d;
			n/=prime[i];
		}
		if(n==1)break;
	}
}
void add_factorial(int n,int d){
	for(int i=1;i<=n;i++)add_integer(i,d);
}
int main(){
	init();
	  for(int i = 2; i <= 10000; i++)
    if(!is_prime[i]) prime.push_back(i);
	int p,q,r,s;
	while(cin>>p>>q>>r>>s){
		memset(e,0,sizeof(e));
      add_factorial(p,1);
	   add_factorial(q,-1);
	   add_factorial(p-q,-1);
	    add_factorial(r,-1);
		 add_factorial(s,1);
		  add_factorial(r-s,1);
		  double ans=1;
		  for(int i=0;i<prime.size();i++){
			  ans*=pow(prime[i],e[i]);
		  }
		  printf("%.5lf\n",ans);
	}
	return 0;
}