标签：UVA306 int 题解 rep cdots Cipher pmatrix ans Array

本文同步发布于 洛谷博客。

题意：多测，给你一个置换 \(A\)，分别对几个字符串 \(s_i\) 执行 \(k_i\) 次，问最后得到的字符串是啥。

我们知道两个置换的乘法：假设有置换 \(f=\begin{pmatrix}a_1 & a_2 & \cdots & a_n\\a_{p_1} & a_{p_2} & \cdots & a_{p_n}\end{pmatrix}\) 和 \(g=\begin{pmatrix}a_{p_1} & a_{p_2} & \cdots & a_{p_n}\\a_{q_1} & a_{q_2} & \cdots & a_{q_n}\end{pmatrix}\)，有 \(f\circ g=\begin{pmatrix}a_1 & a_2 & \cdots & a_n\\a_{q_1} & a_{q_2} & \cdots & a_{q_n}\end{pmatrix}\)。

我们可以用这个乘法运算代替整数乘法实现快速幂，即可在 \(O(n\log{k})\) 的时间复杂度内求出 \(A^k\)。然后注意这个结果是一个置换，要知道最后第 \(i\) 个位置是几，就要看第二行为 \(i\) 的是第几个位置，输出对应字符即可。

//By: Luogu@rui_er(122461)
#include <bits/stdc++.h>
#define rep(x,y,z) for(int x=y;x<=z;x++)
#define per(x,y,z) for(int x=y;x>=z;x--)
#define debug printf("Running %s on line %d...\n",__FUNCTION__,__LINE__)
using namespace std;
typedef long long ll;
const int N = 205; 

int n, k;
template<typename T> void chkmin(T &x, T y) {if(x > y) x = y;}
template<typename T> void chkmax(T &x, T y) {if(x < y) x = y;}
struct Array {
	int p[N];
	Array() {memset(p, 0, sizeof(p));}
	Array(int n) {rep(i, 1, n) p[i] = i;}
	~Array() {}
	friend Array operator * (const Array& a, const Array& b) {
		Array ans = Array(n);
		rep(i, 1, n) ans.p[i] = a.p[b.p[i]];
		return ans;
	}
	friend Array operator ^ (Array a, int k) {
		Array ans = Array(n);
		for(;k;k>>=1,a=a*a) if(k & 1) ans = ans * a;
		return ans;
	}
}a, b, c;

int main() {
	while(scanf("%d", &n) == 1 && n) {
		rep(i, 1, n) scanf("%d", &a.p[i]);
		while(scanf("%d", &k) == 1 && k) {
			string s;
			getline(cin, s); // there is a leading ' '
			rep(i, 1, n) s = s + ' ';
			b = a ^ k;
			rep(i, 1, n) c.p[b.p[i]] = i;
			rep(i, 1, n) printf("%c", s[c.p[i]]);
			puts("");
		}
		puts("");
	}
	return 0;
}

标签：UVA306,int,题解,rep,cdots,Cipher,pmatrix,ans,Array
来源： https://www.cnblogs.com/ruierqwq/p/UVA-306.html

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