在非模意义下直接用公式阶乘求无疑非常容易炸,那么我们就需要更稳定的求组合数的方法
1、用组合数性质C(n,m)=C(n,m-1)+C(n-1,m-1),写递归程序
LL C(int m,int n)
{
if(m==n||!n) return 1;
else if(c[m][n]) return c[m][n];
else return c[m][n]=C(m-1,n)+C(m-1,n-1);
}
注意:C数组大小一定要开够
2、化简C(m,m)=m!/(n!*(m-n)!)=(m-n+1)…m/n!所以我们枚举一个i,计算phi i=1~n (m-n+i)/i即可;
LL C(int n,int m){
long long res=1;
for(int i=1;i<=m;i++)
res=res*(LL)(n-m+i)/i;
return res;
}
最后附一道有关组合数的题目:洛谷 P3166 [CQOI2014]数三角形
发现正着算十分不好算,我们可以考虑用荣斥倒着算,显然取三个点的取法有C((m+1)*(n+1),3)个,在竖直水平的线上被算多的三角形显然有(n+1)C(m+1,3)+(m+1)C(n+1,3)个,然后考虑在斜线上被算多的三点。我们可以枚举a、b表示点(0,0)到(a,b)的线,那么经过平移对称我们可以很明显的发现共有(m-a+1)(n-b+1)*2条线,然后我们只要计算出其上的整点个数就可以轻松地算出此形态的算多情况。直线方程为y=(b/a)x,即使要求出x使左式值为整,令gcd(a,b)=g,则y=((b/g)/(a/g))x,由于要求x小于a,那么显然有g个点满足整点,由于还算上了0点,所以整点个数就是g-1,那么我们就可以轻松计算了。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=1e3+10;
int m,n;
LL ans,c[maxn*maxn][5];
//int exi[SIZE][SIZE];
/*LL C(int n,int m){//C(n,m)=C(n-1,m)+C(n-1,m-1)
if(m==0 || m==n) return 1;
if(c[n][m]) return c[n][m];
else return c[n][m]=C(n-1,m)+C(n-1,m-1);
}*/
LL C(int m,int n)
{
if(m==n||!n) return 1;
else if(c[m][n]) return c[m][n];
else return c[m][n]=C(m-1,n)+C(m-1,n-1);
}
/*LL C(int n,int m){
long long res=1;
for(int i=1;i<=m;i++)
res=res*(LL)(n-m+i)/i;
return res;
}*/
LL gcd(int a,int b){return !b?a:gcd(b,a%b);}
int main()
{
scanf("%d%d",&n,&m);
ans=1LL*C((m+1)*(n+1),3);
ans-=1LL*(m+1)*C(n+1,3);
ans-=1LL*(n+1)*C(m+1,3);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
ans=ans-1LL*(n-i+1)*(m-j+1)*2*(gcd(i,j)-1);
}
printf("%lld",ans);
return 0;
}