cf711E. ZS and The Birthday Paradox

最新推荐文章于 2021-06-04 14:25:18 发布

*ACoder*

最新推荐文章于 2021-06-04 14:25:18 发布

阅读量161

点赞数

分类专栏： # 数学杂题

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/102009600

版权

数学杂题专栏收录该内容

93 篇文章 0 订阅

订阅专栏

链接

点击跳转

题解

$\frac{2^n(2^n-1)(2^n-2)\dots(2^n-(k-1))}{2^{nk}}$

最头疼的地方就在于约分

然而幸好 $2^n-x$ 中包含 $2$ 的幂次和 $x$ 中包含 $2$ 的幂次是相等的，所以就可以 $O (l o g k)$ 求出这个幂次了

为啥相等？简单口胡一下：

对于 $x < 2^n$ ，要使 $2^k|(2^n-x)$ ，已经有 $2^k|2^n$ ，则必有 $2^k|x$

如果 $2^(k+1)\not|(2^n-x)$ ，因为 $2^k|2^n$ ，所以必定 $2^k\not|x$

代码

import math
n, k = [int(x) for x in input().split()]
if n<70 and k>2**n:
    print(1,1)
    exit(0)
mod = int(1e6)+3

def fastpow(a,b):
    t, ans = a, 1
    while b:
        if(b&1):
            ans = ans*t%mod
        t = t*t %mod
        b>>=1
    return ans

t=k-1
cnt=0
while t:
    cnt += t>>1
    t>>=1

x=0
t=fastpow(2,n)
if k<mod:
    x=1
    for i in range(1,k):
        x = x*(t-i)%mod
y=fastpow(2,n*(k-1))

inv = fastpow(2,mod-2)
inv = fastpow(inv,cnt)

x=(x*inv%mod+mod)%mod
y=(y*inv%mod+mod)%mod

x=(y-x+mod)%mod

print(x,y)