记录编号 |
271330 |
评测结果 |
AA |
题目名称 |
[SHOI 2015] 超能粒子炮·改 |
最终得分 |
100 |
用户昵称 |
0 |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.431 s |
提交时间 |
2016-06-15 20:28:38 |
内存使用 |
41.85 MiB |
显示代码纯文本
#include<cstdio>
#include<map>
using namespace std;
typedef long long ll;
const int mod=2333;
ll n,m;
int C[mod+1][mod+1],sum[mod+1][mod+1];
int Lucas(ll n,ll m)
{
if(n<mod&&m<mod) return C[n][m];
return C[n%mod][m%mod]*Lucas(n/mod,m/mod)%mod;
}
int work(ll n,ll m)
{
if(n<mod) return sum[n][m];
ll ret=0,s=(m+1)/mod;
ret=sum[n%mod][mod-1]*work(n/mod,s-1)%mod;
if((m+1)%mod) ret=(ret+sum[n%mod][m%mod]*Lucas(n/mod,s))%mod;
return ret;
}
int main()
{
freopen("cnlzb.in","r",stdin);
freopen("cnlzb.out","w",stdout);
int T=0;
for(int i=0;i<mod;i++){
sum[i][0]=C[i][0]=1;
for(int j=1;j<=i;j++){
C[i][j]=C[i-1][j-1]+C[i-1][j];
if(C[i][j]>=mod) C[i][j]-=mod;
}
sum[i][0]=C[i][0];
for(int j=1;j<mod;j++){
sum[i][j]=sum[i][j-1]+C[i][j];
if(sum[i][j]>=mod) sum[i][j]-=mod;
}
}
scanf("%d",&T);
while(T--){
scanf("%lld%lld",&n,&m);
printf("%d\n",work(n,m));
}
return 0;
}