超能粒子炮·改 HYSBZ - 4591

题目链接：

http://www.lydsy.com/JudgeOnline/problem.php?id=4591

就是求：

                   C(n,0)+C(n,1)+...+C(n,k)mod2

lucas  定理： C(n,m)%P=C(n/p,m/p)*C(n%p,m%p)%p；

公式推了半天，结果思路混乱，不知道推出了个啥，哎，果断搜了一波代码，看了半天也不是很懂，先记录一下吧！

这个博客写的可以的：点击打开链接

果然还是很菜，先贴代码吧：

#include<cstdio>
#include<cmath>
#include<iostream>
#define maxn 2500
using namespace std;

const int MOD=2333;
typedef long long ll;
ll c[maxn][maxn],s[maxn][maxn];

void init()
{
   c[0][0]=1;
   for(int i=1;i<=MOD;i++)
   {
       c[i][0]=c[i][i]=1;
       for(int j=1;j<=i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%MOD;
   }
   for(int i=0;i<=MOD;i++)
   {
       s[i][0]=c[i][0];
       for(int j=1;j<=MOD;j++)
       {
            s[i][j]=(s[i][j-1]+c[i][j])%MOD;
       }
    }
}
ll luckc(ll n,ll k)
{
    if(k==0)  return 1;
    return (c[n%MOD][k%MOD]*luckc(n/MOD,k/MOD))%MOD;
}
ll lucks(ll n,ll k)
{
    if(k<0)  return 0;
    return (lucks(n/MOD,k/MOD-1)*s[n%MOD][MOD-1]+luckc(n/MOD,k/MOD)*s[n%MOD][k%MOD])%MOD;
}
int main()
{
   int t;
   ll n,k;
   init();
   scanf("%d",&t);
   while(t--)
   {
       scanf("%lld%lld",&n,&k);
       printf("%lld\n",lucks(n,k));
   }
   return 0;
}

keep driving