Product it again

最新推荐文章于 2024-01-29 12:01:17 发布

weixin_30871905

最新推荐文章于 2024-01-29 12:01:17 发布

阅读量90

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/lawyer/p/6761475.html

本文介绍了一种求解特定范围内所有整数的质数幂乘积的算法，并提供了具体的C++实现代码。该算法利用了质数分布的特性，通过预处理得到所有质数，然后计算每个质数对该乘积的贡献。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：求解 $$\prod_{1 \leq i \leq n} \prod_{1 \leq j \leq m} {(i,j)}$$

解法：

满脑子的反演

考虑对于第一个质数 $p$ 的贡献为 $p^{[\frac{n}{p}][\frac{m}{p}] + [\frac{n}{p^2}][\frac{m}{p^2}] ... }$

这样1~n的质数大概有$O(\frac{n}{logn})$，对于每一个质数$O(logn)$，总效率大概为 $O(n)$

#include <iostream>
#include <cstdio>
#include <cstring>

#define LL long long
#define N 10000010
#define P 1000000007LL

using namespace std;

int n, m, tot, prime[N];
bool v[N];

LL qpow(LL x, LL n)
{
    LL ans = 1;
    for(; n; n >>= 1, x = x * x % P)
        if(n & 1)
            ans = ans * x % P;
    return ans;
}

int main()
{
    int T;
    for(int i = 2; i < N; i++)
        if(!v[i])
        {
            prime[++tot] = i;
            for(int j = i+i; j < N; j += i)
                v[j] = 1;
        }
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d %d", &n, &m);
        if(n > m) swap(n, m);
        LL ans = 1;
        for(int i = 1; i <= tot; i++)
            if(prime[i] <= n)
            {
                LL tmp = prime[i];
                LL cnt = 0;
                while(tmp <= n)
                {
                    cnt += (n/tmp) * (m/tmp);
                    tmp *= prime[i];
                }
                ans = ans * qpow(prime[i], cnt) % P;
            }
        printf("%lld\n", ans);
    }
}