CF451E-Devu and Flowers【组合计数,容斥】

正题

题目链接:https://www.luogu.org/problem/CF451E


题目大意

长度为nnn的序列aia_iai,有如下要求

  1. ai∈[0...xi]a_i\in[0...x_i]ai[0...xi]
  2. ∑i=1nai=s\sum_{i=1}^na_i=si=1nai=s

求序列个数。


解题思路

考虑容斥,这样我们就可以将条件转换为ai>xia_i>x_iai>xi,然后我们可以让n−xi−1n-x_i-1nxi1这样就将条件变为了ai≥0a_i\geq 0ai0,然后可以用组合数插板法计算,也就是答案为
∑SCs−∑xi∈S(xi+1)−1n−1∗(−1)∣S∣\sum_{S}C_{s-\sum_{x_i\in S}(x_i+1)-1}^{n-1}*(-1)^{|S|}SCsxiS(xi+1)1n1(1)S


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll XJQ=1e9+7;
ll n,s,inv[50],x[50],ans;
bool v[50];
ll power(ll x,ll b)
{ll ans=1;while(b){if(b&1) ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans;
}
ll C(ll n,ll m)
{if(n<0||m<0||n<m)return 0;if(n%XJQ==0||m==0)return 1;ll x=1;for(ll i=n-m+1;i<=n;i++)x=x*(i%XJQ)%XJQ;return x*inv[m]%XJQ;
}
void dfs(ll dep,ll _,bool flag)
{ll _n=s+n;for(ll i=1;i<=n;i++)if(v[i]) _n-=x[i]+1;if(flag)ans=((ans+_*C(_n-1,n-1))%XJQ+XJQ)%XJQ;if(dep>n) return;v[dep]=1;dfs(dep+1,-_,1);v[dep]=0;dfs(dep+1,_,0);
}
int main()
{scanf("%lld%lld",&n,&s);inv[1]=1;for(ll i=2;i<=20;i++)inv[i]=inv[i-1]*power(i,XJQ-2)%XJQ;for(ll i=1;i<=n;i++)scanf("%lld",&x[i]);ans=C(s+n-1,n-1);dfs(1,1,0);printf("%lld",ans);
}

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

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

相关文章

SpringCloud Ribbon(四)之全局配置与服务配置

一、Ribbon全局配置 RibbonClients注解配置导入全局Ribbon客户端默认设置 public class GlobalRibbonConfig {//负载均衡策略Beanpublic IRule rule() {return new MicroServiceLoadBalanceRule();}//负载均衡器Beanpublic <E extends Server> ILoadBalancer loadBalan…

.NET下使用socket.io随笔记录

一、问题背景目前公司在互联网产品上需要程序与前端部分要进行一个实时交互&#xff0c;在进行一定程度上的选型后&#xff0c;决定使用socket.io框架进行一个实践&#xff0c;算是公司的一个新的 尝试&#xff0c;也算是给自己增加增长见闻&#xff0c;由于我是做后端以及桌面…

bzoj2111,P2606-[ZJOI2010]排列计数【Lucas,组合计数,dp】

正题 题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id2111 https://www.luogu.org/problem/P2606 题目大意 长为nnn的序列PPP&#xff0c;然后要求Pi>P⌊i2⌋P_i>P_{\lfloor \frac{i}{2}\rfloor}Pi​>P⌊2i​⌋​。求排列个数。 解题思路 若用iii连向…

SpringCloud Ribbon(五)之服务实例列表ServerList

一、服务实例列表ServerList 服务实例列表&#xff08;ServerList&#xff09;为负载均衡器&#xff08;Loadbalancer&#xff09;提供服务的可用实例列表。 负载均衡器&#xff08;Loadbalancer&#xff09;通过服务实例列表&#xff08;ServerList&#xff09;从注册中心&a…

超级简便的容器化部署工具(使用 ASP.NET Core 演示)

Docker 改变了我们部署网站的方式&#xff0c;从原先的手动编译打包上传&#xff0c;到现在的构建镜像然后推送部署&#xff0c;让我们在配置环境上所花费的时间大大减少了。不仅如此&#xff0c;通过一系列相关的工具配合&#xff0c;可以很轻松的实现 CI、CD。本文即将介绍的…

学习手记(2019/7/05~2019/8/31)——快乐暑假

文章目录二分答案的作用堆和区间很糙ddp线段树合并网络流结论の1树上莫队对角线与GCD区间与扫描线与方案数欧拉欧拉*1斯坦纳树切比雪夫距离二分匹配结论の1min-max容斥计算几何の -1二分答案的作用 求最大值最小ororor最小值最大将求值问题转换为判断问题在判断问题之间相互转…

SpringCloud Ribbon(六)之服务实例过滤器ServerListFilter

一、服务实例过滤器ServerListFilter 服务实例过滤器&#xff08;ServerListFilter&#xff09;为负载均衡器&#xff08;Loadbalancer&#xff09;提供从服务实例列表&#xff08;ServerList&#xff09;获取的服务实例过滤出符合要求的服务实例。 负载均衡器&#xff08;Lo…

听说你开发.NET还在用VS,小哥哥给你推荐全平台的Rider

前言.NET平台的开发一直都只能使用Visual Studio来开发&#xff0c;自从dotnet core 发布后不久&#xff0c;jetbrains 发布了Rider预览版&#xff0c;到目前为止的正式版2017.3.1。博主都使用过&#xff0c;因为博主的主力开发语言是C#&#xff0c;所以一直以来被捆绑到Window…

jzoj6293-迷宫【ddp,线段树,矩阵乘法】

正题 题目大意 一个n∗mn*mn∗m的迷宫&#xff0c;不能往左走&#xff0c;有墙&#xff0c;每次修改一个点或询问两个点之间的最短距离。 解题思路 考虑到nnn的值很小&#xff0c;所以我们可以用矩阵转移&#xff0c;然后要求支持修改和查询所以我们考虑ddpddpddp。 我们可以…

Scala与Java差异(六)之类定义

一、类定义 &#xff08;1&#xff09;定义类&#xff0c;包含field以及方法 class ScalaClass {private var field "one"def aMethod() { print("field : " field ) } def getField field } &#xff08;2&#xff09;调用对象方法 val scalaCla…

jzoj6297-世界第一的猛汉王【切比雪夫距离,扫描线】

正题 题目大意 有若干个红点和蓝点&#xff0c;对于每一对红点和蓝点&#xff0c;若距离大于DDD则蓝点压制红点&#xff0c;否则红点压制蓝点。然后红点和蓝点之间也有不定的压制关系。 求有多少个三角要求AAA压制BBB&#xff0c;BBB压制CCC&#xff0c;CCC压制AAA且至少包含…

给Ocelot做一个Docker 镜像

写在前面在微服务架构中&#xff0c;ApiGateway起到了承前启后&#xff0c;不仅可以根据客户端进行分类&#xff0c;也可以根据功能业务进行分类&#xff0c;而且对于服务调用服务也起到了很好的接口作用。目前在各个云端中&#xff0c;基本上都提供了ApiGateway的功能&#xf…

Spark SQL(一)之简介

Spark SQL是用于结构化数据处理的Spark模块。与基本的Spark RDD API不同&#xff0c;Spark SQL提供的接口为Spark提供了有关数据结构和正在执行的计算的更多信息。在内部&#xff0c;Spark SQL使用这些额外的信息来执行额外的优化。 与Spark SQL交互的方法有多种&#xff0c;包…

hdu4336-Card Collector【min-max容斥,期望概率】

正题 题目链接:http://acm.hdu.edu.cn/showproblem.php?pid4336 题目大意 每次获得第iii张牌的概率为pip_ipi​(每次只能获得一张)&#xff0c;期望多少回合后获得所有牌。 解题思路 min−maxmin-maxmin−max容斥:max{S}∑T⊆S(−1)∣T∣−1min{T}max\{S\}\sum_{T\subseteq …

.NET Core UI框架Avalonia

.NET Core UI框架Avalonia&#xff0c;Avalonia是一个基于WPF XAML的跨平台UI框架&#xff0c;并支持多种操作系统&#xff1a;Windows&#xff08;.NET Framework&#xff0c;.NET Core&#xff09;&#xff0c;Linux&#xff08;GTK&#xff09;&#xff0c;MacOS&#xff0c…

Spark SQL(二)之DataSet操作

一、创建DataSet 使用SparkSession&#xff0c;应用程序可以从现有的RDD&#xff0c;Hive表的或Spark数据源创建DataFrame 。 &#xff08;1&#xff09;基于JSON的内容创建一个DataFrame //hdfs Dataset<Row> df spark.read().json("hdfs://master:9000/test.j…

jzoj6296-投票【期望dp,贪心】

正题 题目大意 nnn个人&#xff0c;第iii投票的概率是pip_ipi​&#xff0c;选择kkk个求最大的平票概率。 解题思路 我们显然要让kkk人中一半投票的概率大&#xff0c;一半投票的概率小。 所以我们可以先进行排序&#xff0c;这样我们发现答案一定是选取一段前缀和一段后缀。…

揽货最短路径解决方案算法 - C# 蚁群优化算法实现

需求为&#xff08;自己编的&#xff0c;非实际项目&#xff09;&#xff1a;某配送中心进行揽货&#xff0c;目标客户数为50个客户&#xff0c;配送中心目前的运力资源如下&#xff1a;现有车辆5台单台运力最大行驶距离200千米单台运力最大载重公斤1吨问&#xff1a;运力怎样走…

Spark SQL(三)之视图与执行SQL

一、视图与sql执行 SparkSession能够以编程方式运行SQL查询并返回结果Dataset<Row> Dataset<Row> df spark.read().json("hdfs://master:9000/test.json"); df.createOrReplaceTempView("man");Dataset<Row> sqlDF spark.sql("…

hdu4965-Fast Matrix Calculation【矩阵乘法】

正题 题目链接:http://acm.hdu.edu.cn/showproblem.php?pid4965 题目大意 给出矩阵A,BA,BA,B&#xff0c;求(AB)n(AB)^n(AB)n。然后对于每个元素%6\% 6%6后取和。 题目大意 我们发现如果直接让AB∗ABAB*ABAB∗AB这样的时间复杂度是n3n^3n3&#xff0c;显然不可过。但是我们…