Choose and divide uva 10375

求C(i,j)/ C(u,v)

公式化简

 

筛出所有质数

除法的时候： 因数分解， 求对每个质数的个数的贡献

  #include <iostream>
  #include <cstring>
  #include <cmath>
  #include <iomanip>
  #include <vector>
  using namespace std;
   const int M =1e4+3;
 int vis[M],p,q,s,r;
 vector<int> prm;
 int A[M];
 
 void init(int top){
 	int i,j;
 	for(i=2;i<=top;i++)
 	 if(vis[i]==0){
 	 	prm.push_back(i);
 	  	for(j=i*2;j<=top;j+=i) vis[j]=1;
 	 } 
 }
 
 void ADD(int x,int v){
 	for(int i=0;i<prm.size();i++){
 		while(x%prm[i]==0){
 			A[i]+=v;
 			x/=prm[i];
 		}
 		if(x==1) break;
 	}
 }
 void help(int x,int v){
 	for(int i=1;i<=x;i++) ADD(i,v);
 }
 void solve(){
 	help(p,1),help(r-s,1),help(s,1),
	help(p-q,-1),help(q,-1),help(r,-1);
	double ans=1;
	for(int i=0;i<prm.size();i++){
		ans*=pow(double(prm[i]),A[i]);
	}
	printf("%.5lf\n",ans);
 } 
 signed main(){
	init(1e4);
	while(scanf("%d%d%d%d",&p,&q,&r,&s)==4){
		memset(A,0,sizeof A);
		solve();
	}
 }