当前位置: 首页 > 编程日记 > 正文

【数学 技巧】2.14计数

有趣的组合数学题;考试时候打满确实挺不容易的……

题目描述

对于一个 $n$ 阶排列 $p$,我们建立一张无向简单图 $G(p)$,有 $n$ 个节点,标号从 $1$ 到 $n$,每个点向左右两侧最近的比它大的点以及比它小的点连边。 形式化地,在 $G(p)$ 中,$\forall u<v$,边 $(u,v)$ 存在当且仅当以下四个条件至少一个成立:

  • $p_u<p_v$,且不存在 $u<i<v$ 满足 $p_u<p_i$;

  • $p_u>p_v$,且不存在 $u<i<v$ 满足 $p_u>p_i$;

  • $p_u<p_v$,且不存在 $u<i<v$ 满足 $p_i<p_v$;

  • $p_u>p_v$,且不存在 $u<i<v$ 满足 $p_i>p_v$。

现在在所有的 $n$ 阶排列中随机选择一个排列 $p$,请求出 $G(p)$ 中三元简单环的期望个数,答案对 $998244353$ 取模。

数据规模与约定

对于所有数据,$1\le n<998244353$。


题目分析

基础的暴力

首先会有个$(C^3_n)^2$的想法:枚举所有的三元环位置以及三个位置各是什么数,再根据这来统计每个三元环的合法贡献。

观察性质:从零开始推式子

提前透露一个模型:

$a+b+1=k, \, {n\choose k}=\sum_i{i-1\choose a}{n-i\choose b}$

感性理解来看相当于枚举一个中间点$i$并将它强制选取,再在它的两端钦定各选$a$,$b$个。这样就能有效地避免重复。

由于这题有非常特殊的建边性质,所以可以先来考虑怎么样形态的三元组才是合法的。

容易发现,三元环的极值一定是在中间点,那么就只存在四种形态的三元组。同时,这四种形态的三元环具有很好的对称性,所以它们方案数是相同的。

下面只考虑一种情况如上图所示。在这种情况下,想要三点间都存在边只能是能让$\max a_{i...k} < a3$.

接下去从最暴力的式子开始推起:

$\sum_{i=1}^{n-2}{n-i\choose2}\sum_{k_1+k_2<i}{i-1\choose k_1}{i-1-k_1\choose k_2}k_1!k_2!(n-k_1-k_2-3)!(n-k_1-k_2-2)$

其中i表示枚举a1的取值;k1,k2表示a1...a3,a3...a2之间夹的数的个数。最后的一项$(n-k_1-k_2-2)$表示将i...k这些数看作一个整体,插空在其他数里的方案数。

然后将后面的组合数拆开:

$\sum_{i=1}^{n-2}{n-i\choose2}\sum_{k_1+k_2<i}{(i-1)!(n-k_1-k_2-2)!\over (i-k_1-k_2-1)!}$ 

注意到这里变量$k_1+k_2$作为一个整体出现,所以可用$k$来代替:

$\sum_{i=1}^{n-2}{n-i\choose2}\sum_{k=0}^{i-1}{(i-1)!(n-k-2)!\over (i-k-1)!}(k+1)$

最后项$k+1$是因为有$k+1$种$k_1+k_2=k$的方案。这一点比较细节

接下去是换循环顺序和考虑组合数意义的环节:

$\sum_{k=0}^{n-3}\sum_{i=k+1}^{n-2}{(i-1)!\over(i-k-1)!}(n-k-2)!(k+1){n-i\choose2}$

将无关$i$的变量提出,同时在外部乘$k!$来构造内部的组合数:

$\sum_{k=0}^{n-3}(n-k-2)!(k+1)k!\sum_{i=k+1}^{n-2}{(i-1)!\over(i-k-1)!k!}{n-i\choose2}$

$\sum_{k=0}^{n-3}(n-k-2)!(k+1)!\sum_{i=k+1}^{n-2}{i-1\choose k}{n-i\choose2}$

考虑这两个组合数相乘的实际意义。第一个相当于钦定选第$i$个数,再在前$i-1$个数里选出$k$个数;第二个相当于在剩下的$n-i$个数里再选出2个数。那么这一种限制就是上面提到的模型,即将$i$从$k+1$枚举到$n$的总和,相当于是在$n$个数中取出$k+3$个数的方案数。

因此式子可化成这个样子:

$\sum_{k=0}^{n-3}(n-k-2)!(k+1)!{n\choose k+3}$

$n!\sum_{k=0}^{n-3}{(n-k-2)!(k+1)!\over(n-k-3)!(k+3)!}$

由于我们要求的是期望,那么就能把$n!$直接给去了。

$\sum_{k=0}^{n-3}{n-(k+2)\over(k+2)(k+3)}$

$\sum_{k=2}^{n-1}{n-k\over k(k+1)}$

$n\sum_{k=2}^{n-1}{1\over k(k+1)}-\sum_{k=2}^{n-1}{1\over k+1}$

对于前一项式子有:

$n\sum_{k=2}^{n-1}{1\over k(k+1)}=n\sum_{k=2}^{n-1}[{1\over k}-{1\over (k+1)}]=n({1\over 2}-{1\over n})$

而后一项式子是经典的求自然数倒数和,解决方法是分段打表。

节略后的代码如下:

 1 #include<bits/stdc++.h>
 2 typedef long long ll;
 3 const int MO = 998244353;
 4 const int seg[] = {0,227146302,963788543,596373246...};      //此处省略……65K
 5 
 6 int n,fac,size,pos;
 7 ll ans,cnt;
 8 
 9 int qmi(int a, int b)
10 {
11     int ret = 1;
12     for (; b; b>>=1, a=1ll*a*a%MO)
13         if (b&1) ret = 1ll*ret*a%MO;
14     return ret;
15 }
16 int main()
17 {
18     scanf("%d",&n);
19     ans = (2*n-4)%MO, size = 150000;
20     pos = (n-2)/size, cnt = seg[pos];
21     for (int i=1, mx=n-2-pos*size, sta=pos*size+2; i<=mx; i++)
22         cnt = 1ll*(cnt+qmi(sta+i, MO-2))%MO;
23     printf("%lld\n",(ans-(cnt<<2)%MO+MO)%MO);
24     return 0;
25 }

END

转载于:https://www.cnblogs.com/antiquality/p/10385163.html

相关文章:

冒泡排序 算法

算法思路&#xff1a; 从第一个元素开始遍历&#xff0c;比较当前元素和下一个元素的大小不符合&#xff0c;则交换结束最后一个元素&#xff0c;则重新遍历 实现&#xff1a; void bubble_sort(vector<int> &arr) {for (int i 0;i < arr.size() - 1; i) {//交…

5单个编译总会编译全部_5分钟读懂JavaScript预编译

大家都知道JavaScript是解释型语言&#xff0c;既然是解释型语言&#xff0c;就是编译一行&#xff0c;执行一行&#xff0c;那又何来预编译一说呢?脚本执行js引擎都做了什么呢&#xff1f;今天我们就来看看吧。1-JavaScript运行三部曲语法分析预编译解释执行语法分析很简单&a…

2014百度面试题目---“求比指定整数大且最小的不重复数”解答

题目&#xff1a;给定任意一个正整数&#xff0c;求比这个数大且最小的“不重复数”&#xff0c;“不重复数”的含义是相邻两位不相同&#xff0c;例如1101是重复数&#xff0c;而1201是不重复数。 代码&#xff1a; #include <iostream> using namespace std;bool istha…

3DMAX 批量 场景 对象 导出 .X格式 脚本

一、首先你需要下载一个 Total Commader文件管理软件。利用这个软件你可以收集文件夹下包含子文件夹下的max文件&#xff08;或完整路径&#xff09;打开TotalCMD后使用查找文件&#xff1a;&#xff08;如图红框中的操作&#xff09;1.2.3. 复制文件名和完整路径后粘贴到文本文…

C++复数类面向对象的参考

#include <bits/stdc.h> #include <future> #include <thread>using namespace std;class Complex { public:Complex (double r 0, double i 0): re (r), im (i) {} ///冒号后面是初始化的过程&#xff0c;注意分清初始化和赋值的区别Comp…

快速排序 算法

算法思路 序列终任意选择一个数&#xff0c;把序列分为比这个数大和比这个数小的两个子序列不断重复以上步骤(递归) 代码实现 int partition1(vector<int> &arr, int begin , int end) {int ret arr[begin];int index begin 1;for (int i index;i < end; i…

利用Spring AOP与JAVA注解为系统增加日志功能

Spring AOP一直是Spring的一个比较有特色的功能&#xff0c;利用它可以在现有的代码的任何地方&#xff0c;嵌入我们所想的逻辑功能&#xff0c;并且不需要改变我们现有的代码结构。 鉴于此&#xff0c;现在的系统已经完成了所有的功能的开发&#xff0c;我们需要把系统的操作日…

maven引入hadoop_如何添加Hadoop依赖通过Maven

匿名用户1级2017-09-09 回答Hadoop开发中需要用到至少不下10个的依赖包&#xff0c;它们相互间的依赖关系比较复杂&#xff0c;不同版本的依赖关系也有所不同&#xff0c;而间接依赖导致的程序错误并不会在运行之前报错&#xff0c;因此确定适合一个版本的依赖包&#xff0c;会…

js 闭包作用

2019独角兽企业重金招聘Python工程师标准>>> 一、变量的作用域 要理解闭包&#xff0c;首先必须理解Javascript特殊的变量作用域。 变量的作用域无非就是两种&#xff1a;全局变量和局部变量。 Javascript语言的特殊之处&#xff0c;就在于函数内部可以直接读取全…

vue实用组件——页面公共头部

可伸缩自适应的页面头部&#xff0c;屏幕适应范围更广泛 效果如下&#xff1a; 代码如下&#xff1a; <template> <div class"site-header"> <div class"logo"><img src"/assets/icons/logo.png" alt"">&…

插入排序 算法

算法思路 维护一段有序数列同时遍历待排序数列&#xff0c;在有序数列中找到合适的位置插入元素 基本代码 实现如下: void insertion(vector<int>& arr){for(int i1;i<arr.size();i){int tempi;for(int ji-1;j>0;j--){//有序序列不断得增加if(arr[temp]<…

线段树入门【转】

文章来自 &#xff1a; http://blog.csdn.net/x314542916/article/details/7837276 学习算法&#xff0c;自己收藏着。 线段树的入门级 总结 线段树是一种二叉搜索树&#xff0c;与区间树相似&#xff0c;它将一个区间划分成一些单元区间&#xff0c;每个单元区间对应线段树中的…

python自动化框架pytest pdf_pytest+python下的UI自动化基础框架

整体设计模式&#xff1a;config目录&#xff1a;存放一些公共的静态文件&#xff0c;如项目名称&#xff0c;配置文件等这些环境变量(可以用其他组件替换&#xff0c;如sql&#xff0c;主要能把配置文件的内容被程序识别)。httptrquest目录&#xff1a;存放接口代码&#xff0…

ny520 最大素因子 筛选法求素数

最大素因子时间限制&#xff1a;1000 ms | 内存限制&#xff1a;65535 KB难度&#xff1a;2 描述 GreyAnts最近正在学习数论中的素数&#xff0c;但是现在他遇到了一个难题&#xff1a;给定一个整数n,要求我们求出n的最大素因子的序数,例如&#xff1a;2的序数是1,3的序数是2…

JAVA_SE之内部类

内部类分类&#xff1a; 1. 成员内部类 静态内部类 非静态内部类 2. 局部内部类 3. 匿名内部类 1. 成员内部类&#xff1a; package com.atguigu.java; /** 类的第5个成员&#xff1a;内部类* 1.相当于说&#xff0c;我们可以在类的内部再定义类。外面的类&#xff1a;外部类。…

希尔排序 算法

算法思路 插入排序的改进版&#xff0c;选择插入距离远的元素选择一个间距&#xff0c;将序列分成很多子序列并行插入排序降低间距&#xff0c;并重复插入元素&#xff0c;直到间距将为1&#xff0c;完成排序。 算法实现 void shell_sort(vector<int> &arr, int b…

解决Apache CXF 不支持传递java.sql.Timestamp和java.util.HashMap类型问题

在项目中使用Apache开源的Services Framework CXF来发布WebService&#xff0c;CXF能够很简洁与Spring Framework 集成在一起&#xff0c;在发布WebService的过程中&#xff0c;发布的接口的入参有些类型支持不是很好&#xff0c;比如Timestamp和Map。这个时候我们就需要编写一…

python教学上机实验报告怎么写_Python基础(下)

不要忘了冒号啊&#xff01;&#xff01;&#xff01;&#xff01;&#xff01;对于基本数据类型的变量&#xff0c;变量传递给函数后&#xff0c;函数会在内存中复制一个新的变量&#xff0c;从而不影响原来的变量。(我们称此为值传递)但是对于表来说&#xff0c;表传递给函数…

比较有用的样式

背景图水平垂直居中 background:#ebebeb url(/Images/BlogHTImg/bkht_jia.jpg) center center no-repeat; 背景图居左垂直居中 background:#ebebeb url(/Images/BlogHTImg/bkht_jia.jpg) left center no-repeat; background:#ebebeb url(/Images/BlogHTImg/bkht_jia.jpg) 5px…

Python:线程之定位与销毁

背景 开工前我就觉得有什么不太对劲&#xff0c;感觉要背锅。这可不&#xff0c;上班第三天就捅锅了。 我们有个了不起的后台程序&#xff0c;可以动态加载模块&#xff0c;并以线程方式运行&#xff0c;通过这种形式实现插件的功能。而模块更新时候&#xff0c;后台程序自身不…

选择排序 算法

算法思路 维护一段有序数列&#xff0c;同时遍历待排序数列&#xff0c;找到最小的元素插入有序数列中重复&#xff0c;直到待排序数列没有剩余元素 代码实现 void select_sort(vector<int> &arr) {for (int i 0;i < arr.size(); i) {int temp arr[i];int in…

hdu2236 无题II 最大匹配 + 二分搜索

中文题目&#xff0c;题意大家都明白。 看到“不同的行和列”就觉得要用二分匹配来做。要求最大值与最小值的差值最小&#xff0c;是通过枚举边的下限和上限来完成。 枚举过程是这样的&#xff0c;在输入的过程可以记录下边权的最大值MAX和最小值MIN。那么他们的边权的差值的最…

python十大标准_python对标准类型的分类

python的标准类型可以按照三种方式分类。一、按存储模型分类按存储模型分可以分为原子(标量)类型和容器类型。原子(标量)类型指对象(这里的对象不是对象数据类型&#xff0c;而是任何可能的值)的值只能含有一种数据类型&#xff0c;比如数值和字符串。容器类型指它们的值可以含…

mysql慢查询开启及分析方法

最近服务维护的公司的DB服务器&#xff0c;总是会出现问题&#xff0c;感觉需要优化一下了&#xff0c;登陆上去&#xff0c;发现慢查询日志都没有开&#xff0c;真是惭愧&#xff0c; 故果断加上慢查询日志&#xff0c;经过分析sql记录&#xff0c;发现问题很多&#xff0c;开…

如何在调试页面的时候清除页面的缓存?

1.按F12,弹出下图 2.点击右上角的三个点: 3.点击settings 4.找到Network,下面的Disable cache(while DevTools is open) 转载于:https://www.cnblogs.com/studybrother/p/10396990.html

JAVA图片处理--缩放,切割,类型转换

import java.io.*; import java.awt.*; import java.awt.image.*; import java.awt.Graphics; import java.awt.color.ColorSpace; import javax.imageio.ImageIO;public class ChangeImageSize {/** *//*** 缩放图像* param srcImageFile 源图像文件地址* param result …

文本框自动提示_Excel办公小技巧,使用艺术字与文本框,就是那么的简单

Excel中的艺术字同时拥有文字和图形两种对象的属性&#xff0c;不仅可以修改其中的内容&#xff0c;还可以调整形状的大小、设置边框以及内部填充等效果&#xff0c;常在编辑表格标题或者输入一些比较有提示性的文本时使用&#xff0c;在突出关键内容的同时美化表格效果添加艺术…

Linux之父盟友分道扬镳 直言开源模式软肋

Linux之父盟友分道扬镳 直言开源模式软肋2005-09-06 12:53:00标签&#xff1a;linux职场开源休闲从1993年起&#xff0c;Larry McVoy就一直是Linux之父Linus Torvalds最忠实的盟友之一。 然而经历了这些年后&#xff0c;McVoy开始相信&#xff0c;开源这种风靡一时、纷纷被…

身份证第18位计算

本文计算方式源自 百度百科&#xff0c;根据计算方式&#xff0c;Java计算代码如下文所示。 计算方法 1、将前面的身份证号码17位数分别乘以不同的系数。从第一位到第十七位的系数分别为&#xff1a;7&#xff0d;9&#xff0d;10&#xff0d;5&#xff0d;8&#xff0d;4&…

归并排序 算法

算法思路 将一个数列不断拆分为子序列&#xff0c;直到只剩下0或者1个元素再将子序列按顺序合并为原来数列的大小&#xff0c;完成排序 代码实现 //合并两个有序数组 vector<int> merge_two_sort(vector<int> &arr1, vector<int> &arr2) {vector&…