uva 10375 Choose and divide

数学题：（组合数学，分解质因子）

题意就是C(m,n)/C(s,t)，凡是组合数就提醒我们要注意数据范围，一般处理组合数学问题的策略就是能约分的先约分，能化简的先化简，能除的先除，不要全部乘起来再除掉

1.转化公式

C(m,n) = [ m*(m-1)*(m-2)...(m-n+1) ] / [ n! ]

2.两者相除

C(p,q)/C(r,s) = [  p*(p-1)...(p-q+1) * s! ] / [ r*(r-1)....(r-s+1)*q! ]

但是不能全部乘起来再除（爆），也不能一边乘一边除（精度），把每个数字分解质因子，分子中分解出来的质因子，个数增加，分母分解出来的质因子，个数减少

最后查看所有的质因子，个数为正的要乘，个数为负的要除，要一边乘一边除否则又会爆

 

代码写得不好，可以优化一下提高时间的，但是没心机写了就不优化了，跑了0.980s

#include <cstdio>

#include <cstring>

#define MAX 10000

int c[MAX+100];  //质因子

//将i分解为质因子并使质因子数增加



void div(int n ,int f)

{

    for(int i=2; n>1; i++) if(!(n%i)) //能分解因子

    {

        while(n>1 & !(n%i))

        {

            c[i]+=f;

            n/=i;

        }

    }

}

void solve(int p,int q,int s,int t)

{

    //分子:  p……p-q+1 * t!

    //分母： s……s-t+1 * m!



    double a,b,ans;

    memset(c,0,sizeof(c));

    for(int i=p; i>=p-q+1; i--) div(i,1);   

    for(int i=t; i>=1; i--)     div(i,1);

    //将i分解为质因子并使质因子数增加



    for(int i=s; i>=s-t+1; i--) div(i,-1);

    for(int i=q; i>=1; i--)     div(i,-1);

    //将i分解质因子并使质因子数减少



    ans=a=b=1;

    for(int i=1; i<=MAX; i++)  //扫描所有质因子

        if(c[i]>0)

            for(int k=1; k<=c[i]; k++) ans*=i;

        else if(c[i]<0)

            for(int k=c[i]; k<0; k++)  ans/=i;



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

}



int main()

{

    int p,q,s,t;

    while(scanf("%d%d%d%d",&p,&q,&s,&t)!=EOF)

    {

        solve(p,q,s,t);

    }

    return 0;

}