UVA 10375 Choose and Divide +計算組合數

今天做了幾個數論的題,感覺比DP容易一點,但同樣也有點坑。

題意:計算C(a,b)/C(n,m) ,a,b,n,m<=INT_MAX;
分析:令c=a-b,k=m-n; 列式得ans=(a!m!k!)/(n!b!c!)。用上面一行與下面一行相消後就可以直接計算了。
具體實現見代碼:

#include<bits/stdc++.h>
#define LL long long
using namespace std;
double ans;
double arr[2][30000];  ///儲存上下兩行的數據
int cnt[2];  ///上下兩行數的個數
void add(int a,int b,int j){ ///向兩行中添加相消後剩餘的數
    for(int i=b+1;i<=a;i++) arr[j][cnt[j]++]=i;  
    
}
void slove(int j){ ///計算答案
    ans=1.0;
    for(int i=0;i<cnt[j];i++) ans*=(arr[0][i]/arr[1][i]);
    if(!j) for(int i=cnt[0];i<cnt[1];i++) ans/=arr[1][i];
    else for(int i=cnt[1];i<cnt[0];i++) ans*=arr[0][i];
}
int main(){
    int a[3],b[3];
    //freopen("in.txt","r",stdin);
    while(scanf("%d%d%d%d",&a[0],&b[0],&b[1],&a[1])!=EOF){ 
        a[2]=b[1]-a[1]; ///爲了方便後面的排序,先把上下兩行中同行的數記錄在一行
        b[2]=a[0]-b[0];
        sort(a,a+3); ///排序後兩行中較大的數與較大的數相消,能有效減少計算量
        sort(b,b+3);
        cnt[1]=cnt[0]=0;
        for(int i=0;i<3;i++){
            if(a[i]>b[i]) add(a[i],b[i],0);
            else if(a[i]<b[i]) add(b[i],a[i],1);
            else continue;
        }
        cnt[0]>cnt[1]?slove(1):slove(0);
        printf("%.5f\n",ans);
    }
    return 0;
}

寫完後發現有大佬幾行代碼就寫完了
這就是大佬與萌新的區別嗎
我永遠都是萌新QAQ

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章