【HDU 3430】中国剩余定理

题意：洗牌，1-n的序列，按照序列1的形式变化，问需要多少次操作可以变成序列2，若无输出-1；

题解：先找到循环规律，发现每位的变化情况若是可以到达，则存在循环节，所以使得第i位成功转换的操作次数是若干个循环节后再操作几次，故可以使用不互质的中国剩余定理求解，套在实际问题中需处理该问题的一些特殊输入。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<algorithm>
#define maxn 550
#define mod 1000000007
using namespace std;
typedef long long ll;
int n;
int a[maxn],op[maxn];
ll m[maxn],r[maxn];
int tmp[maxn],nxt[maxn];
int vis[maxn];

void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y){
    if(b==0){
        x=1;y=0;d=a;
        return;
    }
    ex_gcd(b,a%b,d,y,x);  
    y-=x*(a/b);
}

ll gcd(ll a,ll b){
    return b==0?a:gcd(b,a%b);
}

ll ex_CRT(){
    ll a,b,c,c1,c2,x,y,d;
    a=m[1];  
    c1=r[1];
    for(int i=2;i<=n;i++){
        b=m[i];
        c2=r[i];
        c=c2-c1;
        ex_gcd(a,b,d,x,y);
        if(c%d) return -1;
        ll b1=b/d;
        x=((c/d*x)%b1+b1)%b1;
        c1=a*x+c1;
        a=a*b1;
    }
    if(c1==0){
        c1=1;
        for(int i=1;i<=n;i++){
            c1=c1*m[i]/gcd(c1,m[i]);
        }    
    }
    return c1;
}



int main(){
	scanf("%d",&n);
	while(n!=0){
		bool f=true;
		for(int i=1;i<=n;i++) scanf("%d",&a[i]);
		for(int i=1;i<=n;i++){
			scanf("%d",&op[i]);
			if(op[i]!=i) f=false;
		} 
		if(f){
			printf("0\n");
			scanf("%d",&n);
			continue;
		}
		
		memset(r,-1,sizeof(r));
		memset(vis,0,sizeof(vis));
		for(int i=1;i<=n;i++){
			tmp[i]=i;
			m[i]=0;
		} 
		int cnt=0;
		while(cnt!=n){
			for(int i=1;i<=n;i++){
				nxt[i]=tmp[a[i]];
				if(vis[i]) continue;
				m[i]++;
				if(nxt[i]==i){
					cnt++;
					vis[i]=1;
				} 
				if(nxt[i]==op[i]) r[i]=m[i];
			}
			for(int i=1;i<=n;i++) tmp[i]=nxt[i];
		}
		
		f=false;
		for(int i=1;i<=n;i++){
			if(r[i]<0) {
				f=true;
				break;
			}
			r[i]=r[i]%m[i];
		} 
		if(f){
			printf("-1\n");
			scanf("%d",&n);
			continue;
		} 
	/*	for(int i=1;i<=n;i++){
			printf("m=%d r=%d\n",m[i],r[i]);
		}*/
		printf("%lld\n",ex_CRT());
		
		scanf("%d",&n);
	}

	return 0;
}