Chat Group(快速幂)

 还有一个规律 :C(n,i+1)=C(n,i)*((n-i+1)/i);

杨辉三角 可以直接有2^n的快速幂算出来

#include<iostream>
#include<string>
#include<cstdio>
#define mod 1000000007
using namespace std;
long long extend_gcd(long long a,long long b,long long &x,long long &y)
{
    long long d;
    if(b==0)
    {x=1; y=0;  return a;}
    d=extend_gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
long long mod_reverse(long long a, long long p)
{
    long long x,y;
    long long d=extend_gcd(a,p,x,y);
    if(d==1)
    return (x%mod+mod)%mod;
    else
    return -1;
}
long long pow_quick(long long a,long long p)
{
	long long ans=1;
	while(p>0)
	{
		if(p&1)
        ans=ans*a%mod;
		p>>=1;
		a=a*a%mod;
	}
	return ans;
}
int main()
{
    int t;
    long long n,k;
    long long i;
    int Cas=1;
    scanf("%d",&t);
    while(t--)
    {
       scanf("%lld%lld",&n,&k);
       long long sum1=1;
       long long temp1=1;
	   for(i=1;i<=k-1;i++)
       {
         temp1=(((((temp1%mod)*(n-i+1)%mod)%mod)%mod)*(mod_reverse(i,mod)%mod))%mod;
         temp1%=mod;
         sum1+=temp1;
         sum1%=mod;
       }
       long long s=pow_quick(2,n);
       long long ans1=(s+mod-sum1)%mod;
	   printf("Case #%d: %lld\n",Cas++,ans1);
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值