BZOJ4591 [Shoi2015]超能粒子炮·改

2019-04-14 15:33发布

由于模数较小,考虑一下卢卡斯定理 在模p意义下,将式子化简如图 化简过程中由于带模的项可以预处理,所以把带除的项合并到一起

n,k小于p的S可以预处理,C可以卢卡斯算 #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define MAXN 2500 #define MAXM 1010 #define ll long long #define INF 1000000000 #define MOD 2333 #define eps 1e-8 ll fac[MAXN],ine[MAXN]; ll sum[MAXN][MAXN]; ll C(ll n,ll m){ if(!m){ return 1; } if(m>n){ return 0; } if(m