bzoj1965: [Ahoi2005]SHUFFLE 洗牌

链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1965

题意：中文题。。。

分析：一个置换的题。。我们对样例进行简单的分析可知经过一次变换有(1->2,2->4,3->6,4->1,5->3,6->5)，然后我们就能发现p[a]=a*2%(n+1)。。然后我们设答案为X即X=p[l]，所以有X*(2^m)=l%(n+1)并且2模n+1的逆元为n/2+1，所以我们将等式变换一下就有X=l*(n/2+1)^m%(n+1)，然后就只要用快速幂处理这个式子就行了，同时要注意数据是10^10可能乘爆，用快速乘即可。

代码：

#include<map>#include<set>#include<cmath>#include<queue>#include<math.h>#include<cstdio>#include<vector>#include<string>#include<cstring>#include<iostream>#include<algorithm>#pragma comment(linker, "/STACK:102400000,102400000")using namespace std;const int N=100010;const int MAX=151;const int MOD1=100000007;const int MOD2=100000009;const double EPS=0.00000001;typedef long long ll;const ll MOD=5000011;const ll INF=10000000010;typedef unsigned long long ull;ll mul(ll a,ll b,ll M) {    ll ret=0;    for (ll i=b;i;i>>=1,a=(a+a)%M)    if (i&1) ret=(ret+a)%M;    return ret;}ll qpow(ll a,ll b,ll M) {    ll ret=1;    for (ll i=b;i;i>>=1,a=mul(a,a,M))    if (i&1) ret=mul(ret,a,M);    return ret;}int main(){    ll n,m,l,ans;    scanf("%lld%lld%lld", &n, &m, &l);    ans=l*qpow(n/2+1,m,n+1);    printf("%lld\n", (ans+n+1)%(n+1));    return 0;}/*6 2 3*/