Luogu P5368 [PKUSC2018]真实排名
时间:2019-12-21
本文章向大家介绍Luogu P5368 [PKUSC2018]真实排名,主要包括Luogu P5368 [PKUSC2018]真实排名使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
老年选手只会做SB题了(还调了好久)
很容易想到分类讨论,按第\(i\)个人有没有翻倍来算
若\(a_i\)未翻倍,显然此时将\([0,\lceil \frac{a_i}{2}\rceil)\)的数和\([a_i,\infty)\)的数翻倍都可以,记它们的个数为\(x\),则贡献为\(C_x^k\)
若\(a_i\)翻倍了,此时我们要算出\(i\)的排名变化了多少,记为\(dlt\)。然后在\([a_i,2a_i)\)之间的数翻倍之后都是会超过\(2a_i\)的,记为\(x\),因此这部分就是\(C_x^{dlt}\)。
然后此时还剩下\(k-1-dlt\)次操作,显然在\(a_i\)翻倍成\(2a_i\)之后,\([0,a_i)\)与\([2a_i,\infty]\)的数翻不翻倍都不影响答案,记个数为\(y\),贡献还要乘上\(C_y^{k-1-dlt}\)
具体实现的时候可以离散化之后对值等于某个数的所有位置一起讨论
#include<cstdio>
#include<algorithm>
#define RI int
#define CI const int&
using namespace std;
const int N=100005,mod=998244353;
int n,k,a[N],rst[N],pfx[N],suf[N],id[N],ans[N],num,fact[N],inv[N];
inline void inc(int& x,CI y)
{
if ((x+=y)>=mod) x-=mod;
}
inline int quick_pow(int x,int p=mod-2,int mul=1)
{
for (;p;p>>=1,x=1LL*x*x%mod) if (p&1) mul=1LL*mul*x%mod; return mul;
}
inline void init(CI n)
{
RI i; for (fact[0]=i=1;i<=n;++i) fact[i]=1LL*fact[i-1]*i%mod;
for (inv[n]=quick_pow(fact[n]),i=n-1;~i;--i) inv[i]=1LL*inv[i+1]*(i+1)%mod;
}
inline int C(CI n,CI m)
{
if (n<0||m<0) return 0; if (m==0) return 1; if (n<m) return 0;
return 1LL*fact[n]*inv[m]%mod*inv[n-m]%mod;
}
inline int GB(CI x) //>=x
{
return lower_bound(rst+1,rst+num+1,x)-rst;
}
inline int LB(CI x) //<=x
{
return upper_bound(rst+1,rst+num+1,x)-rst-1;
}
int main()
{
RI i; for (scanf("%d%d",&n,&k),i=1;i<=n;++i) scanf("%d",&a[i]),rst[i]=a[i];
rst[n+2]=1e9+1; sort(rst+1,rst+n+3); num=unique(rst+1,rst+n+3)-rst-1;
for (i=1;i<=n;++i) ++pfx[id[i]=GB(a[i])],++suf[id[i]];
for (i=num-1;i;--i) suf[i]+=suf[i+1]; for (i=2;i<=num;++i) pfx[i]+=pfx[i-1];
for (init(n),i=1;i<num;++i)
{
if (!rst[i]) { ans[i]=C(n,k); continue; }
int ls=pfx[LB(rst[i]-1>>1)]; ans[i]=C(suf[i]-1+ls,k);
int c=suf[GB(rst[i]<<1)],dlt=suf[i]-c-1,lt=pfx[LB(rst[i]-1)];
inc(ans[i],1LL*C(pfx[LB((rst[i]<<1)-1)]-lt-1,dlt)*C(c+lt,k-dlt-1)%mod);
}
for (i=1;i<=n;++i) printf("%d\n",ans[id[i]]); return 0;
}
原文地址:https://www.cnblogs.com/cjjsb/p/12077911.html
- 曾六位数被秒的yadea.com 终端真是雅迪!
- 2018年IT行业必须面对的7个改变
- SqlTransaction 数据库编程事务使用示例
- 再谈web.config/app.config敏感数据加/解密的二种方法
- NVIDIA颁布最新禁令,GeForce不能再行深度学习了!
- 2017年终盘点:请留意这七项科学技术的最新进展
- 什么叫做真正的无人驾驶汽车?难道真的没有人么?
- 零基础学习人工智能之Python篇1-Python定义
- aspx超时问题容易忽略的细节
- 给最后一周下个猛料,JavaScript 2017 使用调查!
- MongoDB快速入门
- AI:被过度夸大
- setInterval()与clearInterval()的用法
- Jquery开发插件的方法
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- flutter 轮播图动态加载网络图片的方法
- Kotlin如何捕获上下文中的变量与常量详解
- Android使用webView长按保存下载网络图片
- Android实现万能自定义阴影控件实例代码
- 浅析android studio3.5中使用recycleview的包
- Android自定义View实现拼图小游戏
- 解决android设备断电重启后WIFI不能自动重连的BUG(收藏)
- Android gradle配置抽取合并的操作步骤
- Android BottomNavigationBar底部导航的使用方法
- Android超清晰6.0权限申请AndPermission
- Android仿微信录制语音功能
- Android仿微信语音对讲录音功能
- flutter 屏幕尺寸适配和字体大小适配的实现
- flutter传递值到任意widget(当需要widget嵌套使用需要传递值的时候)
- android 9.0 launcher3 去掉抽屉式显示所有 app(代码详解)