HDU 6116 路径计数
时间:2019-12-21
本文章向大家介绍HDU 6116 路径计数,主要包括HDU 6116 路径计数使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
HDU 6116 路径计数
普通生成函数常用于处理组合问题,指数生成函数常用于处理排列问题。
考虑 对于 $ a $ 个 $ A $ 分为很多堆,这么分的方案数是 $ C_{a-1}^{i-1} $
然后对于每一堆我们看成一个数来放,并且所有堆都这样做,这样的话总的方案数量是 $ \frac{(i+j+k+l)!}{i!j!k!l!} $
就算所有一堆看成的数的排列是不存在相邻相等的,至少都有 $ n-i-j-k-l $ 对相邻的相同的数。
然后就可以容斥了,枚举 $ i+j+k+l $ 直接计算就好了。
$ ans = \displaystyle \sum_{x=1}^{n} (-1)^{n-x} x! \sum_{i+j+k+l=x} \frac{C_{a-1}^{i-1} C_{b-1}^{j-1} C_{c-1}^{k-1} C_{d-1}^{l-1}}{i ! j ! k ! l !} $
后面其实就是四个指数生成函数乘积了。
#include<bits/stdc++.h>
using namespace std;
#define MAXN 200010
#define P 998244353
#define clr( a ) memset( a , 0 , sizeof a )
typedef long long ll;
int wn[2][MAXN];
int Pow( int x , int y ) {
int res=1;
while(y) {
if(y&1) res=res*(ll)x%P;
x=x*(ll)x%P,y>>=1;
}
return res;
}
void getwn(int l) {
for(int i=1;i<(1<<l);i<<=1) {
int w0=Pow(3,(P-1)/(i<<1)),w1=Pow(3,P-1-(P-1)/(i<<1));
wn[0][i]=wn[1][i]=1;
for(int j=1;j<i;++j)
wn[0][i+j]=wn[0][i+j-1]*(ll)w0%P,
wn[1][i+j]=wn[1][i+j-1]*(ll)w1%P;
}
}
int rev[MAXN];
void getr(int l) { for(int i=1;i<(1<<l);++i) rev[i]=(rev[i>>1]>>1)|((i&1)<<l-1); }
void NTT(int *A,int len,int f) {
for(int i=0;i<len;++i) if(rev[i]<i) swap(A[i],A[rev[i]]);
for(int l=1;l<len;l<<=1)
for(int i=0;i<len;i+=(l<<1))
for(int k=0;k<l;++k) {
int t1=A[i+k],t2=A[i+l+k]*(ll)wn[f][l+k]%P;
A[i+k]=(t1+t2)%P;
A[i+l+k]=(t1-t2+P)%P;
}
if( f == 1 ) for(int inv=Pow(len,P-2),i=0;i<len;++i) A[i]=A[i]*(ll)inv%P;
}
int a , b , c , d;
int J[MAXN] , invJ[MAXN] , inv[MAXN];
int cc( int a , int b ) {
if( b > a ) return 0;
return 1ll * J[a] * invJ[b] % P * invJ[a - b] % P;
}
int A[MAXN] , B[MAXN] , C[MAXN] , D[MAXN];
int main() {
J[0] = inv[1] = invJ[0] = J[1] = invJ[1] = 1;
for( int i = 2 ; i < MAXN ; ++ i ) inv[i] = 1ll * ( P - P / i ) * inv[P % i] % P , J[i] = 1ll * J[i - 1] * i % P , invJ[i] = 1ll * invJ[i - 1] * inv[i] % P;
while( cin >> a >> b >> c >> d ) {
clr( A ) , clr( B ) , clr( C ) , clr( D );
int n = a + b + c + d;
for( int i = 1 ; i <= a ; ++ i ) A[i] = 1ll * cc( a - 1 , i - 1 ) * invJ[i] % P;
for( int i = 1 ; i <= b ; ++ i ) B[i] = 1ll * cc( b - 1 , i - 1 ) * invJ[i] % P;
for( int i = 1 ; i <= c ; ++ i ) C[i] = 1ll * cc( c - 1 , i - 1 ) * invJ[i] % P;
for( int i = 1 ; i <= d ; ++ i ) D[i] = 1ll * cc( d - 1 , i - 1 ) * invJ[i] % P;
int len = 1 , l = 0;
while( len <= n ) len <<= 1 , ++ l;
getwn( l ) , getr( l );
NTT( A , len , 0 ) , NTT( B , len , 0 ) , NTT( C , len , 0 ) , NTT( D , len , 0 );
for( int i = 0 ; i < len ; ++ i ) A[i] = 1ll * A[i] * B[i] % P * C[i] % P * D[i] % P;
NTT( A , len , 1 );
ll res = 0;
for( int i = 1 ; i <= n ; ++ i )
res += ( ( n - i & 1 ) ? -1ll : 1ll ) * J[i] * A[i] % P , res += P , res %= P;
cout << res << endl;
}
}
原文地址:https://www.cnblogs.com/yijan/p/hdu6116.html
- PHP数据结构(二十六) ——基数排序实现36进制数排序
- Apache配置
- jquery事件
- 设计模式专题(二)——策略模式
- ASP.NET AJAX(10)__Authentication ServiceAuthentication ServiceAuthentication Service属性Authentication
- 高效开发 MVVM 和 databinding 你需要使用的工具
- ASP.NET AJAX(9)__Profile Service什么是ASP.NET Profile如何使用ASP.NET ProfileProfile ServiceProfile Service预
- 设计模式专题(三)——装饰模式
- ASP.NET AJAX(8)__Microsoft AJAX Library中异步通信层的使用什么是异步通信层Micorsoft AJAX Library异步通信层的组成WebRequestExec
- ASP.NET AJAX(7)_Microsoft AJAX Library扩展客户端组件继承时需要注意的问题扩展类型如何修改已有类型
- ASP.NET AJAX(6)__Microsoft AJAX Library中的面向对象类型系统命名空间类类——构造函数类——定义方法类——定义属性类——注册类类——抽象类类——继承类——调用父类方
- 设计模式专题(四)——代理模式
- Array数组函数(一)
- ASP.NET AJAX(5)__JavaScript原生类型以及Microsoft AJAX Library什么是Microsoft AJAX LibraryObject原生类型Object.pro
- 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 数组属性和方法
- 常见编程模式之快慢指针
- python pywifi模块——暴力破解wifi
- 面试题系列第3篇:Integer等号判断的内幕,你可能不知道?
- Go by Example 中文:工作池
- 推荐一款万能抓包神器:Fiddler Everywhere
- 猿实战04——el-upload结合nginx之通用图片处理
- 30 多个有内味道且笑死的人代码注释
- Logstash-input-jdbc 同步 mysql 准实时数据至 ElasticSearch 搜索引擎
- 总结一些,我在书写 CSS 的时候,经常犯的错误!
- 通俗理解 set,dict 背后的哈希表
- K8S 生态周报| Google 选择 Cilium 作为 GKE 下一代数据面
- [Introduction]万字手撕Go http源码server.go
- Python 3.9 值得关注的更新点
- 趣味算法:JS实现红绳算法(匹配合适的另一半)
- 2019夏令营笔试编程题