x*y的所有因数=x的所有因数分别乘以y的所有因数

例题C-小红的因子幂和_牛客周赛 Round 143(重现赛)

思路：

我们先分别求出x，y的所有因子，接着用一个集合来存储x*y的所有因数，也就是x的所有因数分别乘以y的所有因数，然后用ksm求值便可

代码如下：

#include <bits/stdc++.h>
using namespace std;
#define endl "\n"
#define int  long long
const int mod = 1e9 + 7;
int q_mul(int n, int x) {
    int res = 1;//累乘变量
    n %= mod;
    while (x != 0) {//只要x不为0就始终执行下去
        if (x % 2 == 1)res =res*n%mod;//当前的2进制位为1，乘上当前位置的n，即我们的任务2
        x /= 2;//x除以2，和短除法一样
        n =n*n%mod;//n每次乘以自身，即我们的任务1
    }
    return res;
}
vector<int> isp(int n) {
    vector<int>v;
    for (int i=1;i*i<=n;i++) {
        if (n%i==0) {//判断i是不是n的因子（可以整除就是）
            v.push_back(i);
            if (i!=n/i) {//与i成对的因子
                v.push_back(n/i);
            }
        }
    }
    return v;
}
void solve() {
    int x,y;
    cin >> x >> y;
    vector<int>a=isp(x);
    vector<int>b=isp(y);
    set<int>s;
    for (int i=0;i<a.size();i++) {
        for (int j=0;j<b.size();j++) {
            s.insert(a[i]*b[j]);
        }
    }
    int ans=0;
    for (auto it=s.begin();it!=s.end();it++) {
        ans=(ans+q_mul(*it,*it))%mod;
    }
    cout << ans << endl;
}

signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int T=1;
    //cin >> T;
    while (T--)
        solve();
    return 0;
}