[51Nod 1189阶乘分数]数学

最新推荐文章于 2022-05-03 20:20:10 发布

Xingw-Xiong

最新推荐文章于 2022-05-03 20:20:10 发布

阅读量807

点赞数 1

分类专栏： ACM____数学文章标签： math

本文链接：https://blog.csdn.net/ACMore_Xiong/article/details/54318002

版权

ACM____数学专栏收录该内容

41 篇文章 0 订阅

订阅专栏

[51Nod 1189阶乘分数]数学

知识点：math prime factor

1. 题目链接

[51Nod 1189阶乘分数]

2. 题意描述

$\frac{1}{N!} = \frac{1}{X} + \frac{1}{Y}（0<x<=y）$ ，给出N，求满足条件的整数解的数量。例如： $N = 2，1/2 = 1/3 + 1/6，1/2 = 1/4 + 1/4$ 。由于数量可能很大，输出 $\mod 10^9 + 7$ 。 $（1 <= N <= 1000000)$

3. 解题思路

贴一下Wannafly Union的题解。
通常这种题目，x和y都具有轮换性时，就要特别注意了
通常最后式子能化简成一个具有轮换性的式子，比如这题，化简后得到

(X - N!) * (Y - N!) = (N!) 2

$(X-N!)*(Y-N!)=(N!)^2$
步骤如下：

(N!) 2 = (X Y X + Y) 2 = X X + Y * Y X + Y = X X + Y * Y X + Y = (X - X Y X + Y) * (Y - X Y X + Y) = (X - N!) * (Y - N!)

$\begin{align} (N!)^2 & =(\frac{XY}{X+Y})^2\\ & =\frac{X}{X+Y}*\frac{Y}{X+Y}\\ & =\frac{X}{X+Y}*\frac{Y}{X+Y}\\ & =(X-\frac{XY}{X+Y})*(Y-\frac{XY}{X+Y})\\ & =(X-N!)*(Y-N!) \end{align}$
所以我们只需要对阶乘求约数个数，就知道方案数了！
到了这一步后，大家就很熟练了。只要对质因数分解，由于有一个平方，把每个质数的右上标全部乘以2。
最后得到

(N!)=pa11∗pa22∗⋯∗pamm $(N!)=p_1^{a_1}*p_2^{a_2}*\cdots*p_m^{a_m}$
其中

pi $p_i$ 是质数，表示分解出的质因子的个数。
如果x与y的大小没有限制，那么答案就是

∏(ai+1) $\prod{(a_i+1)}$
但是本题规定了

x≤y $x\le y$ , 那么答案就是

⌊1+∏(ai+1)2⌋ $\lfloor\frac{1+\prod{(a_i+1)}}{2}\rfloor$ （逆元）

4. 实现代码

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
typedef long double LB;
typedef pair<int, int> PII;
typedef pair<LL, LL> PLL;

const int INF = 0x3f3f3f3f;
const LL INFL = 0x3f3f3f3f3f3f3f3fLL;
const LB eps = 1e-8;
const int MAXN = 1000000 + 5;
const int MOD = 1e9 + 7;

int n;
int prime[MAXN + 1];
int cnt[MAXN];
int factor[100][2];
int fatCnt;
LL inv2, ans;

template<typename T> T umax(T& a, T b) { a = max(a, b); }

LL qpow(LL a, LL b) {
    LL ret = 1;
    while(b > 0) {
        if(b & 1) ret = ret * a % MOD;
        a = a * a % MOD;
        b >>= 1;
    }
    return ret;
}

void getPrime() {
    memset(prime, 0, sizeof(prime));
    for(int i = 2; i <= MAXN; i++) {
        if(!prime[i])prime[++prime[0]] = i;
        for(int j = 1; j <= prime[0] && prime[j] <= MAXN / i; j++) {
            prime[prime[j]*i] = 1;
            if(i % prime[j] == 0) break;
        }
    }
}

int getFactors(int x) {
    fatCnt = 0;
    int tmp = x;
    for(int i = 1; prime[i] <= tmp / prime[i]; i++) {
        factor[fatCnt][1] = 0;
        if(tmp % prime[i] == 0) {
            factor[fatCnt][0] = prime[i];
            while(tmp % prime[i] == 0) {
                factor[fatCnt][1]++;
                tmp /= prime[i];
            }
            fatCnt++;
        }
    }
    if(tmp != 1) {
        factor[fatCnt][0] = tmp;
        factor[fatCnt++][1] = 1;
    }
    return fatCnt;
}

int main() {
#ifdef ___LOCAL_WONZY___
    freopen("input.txt", "r", stdin);
#endif // ___LOCAL_WONZY___
    getPrime();
    inv2 = qpow(2, MOD - 2);
    while(~scanf("%d", &n)) {
        memset(cnt, 0, sizeof(cnt));
        int maxv = -1;
        for(int i = 1; i <= n; ++i) {
            getFactors(i);
            for(int j = 0; j < fatCnt; ++j) cnt[factor[j][0]] += factor[j][1], umax(maxv, factor[j][0]);
        }
        ans = 1;
        for(int i = 2; i <= maxv; ++i) {
            if(cnt[i] == 0) continue;
            ans = ans * ((cnt[i] << 1) + 1) % MOD;
        }
        ans = (ans + 1) * inv2 % MOD;
        printf("%lld\n", ans);
    }
    return 0;
}