中位数
来源:互联网 发布:数据库进销存管理系统 编辑:程序博客网 时间:2024/06/07 13:46
题目大意
定义一个序列a的特征值:
如果n=1则a[1]即为特征值
否则构造序列b,bi为ai、ai-1、ai+1的中位数。
a的特征值等于b的特征值。
求特征值为k的长度为n的排列a有多少种。
做
转化为>=k
那么只有0和1
从中间0101交替往两边对称,直到不对称时,两个数是多少决定了特征值是0还是1。
因此可以枚举对称长度,剩余组合数。
#include<cstdio>#include<algorithm>#define fo(i,a,b) for(i=a;i<=b;i++)#define fd(i,a,b) for(i=a;i>=b;i--)using namespace std;typedef long long ll;const int maxn=1000000+10,mo=998244353;int fac[maxn],inv[maxn];int i,j,k,l,t,n,m,ans;int qsm(int x,int y){ if (!y) return 1; int t=qsm(x,y/2); t=(ll)t*t%mo; if (y%2) t=(ll)t*x%mo; return t;}int C(int n,int m){ if (n<m||m<0) return 0; return (ll)fac[n]*inv[m]%mo*inv[n-m]%mo;}int calc(int k){ if (k==0) return 0; if (k==1) return n==1; int t=(k==n/2+1); int i; fo(i,1,n/2){ (t+=C(n-2*i-1,k-i-1)*2%mo)%=mo; (t+=C(n-2*i-1,k-i-2))%=mo; } return (ll)t*fac[k]%mo*fac[n-k]%mo;}int main(){ freopen("median.in","r",stdin);freopen("median.out","w",stdout); scanf("%d%d",&n,&k); fac[0]=1; fo(i,1,n) fac[i]=(ll)fac[i-1]*i%mo; inv[n]=qsm(fac[n],mo-2); fd(i,n-1,0) inv[i]=(ll)inv[i+1]*(i+1)%mo; ans=(calc(n-k+1)-calc(n-k))%mo; (ans+=mo)%=mo; printf("%d\n",ans);}
0 0
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 中位数
- 深层学习为何要“Deep”(下)
- 腾讯2017实习生招聘在线笔试编程题(技术类-PC客户端开发)
- UML建模、数据库设计和UI设计工具
- 数学之美番外篇:平凡而又神奇的贝叶斯方法
- 并查集
- 中位数
- POJ 2891 Strange Way to Express Integers (扩展欧几里得)
- java字符串之split
- Ubuntu16.04配置apache2+CGI
- 题目1124:Digital Roots 九度OJ
- 数据权限实现(Mybatis拦截器+JSqlParser)
- raise IOError, 'Not a gzipped file'
- 常用计算机远程控制软件介绍
- JS 调取摄像头