[bzoj1965] 洗牌
题意:给你一副\(10^{10}\)级别的扑克牌,一开始按1-n排列,然后每次重排,将原排列分成两份,取第二份的第一个为新排列的第一个,取第一份的第一个为新排列的第二个,依次类推,问你重排m次后,第l位上的牌是多少。
题解:
数论
发现规律:重排一次后,\(x\)的位置变为\(2*x\pmod{n+1}\)
然后设x为l位上的数字,有\(x*2^m{\equiv}l\pmod{n+1}\)
\(2\)在\(n+1\)下的逆元是\(n/2+1\),则有\(x\equiv(n/2+1)^m*l\pmod{n+1}\)
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
ll gi() {
ll x=0,o=1; char ch=getchar();
while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
if(ch=='-') o=-1,ch=getchar();
while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
return o*x;
}
ll mul(ll x, ll y, ll p) {
ll ret=0;
while(y) {
if(y&1) (ret+=x)%=p;
(x+=x)%=p,y>>=1;
}
return ret;
}
ll qpow(ll x, ll y, ll p) {
ll ret=1;
while(y) {
if(y&1) ret=mul(ret,x,p)%p;
x=mul(x,x,p)%p,y>>=1;
}
return ret;
}
int main() {
ll n=gi(),m=gi(),l=gi();
ll inv=qpow(n/2+1,m,n+1);
printf("%lld\n", mul(inv,l,n+1));
return 0;
}