洛谷-4345 [SHOI2015]超能粒子炮·改

最新推荐文章于 2022-08-05 18:21:52 发布

mkopvec

最新推荐文章于 2022-08-05 18:21:52 发布

阅读量349

点赞数

分类专栏：算法题

本文链接：https://blog.csdn.net/mkopvec/article/details/102766950

版权

算法题专栏收录该内容

311 篇文章 3 订阅

订阅专栏

题目描述
曾经发明了脑洞治疗仪与超能粒子炮的发明家 SHTSC 又公开了他的新发明：超能粒子炮・改——一种可以发射威力更加强大的粒子流的神秘装置。
超能粒子炮・改相比超能粒子炮，在威力上有了本质的提升。它有两个参数n,k，它会向每个编号为0到k（包含两端）的位置iii发射威力为 $C_{n}^{i} mod 2333$ 的粒子流。
现在 SHTSC 给出了他的超能粒子炮・改的参数，让你求出其发射的粒子流的威力之和除以2333所得的余数。
输入格式
第一行一个整数 t表示数据组数。之后 t 行，每行两个整数 n、k，含义如题面描述。
输出格式
t 行，每行一个整数，表示其粒子流的威力之和模 2333 的值。

输入输出样例
输入 #1
3
5 5
10 7
1145 14

输出 #1
32
968
763

说明/提示

解释：卢卡斯定理， $C_n^k \mod p=C_{n/p}^{k/p}*C^{k\mod p}_{n\mod p}$
我们设 $f(n,k)=\sum_{i=0}^kC_{n}^i=\sum_{i=0}^{\frac{k}{k}}C_{\frac{n}{p}}^i*(\sum_{j=0}^{p-1}C_{n \mod p}^j)$
$=\sum_{i=0}^{\frac{k}{p}-1}C_{\frac{n}{p}}^i*(\sum_{j=0}^{p-1}C_{n \mod p}^j)+C_{\frac{n}{p}}^{\frac{k}{p}}*(\sum_{\frac{n}{p}}^{\frac{k}{p}}C_{n \mod p}^j)$
$\mod p,p-1)*f(\frac{n}{p},\frac{k}{p}-1)+C_{n \mod p}^{k \mod p}*f(n \mod p,k \mod p)$
然后我们分块就可以啦

#include<iostream>
#define mod 2333
using namespace std;
typedef long long ll;
ll C[2335][2335]={0};
ll F[2335][2335]={0};
ll Lucas(ll n,ll m){
    if(!m) return 1;
    if(n==m) return 1;
    if(n<m) return 0;
    return C[n%mod][m%mod]*Lucas(n/mod,m/mod)%mod;
}
ll f(ll n,ll k){
    if(k<0) return 0LL;
    if(k==0) return 1LL;
    if(n==0) return 1LL;
    if(n<=mod&&k<=mod) return F[n][k];
    return (F[n%mod][mod-1]*f(n/mod,k/mod-1)%mod+Lucas(n/mod,k/mod)*F[n%mod][k%mod]%mod)%mod;
}
int main(){
    ios::sync_with_stdio(false);
    C[0][0]=1;
    for(int i=1;i<=2333;i++){
        C[i][0]=1;
        for(int j=1;j<=i;j++){
            C[i][j]=C[i-1][j-1]+C[i-1][j];
            C[i][j]%=mod;
        }
    }
    F[0][0]=1;
    for(int i=0;i<=2333;i++){
        F[i][0]=1;
        for(int j=1;j<=2333;j++){
            F[i][j]+=F[i][j-1]+C[i][j];
            F[i][j]%=mod;
        }
    }
    int T=0;cin>>T;
    while(T--){
        ll n,k;
        cin>>n>>k;
        cout<<f(n,k)<<endl;
    }
    return 0;
}