CRB and Candies （组合的LCM）

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

Lazines_by

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

阅读量378

点赞数

分类专栏：排列组合

本文链接：https://blog.csdn.net/xiaolonggezte/article/details/77505531

版权

排列组合专栏收录该内容

12 篇文章 0 订阅

订阅专栏

组合数的lcm:

C(n,0) ,C(n,1),C(n,2)……C(n,n)的lcm为1,2,……，n,n+1的lcm除以n+1,其中1--(n + 1)的lcm可以考虑是每一个质数的最大幂，只有当单独一个数形成数的时候，幂最大。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x)
#define PI acos(-1.0)
#define ITER set<int>::iterator
const int Mod = 1e9 + 7;
const int maxn = 1000000 + 10;

bool is_[maxn];
int prime[maxn];
int len = 0;
void Init()
{
    clr(is_,true);
    for(int i = 2; i < maxn; i ++)
    {
        if(is_[i])
        {
            prime[len ++] = i;
            for(int j = i * 2; j < maxn; j += i)
                is_[j] = false;
        }
    }
}
ll pows(ll x,ll n)
{
    ll ret = 1;
    while(n)
    {
        if(n & 1)ret = ret * x % Mod;
        x = x * x % Mod;n >>= 1;
    }
    return ret;
}
int main()
{
    Init();int n;int Tcase;scanf("%d",&Tcase);
    while(Tcase --)
    {
        scanf("%d",&n);ll ans = 1;n ++;
        for(int i = 0; i < len && prime[i] <= n; i ++)
        {
            int t = 0;ll temp = 1;
            while(temp * prime[i] <= n){temp *= prime[i];t ++;}
            ans = ans * temp % Mod;
        }
        printf("%I64d\n",ans * pows(n,Mod - 2) % Mod);
    }
    return 0;
}