置换群快速幂运算

强烈推荐阅读潘震皓的《置换群快速幂运算》，阅读完之后可以尝试做下面的三道题目~~

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