每个人都知道有一个电脑游戏名称“炉石”,最近xzz非常喜欢这个游戏。如果你想赢,你需要能力和好运。
有一个咒语名称为“奥术飞弹”,可以在所有敌人之间随意分配3种伤害。
xzz有一个“奥术飞弹加号”可以在所有敌人中随机分配n点伤害。敌方战场只有两个角色:敌方英雄和敌方奴才。
现在我们假设敌人的英雄有无限的健康,敌人的仆从有健康。
xzz想要知道“奥术飞弹加号”杀死敌人小兵的概率(在“奥术飞弹加”之后,小兵的生命值不超过0),你能帮助他吗?
输入
输入的第一行包含一个整数T,即测试用例的数量。T测试案例如下。(1≤T≤100000)
每个测试用例由一行包含两个整数n和m组成。(0≤m≤n≤1000)
产量
对于每个测试用例,因为概率可以写成x / y,请输出x * y ^ -1 mod 1000000007.(y * y ^-1≡1mod 10 ^ 9 + 7)
样本输入
2 3 3 2 1
样本输出
125000001 750000006
题解:
分子为求c[n][m]+c[n][m+1]+···+c[n][n]的和
分母为2^n,(但要求逆元,用费马小定理)
组合数利用杨辉三角打表,注意取余
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<math.h>
#include<queue>
#include<bitset>
typedef long long ll;
const int xmax=1e6+100;
const int mod=1000000007;
using namespace std;
/**李洛克:只要能够靠努力的话,天才也是可以被打败的,那个时候的你曾经这样的跟我说过!
我想证明,即便是不会忍术,不会幻术,也能成为一名优秀的忍者,这是我生命的全部,这就是我的忍道! ——岸本齐史**/
ll qpow(ll a,ll b)
{
ll ans=1;
while(b)
{
if(b&1) ans=(ans*a)%mod;
a=(a*a)%mod;
b>>=1;
}
return ans%mod;
}
ll c[1037][1030];
ll res[1030];
void db()
{
c[0][0]=1;
c[1][0]=1;c[1][1]=1;c[0][1]=1;
for(int i=2;i<=1021;i++)
for(int j=0;j<=i;j++)
{
if(j==0)
c[i][j]=1;
else
c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
}
for(int i=1;i<=1027;i++)
for(int j=1;j<=i;j++)
c[i][j]=(c[i][j]+c[i][j-1])%mod;
for(int i=1;i<=1027;i++)
{
ll pr=qpow(2,i);
res[i]=qpow(pr,mod-2);
}
}
int main()
{
db();
int t;
int n,m;
cin>>t;
while(t--)
{
ll ans=0;
cin>>n>>m;
if(m==0)
printf("1\n");
else if(n<m)
printf("0\n");
else
{
ans=(c[n][n-m]*res[n])%mod; ///因为前面打表时求的是1-i的和,所以这里输出c[n][n-m]
printf("%lld\n",ans);
}
}
return 0;
}