@OceanEye7年前
05/24
21:29
代码
无脑用Lucas化简 分类后分段递归处理
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 |
#include <cstdio> #include <algorithm> #include <cstring> #include <cstdlib> #include <cmath> #define ll long long #define p(x) ('0'<=x&&x<='9') #define For(i,l,r) for(int i=l;i<=r;i++) char cc; int C; template <class T> void read( T &x ) { x=0; cc=getchar(); C=1; while(!p(cc)) { if(cc=='-') C=-1; cc=getchar(); } while(p(cc)) { x=x*10+cc-48; cc=getchar(); } x*=C; } using namespace std; // lucas C(n/m)%p = C(n/p,m/p)*C(n%p,m%p)%p #define MOD 2333 #define M MOD+100 int sum[M][M],c[M][M]; void Dec(int &x) { while(x>=MOD) x-=MOD; } void pre() { c[0][0]=sum[0][0]=1; For(i,1,MOD-1) sum[0][i]=1; For(i,1,MOD-1) { c[i][0]=sum[i][0]=1; for(int j=1;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1]),Dec(c[i][j]); for(int j=1;j<MOD;j++) sum[i][j]=(sum[i][j-1]+c[i][j]),Dec(sum[i][j]); } } int calc_single(ll n,ll m) { if(m<0||n<m) return 0; if(m<MOD&&n<MOD) return c[n][m]; int ret=1,_,__; for(;n&&m&&ret;n/=MOD,m/=MOD) { _=n%MOD; __=m%MOD; if(_<__) return 0; ret*=c[_][__]; ret%=MOD; } return ret; // return calc_single(n/MOD,m/MOD)*c[n%MOD][m%MOD]%MOD; } int calc(ll n,ll k) { if(k<0) return 0; return (calc(n/MOD,k/MOD-1)*sum[n%MOD][MOD-1]%MOD+sum[n%MOD][k%MOD]*calc_single(n/MOD,k/MOD)%MOD)%MOD; } int main() { ll t,n,k; pre(); for(read(t);t;t--) { read(n); read(k); printf("%d\n",calc(n,k)); } return 0; } |