置换群的题

题目链接

Leonardo's Notebook

题目大意:对于一个给定的置换B,是否存在置换A^2=B。

思路:

        对于一个置换t,它可以由多个无交集的轮换相乘得到。比如(3245176),可以看成三个轮换,即1->3->4->5->1,  2->2,  6->7->6,也可以看成(1345)(2)(67)。可以容易的知道,对于每个轮换都不会影响到其他轮换,是彼此独立的。

        对于每个轮换,有几个性质:

        设轮换为t,它的长度为n,则:

        1.t^k的轮换乘个数为gcd(k,n),且每个轮换分别是轮换t中下标i mod gcd(k,n)=0,1,2…的元素的连接,即以未选的数开始,每隔gcd(k,n)作为该轮换的值。比如(234561)^2是

(246)(351)。

        2.t^n是一次循环,也就是说t^(n+1)=t。

        对于这题幂是2,容易知道对于奇数长度l的轮换,gcd(2,l)等于1,也就是说对于一个奇数长度的轮换的平方也是一个奇数轮换,且元素不变,也就是说对于任意的奇数长度轮换一定可以由某个同长度的奇数轮换的平方得到。

        对于偶数长度l的轮换,gcd(2,l)等于2,也就是说对于一个偶数长度的轮换,它的平方一定是拆成另外两个长度相等的轮换,也就是说对于任意一个偶数长度的轮换只要成对(某个偶数长度的轮换个数也要都是偶数个)才能由A^2分裂出来。

代码:

#include<bits/stdc++.h>
using namespace std;
const int N=3e6+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
//const ll P=2281701377;
const ll P=998244353;


void solve(){
	string s;
	cin>>s;
	s=' '+s;
	int cnt[300]={0};
	int vis[300]={0};
	int b[300];
	for(int i=1;i<s.size();i++){
		b[i]=s[i]-'A'+1;
	}
	for(int i=1;i<s.size();i++){
		if(!vis[i]){
			int n=0;
			int p=i;
			while(!vis[p]){
				n++;
				vis[p]=1;
				p=b[p];
			}
			cnt[n]++;
		}
	}
	for(int i=1;i<=30;i++){
		if((i&1)==0){
			if(cnt[i]&1){
				cout<<"No"<<endl;
				return;
			}
		}
	}
	cout<<"Yes"<<endl;
}
int main(){
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int t=1;
    cin>>t;
    while(t--){
        solve();
    }

}


题目链接

Permutations POJ - 2369

题目大意:问最少几次轮换变回原位

思路:就是将轮换分成几个不相交的轮换乘积,然后这些乘积长度的最小公倍数,求最小公倍数的

时候,最大公约数是前面数的最小公倍数和当前这个数的最多公约数

代码:

//#include<bits/stdc++.h>
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
const int N=3e6+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
//const ll P=2281701377;
const ll P=998244353;
#define int ll
int gcd(int a,int b){
	if(b==0) return a;
	return gcd(b,a%b);
}
int vis[N];
void solve(){
	int n;cin>>n;
	vector<int> v,vv;
	v.push_back(-1);
	for(int i=1;i<=n;i++){
		int x;cin>>x;
		v.push_back(x);
	}
	for(int i=1;i<=n;i++){
		if(!vis[i]){
			int n=0;
			int p=i;
			while(!vis[p]){
				n++;
				vis[p]=1;
				p=v[p];
			}
			vv.push_back(n);
		}
	}
	ll mul=1,d=1;
	for(int i=0;i<vv.size();i++){
		int x=vv[i];
		d=gcd(x,mul);
		mul=mul*x/d;
	}
	cout<<mul;
}
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int t=1;
    //cin>>t;
    while(t--){
        solve();
    }

}


  • 5
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值