快速幂模板,仅供参考
ll pow_mod(ll a,ll b,ll c)//a^b%c
{
ll ans=1;
ll base=a%c;
while(b)
{
if(b&1)
ans=(ans*base)%c;
base=(base*base)%c;
b>>=1;
}
return ans;
}
很多人可能会对b&1和b>>=1这两步操作疑惑,这两步运用到的是位运算,&表示二进制的每一位数进行与运算,&实际上相当于and,也就是只有1&1=1,其他都等于0。了解到这点后我们可以去看b&1,b的二进制未知,1的八位二进制为0000 0001,那么想要满足b&1>=1的情况下,就是看b的二进制最后一位是否为1,也就是奇数,所以这步相当于b%2。接着b>>=1,是将b的二进制数向右移一位也就是进行除以2的操作,所以这步相当于b/=2。
组合数模板,n中选m个,取余p
ll c(int n,int m)
{
ll x=1,y=1;
for (int i=n;i>=n-m+1;i--) x=x*i;
for (int i=1;i<=m;i++) y=y*i;
return x/y;
}