强烈推荐阅读潘震皓的《置换群快速幂运算》,阅读完之后可以尝试做下面的三道题目~~
POJ 1721 CARDS
这道题目在上面的论文中已经讲述了,因为题目的每一次洗牌便是当前的置换的一次的平方运算,因此经过S次洗牌之后,相当于做了(T)^(2^S)的运算,由于题目中已经假设原来的牌的数目为奇数,因此对于置换的平方运算gcd(n,2)=1,因此对于每次平方运算置换内的循环是不会产生分裂的,仍旧为一个大循环。我们知道对于一个长度为n的循环T^n整数幂运算后,置换变成了恒等置换,因此在计算(T)^(2^S)的运算时,只需要计算(T)^((2^S)%(2^A))即可(2^A=1(mod n))。题目给出的是目标序列和洗牌的次数S,求解原先的序列即(T)^(1/(2^S)),即置换的分数幂运算,由上面我们知道((T)^(1/2^S))^(2^A)仍为该置换,因此我们可以将上述的式子变成(T)^(2^(A-S%A))的整数幂运算,这便是原先的牌的序列。对于平方运算我们直接采用模拟的方法求解即可~~
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int elem[1010], t_elem[1010];
int N, S, A;
int power_mod(int a, int b, int c)
{
int ret = 1;
while(b)
{
if(b & 1) ret = ret * a % c;
a = a * a % c;
b >>= 1;
}
return ret;
}
void solve(int cnt)
{
for(int i = 1; i <= cnt; ++i)
{
for(int j = 1; j <= N; ++j)
{
t_elem[j] = elem[elem[j]];
}
for(int j = 1; j <= N; ++j)
{
elem[j] = t_elem[j];
}
}
}
int main()
{
//freopen("aa.in", "r", stdin);
//freopen("bb.out", "w", stdout);
while(scanf("%d %d", &N, &S) != EOF)
{
for(int i = 1; i <= N; ++i)
{
scanf("%d", &elem[i]);
}
for(A = 1; A <= N && power_mo