HDU4390(容斥原理)

题目：Number Sequence

#include <iostream>
#include <string.h>
#include <stdio.h>using namespace std;
typedef long long LL;const int N=1000005;
const LL MOD=1000000007;int p[N];
int a[N],b[N];
bool prime[N];
LL C[1005][1005];
int k;void isprime()
{int i,j;k=0;memset(prime,true,sizeof(prime));for(i=2; i<N; i++){if(prime[i]){p[k++]=i;for(j=i+i; j<N; j+=i){prime[j]=false;}}}
}void Init()
{int i,j;for(i=0; i<1005; i++){C[i][0]=C[i][i]=1;for(j=1; j<i; j++){C[i][j]=(C[i-1][j]%MOD+C[i-1][j-1]%MOD)%MOD;}}
}void Solve(int b[],int n)
{memset(a,0,sizeof(a));for(int i=0; i<n; i++){for(int j=0; p[j]<=b[i]; j++){if(b[i]%p[j]==0){while(b[i]%p[j]==0){a[j]++;b[i]/=p[j];}}}}
}int main()
{isprime();Init();LL ans,n;while(~scanf("%I64d",&n)){ans=0;for(int i=0; i<n; i++)scanf("%d",&b[i]);Solve(b,n);for(int i=0; i<n; i++){LL tmp =C[n][i];for(int j=0; j<k; j++)tmp=(tmp*C[a[j]+n-i-1][a[j]])%MOD;ans+=(i&1)?-tmp:tmp;ans%=MOD;}printf("%d\n",(ans%MOD+MOD)%MOD);}return 0;
}

HDU4390(容斥原理)相关推荐

HDU 4135 Co-prime(容斥原理)
Co-prime 第一发容斥,感觉挺有意思的 →_→ [题目链接]Co-prime [题目类型]容斥 &题意: 求(a,b)区间内,与n互质的数的个数. $a,b\leq 10^{15}$ ...
BZOJ 2440: [中山市选2011]完全平方数 [容斥原理莫比乌斯函数]
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 3028 Solved: 1460 [Submit][S ...
bzoj 4710 [Jsoi2011]分特产组合数学+容斥原理
题面题目传送门解法考虑容斥原理显然,我们可以枚举有多少个人没有收到然后就转化成一个组合问题了假设现在有$x$个物品,$n$个人,可以有人没有被分到,那么分给这$n$个人的方案数 ...
【bzoj1853】[Scoi2010]幸运数字容斥原理+搜索
题目描述在中国,很多人都把6和8视为是幸运数字!lxhgww也这样认为,于是他定义自己的"幸运号码"是十进制表示中只包含数字6和8的那些号码,比如68,666,888都是&quo ...
hdu 4366 Card Collector (容斥原理)
http://acm.hdu.edu.cn/showproblem.php?pid=4336 题意:有 n 张卡片 ,每张卡片出现的概率是 pi 每包至多有一张卡片 ,也有可能没有卡片 . 求 ...
LightOJ 1095 Arrange the Numbers（容斥原理）
题目链接:http://www.lightoj.com/volume_showproblem.php?problem=1095 题意:给出数字n,m,K(0<K<=m<=n).在n个 ...
（每日一题）P1447 [NOI2010] 能量采集（莫反套路 + 欧拉反演 / 容斥原理）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划每日一题(莫反 / 多项式 / 母函数 / 群论) 2021.4.19 莫反 Problem Solu ...
解题报告（二）C、（darkBZOJ 3771）Triple（生成函数 + FFT + 容斥原理）（3）
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
BZOJ 4517: [Sdoi2016]排列计数 [容斥原理]
4517: [Sdoi2016]排列计数题意:多组询问,n的全排列中恰好m个不是错排的有多少个容斥原理强行推♂倒她 $恰好m个不是错排 $ \[ =\ \ge m个不是错排 - \ge m+1个不 ...

HDU4390(容斥原理)

HDU4390(容斥原理)相关推荐

最新文章

热门文章