cf451E. Devu and Flowers(产生不同多重集数量)

cf451E. Devu and Flowers

题意:

有n个箱子,第i个箱子里有ai朵花,同一个箱子里花的颜色一样,不同箱子里的花颜色不一样。现在在这些箱子里选出m朵花组成一束,求一共有多少种方案。要求任意两束花都不一样

题解:

设第i个箱子里花的颜色是Bi,则本题就等价于从多集合S={A1 * B1,A2 * B2 …An * Bn}中选出M个元素能够产生的不同多重集的数量。根据多重集组合数的结论有:
CN+M−1N−1−∑i=1NCN+M−Ai−1N−1+∑1<=i<j<=NCN+M−Ai−Aj−3N−1−....+(−1)NCN+M−∑i=1NAi−(N+1)N−1C_{N+M-1}^{N-1}-\sum_{i=1}^{N}C_{N+M-A_{i}-1}^{N-1}+\sum_{1<=i<j<=N}C_{N+M-A_{i}-A_{j}-3}^{N-1}-....+(-1)^{N}C_{N+M-\sum_{i=1}^{N}A_{i}-(N+1)}^{N-1}CN+M1N1i=1NCN+MAi1N1+1<=i<j<=NCN+MAiAj3N1....+(1)NCN+Mi=1NAi(N+1)N1
具体证明略
对于这种容斥,我们一般用x在二进制表示下,第i位为1则说明选中第i位,x的二进制下共有p个1,就相当于选了p位,就可以用x来枚举上述式子
(−1)p∑N+M−Ai1−Ai2....−Aip−(p+1)N−1(-1)^p\sum_{N+M-A_{i_{1}}-A_{i_{2}}....-A_{i_{p}}-(p+1)}^{N-1}(1)pN+MAi1Ai2....Aip(p+1)N1
这样就可以得到容斥原理计算多重集组合数的公式的每一项
在求组合数过程中我们可以用Lucas进行优化
详细看代码

代码:

// Problem: E. Devu and Flowers
// Contest: Codeforces - Codeforces Round #258 (Div. 2)
// URL: https://codeforces.com/contest/451/problem/E
// Memory Limit: 256 MB
// Time Limit: 4000 ms
// Data:2021-09-02 11:04:13
// By Jozky#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{
#ifdef LOCALstartTime= clock();freopen("in.txt", "r", stdin);
#endif
}
void Time_test()
{
#ifdef LOCALendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
// const int maxn=;
const int mod= 1e9 + 7;
ll a[30];
ll inv[30];
ll poww(ll a, ll b)
{ll ans= 1;while (b) {if (b & 1)ans= ans * a % mod;a= a * a % mod;b>>= 1;}return ans % mod;
}
void init()
{for (int i= 1; i <= 20; i++)inv[i]= poww(i, mod - 2);
}
ll n, m;
ll C(ll n, ll m)
{if (m < 0 || n < 0 || n < m)return 0;n%= mod;if (n == 0 || m == 0)return 1;ll a= 1, b= 1;for (int i= 0; i < m; i++) {a= (a * (n - i)) % mod;b= (b * inv[i + 1]) % mod;}return a * b % mod;
}
ll Lucas(ll n, ll m)
{if (n < m)return 0;if (m == 0)return 1;return Lucas(n / mod, m / mod) * C(n % mod, m % mod) % mod;
}
int main()
{//rd_test();init();cin >> n >> m;for (int i= 1; i <= n; i++)read(a[i]);ll ans= 0;for (int x= 0; x < (1 << n); x++) {if (x == 0) {ans= (ans + Lucas(n + m - 1, n - 1)) % mod;}else {ll t= n + m;int p= 0;for (int i= 0; i < n; i++) {if ((x >> i) & 1) {p++;t-= a[i + 1];}}t-= (p + 1);// cout << "t=" << t << endl;if (p & 1) {ans= (ans - Lucas(t, n - 1)) % mod;}else {ans= (ans + Lucas(t, n - 1)) % mod;}}// cout << "ans=" << ans << endl;}cout << (ans + mod) % mod << endl;//Time_test();
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316191.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.NetCoreLinuxDockerPortainer踩坑历险记

最近有一个云服务器和数据库的迁移任务&#xff0c;踩坑爬坑无数次&#xff0c;觉得必须要记录一下。大家瓜子花生准备好&#xff0c;听我慢慢讲故事#手动笑哭#。故事背景公司是做电商业务的&#xff0c;在天猫有几家旗舰店数据量也很大。阿里有一个称为聚石塔的平台&#xff0…

F - Parenthesis Checking

​​​​​​​1 . ( - > 1 , ) - > -1; 先想前缀和&#xff0c;求出l&#xff0c;r中的最小值&#xff0c;如果最小值小于l-1的前缀和&#xff0c;那么说明有某个点的 ) 大于 ( 数量 不满足条件&#xff1b; 用线段树优化区间修改和区间查询&#xff1b; 2 . …

cf208E. Blood Cousins

cf208E. Blood Cousins 题意&#xff1a; 给你一个森林&#xff0c;m次询问&#xff0c;每次询问(v,p),问v的p-cousin有多少&#xff1f;p-cousin指的是与v在同一层且他们到lca的距离都是p 题解&#xff1a; 对于每次询问(v,p),我们都可以通过其找到v的p距离的父亲节点fa&a…

深入源码理解.NET Core中Startup的注册及运行

开发.NET Core应用&#xff0c;直接映入眼帘的就是Startup类和Program类&#xff0c;它们是.NET Core应用程序的起点。通过使用Startup&#xff0c;可以配置化处理所有向应用程序所做的请求的管道&#xff0c;同时也可以减少.NET应用程序对单一服务器的依赖性&#xff0c;使我们…

并查集虚拟节点

​ HDU 3234 Exclusive-OR AC代码简洁版 AC代码不好看的 Almost Union-Find UVA - 11987 AC代码 #include <iostream> #include <cstdio> #include <algorithm> #include <cstring> #include <sstream> using namespace std; typedef long …

cf246E. Blood Cousins Return

cf246E. Blood Cousins Return 题意&#xff1a; 给你一个森林&#xff0c;每个点都有自己的种类&#xff0c;问以v为根节点的子树中&#xff0c;与v距离为k的节点有多少种 题解&#xff1a; 和cf208E. Blood Cousins这个题差不多&#xff0c;就是多了一个种类&#xff0c;…

Python、Java、TypeScript 和 Perl 作者谈语言设计

Python 作者 Guido van Rossum、Java 作者 James Gosling、Turbo Pascal 和 TypeScript 作者 Anders Hejlsberg&#xff0c;以及 Perl 作者 Larry Wall 本月早些时候齐聚一堂&#xff0c;讨论了(YouTube)编程语言设计的过去和未来。Guido van Rossum 说&#xff0c;设计一种语言…

线段树+树状数组

感觉这东西就是每棵线段树管的区间变了&#xff0c;以前学的时候线段树总是只管一个点或者管&#xff08;1-i&#xff09;这些点&#xff0c;但是这东西如果加上树状数组的思想&#xff0c;每棵线段树管&#xff08; i-(i&-i)1 ~ i &#xff09;这些区间&#xff0c;那么动…

cf1009F. Dominant Indices

cf1009F. Dominant Indices 题意&#xff1a; 1号节点为根&#xff0c;问对于以每个点为根的子树种&#xff0c;求某个深度节点最多的层数。如果多个相等&#xff0c;输出深度小的。 题解&#xff1a; 直接dsu就完事了&#xff0c;相当于是求子树的众数&#xff0c;但是注意…

Asp.Net Core Web应用程序—探索

前言作为一个Windows系统下的开发者&#xff0c;我对于Core的使用机会几乎为0&#xff0c;但是考虑到微软的战略规划&#xff0c;我觉得&#xff0c;Core还是有先了解起来的必要。因为&#xff0c;目前微软已经搞出了两个框架了&#xff0c;一个是Net标准(.NetFramework)&#…

智能优化算法应用:基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用&#xff1a;基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用&#xff1a;基于跳蛛算法3D无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.跳蛛算法4.实验参数设定5.算法结果6.参考文献7.MA…

Transformation HDU - 4578

思路&#xff1a;三种改变的顺序&#xff0c;先变 后乘 再 加。 #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using namespace std; typedef long long LL; const int N 5e5 4, mod 1e4 7; int n, m;struct…

cf375D. Tree and Queries

cf375D. Tree and Queries 题意&#xff1a; 给你一颗有根树&#xff0c;每个点都有一个颜色&#xff0c;有m次询问&#xff0c;问以u为根的子树中&#xff0c;相同颜色数量超过k的有多少种颜色&#xff1f; 题解&#xff1a; 这个题做法很多&#xff0c;有莫队分块&#x…

让你的ASP.NET Core应用程序更安全

对于ASP.NET Core应用程序&#xff0c;除了提供认证和授权机制来保证服务的安全性&#xff0c;还需要考虑下面的一些安全因素&#xff1a;CSRF强制HTTPS安全的HTTP HeadersCSRFASP.NET Core通过AntiForgeryToken来阻止CSRF攻击&#xff0c;一般来说&#xff0c;当用户做表单提交…

201409-5 拼图

七十分超时代码&#xff0c;具体思路可以仿照acwing蒙德里安的梦想 就是枚举第i层放 a &#xff0c;第i-1层放b的合法方案其中b和a都是同一块&#xff1b; #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using nam…

.net core 中间件管道底层剖析

.net core 管道&#xff08;Pipeline&#xff09;是什么&#xff1f;由上图可以看出&#xff0c;.net core 管道是请求抵达服务器到响应结果返回的中间的一系列的处理过程&#xff0c;如果我们简化一下成下图来看的话&#xff0c;.net core 的管道其实就是中间件的部分。微软中…

Codeforces Round #741 (Div. 2)

Codeforces Round #741 (Div. 2) 题号题目知识点AThe Miracle and the SleeperBScenes From a MemoryCRingsD1Two Hundred Twenty One (easy version)D2Two Hundred Twenty One (hard version)ERescue Niwen!FTubular Bells

201609-5 祭坛

超时暴力60分 #include <iostream> #include <algorithm> #include <cstring> #include <sstream> using namespace std; typedef long long LL; const int N 5e5 4, mod 1e9 7;vector<int> v[N]; int y[N], sizy, x[N], sizx, numx[N], a[…

cf1562E. Rescue Niwen!

cf1562E. Rescue Niwen! 题意&#xff1a; 我们定义一个字符串s1s2s3…sn的展开为:s1,s1s2,s1s2…sn,s2,s2s3,s2s3…sn,…,sn 找到字符串s“扩展”的最大递增子序列的大小(根据字典序大小比较) 题解&#xff1a; 第一感觉就是求最长上升子序列的变形 按照字典序大小比较规则…

架构为什么要以领域为核心

很久以前, 人们以为地球是太阳系的中心.然后一位聪明人, 哥白尼, 他改变了我们对太阳系的看法. 他认为太阳是太阳系的中心:这是对太阳更好的一个解释, 更简单也更具说服力. 事实上, 以太阳为中心的模型确实是更优雅的.上面这件事也发生在软件开发里. 下面这个就是很多开发者惯用…