TOJ 4111 Binomial efficient calculate C(n,k)

For Mod 2^n,you can also use &(2^n-1),it will more faster.

M % (2^n) = M & (2^n-1);

It's use a way to get all quality factors of N!.

The portal:http://acm.tju.edu.cn/toj/showp4111.html

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>

const long long  N = 1000000;
long long tag[1000006],p[1000006],numtot;
long long Mod;

long long a1[1000006];

long long Fast_Power(long long n,long long k){
    long long ans = 1;
    while(k){
        if(k & 1) ans = ( ans * n ) % Mod;
        k /= 2;
        n = ( n * n ) % Mod;
    }
    return ans;
}

void Get_prime()
{
    long long  cnt = 0;
    memset(tag,0,sizeof(tag));
    for (long long  i = 2; i <= N; i++)
    {
        if (!tag[i])    p[cnt++] = i;
        for (long long  j = 0; j < cnt && p[j] * i <= N; j++)
        {
            tag[i*p[j]] = 1;
            if (i % p[j] == 0)
                break;
        }
    }
    numtot = cnt;
}

void Deal_with(){
    long long  T,n,k;
    scanf("%lld",&T);
    while(T--){
        scanf("%lld %lld",&n,&k);
        //memset(a1,0,sizeof(a1)); Can't use memset,It's TLE.
        int tempi;
        for(long long  i=0;i<numtot;i++){
            if(p[i] > n)break;
            a1[p[i]] = 0;
            for(long long  j=p[i];j<=n;j *= p[i])a1[p[i]] += n / j;
            for(long long  j=p[i];j<=k;j *= p[i])a1[p[i]] -= k / j;
            for(long long  j=p[i];j<=n-k;j *= p[i])a1[p[i]] -= (n-k) / j;
            tempi = i;
        }
        Mod = 1;
        for(long long  i=1;i<=32;i++){
            Mod *= 2;
        }
        long long ans = 1;
        for(long long  i=0;i<=tempi;i++){//It should use tempi,Because if i > tempi,it's no initialization.
            ans = ( ans * Fast_Power(p[i],a1[p[i]]) ) % Mod;
        }
        printf("%lld\n",ans);
    }
}

int main(void){
    Get_prime();
    Deal_with();
    return 0;
}


1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 、5资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值