Codeforces711 E. ZS and The Birthday Paradox（数论）

题意：

解法：

$$\begin{gathered} 当k>2^n时,显然概率=1,下面只考虑k<=2^n的情况. \\ 容易想到将问题转化为1-(k个人生日互不相同的概率). \\ k个人生日互不相同的概率=\frac{A(2^n,k)}{2^{nk}} \\ =\frac{{\prod }_{i=2^n-k+1}^{2^n}i}{2^{nk}}=\frac{{\prod }_{i=2^n-k+1}^{2^n-1}i}{2^{n(k-1)}} \\ 分母可以直接快速幂. \\ 对于分子部分,由于需要对1e6+3取模, \\ 那么当分子存在1e6+3时,取模就会变成0, \\ 观察到分子的区间为[2^n-1-k+1,2^n-1],区间长度=k-1, \\ 当k-1>=1e6+3时,一定存在一个1e6+3的倍数, \\ 当k-1<1e6+3时,直接暴力计算分子即可. \\ 还有一个问题是分数A/B需要先约分再取模, \\ 由于分母只包含2,因此gcd一定是2的幂次,指数=分子中2的个数, \\ 如何计算分子中2的个数? \\ 由于分子是连乘形式,因此等价于对每个乘积项算出2的个数, \\ 然后将所有个数相加就是分子中2的个数了. \\ 这里有一个重要结论: \\ 如果对于x<2^n,x能提取出2^d,那么2^n-x也能提取出2^d. \\ 因此计算\prod _{i=2^n-k+1}^{2^n-1}i能提取出来多少个2, \\ 等价于计算\prod _{i=1}^{k-1}i=(k-1)!能提取出多少个2. \\ k-1<=1e6+3,(k-1)!的质因子2数量很容易求了. \\ 最后我证明一下上面的结论: \\ x能提取出2^d,那么x=2^d\ast p,其中p与2互质,即p是奇数, \\ {2}^n-x=2^d\ast 2^{n-d}-2^d\ast p=2^d(2^{n-d}-p), \\ 显然2^{n-d}-p也与2互质,即也是奇数, \\ 因此2^n-x=2^d(2^{n-d}-p)提取出的也是2^d. \end{gathered}$$

code：

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxm=2e6+5;
const int mod=1e6+3;
int n,k;
int cal(int n,int x){
    int ans=0;
    while(n){
        ans+=n/x;n/=x;
    }
    return ans;
}
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;
    }
    return ans;
}
void solve(){
    cin>>n>>k;
    if(n<=60&&k>(1LL<<n)){//k>2^n的情况
        cout<<"1 1"<<endl;return ;
    }else{//其他情况一定k<=2^n,(n>60则2^n>1e18>k)
        //分子=mul[2^n-1-k+1,2^n-1]
        int t=ppow(2,n,mod);//2^n
        int up=1;
        for(int i=1;i<=k-1;i++){
            up=up*(t-1-i+1)%mod;
            if(!up)break;
        }
        //分母=(2^n)^(k-1)
        int down=ppow(ppow(2,n,mod),k-1,mod);
        //约分
        int cnt=cal(k-1,2);//(k-1)!中2的个数
        int inv=ppow(ppow(2,cnt,mod),mod-2,mod);//2^cnt的逆元
        up=(up*inv%mod+mod)%mod;
        down=(down*inv%mod+mod)%mod;
        //1-A/B
        up=(down-up+mod)%mod;
        //输出答案
        cout<<up<<' '<<down<<endl;
    }
}
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);
    solve();
    return 0;
}

---