2020牛客第五场E题 lcm+大数模板+一点点思维

E题
题意：给你一段序列p[i]，然后让你在1~n的序列中任意挑出一组使得a[i]，使得经过有限次的循环后令a数组变为严格递增的数组，同时有限次的循环指的是创建一个辅助数组b，另b=a，然后a[i]=b[p[i]],也就是一步一步进行置换，然后最后得到的数组要是1到n，问你最多有多少个这样的数组a。
思路：思维很简单，因为你选择的a无论如何都会变成1234…n,所以可以设你第一个选择的数组就是123…n，然后对这个数组进行循环，可以肯定，经过有限次的循环之后，这个数组有会变成123…n，同时这个循环所有经过的数组都是数组a的一个解，观察可得出，最后的ans就是他们每个环的最小公倍数，但是，注意但是，就是这个大数的lcm的板子，卡了我整场比赛 。

代码如下：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <vector>
#include <cmath>
#include <queue>
#include <set>
#include <map>
#include <stack>
#include <unordered_map>
using namespace std;
typedef unsigned long long ull;
typedef pair<int, int> pii;
#define rep(i, a, n) for(int i = a; i < n; i++)
#define per(i, a, n) for(int i = n-1; i >= a; i--)
#define IOS std::ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
#define INF 1ll<<60
#define fopen freopen("file.in","r",stdin);freopen("file.out","w",stdout);
#define fclose fclose(stdin);fclose(stdout);
const int maxn = 1e5+10;
typedef long long ll;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
    while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
    return x*f;
}
int pre[maxn];
int n;
/*
这里用并查集是为了找他们有多少个环，同时环内有多少个数，应该还有更好的方法，
我觉得这个方法挺傻的。。。之后还用了map去记录，然后用一个数组存起来
*/
int find(int x){
    return x==pre[x]?x:find(pre[x]);
}
/*
	从这里开始是大数+lcm的板子
*/
vector<int> g_primeVector = { 2, 3, 5, 7, 11 }; // 素数表
bool UpdatePrimeVec(int n) // 更新素数表
{
    if (g_primeVector.back() >= n){
        return false;
    }
    for (int num = g_primeVector.back() + 1; num <= n; num++)
    {
        for (auto prime : g_primeVector)
        {
            if (num % prime == 0) {
                break;
            }
            if (prime > sqrt(num)) {
                g_primeVector.push_back(num);
                break;
            }
        }
    }
    return true;
}
vector<int> Mutiple(vector<int> num1, int num2)//大数乘法
{
    vector<int> res;
    int c = 0;
    for (auto i = num1.begin(); i != num1.end(); i++) {
        int tmp = *i * num2 + c;
        c = tmp / 10;
        tmp = tmp % 10;
        res.push_back(tmp);
    }
    while (c) {
        res.push_back(c % 10);
        c /= 10;
    }
    return res;
}
vector<int> GetLeastComMult(int arr[], int m) //获取最小公倍数
{
    vector<int> leastComMult(1, 1);
    UpdatePrimeVec(n);
    vector<int> multTimes(g_primeVector.size(), 0);
    for (int i = 1; i <= m; i++)
    {
        int num = arr[i];
        for (int j = 0; j < g_primeVector.size() && num > 0; j++)
        {
            int time = 0;
            while ((num % g_primeVector[j] == 0) && (num > 0))
            {
                time++;
                num /= g_primeVector[j];
            }
            if (multTimes[j] < time) {
                multTimes[j] = time;
            }
        }
    }
    for (int m = 0; m < multTimes.size(); m++)
    {
        while (multTimes[m]--)
        {
            leastComMult = Mutiple(leastComMult, g_primeVector[m]);
        }
    }
    return leastComMult;
}
 
void CoutVec(vector<int> numPlus)// 输出大数
{
    for (int i = numPlus.size() - 1; i >= 0; i--)
    {
        cout << numPlus[i];
    }cout << ::endl;
}
/*
	到这里结束，同时输出lcm之后的数
*/
int a[maxn];
int arr[maxn];
int main(){
    n=read();
    rep(i,1,n+1) pre[i]=i;
    rep(i,1,n+1){
        a[i]=read();
        int r1=find(a[i]), r2=find(i);
        if(r1!=r2){
            if(r1>r2) pre[r2] = r1;
            else pre[r1] = r2;
        }
    }
    std::map<int, int> m;
    rep(i,1,n+1){
        m[find(i)]++;
    }
    std::map<int, int>::iterator it;
    int pos = 1;
    for(it=m.begin(); it != m.end(); it++){
        arr[pos++] = it->second;
    }
    CoutVec(GetLeastComMult(arr, pos));
    return 0;
}