清华大学生命科学博士就业_已拥有的是全部的生命

清华大学生命科学博士就业_已拥有的是全部的生命不错的组合数学题。同时这也驱使我去看积灰好久的《具体数学》(看了yu大的blog后)。然后看得头秃……得到一个不等式前缀和大于等于取了的个数。所以如果把每个卡的值减一,问题就变成了求一个排列,使得前

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

不错的组合数学题。同时这也驱使我去看积灰好久的《具体数学》(看了yu大的blog后)。然后看得头秃……

得到一个不等式前缀和大于等于取了的个数。所以如果把每个卡的值减一,问题就变成了求一个排列,使得前缀和都非负。

\[\forall i, \mathbf{s.t.} a_1 + a_2 + \dots + a_i \geq 0 \]

显然 \(\sum a_i = 0\)

看样子很难搞。

可以发现,如果把序列循环位移,可能会得到新的方案,那么对每个环进行计算。

但是环可能会贡献多个值,如果这么做,答案就是 \(\sum_{C} f\left(C\right)\),其中 \(C\) 是圆排列, \(f\left(C\right)\) 是断环的方式数,使得前缀和大于等于 \(0\),也等价于找到一个合法的断环方式,此时 \(f\left(C\right)\) 就是使得部分和为 \(0\) 的位置数。

这不就是定义吗?显然一个环的多个贡献给我们的计数造成了麻烦。

考虑转换一下,使环的贡献只能是 \(1\)。

如果见过 \(\texttt{Raney}\) 引理,就会发现一个类似的地方:

下面摘自《具体数学》301页 原话:

如果 \(\left< x_1, x_2, \dots , x_m \right>\) 是任何一个其和为 \(+1\) 的整数序列,那么它的循环移位

\[\left< x_1, x_2, \dots , x_m \right>,\left< x_2, \dots ,x_m,x_1 \right>,\dots,\left< x_m, x_1,\dots ,x_{m-1} \right> \]

中恰好有一个满足所有的部分和都是正数。

可以感性理解一下

然后提供了一个例子,就是往卡特兰数计数时,序列最前面加一个 \(1\)。然后容易计数。

这道题的部分和非负,如果在序列前面加个 \(1\),不就是 \(\texttt{Raney}\) 引理了?

我们考虑像这个引理一样,在序列头加一个 \(1\),可以发现,并不是所有地方都能加,例如 \(2,-1,-1\),如果在 \(2\) 后面加 \(1\) 显然是不行的。然后发现加的方案数正好和上面的 \(f\left( C\right)\) 一样——这不是加了跟没加一样?

但是不用担心,我们考虑黑科技,不加 \(1\),改成加 \(-1\)。

我们往原序列里插入一个 \({-1}^*\) ,此时问题转换为对圆排列找出一种断环方式,使得前 \(m\) 个的部分和非负,并且元素和为 \(-1\)。

易证最后一个必定是 \(-1\),并且用类似方法可以证明一定存在且仅存在一个满足条件的断环方式。(我不会证,逃

由于有标号,有 \(m – n + 1\) 个 \(-1\),于是最后一个 \(-1\) 就有 \(m – n + 1\) 种标号。我们要钦定最后一个 \(-1\) 为 \({-1}^*\),所以答案要除以 \(m – n + 1\)。

\(m + 1\) 个元素的圆排列方案数为 \(m!\),于是答案就是 \(\frac{m!}{m – n + 1}\)。

复杂度 \(O\left(m\right)\),可以通过此题。

使用快速阶乘可以得到更低的复杂度以及更大的常数。

#include <bits/stdc++.h>

const int mod = 998244353;
typedef long long LL;
int n, m;
int main() {
	std::cin >> n;
	for (int i = 1, t; i <= n; ++i) std::cin >> t, m += t;
	int ans = 1;
	for (int i = 2; i <= m; ++i)
		if (i != m - n + 1)
			ans = (LL) ans * i % mod;
	std::cout << ans << std::endl;
	return 0;
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/167479.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • oracle错误904解决方法_遇到oracle错误 12154

    oracle错误904解决方法_遇到oracle错误 12154案例情景–在一次Oracle数据库导出时:C:\DocumentsandSettings\Administrator>exp[emailprotected]_dbfile=E:\lsxy.dmpowner=lsxyExport:Release11.2.0.1.0-Productionon星期一11月2614:07:182012Copyright(c)…

    2022年9月20日
  • FFmpeg从入门到精通(一)

    FFmpeg从入门到精通(一)FFmpeg介绍FFmpeg既是一种音视频编解码工具,也是一组音视频解码开发套件。作为编解码开发套件,它为开发者提供了丰富的音视频处理的调用接口。FFmpeg基本组成先看图,再看字。FFmpeg的封装模块AVFormatAVFormat中实现了目前多媒体领域中大多数媒体封装格式,包括封装和解封装,如:MP4,FLV等文件封装格式、HTTP,RTMP,HLS等网络协议封装格式。FFm…

    2022年6月26日
  • 提升网页打开速度的实用方法

    提升网页打开速度的实用方法

    2021年10月16日
  • 2021韩顺平图解linux_狗剩学习笔记

    2021韩顺平图解linux_狗剩学习笔记韩顺平图解Linux全面升级https://www.bilibili.com/video/BV1Sv411r7vdP001_韩顺平图解Linux全面升级_课程内容20:37P002_韩顺平图解Linux全面升级_应用领域05:05P003_韩顺平图解Linux全面升级_概述16:37P004_韩顺平图解Linux全面升级_Linux与Unix18:09P005_韩顺平图解Linux全面升级_vmware15.5安装17:36P006_韩顺平图解

    2022年5月18日
  • Kafka集群扩容「建议收藏」

    Kafka集群扩容「建议收藏」Kafka集群的服务发现是有zk实现的,因此想kafka集群添加新的broker就非常容易。我们只需要为新的broker设置一个唯一的broker.id,然后启动新增的kafka就行。Kafka集群会自动发现新增的broker并同步原数据,包括当前集群有哪些topics以及topics的分区信息等。但是我们需要的注意的一点,新增的broker不能分担之前已经创建的topic的负载,必须手动执行…

    2022年5月16日
  • jrtplib java,jrtplib 分包处理

    转载自:http://blog.csdn.net/sxcong/article/details/3736354听说jrtplib写的不错,终于找到时间下来看看。下载,直接用VC6编译,很容易。然后打开VC,建立工程,测试examples下那几个收发程序,的确用起来很简单。想想以前都是自己封装UDP,现在的程序员真幸福。不过,在发送视频数据时出了问题,跟踪进去看了一下,里面设置最大帧数据长度为140…

    2022年4月17日

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号