记录编号 271330 评测结果 AA
题目名称 [SHOI 2015] 超能粒子炮·改 最终得分 100
用户昵称 Gravatar0 是否通过 通过
代码语言 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;
}