UVa 10375 (唯一分解定理) Choose and divide

题意：

求组合数C(p, q) / C(r, s)结果保留5为小数。

分析：

先用筛法求出10000以内的质数，然后计算每个素数对应的指数，最后再根据指数计算答案。

 #include <cstdio>
 #include <cmath>
 #include <cstring>

 const int maxn = ;
 int pri[maxn], cnt, e[maxn];    //e记录每个素数的质数
 bool vis[maxn + ];

 void add_interger(int n, int d)    //乘以n的d次幂
 {
     for(int i = ; i < cnt; ++i)
     {
         while(n % pri[i] == )
         {
             n /= pri[i];
             e[i] += d;
         }
         if(n == ) return;
     }
 }

 void add_factorial(int n, int d)    //乘以(n!)的d次幂
 {
     for(int i = ; i <= n; ++i)
         add_interger(i, d);
 }

 int main()
 {
     //freopen("10375in.txt", "r", stdin);

     int m = sqrt(maxn + 0.5);
     for(int i = ; i <= m; ++i) if(!vis[i])
         for(int j = i*i; j <= maxn; j += i) vis[j] = true;
     cnt = ;
     for(int i = ; i < maxn; ++i) if(!vis[i]) pri[cnt++] = i;
     //for(int i = 0; i < 10; ++i) printf("%d\n", pri[i]);

     int p, q, r, s;
     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, );

         double ans = 1.0;
         for(int i = ; i < cnt; ++i)
             ans *= pow(pri[i], e[i]);

         printf("%.5f\n", ans);
     }

     return ;
 }

代码君