Solution -「LOCAL」「cov. HDU 6864」找朋友
\(\mathcal{Description}\)
Link.(几乎一致)
给定 \(n\) 个点 \(m\) 条边的仙人掌和起点 \(s\),边长度均为 \(1\)。令 \(d(u)\) 表示 \(u\) 到 \(s\) 的最短距离。对于任意一个结点的排列 \(\{p_1,p_2,\cdots,p_n\}\),记 \(t_i\) 满足 \(p_{t_i}=i\),称排列合法,当且仅当:
\]
求合法排列数,对 \(998244353\) 取模。
\(n\le10^4\),\(m\le2\times10^4\),保证不存在 \((u,v)\in E\),使得 \(d(u)=d(v)\)。
\(\mathcal{Solution}\)
考虑一个偶环(题目保证无奇环),起点终点在左右两端,上下各有 \(l\) 个结点相连。可见上下的点间是互不影响的,我们只需要分别保证上方和下方结点的相对位置。
再考虑一棵树,每个结点必须先于其子树内的点出现,所有方案为 \(n!\),每个结点 \(u\) 就会使其 \(\times\frac{1}{siz_u}\)。
对于仙人掌,处理出一棵 BFS 树,并得到环的信息。对于非环上的点,直接按树上的点来贡献系数。否则,对于一个环,如图:
DP 求解,当前子树顺序已确定,令 \(f(i,j)\) 表示用左边前 \(i\) 个和右边前 \(j\) 个时对答案贡献的系数。转移比较显:
\]
其中 \(siz_i\) 表示 \(i\) 在 BFS 树上的子树大小,需要特殊处理 \(i=0\) 或 \(j=0\) 的情况。
\(\mathcal{Code}\)
#include <queue>
#include <cstdio>
#include <vector>
typedef std::pair<int, int> pii;
const int MAXN = 1e4, MAXM = 2e4, MOD = 998244353;
int n, m, s, ecnt = 1, inv[MAXN + 5], head[MAXN + 5], dist[MAXN + 5];
int fa[MAXN + 5], siz[MAXN + 5], sL[MAXN + 5], sR[MAXN + 5], f[MAXN + 5][MAXN + 5];
bool cut[MAXM + 5], vis[MAXN + 5];
std::vector<pii> cir;
struct Edge { int to, nxt; } graph[MAXM * 2 + 5];
inline void link ( const int s, const int t ) {
graph[++ ecnt] = { t, head[s] };
head[s] = ecnt;
}
inline void initBFTree () {
std::queue<int> que;
que.push ( s ), dist[s] = 1;
while ( !que.empty () ) {
int u = que.front (); que.pop ();
for ( int i = head[u], v; i; i = graph[i].nxt ) {
if ( !dist[v = graph[i].to] ) {
fa[v] = u, dist[v] = dist[u] + 1;
que.push ( v );
} else if ( dist[v] > dist[u] ) {
cut[i >> 1] = true;
cir.push_back ( pii ( u, v ) );
}
}
}
}
inline void initSize ( const int u ) {
siz[u] = 1;
for ( int i = head[u], v; i; i = graph[i].nxt ) {
if ( !cut[i >> 1] && ( v = graph[i].to ) ^ fa[u] ) {
initSize ( v ), siz[u] += siz[v];
}
}
}
int main () {
freopen ( "abgfriend.in", "r", stdin );
freopen ( "abgfriend.out", "w", stdout );
scanf ( "%d %d %d", &n, &m, &s );
int ans = inv[1] = 1;
for ( int i = 2; i <= n; ++ i ) {
ans = 1ll * i * ans % MOD;
inv[i] = 1ll * ( MOD - MOD / i ) * inv[MOD % i] % MOD;
}
for ( int i = 1, u, v; i <= m; ++ i ) {
scanf ( "%d %d", &u, &v );
link ( u, v ), link ( v, u );
}
initBFTree ();
initSize ( s );
for ( int i = 0; i ^ cir.size (); ++ i ) {
int u = cir[i].first, v = cir[i].second, cnt = 0;
for ( int p = u, q = fa[v]; p ^ q; p = fa[p], q = fa[q] ) {
vis[p] = vis[q] = true, ++ cnt;
sL[cnt] = siz[p], sR[cnt] = siz[q];
}
for ( int i = 0; i <= cnt; ++ i ) {
for ( int j = 0; j <= cnt; ++ j ) {
if ( !i && !j ) f[i][j] = 1;
else if ( !i ) f[i][j] = 1ll * f[i][j - 1] * inv[sR[j]] % MOD;
else if ( !j ) f[i][j] = 1ll * f[i - 1][j] * inv[sL[i] + siz[v]] % MOD;
else f[i][j] = 1ll * ( f[i - 1][j] + f[i][j - 1] ) * inv[sL[i] + sR[j]] % MOD;
}
}
ans = 1ll * ans * f[cnt][cnt] % MOD;
}
for ( int i = 1; i <= n; ++ i ) {
if ( !vis[i] ) {
ans = 1ll * ans * inv[siz[i]] % MOD;
}
}
printf ( "%d\n", ans );
return 0;
}
\(\mathcal{Details}\)
一开始局部变量 cnt
没赋初值,Windows 贴心地帮助兔子清了零,然后在 Lemon 上测 RE 一大片 qwq……
Solution -「LOCAL」「cov. HDU 6864」找朋友的更多相关文章
- Solution -「LOCAL」「cov. HDU 6816」折纸游戏
\(\mathcal{Description}\) Link(削弱版). \(n\) 张纸叠在一起对折 \(k\) 次,然后从上到下为每层的正反两面写上数字,求把纸重新摊平后每张纸上的数字序列 ...
- Solution -「CTS 2019」「洛谷 P5404」氪金手游
\(\mathcal{Description}\) Link. 有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...
- 「题解」「美团 CodeM 资格赛」跳格子
目录 「题解」「美团 CodeM 资格赛」跳格子 题目描述 考场思路 思路分析及正解代码 「题解」「美团 CodeM 资格赛」跳格子 今天真的考自闭了... \(T1\) 花了 \(2h\) 都没有搞 ...
- 【翻译】西川善司的「实验做出的游戏图形」「GUILTY GEAR Xrd -SIGN-」中实现的「纯卡通动画的实时3D图形」的秘密,后篇
http://www.4gamer.net/games/216/G021678/20140714079/ 连载第2回的本回, Arc System Works开发的格斗游戏「GUILTY G ...
- Android内存管理(4)*官方教程 含「高效内存的16条策略」 Managing Your App's Memory
Managing Your App's Memory In this document How Android Manages Memory Sharing Memory Allocating and ...
- SSH连接时出现「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」解决办法
用ssh來操控github,沒想到連線時,出現「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」,後面還有一大串英文,這時當然要向Google大神求助 ...
- 「Windows MFC 」「Edit Control」 控件
「Windows MFC 」「Edit Control」 控件
- 「ZJOI2019」&「十二省联考 2019」题解索引
「ZJOI2019」&「十二省联考 2019」题解索引 「ZJOI2019」 「ZJOI2019」线段树 「ZJOI2019」Minimax 搜索 「十二省联考 2019」 「十二省联考 20 ...
- Loj #6069. 「2017 山东一轮集训 Day4」塔
Loj #6069. 「2017 山东一轮集训 Day4」塔 题目描述 现在有一条 $ [1, l] $ 的数轴,要在上面造 $ n $ 座塔,每座塔的坐标要两两不同,且为整点. 塔有编号,且每座塔都 ...
随机推荐
- Cannot uninstall 'pyparsing'. It is a distutils installed project
我的环境: [root@ansible ~]# python -V Python 2.7.5 [root@ansible ~]# cat /etc/redhat-release CentOS Linu ...
- Go语言系列之标准库fmt
fmt包实现了类似C语言printf和scanf的格式化I/O.主要分为向外输出内容和获取输入内容两大部分. 向外输出 标准库fmt提供了以下几种输出相关函数. Print Print系列函数会将内容 ...
- 前端 vue 等刷新清浏览器缓存的方法
vue开发过程中发现内存一直往上飙,360时不时的提示下清理内存,测试发现每刷新一次页面内存就增加1%. 清掉浏览器的缓存等各种数据以后内存明显变小,刷新页面导致内存爆表. 这种问题着实不应该存在,记 ...
- 51 Nod 1133 不重叠的线段 (贪心算法)
原题链接:https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1133 题目分析:感觉这到第不应该被分到二级算法题,比 109 ...
- JAVA-JDK1.8-ConCurrentHashMap-源码并且debug说明
概述 在上述的随笔中已经介绍了JDK1.7版本的ConCurrentHashMap源码和测试了,现在这篇随笔主要介绍JDK1.8版本的ConCurrentHashMap,这个版本抛弃了分段锁的实现,直 ...
- 从AAB文件生成APK文件
开头先说方法 1.下载 bundletool 2.开始生成 java -jar bundletool.jar build-apks --bundle=[aab文件路径(例如:C:\Users\Admi ...
- 微服务架构 | 3.2 Alibaba Nacos 注册中心
目录 前言 1. Nacos 基础知识 1.1 Nacos 命名方式 1.2 Nasoc 是什么 1.3 Nacos 的 4 个关键特性 1.4 Nacos 生态图 1.5 Nacos 架构图 1.6 ...
- Rust 使用 dotenv 来设置环境变量
在项目中,我们通常需要设置一些环境变量,用来保存一些凭证或其它数据,这时我们可以使用 dotenv 这个 crate. 首先在项目中添加 dotenv 这个依赖: 例如在下面这个项目中,需要设置数据库 ...
- golang中闭包讲解
1. 非闭包时的结果 package main import "fmt" func main() { var funcSlice []func() for i := 0; i &l ...
- java-异常-异常应用
1 package p1.exception; 2 3 4 /* 5 * 老师用电脑上课. 6 * 7 * 问题领域中涉及两个对象. 8 * 老师,电脑. 9 * 10 * 分析其中的问题. 11 * ...