洛谷 P2054 [AHOI2005]洗牌

题目描述

为了表彰小联为Samuel星球的探险所做出的贡献,小联被邀请参加Samuel星球近距离载人探险活动。

由于Samuel星球相当遥远,科学家们要在飞船中度过相当长的一段时间,小联提议用扑克牌打发长途旅行中的无聊时间。玩了几局之后,大家觉得单纯玩扑克牌对于像他们这样的高智商人才来说太简单了。有人提出了扑克牌的一种新的玩法。

对于扑克牌的一次洗牌是这样定义的,将一叠N(N为偶数)张扑克牌平均分成上下两叠,取下面一叠的第一张作为新的一叠的第一张,然后取上面一叠的第一张作为新的一叠的第二张,再取下面一叠的第二张作为新的一叠的第三张……如此交替直到所有的牌取完。

如果对一叠6张的扑克牌1 2 3 4 5 6,进行一次洗牌的过程如下图所示:

从图中可以看出经过一次洗牌,序列1 2 3 4 5 6变为4 1 5 2 6 3。当然,再对得到的序列进行一次洗牌,又会变为2 4 6 1 3 5。

游戏是这样的,如果给定长度为N的一叠扑克牌,并且牌面大小从1开始连续增加到N(不考虑花色),对这样的一叠扑克牌,进行M次洗牌。最先说出经过洗牌后的扑克牌序列中第L张扑克牌的牌面大小是多少的科学家得胜。小联想赢取游戏的胜利,你能帮助他吗?

输入输出格式

输入格式:

 

输入文件中有三个用空格间隔的整数,分别表示N,M,L

(其中0<N≤10^10 ,0 ≤M≤10^10,且N为偶数)。

 

输出格式:

 

单行输出指定的扑克牌的牌面大小。

 

输入输出样例

输入样例#1:
6 2 3
输出样例#1:
6

说明

0<N≤10^10 ,0 ≤M≤10^10,且N为偶数

题解:扩展欧几里得+快速幂

发现经过m次洗牌,原来标号为i的牌在(i*2^m)%(n+1)的位置。

那么问题就是求 (i*2^m)≡l(mod n+1),用扩展欧几里得求。

代码:

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;

LL n,m,l,x,y;

LL ksm(LL x,LL m){
    LL ret=1;
    while(m){
        if(m&1)ret=ret*x%n;
        x=x*x%n;
        m>>=1;
    }
    return ret;
}

LL exgcd(LL a,LL b,LL &x,LL &y){
    if(b==0){
        x=1;y=0;
        return a;
    }
    int t,r=exgcd(b,a%b,x,y);
    t=x;x=y;y=t-a/b*y;
    return r;
}

int main(){
    scanf("%lld%lld%lld",&n,&m,&l);n++;LL k=ksm(2,m);
    LL gcd=exgcd(k,n,x,y);
    x=(x%n+n)%n,gcd=l/gcd*x;
    cout<<gcd%n<<endl;
    return 0;
}
原文地址:https://www.cnblogs.com/zzyh/p/7665694.html