P5520 【[yLOI2019] 青原樱】题解

整理博客的时候改了下分类标签,重新审一下

题目传送门

翻了翻题解区,发现基本没和我写的一样的(主要是都比我的写的简单

看题目:

第一眼,数学题;第二眼:组合数

接着想起来那道放苹果

n个位置,m棵树,就有n-m个空位,记space=n-m

转化问题为:

  1. space个空位插入m-1个位置(即左右两边都不留空)
  2. sapce个空位插入m个位置(即左边或右边留空,这种情况的答案要乘2)
  3. space个空位插入m+1个位置(即左右两边都留空位)

现在将space个空位看作space个苹果,m或m-1或m+1个位置看成盘子,因为每个位置(盘子)都要有至少一个空位(苹果),所以,这和放苹果就没什么区别了

运用一点隔板法解决:

以情况1为例: space个苹果间有space-1个间隔,因为要放进m-1个盘子,所以只需在space-1个间隔中选m-2个插入隔板,不重复不考虑顺序,所以:\(\tbinom {space-1} {m-2}\)

那么另外两种也就简单了,分别为:\(\tbinom {space-1} {m-1}\),\(\tbinom{space-1} {m}\)

考虑如何算组合数,发现p不一定为质数,所以对组合数计算的每个数分解,再约分

我用he数组表示一个数的最小质因子,为0说明不为合数

用sum数组记录每个质数的指数,是分子就加一,分母就减一

之后再快速幂乘起来

因为那m棵树是不同的,所以有m!种排列方式,也就是ans\(\times {m!}\)

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<iomanip>
#include<cstring>
#define R register
#define EN printf("\n")
#define LL long long
inline LL read(){
LL x=0,y=1;
char c=getchar();
while(c<'0'||c>'9'){if(c=='-') y=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+(c^48);c=getchar();}
return x*y;
}
int he[2000006],prime[2000006],cnt;//he[i]为i的最小质因子
LL sum[2000006];
inline void getprime(LL n){
for(R int i=2;i<=n;i++){
if(!he[i]) prime[++cnt]=i;
for(R int j=1;j<=cnt&&i*prime[j]<=n;j++){
he[prime[j]*i]=prime[j];
if(!(i%prime[j])) break;
}
}
}
inline LL pow(LL a,LL b,LL p){
LL ret=1;
while(b){
if(b&1) ret=(ret*a)%p;
a=(a*a)%p;
b>>=1;
}
return ret;
}
inline void fenjie(int x,int v){
while(he[x]){
sum[he[x]]+=v;
x/=he[x];
}
if(x>1) sum[x]+=v;
}
inline LL C(LL n,LL k,LL p){
memset(sum,0,sizeof sum);
if(k>n) return 0;
if(k==n) return 1%p;
if(k==0) return 1%p;
if(k==1) return n%p;
LL ret=1;
for(R int i=n-k+1;i<=n;i++) fenjie(i,1);
for(R int i=1;i<=k;i++) fenjie(i,-1);
for(R int i=2;i<=n;i++){
if(sum[i]) ret=(ret*pow(i,sum[i],p))%p;
}
return ret;
}
int main(){
LL ty=read(),n=read(),m=read(),p=read();
LL space=n-m;
getprime(n);
LL jc=1;
for(R int i=2;i<=m;i++) jc=(jc*i)%p;
LL ans=C(space-1,m-2,p);
ans=(ans+C(space-1,m-1,p)*2)%p;
ans=(ans+C(space-1,m,p))%p;
ans=(ans*jc)%p;
printf("%lld",ans);
return 0;
}

P5520 【[yLOI2019] 青原樱】的更多相关文章

  1. 【组合数学】【P5520】[yLOI2019] 青原樱

    A [yLOI2019] 青原樱 Background 星川之下皆萤火尘埃 我独行在人潮你天真而待 相遇若是借丹青着色 青原上 绯樱如海 --银临<青原樱>(Cover 人衣大人) Des ...

  2. 洛谷P5520 【[yLOI2019] 青原樱】

    这题是小学奥数啊. 题意:求\(m\)个不同物品两两不相邻的方案数. 直接排列组合. 我们可以减掉他们之间最少需要空出来的位数--\(m-1\)个空位 像这样,我们只用留\(m-1\)个空位放在每两个 ...

  3. asp.net MVC4——省市三级联动数据库

    数据库设计

  4. 通用js地址选择器

    用js实现通用的地址选择器,省份,城市,地区自动关联更新 点击下面查看详细代码: http://runjs.cn/code/s8sqkhcv 关键地址库代码: var addr_arr = new A ...

  5. java 随机生成身份证代码

    import java.util.Calendar; import java.util.Collection; import java.util.HashMap; import java.util.I ...

  6. 全国城市三级联动 html+js

    全国城市三级联动,没有css,所以屏幕的自适应必须自己想办法,手机端慎用(最好不要用,因为有些我也说不出的展示问题). html页面 <!DOCTYPE html> <html> ...

  7. 省市级联.net

    初学javascript,编译省市级联,使用json在一般处理程序中编译,利用ajax传递数据到web前台 <html xmlns="http://www.w3.org/1999/xh ...

  8. Android 三级联动选择城市+后台服务加载数据库

    技术渣,大家将就着看 首先我们需要一个xml数据保存到数据库,这里我从QQ下面找到一个loclist.xml文件 <CountryRegion Name="中国" Code= ...

  9. react 写的省市三级联动

    <!DOCTYPE html><html><head> <meta charset="utf-8"> <title>Ba ...

随机推荐

  1. javascript入门 之 zTree(十一 托拽事件(一))

    <!DOCTYPE html> <HTML> <HEAD> <TITLE> ZTREE DEMO - drag & drop</TITLE ...

  2. alg-查找只出现一次的数

    //只有2个数出现1次,其余的数都出现2次 class Solution { public: vector<int> singleNumber(const vector<int> ...

  3. flex实现三列布局

    css3新引入的flex在某些情况下布局非常实用 因为它是弹性盒子所以自适应效果会很棒 不过各项布局方案还是各有优劣 <!DOCTYPE html> <html lang=" ...

  4. 这本最适合夯实基础的经典 Java 书籍,可能80% 的 Java 程序员没有认真看过!

    公众号[程序员书单]出品,转载请注明出处 作者:黄小斜 今天要给大家带来的一本书,是大名鼎鼎的head first系列丛书的一本<head first Java>相信很多学习Java的朋友 ...

  5. stand up meeting 1/7/2016

    part 组员                今日工作              工作耗时/h 明日计划 工作耗时/h    UI 冯晓云 调研下滑条的存在问题,尝试替换方案     6 全面实行替换 ...

  6. 再接再厉,JSONViewer现已支持Firefox、Microsoft Edge、360浏览器,可能是最好用的JSON格式化工具

    之前写的JSONViewer,截至目前在谷歌商店里已经有1000+的自然下载量了 为什么开发JSONViewer? 日常开发中,拿到接口输出的JSON一般会去在线的JSON格式化网站查看,但是在线格式 ...

  7. DOS事件

    1 onblur 失去焦点 2 onchange   改变 3 onkeyup      按键弹起 4 onmouse over 鼠标移上去 5 onmouse leave   鼠标离开 6 onmo ...

  8. 【翻译】借助 NeoCPU 在 CPU 上进行 CNN 模型推理优化

    本文翻译自 Yizhi Liu, Yao Wang, Ruofei Yu.. 的  "Optimizing CNN Model Inference on CPUs" 原文链接: h ...

  9. Springboot:员工管理之首页(十(2))

    访问首页可以通过两种方式: 1:编写controller 2:自定义扩展视图解析器(推荐使用) 1:编写Controller com\springboot\controller\IndexContro ...

  10. 移植seetafaceengine-master、opencv到ARM板

    0.前言 在要移植opecv和SeetaFaceEngine-master到ARM板子上运行的所有步骤之前,有几点需要注意的: 查看板子运行的Kernel版本 交叉编译工具链的gcc版本,关键就是工具 ...