将要求的数离散为素数的指数，然后遍历一遍素数表。。。。
#include <cstdio> #include <iostream> #include <cstring> #include <queue> #include <vector> #include <cmath> using namespace std; typedef long long ll; int prime[10001]; int is_prime[10001]; int ret; void init() { ret=0; memset(is_prime,0,sizeof(is_prime)); for(int i=2;i<=10000;i++) { if(!is_prime[i]) { prime[ret++]=i; for(int j=i*2;j<=10000;j+=i) is_prime[i]=1; } } } int e[10001]; void add2(int n,int d) { for(int i=0;i<ret;i++) { int temp=prime[i]; while(n%temp==0) { n/=temp; e[i]+=d; } if(n==1) break; } } void add(int x,int d) { for(int i=1;i<=x;i++) add2(i,d); } int main() { int p,q,r,s; init(); while(cin>>p>>q>>r>>s) { memset(e,0,sizeof(e)); add(p,1); add(q,-1); add(p-q,-1); add(r,-1); add(s,1); add(r-s,1); double ans=1; for(int i=0;i<ret;i++) { ans*=pow(prime[i],e[i]); } printf("%.5lf ",ans); } return 0; }