codeforce 954 div3 G2题

思路:

质因子分解可以顺着分解,也可以逆着分解

即找到每一个数字的倍数,再找到每一个数字的因数

cpp 复制代码
const int N = 5e5+10;
vector<int> ff[N];
vector<int> f[N];
vector<int> g[N];

void solve(){
    int n;
    cin>>n;
    vector<int> nums(n+1,0);
    ll ans = 0;
    for(int i = 1;i<=n;i++){
        cin>>nums[i];
    }
    for(int i = 1;i<=n;i++){
        int w = gcd(i,nums[i]);
        int x = i/w,y = nums[i]/w;
        //cout<<x<<" "<<y<<" eeeeee"<<endl;
        if(nums[i]%i == 0)  ans--;
        f[x].emplace_back(y);
        g[y].emplace_back(x);
    }
    vector<int>cnt(n+1,0);
    for(int i = 1;i<=n;i++){
        int m = f[i].size();
        if(m == 0)  continue;
        sort(f[i].begin(),f[i].end());
        for(int j = i;j<=n;j+=i){
            for(auto x : g[j]){
                cnt[x]++;
            }
        }
        for(int j = 0;j<m;j++){
            int k = 1;
            while(j<m-1 && f[i][j] == f[i][j+1]){
                j++;
                k++;
            }
            int y = f[i][j];
            ll sum = 0;
            for(auto son : ff[y]){
                sum += cnt[son];
            }
            sum *= k;
            ans += sum;
        }


        for(int j = i;j<=n;j+=i){
            for(auto x : g[j]){
                cnt[x]--;
            }
        }
    }
    for(int i = 1;i<=n;i++){
        f[i].clear();
        g[i].clear();
    }
    cout<<ans/2<<endl;
}
相关推荐
2401_841495641 小时前
【计算机视觉】基于复杂环境下的车牌识别
人工智能·python·算法·计算机视觉·去噪·车牌识别·字符识别
Jonkin-Ma1 小时前
每日算法(1)之单链表
算法
晚风残2 小时前
【C++ Primer】第六章:函数
开发语言·c++·算法·c++ primer
杨云强2 小时前
离散积分,相同表达式数组和公式
算法
地平线开发者2 小时前
征程 6 | BPU trace 简介与实操
算法·自动驾驶
满天星83035772 小时前
【C++】AVL树的模拟实现
开发语言·c++·算法·stl
Lris-KK2 小时前
力扣Hot100--94.二叉树的中序遍历、144.二叉树的前序遍历、145.二叉树的后序遍历
python·算法·leetcode
Mr_WangAndy3 小时前
C++设计模式_行为型模式_责任链模式Chain of Responsibility
c++·设计模式·责任链模式·行为型模式
麦麦鸡腿堡3 小时前
Java的动态绑定机制(重要)
java·开发语言·算法
时间之里3 小时前
【c++】:Lambda 表达式介绍和使用
开发语言·c++