【t013】无聊的军官

Time Limit: 1 second
Memory Limit: 32 MB

【问题描述】

每个学年的开始，高一新生们都要进行传统的军训。今年有一个军训教官十分奇怪，他为了测试学员们的反应能力，每次吹哨后学员们都会变换位置。每次左数第I位学员都会站到第ai个位置，经过若干次之后，队伍又会回到原来的样子。你的任务是计算n个人的队伍至少经过多少次之后，队伍恢复到原来样子。

【输入格式】

输入文件的第一位包含一个整数N(0<N10000)，表示队伍的人数。接下来N行，每行一个正整数ai表示左起第i个人接下来出现在左起第ai个位置上。

【输出格式】

仅包括一行，一个正整数M，表示军官居最少的吹哨次数。

【输入样例1】

523451

【输出样例1】

5

【数据规模】

对于30%的数据，有N ≤100；对于100%的数据，有N ≤10000；对于全部数据，答案均在64位整数范围之内。

【题目链接】:http://noi.qz5z.com/viewtask.asp?id=t013

【题意】

中文题..

【题解】

那个变化位置的序列最后肯定能够形成若干个环的;
每个环的大小就是每个人回到原来位置所需要的步骤数;
那么求这几个环的大小的最小公倍数就好了;
这里我先把环的大小这些数字质因数分解一下;
然后按照每个因数出现次数最多的,作为这个因数的出现次数;
然后每个因数^(因数的出现次数)全部乘起来就好;
先求gcd再相乘除它也能求最小公倍数;
但好像会超时.
(没看到题目下面的数据范围，写了个高精度…)

【完整代码】

#include <cstdio>#include <algorithm>#include <cstdio>#include <vector>#include <cmath>using namespace std;#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1#define LL long long#define rep1(i,a,b) for (int i = a;i <= b;i++)#define rep2(i,a,b) for (int i = a;i >= b;i--)#define mp make_pair#define pb push_back#define fi first#define se second#define rei(x) scanf("%d",&x)#define rel(x) scanf("%I64d",&x)typedef pair<int,int> pii;typedef pair<LL,LL> pll;const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};const double pi = acos(-1.0);const int MAXN = 1e4+100;int a[MAXN],n;bool bo[MAXN];vector <int> v,what,zhi;int num[MAXN],ans[MAXN],al;bool is(int x){    double t = x;    int len = int(sqrt(t));    rep1(i,2,len)        if (x%i==0)            return false;    return true;}int main(){    //freopen("F:\\rush.txt","r",stdin);    rei(n);    rep1(i,1,n)        rei(a[i]);    rep1(i,2,n)        if (is(i))            zhi.pb(i);    rep1(i,1,n)        if (!bo[i])        {            int j = i,cnt = 0;            while (!bo[j])            {                bo[j] = true;                cnt++;                j = a[j];            }            v.pb(cnt);        }    int len = v.size();    rep1(i,0,len-1)    {        int now = 2,x = v[i],ii = 0,len2 = zhi.size();        while (ii<=len2-1 && now<=x)        {            now = zhi[ii];            if (x%now==0)            {                int cnt = 0;                while (x%now==0)                {                    cnt++;                    x/=now;                }                if (num[now]==0)                {                    num[now] = cnt;                    what.pb(now);                }                else                    num[now] = max(num[now],cnt);            }            ii++;        }    }    al = 1;    ans[1] = 1;    len = what.size();    rep1(i,0,len-1)    {        int len2 = num[what[i]];        rep1(j,1,len2)        {            int x = 0;            rep1(k,1,al)            {                ans[k]= ans[k]*what[i]+x;                x = ans[k]/10;                ans[k] = ans[k]%10;            }            while (x>0)            {                al++;                ans[al] = x;                x = ans[al]/10;                ans[al] = ans[al]%10;            }        }    }    rep2(i,al,1)        printf("%d",ans[i]);    return 0;}