根据生日攻击理论的算法,我们很容易得到计算答案的公式,用1减去没有人生日相同的概率
首先我们要考虑约分
因为分母只有2这一个因子,所以我们只要考虑分子里面有多少个2就可以了
显然 所含的2的数目与i所含的2的数目相同
所以我们只要求1~k-1中有多少个2,我们只要枚举2,4,8…一个一个统计一遍就好
然后我们发现了一个重要的性质
如果 ,因为分子是若干个连续的数相乘,所以其中必然有一个是1e6+3的倍数,又因为gcd只有因子2,所以约分完之后分子仍然是1e6+3的倍数,取模之后必为0
这种情况分子就是分母,把分母用快速幂求一下就好
如果 ,就暴力算一下就好
#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdlib>
#include <utility>
#include <cctype>
#include <algorithm>
#include <bitset>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <deque>
#include <stack>
#include <cmath>
#define LL long long
#define LB long double
#define x first
#define y second
#define Pair pair<int,int>
#define pb push_back
#define pf push_front
#define mp make_pair
#define LOWBIT(x) x & (-x)
using namespace std;
const int MOD=1e6+3;
const LL LINF=2e16;
const int INF=1e9;
const int magic=348;
const double eps=1e-10;
const double pi=3.14159265;
inline int getint()
{
char ch;int res;bool f;
while (!isdigit(ch=getchar()) && ch!='-') {}
if (ch=='-') f=false,res=0; else f=true,res=ch-'0';
while (isdigit(ch=getchar())) res=res*10+ch-'0';
return f?res:-res;
}
inline LL getLL()
{
char ch;LL res;bool f;
while (!isdigit(ch=getchar()) && ch!='-') {}
if (ch=='-') f=false,res=0; else f=true,res=ch-'0';
while (isdigit(ch=getchar())) res=res*10+ch-'0';
return f?res:-res;
}
LL n,k;
LL Gcd,Inv;
LL Up,Down;
inline LL quick_pow(LL x,LL y)
{
x%=MOD;LL res=1;
while (y)
{
if (y&1) res=(res*x)%MOD,y--;
x=(x*x)%MOD;y>>=1;
}
return res;
}
inline void calc_gcd()
{
Gcd=1;Inv=quick_pow(2,MOD-2);
for (register LL i=2;i<=k-1;i<<=1)
Gcd=(Gcd*quick_pow(Inv,(k-1)/i))%MOD;
}
namespace brute
{
inline void solve()
{
LL i,cnt=0,starter=quick_pow(2,n);
for (i=(starter-1+MOD)%MOD,Up=1,cnt=1;cnt<=k-1;i=(i-1+MOD)%MOD,cnt++) Up=(Up*i)%MOD;
Up=(Up*Gcd)%MOD;
Down=quick_pow(quick_pow(2,n),k-1);Down=(Down*Gcd)%MOD;
Up=(Down-Up+MOD)%MOD;
printf("%lld %lld\n",Up,Down);
}
}
namespace big
{
inline void solve()
{
Down=quick_pow(quick_pow(2,n),k-1);Down=(Down*Gcd)%MOD;
printf("%lld %lld\n",Down,Down);
}
}
int main ()
{
n=getLL();k=getLL();
if (n<=63 && (1ll<<n)<k) {printf("1 1\n");return 0;}
calc_gcd();
if (k<=1e6+48) {brute::solve();return 0;}
big::solve();return 0;
}