标签：AtCoder ch qp read Contest printf Mul AGC031 define

原文链接https://www.cnblogs.com/zhouzhendong/p/AGC031D.html

扯淡

比赛的时候看到这题之后在草稿纸上写下的第一个式子就是 

$$f(p,q) = pq^{-1}$$

然后就再也没有改过。

发现了一堆奇奇怪怪的性质可是一直没有用。

直到官方题解出来的时候：

$$\Huge f(p,q) = qp^{-1}$$

题解

我们可以把前面的几个 $a_i$ 写出来。

$$\begin {eqnarray*}a_1 &=& p\\a_2 &=& q\\ a_3 &=& qp^{-1} \\ a_4 &=& qp^{-1} q^{-1}\\ a_5&=&qp^{-1}q^{-1}pq^{-1}\\a_6&=&qp^{-1}q^{-1}ppq^{-1}\\a_7&=&qp^{-1}q^{-1}pqpq^{-1}\\a_8&=&qp^{-1}q^{-1}pqp^{-1}qpq^{-1}\end{eqnarray*}$$

设

$$A = qp^{-1}q^{-1}p$$

则可以归纳证明：

$$\forall i>6, a_i = Aa_{i-6} A^{-1}$$

于是直接算一下置换 $A$ 的复合幂就好了。

算复合幂只要写成轮换的形式就可以做到时间复杂度 $O(n)$ ，我偷懒写了 $O(n\log n)$ 的快速幂。

代码

#pragma GCC optimize("Ofast","inline")
#include <bits/stdc++.h>
#define clr(x) memset(x,0,sizeof (x))
#define For(i,a,b) for (int i=a;i<=b;i++)
#define Fod(i,b,a) for (int i=b;i>=a;i--)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define _SEED_ ('C'+'L'+'Y'+'A'+'K'+'I'+'O'+'I')
#define outval(x) printf(#x" = %d\n",x)
#define outvec(x) printf("vec "#x" = ");for (auto _v : x)printf("%d ",_v);puts("")
#define outtag(x) puts("----------"#x"----------")
#define outarr(a,L,R) printf(#a"[%d...%d] = ",L,R);\
						For(_v2,L,R)printf("%d ",a[_v2]);puts("");
using namespace std;
typedef long long LL;
typedef vector <int> vi;
LL read(){
	LL x=0,f=0;
	char ch=getchar();
	while (!isdigit(ch))
		f|=ch=='-',ch=getchar();
	while (isdigit(ch))
		x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	return f?-x:x;
}
const int N=100005;
int n,k;
vi p,q,ip,iq,a[10];
vi Inv(vi A){
	vi B(n);
	For(i,0,n-1)
		B[A[i]]=i;
	return B;
}
vi Mul(vi A,vi B){
	vi C(n);
	For(i,0,n-1)
		C[i]=A[B[i]];
	return C;
}
vi Pow(vi x,int y){
	vi ans;
	For(i,0,n-1)
		ans.pb(i);
	for (;y;y>>=1,x=Mul(x,x))
		if (y&1)
			ans=Mul(ans,x);
	return ans;
}
int main(){
	n=read(),k=read();
	For(i,1,n)
		p.pb(read()-1);
	For(i,1,n)
		q.pb(read()-1);
	ip=Inv(p),iq=Inv(q);
	a[1]=p,a[2]=q;
	For(i,3,6)
		a[i]=Mul(a[i-1],Inv(a[i-2]));
	int len=(k-1)/6;
	vi cir=Pow(Mul(q,Mul(ip,Mul(iq,p))),len);
	vi rem=a[k-len*6];
	vi res=Mul(cir,Mul(rem,Inv(cir)));
	For(i,0,n-1)
		printf("%d ",res[i]+1);
	return 0;
}

标签：AtCoder,ch,qp,read,Contest,printf,Mul,AGC031,define
来源： https://www.cnblogs.com/zhouzhendong/p/AGC031D.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。