Atcoder agc031D

漏了一题可还行。。。这个显然是我不会的群论题，虽然也不用太高深的知识，可能说是找规律题更合适。
用置换的符号，可以发现$f(p,q)=q\cdot p^{-1}$，这里$\cdot$指的是置换的复合。
于是有：

• $a_1=p$
• $a_2=q$
• $a_3=q{p}^{-1}$
• $a_4=q{p}^{-1}{q}^{-1}$
• $a_5=q{p}^{-1}{q}^{-1}p{q}^{-1}$
• $a_6=q{p}^{-1}{q}^{-1}pp{q}^{-1}$
• $a_7=q{p}^{-1}{q}^{-1}pqp{q}^{-1}=(q{p}^{-1}{q}^{-1}p)p(q{p}^{-1}{q}^{-1}p{)}^{-1}$
• $a_8=q{p}^{-1}{q}^{-1}pq{p}^{-1}qp{q}^{-1}=(q{p}^{-1}{q}^{-1}p)q(q{p}^{-1}{q}^{-1}p{)}^{-1}$
  令$A=q{p}^{-1}{q}^{-1}p$，那么有$a_7=A{a}_1{A}^{-1}$，$a_8=A{a}_2{A}^{-1}$，于是容易用归纳法证明$a_n=A{a}_{n-6}{A}^{-1}(n>6)$。这样可以快速知道$a_k=A^{\lfloor \frac{k-1}{6}\rfloor }{a}_{(k-1)\mathrm{m}\mathrm{o}\mathrm{d}6+1}(A^{-1}{)}^{\lfloor \frac{k-1}{6}\rfloor }$。
  置换的复合，求逆，$k$次幂都是$\mathcal{O}(N)$的，于是总时间复杂度仍为$\mathcal{O}(N)$。

#include <bits/stdc++.h>

using namespace std;

bool vis[100005];
int now[100005];

struct Per {
  int *num,n;
  Per() {}
  Per(int a):n(a) {num=new int[n+1];}
  Per inv() {
  	Per b(n);
  	for(int i=1;i<=n;i++) b.num[num[i]]=i;
  	return b;
  }
  Per operator * (Per b) {
  	Per c(n);
  	for(int i=1;i<=n;i++) c.num[i]=num[b.num[i]];
  	return c;
  }
  Per operator ^ (int k) {
  	Per b(n);
  	for(int i=1;i<=n;i++) vis[i]=0;
  	for(int i=1;i<=n;i++)
  	  if (!vis[i]) {
  	  	  int x=i,cnt=0;
  	  	  do {
  	  	  	now[++cnt]=x;
  	  	  	vis[x]=1;
  	  	  	x=num[x];
		  } while (x!=i);
		  for(int j=1;j<=cnt;j++) b.num[now[j]]=now[(j+k-1)%cnt+1];
		}
	return b;
  }
};

Per fir[7];

int main() {
  int n,k;
  scanf("%d%d",&n,&k);
  fir[1]=Per(n);
  fir[2]=Per(n);
  for(int i=1;i<=n;i++) scanf("%d",&fir[1].num[i]);
  for(int i=1;i<=n;i++) scanf("%d",&fir[2].num[i]);
  for(int i=3;i<=6;i++) fir[i]=fir[i-1]*fir[i-2].inv();
  Per a=fir[2]*fir[1].inv()*fir[2].inv()*fir[1];
  Per ans=(a^((k-1)/6))*fir[(k-1)%6+1]*(a.inv()^((k-1)/6));
  for(int i=1;i<=n;i++) printf("%d ",ans.num[i]);
  printf("\n");
  return 0;
}