codeforces Round#369 div2-E ZS and The Birthday Paradox

原创 2016年08月30日 00:04:07

题意: 一年中有2^n天,问有k个人至少有两个人生日在同一天的概率,用分数表示。其中n,k<=10^18 

题解:1e6+3为素数。直接求不好求,我们可以总数-不成立的数。不成立的即为k个人中没有一个人的生日在同一天,即。而当k-1>=mod时,分子中必然含有mod的倍数,取余后分子为0;当k-1<mod时,我们可以直接求解。

由于分母中只含素数2,化简时只需求出分子中2的个数。

注意:k>2^n的情况,直接输出1 1即可

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>

using namespace std;

typedef long long LL;

const int M=1000003;

LL n,K;
LL up,down;
LL inv;
LL Calc_bit(LL K)
{
    if (K==1)return 0;
    return (K/2)+Calc_bit(K/2);
}
LL RP(LL a,LL b)
{
    LL Ans=1;
    for (;b;b>>=1){
        if (b&1)Ans=Ans*a%M;
        a=a*a%M;
    }
    return Ans;
}
void Calc_up()
{
    if (K>=2*M){
        up=0;return ;
    }
    LL simp=Calc_bit(K-1);

    up=1;
    LL bg=RP(2,n);
    for (int i=1;i<=K-1;i++)up=up*(bg-i+M)%M;
    up=up*RP(inv,simp)%M;
}
void Calc_down()
{
    LL simp=Calc_bit(K-1);
    down=RP(2,n);
    down=RP(down,K-1);
    down=down*RP(inv,simp)%M;
}
int main()
{
    inv=RP(2,M-2);
    cin>>n>>K;
    if (K> (1LL<<min(n,60LL)) ){
        puts("1 1");return 0;
    }
    Calc_up();
    Calc_down();
    up=(down-up+M)%M;
    cout<<up<<" "<<down<<endl;
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

Codeforces 711 E. ZS and The Birthday Paradox(数学)——Codeforces Round #369 (Div. 2)

[传送门](http://codeforces.com/contest/711/problem/E) E. ZS and The Birthday Paradoxtime limit per test...

CodeForces 369 div2 E. ZS and The Birthday Paradox 数论

E. ZS and The Birthday Paradox time limit per test 2 seconds memory limit per test 256 m...

Codeforces #369 (Div. 2) E. ZS and The Birthday Paradox (勒让德定理+逆元)

E. ZS and The Birthday Paradox time limit per test 2 seconds memory limit per test 256 megabytes...

CodeForces-711E ZS and The Birthday Paradox(勒让德定理+抽屉原理)

E. ZS and The Birthday Paradox time limit per test 2 seconds memory limit per test 256 m...

Codeforces Problem 711E ZS and The Birthday Paradox(抽屉原理+乘法逆元+费马小定理+组合数+快速幂+概率论)

Codeforces Problem 711E ZS and The Birthday Paradox(抽屉原理+乘法逆元+费马小定理+组合数+快速幂+概率论)...

【28.57%】【codeforces 711E】ZS and The Birthday Paradox

time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandar...

【数论】codeforce711E ZS and The Birthday Paradox

题目大意: 给n和k。 求2^n天,k个人,其中至少两人同一天生日的概率。 1 ≤ n ≤ 1018, 2 ≤ k ≤ 1018 输出分子和分母 MOD 1e6+3 求MOD前,...

Codeforces Round #369 (Div. 2) B. Chris and Magic Square

B. Chris and Magic Square time limit per test 2 seconds memory limit per test 256 megaby...

Codeforces Round #369 (Div. 2) B. Chris and Magic Square【数学,模拟】

题意:输入一个n*n的矩阵,有一个位置为0,问你能否在0处填上一个正整数,使得该矩阵的每一行,每一列,主对角,副对角线上的和都相等,即构成一个幻方。 由于要构成一个幻方,所以只用根据行的和求出要填的数...

Codeforces Round #343 (Div. 2) D. Babaei and Birthday Cake(dp、BIT)

题意: 给定N≤105个蛋糕,编号为1∼N,每个都有体积Vi给定N\le 10^5个蛋糕, 编号为1\sim N, 每个都有体积V_i 任意一个蛋糕都可以放在桌子上,对于蛋糕对(i,j),i...
  • lwt36
  • lwt36
  • 2016年02月26日 23:18
  • 223
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:codeforces Round#369 div2-E ZS and The Birthday Paradox
举报原因:
原因补充:

(最多只允许输入30个字)