【bzoj1965】: [Ahoi2005]SHUFFLE 洗牌 数论-快速幂-扩展欧几里得

【bzoj1965】: [Ahoi2005]SHUFFLE 洗牌

观察发现第x张牌

当x<=n/2 x=2x

当x>n/2 x=2x-n-1

好像就是 x=2x mod (n+1)

 就好了

/* http://www.cnblogs.com/karl07/ */
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL long long

LL n,m,l,g,x,y;

LL gcd(LL a,LL b){return !b ? a : gcd(b,a%b);}

void ex_gcd(LL a,LL b,LL &x,LL &y){
    if (b==0){x=1; y=0; return;}
    ex_gcd(b,a%b,y,x);
    y-=x*(a/b);
}

LL Q_pow(LL x,LL y){
    LL ans=1;
    while(y){
        if (y&1) ans=ans*x%(n+1);
        x=x*x%(n+1);
        y=(y>>1);
    }
    return ans;
}

int main(){
    scanf("%lld%lld%lld",&n,&m,&l);
    m=Q_pow(2,m);
    g=gcd(n+1,m);
    ex_gcd(m/g,(n+1)/g,x,y);
    x=x*(l/g)%(n+1);
    printf("%lld
",(x+(n+1))%(n+1));
    return 0;
}