[HAOI 2018]染色
传送门
Description
一个长度为\(N\)的序列, 每个位置都可以被染成 \(M\)种颜色中的某一种.
出现次数恰好为 \(S\)的颜色种数有\(i\)种, 会产生\(w_i\)的愉悦度.
对于所有染色方案, 能获得的愉悦度的和对\(1004535809\)取模的结果.
Solution
\[ ans=\sum_{i=0}^{lim} w_i\cdot num_i \]
how to get \(num_i\)?
\(f_i\) : the number of occurrences of at least i colors is exactly the number of S
so \(f_i=\binom{m}{i}\cdot \frac{n!}{(s!)^i(n-iS)!}\cdot(m-i)^{n-iS}\)
According to the binomial inversion, we can know that:
\[ num_i=\sum_{j=i}^{lim}(-1)^{j-i}\binom{j}{i}f[j] \]
so
\[ num_i=\frac{1}{i!}\sum_{j=i}^{lim} \frac{(-1)^{j-i}}{(j-i)!}\cdot(f[j]\cdot j!) \]
we can use NTT.
Code
#include<bits/stdc++.h>
#define reg register
#define ll long long
#define db double
using namespace std;
inline int read()
{int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return x*f;
}
const int P=1004535809,G[2]={3,334845270},NN=5e5+5;
int Mul(int x,int y){return 1ll*x*y%P;}
int Add(int x,int y){return (x+y)%P;}
const int MN=1e7+5,MM=1e5+5,MS=155;
int N,M,S,W[MM],f[MM],fac[MN],inv[MN];
int fpow(int x,int y){int r=1;for(;y;y>>=1,x=Mul(x,x))if(y&1)r=Mul(r,x);return r;}
int C(int x,int y){if(x<0||y<0||x<y)return 0;return Mul(fac[x],Mul(inv[y],inv[x-y]));}
int a[NN],b[NN],pos[NN];
void NTT(int *a,bool ty,int L)
{reg int i,j,k,w,wn,x,y;for(i=0;i<L;++i) if(pos[i]<i) swap(a[i],a[pos[i]]);for(i=1;i<L;i<<=1){wn=fpow(G[ty],(P-1)/(i<<1));for(j=0;j<L;j+=(i<<1))for(w=1,k=0;k<i;++k,w=Mul(w,wn)){x=a[j+k],y=Mul(a[j+i+k],w);a[j+k]=Add(x,y);a[j+i+k]=Add(x,P-y);}}if(ty)for(j=fpow(L,P-2),i=0;i<L;++i)a[i]=Mul(a[i],j);
}
int ans;
int main()
{N=read(),M=read(),S=read();reg int i,lim;for(i=0;i<=M;++i) W[i]=read();lim=max(N,M);for(fac[0]=i=1;i<=lim;++i) fac[i]=Mul(fac[i-1],i);for(inv[0]=inv[1]=1,i=2;i<=lim;++i) inv[i]=Mul(inv[P%i],(P-P/i));for(i=1;i<=lim;++i) inv[i]=Mul(inv[i],inv[i-1]);lim=min(M,N/S);for(i=0;i<=lim;++i)f[i]=Mul(Mul(C(M,i),Mul(fac[N],fpow(inv[S],i))),Mul(inv[N-i*S],fpow(M-i,N-i*S)));for(i=0;i<=lim;++i) b[lim-i+1]=Mul(f[i],fac[i]);for(i=0;i<=lim;++i) a[i]=(i&1)?(P-inv[i]):inv[i];int MA;for(MA=1;MA<=(lim<<1);MA<<=1);for(i=0;i<MA;++i) pos[i]=(pos[i>>1]>>1)|((i&1)*(MA>>1));NTT(a,0,MA);NTT(b,0,MA);for(i=0;i<MA;++i) a[i]=Mul(a[i],b[i]);NTT(a,1,MA);for(i=0;i<=lim;++i) ans=Add(ans,Mul(W[i],Mul(a[lim-i+1],inv[i])));return 0*printf("%d\n",ans);
}
Blog来自PaperCloud,未经允许,请勿转载,TKS!
转载于:https://www.cnblogs.com/PaperCloud/p/10925578.html
[HAOI 2018]染色相关推荐
- 2018 CodeM复赛:C. 边的染色
链接:https://www.nowcoder.com/acm/contest/152/C 来源:牛客网 题目描述 小团有一张n个点,m条边的无向图G,有些边上已经被标记了0或1,表示它的边权. ...
- [2018 徐州 网络赛|Hard to prepare ] 环形染色问题的公式解法
题目来源 After Incident, a feast is usually held in Hakurei Shrine. This time Reimu asked Kokoro to deli ...
- 2018 ACM 国际大学生程序设计竞赛上海大都会赛重现赛 F Color it
链接:https://www.nowcoder.com/acm/contest/163/F 来源:牛客网 2018 ACM 国际大学生程序设计竞赛上海大都会赛重现赛 F Color it 时间限制:C ...
- 2018冬令营模拟测试赛(九)
2018冬令营模拟测试赛(九) [Problem A]王子 试题描述 不是所有王子都会遇见自己的中关村,主公,公主. 从前有个王子姓王,王王子遇到了一位美丽的公主,她的名字当然是公公主啦. 王王子对公 ...
- UOJ#449. 【集训队作业2018】喂鸽子
#449. [集训队作业2018]喂鸽子 DP好题 法一:min-max容斥 处理前m个,最快吃饱的鸽子期望的时间 根据期望的定义 考虑每个方案数的概率*期望次数 枚举前m个用了x个,概率都是(1/m ...
- java 区间树_线段树(区间树)之区间染色和4n推导过程
前言 线段树(区间树)是什么呢?有了二叉树.二分搜索树,线段树又是干什么的呢?最经典的线段树问题:区间染色:正如它的名字而言,主要解决区间的问题 一.线段树说明 1.什么是线段树? 线段树首先是二叉树 ...
- 训练日志 2018.10.24
图的连通性看了近一周,自闭了... 求强连通分量的 Kosaraju 算法和基础的 Tarjan 算法理解的较为深刻,但魔改后的求割点与桥以及求点/边双连通分量的 Tarjan 算法着实神奇,越看越懵 ...
- (寒假开黑gym)2018 ACM-ICPC, Syrian Collegiate Programming Contest
layout: post title: (寒假开黑gym)2018 ACM-ICPC, Syrian Collegiate Programming Contest author: "luow ...
- CNCC 2018 经典计算机算法技术论坛全解读 | CNCC 2018
雷锋网 AI 科技评论按:计算机算法是计算机科学的灵魂和基石.在历史长河的检验中,一些计算机算法成为经典留存至今,随着计算机科学的日益发展而不断推陈出新,并对现今的计算机科学产生着重要的影响.然而,当 ...
最新文章
- 以为精通Java 线程池,看到这些误区,还是年轻了...
- ESP8266串口处理
- C#中使用jieba.NET、WordCloudSharp制作词云图
- 【译】用 JavaScript 和 Emoji 做地址栏动画
- Solidity语言基础 和 Etherum ERC20合约基础
- java中构造器的调用顺序
- CodeForces - 976C(结构体排序思维)
- 内置式永磁同步电机IPMSM,基于虚拟信号注入法最大转矩电流比MTPA控制仿真模型
- Python学习笔记(4)-文件读取报错
- 进制转换: 十进制转换成其他进制(一)
- DedeCMS 后台地址泄露漏洞
- CentOS 7.3安装详解
- BUU-CRYPTO1密码学小白 25道入门题 详细解题思路
- 计算几何基础--线段的性质
- python怎么使用int四舍五入_python浮点数舍入(ROUND)方式总结
- DNS无法解析IP_DNS之基本原理
- idea创建SpringBoot工程
- Centos7安装Promethus(普罗米修斯)监控系统完整版
- Android Studio教学视频118集(共18.2G)
- CSS3 动画属性 - 逆战班
热门文章
- 正则基础学习1 这俩天补上正则全部内容,让你简单的学会正则,欢迎收藏当文档
- js超出文字个数展示省略号
- ClickHouse分布式集群部署
- js按照指定月份确认该月的开始时间和结束时间
- 利用ES6-Promise()方法封装原始jsonp实现跨域请求公用方法(告别使用JQuery封装好的jsonp)
- ZfNet解卷积:可视化CNN模型( PythonCode可视化Cifar10)
- http的请求方法 GET、HEAD、POST、PUT、DELETE、CONNECT、OPTIONS、TRACE 说明
- set和dict理论和操作
- 收购Roambi,SAP欲领导商务分析云迁移全球市场
- linux下devel包的作用