今天做了几个数论的题,感觉比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