这题好牛逼啊。。。

虽然说也是 DP 套 DP,但是感觉比 TJOI 那题高明到哪里去了(

我们先考虑如何计算期望。如果设 \(f_i\) 为拿到 \(i\) 张牌后胡的方案数,这个并不是很好做,因为你要考虑前面不胡。

所以我们设 \(f_i\) 为手上有 \(13+i\) 张牌但还没有胡的方案数,那么答案就是:

\[\frac{\sum_{i=1}^{4n-13}f_ii!(4n-13-i)!}{(4n-14)!}+1
\]

发现一个手上牌的顺序不影响是否胡,所以有两个阶乘用来打乱排序(因为牌与牌直接互不相同)。

我们考虑外层的 DP 建立一个“胡牌自动机”,所以需要写一个用来判定是否胡牌的 DP。

加入新的牌就是在自动机上面走一步。我们设 \(dp_{i,x,y,0/1}\) 为加入了前 \(i\) 类牌,现在手上有 \(x\) 对 \((i-1,i)\) 的牌和 \(y\) 张 \(i\),\(0/1\) 表示当前无/有对子。表示的是在这个状态下面子最多个数。

我们将发现如果 \(x,y\geq 3\),那么可以自动对出来一个面子,所以有 \(0\leq x,y<3\)。

我们把这两个维度写到矩阵上面。然后交换维度变成 \(dp_{i,0/1,x,y}\)。

考虑加入 \(x\) 张牌时如何从 \(dp_{i-1}\) 转移到 \(dp_i\)。

考虑加入 \(x\) 张牌的时候,先凑到 \((i-2,i-1)\) 里面去,然后自己和自己凑面子,最后加上已经有的面子个数。

如果 \(x\geq 2\),那么我们可以砍掉一对对子再加进去。

也就是:\((f_{i-1,0},x)\to f_{i,0},(f_{i-1,1},x)\to f_{i,1},(f_{i-1,0},x-2)\to f_{i,1}\)。

然后我们再额外记录对子个数即可。

然后我们用一个 std::map 离散化一下就可以建出“胡牌自动机”了。

这个是主要的,在这个 DFA 上面 DP 反而很简单。

直接考虑当前加 \(x\) 张牌,然后从牌堆里面把这 \(x\) 张牌选出来,丢到下一个 DFA 节点里面去就可以了。

我们可以把所有胡掉的节点压成一个节点,然后标记 \(t=-1\),可以节省空间。

#include<cstdio>
#include<map>
typedef int ui;
const ui M=105,mod=998244353;
const ui C[5][5]={{1,0,0,0,0},{1,1,0,0,0},{1,2,1,0,0},{1,3,3,1,0},{1,4,6,4,1}};
ui n,m,a[M],fac[M*4],ifac[M*4],dp[2][M*4][2100];
inline ui min(const ui&a,const ui&b){
return a>b?b:a;
}
inline void max(ui&a,const ui&b){
if(b>a)a=b;
}
struct HuAutomation{
ui tot;
struct Matrix{
#define For for(ui i=0;i^3;++i)for(ui j=0;j^3;++j)
ui p[3][3];
Matrix(){
For p[i][j]=-1;
}
inline ui*operator[](const ui&x){
return p[x];
}
inline bool operator==(Matrix q)const{
For if(p[i][j]!=q[i][j])return false;return true;
}
inline bool operator<(Matrix q)const{
For if(p[i][j]!=q[i][j])return p[i][j]<q[i][j];return false;
}
inline bool CheckHu()const{
For if(~p[i][j]&&p[i][j]>3)return true;return false;
}
inline void Update(Matrix bel,const ui&x){
For if(~bel[i][j])for(ui k=0;k^3&&i+j+k<=x;++k){
max(p[j][k],min(i+bel[i][j]+(x-i-j-k)/3,4));
}
}
#undef For
};
struct Node{
ui t;Matrix P[2];
inline Node(){
t=0;P[0]=P[1]=Matrix();
}
inline bool operator<(Node it)const{
return t==it.t?P[0]==it.P[0]?P[1]==it.P[1]?0:P[1]<it.P[1]:P[0]<it.P[0]:t<it.t;
}
inline Node Hu()const{
Node x;x.t=-1;return x;
}
inline bool CheckHu()const{
return!~t||t>=7||P[1].CheckHu();
}
inline Node operator+(const ui&x){
if(CheckHu())return Hu();Node ans;
ans.P[0].Update(P[0],x);ans.P[1].Update(P[1],x);ans.t=t;
if(x>=2)ans.P[1].Update(P[0],x-2),++ans.t;
if(ans.CheckHu())ans=Hu();
return ans;
}
}tr[2100];
ui trans[2100][5];
inline ui*operator[](const ui&x){
return trans[x];
}
inline Node Hu()const{
return tr[0].Hu();
}
std::map<Node,ui>P;
#define pos(x) (P.count(x)?P[x]:(tr[P[x]=++tot]=x,tot))
inline void Build(){
tr[1].P[0][0][0]=0;P[tr[1]]=++tot;P[tr[2]=Hu()]=++tot;
for(ui i=1;i<=tot;++i)if(i!=2){
for(ui t=0;t<=4;++t)trans[i][t]=pos(tr[i]+t);
}
}
}Hu;
signed main(){
ui ans(0);
Hu.Build();scanf("%u",&n);m=(n<<2)-13;dp[0][0][1]=1;
fac[0]=fac[1]=ifac[0]=ifac[1]=1;
for(ui i=2;i<=m;++i)fac[i]=1ull*fac[i-1]*i%mod,ifac[i]=1ull*(mod-mod/i)*ifac[mod%i]%mod;
for(ui i=2;i<=m;++i)ifac[i]=1ull*ifac[i-1]*ifac[i]%mod;
for(ui i=1;i<=13;++i){
ui w,t;scanf("%u%u",&w,&t);++a[w];
}
for(ui i=0;i<n;++i){
const ui&q=i&1,p=q^1;
for(ui j=0;j<=m;++j)for(ui S=1;S<=Hu.tot;++S)if(dp[q][j][S]){
for(ui x=0;x<=4-a[i+1];++x){
dp[p][j+x][Hu[S][a[i+1]+x]]=(dp[p][j+x][Hu[S][a[i+1]+x]]+1ull*C[4-a[i+1]][x]*dp[q][j][S])%mod;
}
}
for(ui i=0;i<=m;++i)for(ui S=1;S<=Hu.tot;++S)dp[q][i][S]=0;
}
for(ui i=1;i<=m;++i)for(ui S=1;S<=Hu.tot;++S)if(S!=2){
ans=(ans+1ull*fac[i]*fac[m-i]%mod*dp[n&1][i][S])%mod;
}
printf("%u",1ull*ans*ifac[m]%mod+1);
}

LGP5279题解的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

随机推荐

  1. linux内核中的eventfd

    转载请注明来源:https://www.cnblogs.com/hookjc/ eventfd 在内核版本,2.6.22以后有效.查看内核版本可以用命令 uname -r . [cpp] view p ...

  2. 一劳永逸Java环境配置,以及编写我的第一个Java程序

    Java环境配置,以及编写我的第一个Java程序 配置步骤 1.下载jdk 2.安装步骤 3.配置环境 4.我的第一个Java程序 配置步骤 网上的教程有很多,方法也都不尽相同.今天我就分享一下我的配 ...

  3. HTML笔记整理--上节

    一.认识WEB 「网页」主要是由文字.图像和超链接等元素构成,当然除了这些元素,网页中还可以包括音频.视频以及Flash等. 「浏览器」是网页显示.运行的平台. 「浏览器内核」(排版引擎.解释引擎.渲 ...

  4. 32、python并发编程之背景知识

    目录: 一 引子 二 为什么要有操作系统 三 什么是操作系统 四 操作系统与普通软件的区别 五 操作系统发展史 六 总结视频链接: 一 引子 顾名思义,进程即正在执行的一个过程.进程是对正在运行程序的 ...

  5. 阿里云服务器ECS挂载数据盘—linux系统

    参考阿里云官网帮助文档:https://help.aliyun.com/document_detail/25426.html 里面有些步骤说的不是很清楚,初学者可能操作时会遇到问题.通过这篇文档进行进 ...

  6. python控制浏览器上传文件

    自动化爬虫方法和库很多,难点大多数在登录.可以大致分为:普通验证码,扫码登录,QQ一键登录,拖动验证,无痕验证,人工识别(比如12306登录) 万能大法可以破解一切以上需求,自动控制浏览器行为 参考文 ...

  7. 设计模式在实际业务应用中的介绍之3——外观或门面模式Facade对AOP装配业务工厂的应用

    在C#中实现的基于外观或门面模式打造的业务应用案例 以前一直没有想过写一些东西来把项目中用到的知识点及技术实现做一个归纳整理并分享出来.现在打算逐渐的把项目中的一些东西整理并分享出来,与大家共勉! 外 ...

  8. for循环例子

    代码 点击查看[ForTest.java]代码 //package com.d; import java.util.Scanner; /** * For循环例子 * @date: 2022.2.24 ...

  9. 个人站bitlove.cn

    新启用了个人网站 微爱博客 http://bitlove.cn 有兴趣的可以移步交流

  10. 从数据分析系统总架构理解BI工具的价值所在

    ​现如今,应用商业智能BI工具的企业是越来越多了,由此也可见企业对数据分析的重视.因此,掌握一定的数据分析知识对"打工人"来说是非常重要的.现在小编就来跟大家一起来了解一下商业智能 ...