HDU 5407(2015多校10)-CRB and Candies(组合数最小公倍数+乘法逆元)
题目地址: pid=5407">HDU 5407
题意:CRB有n颗不同的糖果,如今他要吃掉k颗(0<=k<=n),问k取0~n的方案数的最小公倍数是多少。
思路:首先做这道题我们须要必备的几个技能点。
1. LCM(C(n,0), C(n,1),…, C(n,n))=LCM(1,2,3,…n+1)/(n+1)。额,这个有一篇证明Kummer定理
2.(1) 乘法逆元定义:
满足a*k≡1 (mod p)的k值就是a关于p的乘法逆元(a,p互质)。
(2)为什么要用乘法逆元:
当我们要求(a/b) mod p的值,且a非常大,无法直接求得a/b的值时,我们就要用到乘法逆元。我们能够通过求b关于p的乘法逆元k,将a乘上k再模p。即(a*k) mod p。其结果与(a/b) mod p等价。
(3)乘法逆元的解法:
A:逆元能够用扩展欧几里德来解最小的正整数就可以:
a*x%p = 1;
a*x = y*p + 1;
a*x -p*y = 1;
B:当p是质数的时候 a/x mod p=a*x^(p-2) mod p
证明:若 a*b mod p = 1 则a和b互为乘法逆元。
由于p是质数。所以由费马小定理得出x^(p-1) mod p = 1 ,所以x*x^(p-2) mod p = 1是成立的,所以x 和 x^(p-2) 互为乘法逆元。
当p不是质数的时候a/x mod p=a*x^(phi(p)-1) mod p。
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <sstream>
#include <algorithm>
#include <set>
#include <queue>
#include <stack>
#include <map>
#include <bitset>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
typedef long long LL;
const int inf=0x3f3f3f3f;
const double pi= acos(-1.0);
const double esp=1e-7;
const int Maxn=1e6+10;
const int mod=1000000007;
bitset<Maxn>pri;
int prime[Maxn];
int vis[Maxn];
int k;
LL x;
LL f[Maxn];
void is_prime()
{
pri.set();
for(int i=2;i<Maxn;i++){
if(pri[i]){
prime[k++]=i;
for(int j=i+i;j<Maxn;j+=i)
pri[j]=0;
}
}
}
LL Mul(LL a,LL b,LL mod)
{
LL res=0;
while(b>0) {
if(b&1)
res=(res+a)%mod;
b>>=1;
a=(a+a)%mod;
}
return res;
}
LL modxp(LL a,LL b,LL mod)
{
LL res=1;
while(b>0) {
if(b&1)
res=Mul(res,a,mod);
b>>=1;
a=Mul(a,a,mod);
}
return res;
}
void get()
{
memset(vis,0,sizeof(vis));
for(int i=0;i<k;i++){//将素数p的k次方p^k标记一下,找出符合的素数
x=prime[i];
while(x<Maxn){
vis[x]=prime[i];
x*=prime[i];
}
}
f[1]=1;
for(int i=2;i<Maxn;i++){
if(vis[i])
f[i]=(f[i-1]*vis[i])%mod;
else
f[i]=f[i-1]%mod;
}
}
int main()
{
int T,n;
is_prime();
get();
scanf("%d",&T);
while(T--){
scanf("%d",&n);
printf("%lld\n",f[n+1]*modxp(n+1,mod-2,mod)%mod);//(f(n+1)/(n+1))%mod
}
return 0;
}
HDU 5407(2015多校10)-CRB and Candies(组合数最小公倍数+乘法逆元)的更多相关文章
- HDU 5410(2015多校10)-CRB and His Birthday(全然背包)
题目地址:HDU 5410 题意:有M元钱,N种礼物,若第i种礼物买x件的话.会有Ai*x+Bi颗糖果,现给出每种礼物的单位价格.Ai值与Bi值.问最多能拿到多少颗糖果. 思路:全然背包问题. dp[ ...
- hdu 5288||2015多校联合第一场1001题
pid=5288">http://acm.hdu.edu.cn/showproblem.php?pid=5288 Problem Description OO has got a ar ...
- CRB and Candies(组合数学+求逆元+lcm)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5407 题目: Problem Description CRB has N different cand ...
- hdu5414(2015多校10)--CRB and String(字符串匹配)
题目链接:pid=5414">点击打开链接 题目大意:有A.B两个字符串.如今有一种操作能够在A的随意一个字符x后面添加一个字符y(x.=y).问能不能将A变为B. 首先假设A能够变成 ...
- HDU 5358(2015多校联合训练赛第六场1006) First One (区间合并+常数优化)
pid=5358">HDU 5358 题意: 求∑i=1n∑j=in(⌊log2S(i,j)⌋+1)∗(i+j). 思路: S(i,j) < 10^10 & ...
- HDU 5373(2015多校7)-The shortest problem(模拟%11)
题目地址:pid=5373">HDU 5373 题意:给你一个数n和操作次数t,每次操作将n的各位数之和求出来放在n的末尾形成新的n,问t次操作后得到的n能否够被11整除. 思路:就是 ...
- HDU 5326(2015多校3)-Work(dfs)
题目地址:pid=5326">HDU 5326 题意:给一张有向图n个点.n - 1(....输入n-1)条边. A指向B代表A管理B.然后能够间接管理,比方A管理B,B管理C.则A管 ...
- hdu 5361 2015多校联合训练赛#6 最短路
In Touch Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) Total ...
- HDU 5371(2015多校7)-Hotaru's problem(Manacher算法求回文串)
题目地址:HDU 5371 题意:给你一个具有n个元素的整数序列,问你是否存在这样一个子序列.该子序列分为三部分,第一部分与第三部分同样,第一部分与第二部分对称.假设存在求最长的符合这样的条件的序列. ...
随机推荐
- 【NOIP2007】字符串展开解题报告
描述 Description 在初赛普及组的“阅读程序写结果”的问题中,我们曾给出一个字符串展开的例子:如果在输入的字符串中,含有类似于“d-h”或“4-8”的子串,我们就把它当作一种简写,输出时,用 ...
- Python3 集合(无序的set)
ayout: post title: Python3 集合(无序的set) author: "luowentaoaa" catalog: true tags: mathjax: t ...
- (寒假开黑gym)2018 ACM-ICPC, Syrian Collegiate Programming Contest
layout: post title: (寒假开黑gym)2018 ACM-ICPC, Syrian Collegiate Programming Contest author: "luow ...
- dutacm.club 1094: 等差区间(RMQ区间最大、最小值,区间GCD)
1094: 等差区间 Time Limit:5000/3000 MS (Java/Others) Memory Limit:163840/131072 KB (Java/Others)Total ...
- 5、Flask实战第5天:页面跳转和重定向
GET和POST请求 在网络请求中有许多请求方式,比如GET, POST, DELETE, PUT等,最常用的就是GET和POST GET 只会在服务器上获取资源,不会更改服务器的状态,这种请求方式推 ...
- Codeforces #480 Tutorial
Problem A,B,C: 简单的模拟,注意A中p mod q时对q=0特殊处理(注意范围) Problem D: Brief Intro: 给定长度为N的数组A,将A中所有连续子序列分成最少的组, ...
- 【bitset】poj2443 Set Operation
模板题.S[i][j]表示i是否存在于第j个集合里.妈蛋poj差点打成poi(波兰无关)是不是没救了. #include<cstdio> #include<bitset> us ...
- Java序列化对象为字符串并将字符串反序列化为对象
对象的序列化与反序列化其实就是将对象的状态保存下来,一般是保存到文件中,但是其实更常用的是将对象序列化为字符串保存到数据库中,然后在需要读取对象的情况下将字符串反序列化为对象. 可以序列化的类必须 ...
- Problem B: 调用函数,输出100到200之间的所有素数
#include <stdio.h> int isPrime(unsigned int n)//定义素数函数 { int i; || n == ) ; ; i * i <= n; i ...
- UILabel设置行间距
UILabel设置行间距: NSMutableAttributedString *attributedString = [[NSMutableAttributedString alloc]initWi ...