【bzoj4591】【Shoi2015超能粒子炮·改】【lucas定理】

Description

曾经发明了脑洞治疗仪&超能粒子炮的发明家SHTSC又公开了他的新发明：超能粒子炮·改--一种可以发射威力更加

强大的粒子流的神秘装置。超能粒子炮·改相比超能粒子炮，在威力上有了本质的提升。它有三个参数n，k。它会

向编号为0到k的位置发射威力为C(n,k) mod 2333的粒子流。现在SHTSC给出了他的超能粒子炮·改的参数，让你求

其发射的粒子流的威力之和模2333。

Input

第一行一个整数t。表示数据组数。

之后t行，每行二个整数n，k。含义如题面描述。

k<=n<=10^18，t<=10^5

Output

t行每行一个整数，表示其粒子流的威力之和模2333的值。

Sample Input

1
5 5

Sample Output

32

题解：根据lucas定理

s[i][j]=lucasS(n/P,k/P-1)*s[n%P][P-1]+lucasC(n/P,k/P)*s[n%P][k%P];

预处理出s[i][j]和c[i][j]即可.

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#define P 2333
#define N 3010
using namespace std;
int c[N][N],s[N][N],T;
long long n,k;
void pre(){
  c[0][0]=1;
  for (int i=1;i<=P;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])%P;
  }
  for (int i=0;i<=P;i++){
    s[i][0]=c[i][0];
    for (int j=1;j<=P;j++)
      s[i][j]=(s[i][j-1]+c[i][j])%P;
  }
}
int lucasC(long long n,long long k){
  if (k==0) return 1;
  return (c[n%P][k%P]*lucasC(n/P,k/P))%P;
}
int lucasS(long long n,long long k){
  if (k<0) return 0; 
  return (lucasS(n/P,k/P-1)*s[n%P][P-1]+lucasC(n/P,k/P)*s[n%P][k%P])%P; 
}
int main(){
  scanf("%d",&T);
  pre();
  while (T--){
  	scanf("%lld%lld",&n,&k);
  	printf("%d\n",lucasS(n,k)); 
  }	
}