hdu 5407 组合数的最小公倍数lcm

最新推荐文章于 2020-06-28 22:24:59 发布

zspqwe

最新推荐文章于 2020-06-28 22:24:59 发布

阅读量433

点赞数 1

分类专栏：数学题文章标签：组合数学数论

本文链接：https://blog.csdn.net/zsp123psz/article/details/47850489

版权

数学题专栏收录该内容

4 篇文章 0 订阅

订阅专栏

hdu 5407 本题求对于给定数n，求n的所有组合数的最小公倍数

如果a整除b， a整除c，那么a/c和a/b的最小公倍数等于a除以c和b的最大公约数，如下所示：

若c|a且b|a，则lcm(a/c, a/b) = a/gcd(c, b); （式1）

因此

lcm(n!/0!n!, n!/1!(n-1)!, n!/2!(n-1)!, n!/i!(n-i)!.....) = n ! / gcd(0!n!, 1!(n-1)!, 2!(n-1)!, i!(n-i)!.....) （式2）

已知（此处不放证明）

n!中个质因子p的次数为∑_(j=1)^∞[n/p^j ]

i!(n-i)!中质因子p的次数为∑_(j=1)^∞[i/p^j +(n-i)/p^j ]

当且仅当(n+1)%(p^j)==0 对于任意的i, n/p^j == i/p^j+ (n-i)/p^j

只要(n+1)%(p^j) !=0 必有一个i使得n/p^j > i/p^j+ (n-i)/p^j 且 n/p^j - i/p^j - (n-i)/p^j = 1；

因此，当出现(n+1)%(p^j) !=0时，等式2 的右边，分子会比分母多一个p

即只要有一个(n+1)%(p^j) !=0 ，答案就多乘一个p

这题技巧性太强，如果缺乏背景知识，是在太难想出来

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include <cmath>
#include <queue>
#define MAXN 1000010
#define ll long long
using namespace std;
const ll mod = 1000000007;
int prim[MAXN];
bool vis[MAXN];
int tot;
void work()
// 预处理获取质因子
{
    memset(vis, 0, sizeof(vis));
    tot = 0;
    prim[tot++] = 2;
    int t = 1000000;
    for(int i = 3; i <= t; i+= 2)
    {
        if(vis[i]) continue;
        for(ll j = (ll)i*i; j <= t; j += 2*i)
            vis[j] = true;
        prim[tot++] = i;
    }
}
int main()
{
    work();
    //freopen("in.txt", "r", stdin);
    int t;
    cin >> t;
    while (t--)
    {
        ll ans = 1;
        int n;
        scanf("%d", &n);
        int tag = 0;
        for (int i=0; i<tot; i++)
        {
            ll tmp = prim[i];
            while (tmp <= n)
            {
                if ((n+1)%tmp!=0)
                    ans = ans*prim[i]%mod;
                tmp *= prim[i];
            }
        }
        printf("%I64d\n", ans);
    }
    return 0;
}