今天做了幾個數論的題,感覺比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