B、Infoplane in Tina Town-----(BestCoder Round #51 (div.2) )

Infoplane in Tina Town

 Time Limit: 7000/3500 MS (Java/Others)

 

 Memory Limit: 524288/524288 K (Java/Others)

问题描述

Tina Town 的镇中有一块表面平整光滑的大石头, 当人面向它走进时, 石头的表面就会亮起, 显示出它的用途. 这块石头是前人留下来的东西,是Tina Town算力和控制系统的核心, 同时它也可以显示Tina Town的大小事和行人关注的内容,也可以作为公共的计算机使用 ,极大地方便了人们的生活(特别是上街忘带终端的人们啦).

Tina和Town在这块大石头上玩一个游戏. 石头上首先显示出了1到n排列的n个数,Town则随机交换了一些数, 然后Town通过宏将这个交换的过程记录下来了. Town 问 Tina: 你知道执行几次这个宏, 这些数会恢复最早的排列嘛? Tina是一个蠢蠢的女孩子, 当然不知道啦, 于是她想向你请教呢.

由于答案可能很大,你只需要输出答案对于3*2^{30}+1=32212254733∗2​30​​+1=3221225473(一个素数)的模就行了.

输入描述

第一行一个数$T$,表示数据组数.$T\le 5$

接下来一组数据两行,第一行一个整数$n$,表示排列的长度.n\le 3\times 10^6n≤3×10​6​​

第二行一个$n$个数表示一个排列A_1... A_nA​1​​...A​n​​,保证数字两两不同且对于所有的A_iA​i​​都有1\le A_i\le n1≤A​i​​≤n.

输出描述

对于每组数据输出一行一个数$ans$,表示答案.

输入样例

2
3
1 3 2
6
2 3 4 5 6 1

输出样例

2
6

Hint

数据规模较大建议使用读入优化.

分析：首先，这道题读题怎么都不理解题意，看了几遍管理员的解释总算明白了，ORZ。

题意：给一个数组，经过数字循环后，变成原来1~n的递增序列。数字循环只能是每一次把区间内的第一个数字移到最后（可能也可以把最后的数字移到第一个），几个区间是同时移动的。

第一个数据 1 3 2 -> 1 2 3 长度为2 第二个数据 2 3 4 5 6 1 -> 3 4 5 6 1 2 -> 4 5 6 1 2 3 -> 5 6 1 2 3 4 -> 6 1 2 3 4 5 -> 1 2 3 4 5 6 长度为6

另外数据2 1 4 5 3的过程是 2 1 4 5 3-> 1 2 5 3 4-> 2 1 3 4 5 -> 1 2 4 5 3 -> 2 1 5 3 4-> 1 2 3 4 5 长度为6

官方题解：

给出一个置换，求它的循环长度。数据范围 3\times 10^63×10​6​​。

其实没有什么好说的，就分解成循环求长度的最小公倍数就好了。对于这个模数要用unsigned int存，这个最小公倍数的求法不能用欧几里得，直接每次分解质因数，用线性筛预处理一下就好了。

时间复杂度:$O(n)$.分析

对于第一部分寻找循环每个点都恰好遍历一次,$O(n)$;

第二部分分解质因数，时间复杂度为

\sum_{i=1}^kf(|p_k|)∑​i=1​k​​f(∣p​k​​∣). $f=O(\log n)$,\sum_{i=1}^k|p_k|=n∑​i=1​k​​∣p​k​​∣=n.

显然有f(|p_k|)\le p_kf(∣p​k​​∣)≤p​k​​,所以这部分时间复杂度$O(n)$.

第三部分快速幂加乘法同样分析.

CODE：（这是我直接求最小公倍数，始终没有发现BUG，但一直WA，ORZ）

#include <iostream>
#include <cstdio>
#include <string.h>
#define ll long long
//#define mod 3221225473
using namespace std;

const ll mod=3221225473;
const int maxn=3000005;
ll n,a[maxn],num[maxn];
bool vis[maxn];

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

int main()
{
    //模下gcd是不正确的
    int t; scanf("%d",&t);
    while(t--){
        memset(vis,0,sizeof(vis));
        int cnt=0;
        scanf("%lld",&n);
        for(int i=1;i<=n;i++)
            scanf("%lld",&a[i]);
        for(int i=1;i<=n;i++){
            if(vis[i]) continue;
            ll k=i,sum=0;
            while(!vis[k]){
                vis[k]=1;
                k=a[k];
                sum++;
            }
            num[cnt++]=sum;
        }
        ll gcc=num[0];
        for(int i=1;i<cnt;i++){
            gcc=gcd(gcc,num[i]);
        }
        ll ans=num[0]/gcc;
        for(int i=1;i<cnt;i++){
            ans=(ans*num[i])%mod;
        }
        if(cnt==1) printf("%lld\n",num[0]);
        else printf("%lld\n",ans);
    }
    return 0;
}

CODE：

#include <iostream>
#include <cstdio>
#include <string.h>
#define ll long long
//#define mod 3221225473
using namespace std;

const ll mod=3221225473;
const int maxn=3000005;
ll n,a[maxn],num[maxn];
bool vis[maxn];

ll qpow_mod(ll x,ll y){
    ll ret=1;
    while(y){
        if(y&1)
            ret=(ret*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return ret;
}

int main()
{
    int t; scanf("%d",&t);
    while(t--){
        memset(vis,0,sizeof(vis));
        memset(num,0,sizeof(num));
        int cnt=0;
        scanf("%lld",&n);
        for(int i=1;i<=n;i++)
            scanf("%lld",&a[i]);
        for(int i=1;i<=n;i++){
            if(vis[i]) continue;
            ll k=i,sum=0;
            while(!vis[k]){
                vis[k]=1;
                k=a[k];
                sum++;
            }
            for(int j=2;j<=sum;j++){
                int tmp=0;
                while(sum%j==0){
                    tmp++;
                    sum/=j;
                }
                if(tmp>num[j])
                    num[j]=tmp;
            }
        }
        ll ans=1;
        for(int i=2;i<=n;i++){
            if(num[i])
                ans=(ans*qpow_mod(i,num[i]))%mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}

CODE：（很纳闷，用了输入外挂，反而超时了）

#include <iostream>
#include <cstdio>
#include <string.h>
#define ll long long
//#define mod 3221225473
using namespace std;

const ll mod=3221225473;
const int maxn=3000005;
ll n,a[maxn],num[maxn];
bool vis[maxn];

ll qpow_mod(ll x,ll y){
    ll ret=1;
    while(y){
        if(y&1)
            ret=(ret*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return ret;
}

inline int read(){
    char ch;
    int num=0;
    while((ch=getchar())<'0'||ch>'9');
    while(ch>='0'&&ch<='9')
        num=num*10+(ch-'0'),ch=getchar();
    return num;
}

int main()
{
    int t; scanf("%d",&t);
    while(t--){
        memset(vis,0,sizeof(vis));
        memset(num,0,sizeof(num));
        int cnt=0;
        scanf("%lld",&n);
        for(int i=1;i<=n;i++)
            a[i]=read();
        for(int i=1;i<=n;i++){
            if(vis[i]) continue;
            ll k=i,sum=0;
            while(!vis[k]){
                vis[k]=1;
                k=a[k];
                sum++;
            }
            for(int j=2;j<=sum;j++){
                int tmp=0;
                while(sum%j==0){
                    tmp++;
                    sum/=j;
                }
                if(tmp>num[j])
                    num[j]=tmp;
            }
        }
        ll ans=1;
        for(int i=2;i<=n;i++){
            if(num[i])
                ans=(ans*qpow_mod(i,num[i]))%mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}