下面是借助逆元求$C_x^y\bmod p$的部分代码，$jc$ 是阶乘数组已经预处理出来了,求大佬查错

inline int ksm(int ds,int zs)//快速幂
{
	int ret=1;
	while(zs)
	{
		if(zs&1) ret=(ret*ds)%MOD;
		zs>>=1;
		ds=(ds*ds)%MOD;
	}
	return ret;
}
int C(int x,int y)
{
	if(x==0) return 1;
	else return (jc[x]%MOD)*(ksm((jc[y]%MOD)*(jc[x-y]%MOD),MOD-2)%MOD);
}