LCM性质 + 组合数 - HDU 5407 CRB and Candies

题目描述

给定一个数n，求LCM(C(n,0),C(n,1),C(n,2)...C(n,n))的值,(n<=1e6)。题目连接

解题思路

颇有趣的一道数论题！

看了下网上别人的作法，什么Kummer定理我还真没据说过，仔细研究一下那个鬼定理真是涨姿式了!

然而这题我并非用Kummer那货搞的(what?).

其实这题真的很简单（不要打我），为何这样说呢？看了下面的解释你就知道我没骗你。

首先咱们看一下这个式子：LCM(C(n,0),C(n,1),C(n,2)...C(n,n))

当时个人第一感受是：晕，仍是打个表吧！结果，打表程序后台打了四个半小时也没打完=.=（时间复杂度算错了）

作这题首先你得知道这个（基本常识）：

求多个数的最小公倍数，有两种方法：

1）分解质因数法

先把这几个数分解质因数，再把它们一切公有的质因数和其中几个数公有的质因数以及每一个数的独有的质因数所有连乘起来，所得的积就是它们的最小公倍数。

例如，求LCM[12，18，20，60]

由于12=（2）×[2]×[3]，18=（2）×[3]×3，20=（2）×[2]×{5}，60=（2）×[2]×[3]×{5}

其中四个数的公有的质因数为2（小括号中的数），

三个数的公有的质因数为2与3[中括号中的数]，

两个数的公有的质因数为5{大括号中的数}，

每一个数独有的质因数为3。

因此，［12，18，20，60］=2×2×3×3×5=180。

2）公式法

因为两个数的乘积等于这两个数的最大公约数与最小公倍数的积。

即（a，b）×[a，b]=a×b。

因此，求两个数的最小公倍数，就能够先求出它们的最大公约数，而后用上述公式求出它们的最小公倍数。

例如，求[18，20]

即得[18，20]=18×20÷（18，20）=18×20÷2=180。

求几个天然数的最小公倍数，能够先求出其中两个数的最小公倍数，

再求这个最小公倍数与第三个数的最小公倍数，依次求下去，直到最后一个为止。

最后所得的那个最小公倍数，就是所求的几个数的最小公倍数。

知道这个后，作这题选择哪一种方法呢？

若是选择第二种方法，恭喜你，你绝壁和我同样想到打表滚粗！

既然第二种方法不行，确定只能是第一种方法了。

那么要怎么作呢？

首先咱们来看，对于组合数C(n,m)，能够有以下变换：

C(n,m)=n!/[(n-m)!*m!]=n*(n-1)*(n-2)*....(m+1) / (n-m)! 

这一步应该没问题吧！

也就是：C(n,m)=n!/[(n-m)!*m!]=n*(n-1)*(n-2)*....(m+1) / (n-m)!  = n*(n-1)*(n-2)*....(m+1)/1/2/3/4/5/..../(n-m)

咱们把先后结合一下，边乘边除：

对于第k步，就至关于*(n+1-k)且/k，k={1,2,...n-m}.

咱们以n=8为例：

C(8,0)=1

C(8,1)=8*7*6*5*4*3*2 /7/6/5/4/3/2/1

C(8,2)=8*7*6*5*4*3 /6/5/4/3/2/1

C(8,3)=8*7*6*5*4 /5/4/3/2/1

C(8,4)=8*7*6*5 /4/3/2/1

C(8,5)=8*7*6 /3/2/1

C(8,6)=8*7 /2/1

C(8,7)=8 /1

C(8,8)=1

结合求n个数的LCM的方法，咱们将问题转换成：

找i个数共有的质数，而后相乘就可，i={1,2,..n}。

好了，你可能会说：*$#@*@，找i个数共有的质数难道不超时，并且你的代码里连一个0~n的for循环都没有，你在逗我？

不急，看下面：

首先咱们明确一点，C(n,k)的最大质因数是不会大于n的。

那么对于一个质数p来讲，他对"n个数的LCM"的贡献在哪？

是否是就是p^1,p^2,p^3...中的一些？

哪些呢？

前面求组合数中，咱们把C(n,m)分红了分子和分母来看。

若是p^x可以整除(n-1+k)，那么他有多是知足的，可是还不够，还要看是否是会被分母抵消掉。

只有p^x知足(n-1+k)%(p^x)==0且知足k%(p^x)!=0，这个p^x才是知足的，也就是对答案才有贡献，此时ans须要乘以p。

最后一步，约约分可能会更方便：把分子分母合一下，变成了:(n-1)%(p^x)!=0，表示(n-1+k)%(p^x)==0和k%(p^x)!=0不是同时出现的，此时才知足。

OK，推导完毕。

最终方法就是：

先筛出1e6之内的全部素数p，而后判断(n-1)%(p^x)是否!=0，是的话，ans*=p。

 

时间复杂度

O(p_num*sqrt(n))

代码

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-08-21-15.17
* Time: 0MS
* Memory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#define mod 1000000007
#define  LL long long
#define  ULL unsigned long long
using namespace std;
const int NN=1000010;
bool v[NN];
int p[NN],num;
void makePrime(){
   int i,j;
   num=-1;
   for(i=2; i<NN; ++i){
         if(!v[i]) p[++num]=i;
         for(j=0; j<=num && i*p[j]<NN; ++j){
               v[i*p[j]]=true;
               if(i%p[j]==0) break;
         }
   }
}
int main(){
   ios_base::sync_with_stdio(false);
   cin.tie(0);
   makePrime();
   int t;
   scanf("%d",&t);
   while(t--){
         int n;
         scanf("%d",&n);
         LL ans=1;
         for(int i=0; i<=num; ++i){
               for(LL t=p[i]; t<=n; t*=p[i]){
                     if((n+1)%t!=0)
                           ans=ans*p[i]%mod;
               }
         }
         printf("%lld\n",ans);
   }
   return 0;
}