\(\mathcal{Description}\)

  Link.

  给定 \(01\) 序列 \(\{A_n\}\) 和 \(\{B_n\}\),其中 \(1\) 的个数均为 \(k\)。记 \(A\) 中 \(1\) 的位置为 \(\{a_k\}\),\(B\) 中的为 \(\{b_k\}\)。现任意排列 \(\{a_k\}\) 和 \(\{b_k\}\),然后依次交换 \(A_{a_i}\) 和 \(A_{b_i}\),\(i=1,2,\dots,k\)。求使操作完成后 \(A=B\) 的排列方案数对 \(998244353\) 取模的结果。

  \(n\le10^4\)。

\(\mathcal{Solution}\)

  首先转化题意——每次操作用 \(A\) 的一个 \(1\) 和 \(B\) 中一个需要 \(1\) 的位置交换,共进行 \(k\) 次操作,求满足条件的方案数。

  我们定义 \(A_i=1\land B_i=0\) 的位置为“富余点”——它们需要为其它位置提供 \(1\);\(A_i=B_i=1\) 的位置为“公共点”——它们已经满足条件,但可以作为传递 \(1\) 的载体;\(A_i=0\land B_i=1\) 的位置为“缺失点”——它们是“富余点”需要提供到的位置。

  那么,只需要每个“富余点”的 \(1\) 都传递给“缺失点”,\(A\) 就与 \(B\) 相等。考虑一次从一个“富余”到“缺失”的传递,如图:

  其中三角形为“缺失点”,圆形为“公共点”(任意多个,亦可不存在),方形为“富余点”。注意若存在“公共点”,传递是不能反向的,这样会导致与三角形相邻的公共点的 \(1\) 变为 \(0\)。(某课件的笔误,注意一下 www。)


  考虑 DP,记 \(f(i,j)\) 表示在传递链中用了 \(i\) 个“公共点”,用了 \(j\) 个“富余点”(即有 \(j\) 条传递链)时,传递链中的方案数。边界为 \(f(0,0)=1\),转移:

\[f(i,j)=f(i-1,j)\cdot ij+f(i,j-1)\cdot j^2
\]

  前一项,将一个公共点加入一条传递链的末尾。有 \(j\) 条链,新的点可以与已有的 \(i-1\) 个点交换(不是交换操作顺序,而是直接交换链中位置),故有系数 \(ij\)。

  后一项,新建一条“缺失”-“富余”链。首先拿出新的一对“缺失点”和“富余点”,仍考虑到原来的“缺失点”或“富余点”可以和新点交换位置,故有系数 \(j^2\)。

  最后统计答案,发现“公共点”没有必要在链中用完,设在链中用 \(a\) 个“公共点”,还剩下 \(b\) 个。那么在所有“公共点”中选出 \(b\) 个,方案数 \(\binom{a+b}{b}\);公共点内部方案 \((b!)^2\);把这 \(b\) 次操作安排进总共 \(k\) 次操作里,方案数 \(\binom{k}{b}\)。所以设“公共点”有 \(s\) 个,“富余点”有 \(t\) 个,答案为(\(i\) 枚举的是不在链中的“公共点”个数):

\[\sum_{i=0}^{s}\binom{s}i\binom{s+t}{i}(i!)^2f(s-i,t)
\]

  复杂度 \(\mathcal O(n^2)\)。

\(\mathcal{Code}\)

#include <cstdio>
#include <cstring> const int MAXL = 1e4, MOD = 998244353;
int n, sur, bal, f[MAXL + 5][MAXL + 5];
int fac[MAXL + 5], ifac[MAXL + 5];
char A[MAXL + 5], B[MAXL + 5]; inline void addeq ( int& a, const int b ) { if ( ( a += b ) >= MOD ) a -= MOD; }
inline int mul ( long long a, const int b ) { return ( a *= b ) < MOD ? a : a % MOD; } inline int qkpow ( int a, int b ) {
int ret = 1;
for ( ; b; a = mul ( a, a ), b >>= 1 ) ret = mul ( ret, b & 1 ? a : 1 );
return ret;
} inline void init () {
fac[0] = 1;
for ( int i = 1; i <= n; ++ i ) fac[i] = mul ( i, fac[i - 1] );
ifac[n] = qkpow ( fac[n], MOD - 2 );
for ( int i = n - 1; ~ i; -- i ) ifac[i] = mul ( i + 1, ifac[i + 1] );
} inline int comb ( const int n, const int m ) {
return n < m ? 0 : mul ( fac[n], mul ( ifac[m], ifac[n - m] ) );
} int main () {
scanf ( "%s %s", A + 1, B + 1 );
n = strlen ( A + 1 ), init ();
for ( int i = 1; i <= n; ++ i ) {
if ( A[i] ^ '0' && B[i] ^ '0' ) ++ bal;
else if ( A[i] > B[i] ) ++ sur;
}
f[0][0] = 1;
for ( int i = 0; i <= bal; ++ i ) {
for ( int j = 0, cur; j <= sur; ++ j ) {
if ( ! ( cur = f[i][j] ) ) continue;
addeq ( f[i + 1][j], mul ( cur, mul ( i + 1, j ) ) );
addeq ( f[i][j + 1], mul ( cur, mul ( j + 1, j + 1 ) ) );
}
}
int ans = 0;
for ( int i = 0; i <= bal; ++ i ) {
int fre = bal - i;
int self = mul ( mul ( fac[fre], fac[fre] ),
mul ( comb ( bal + sur, fre ), comb ( bal, fre ) ) );
addeq ( ans, mul ( f[i][sur], self ) );
}
printf ( "%d\n", ans );
return 0;
}

Solution -「AGC 019E」「AT 2704」Shuffle and Swap的更多相关文章

  1. Solution -「CTS 2019」「洛谷 P5404」氪金手游

    \(\mathcal{Description}\)   Link.   有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...

  2. 「题解」「美团 CodeM 资格赛」跳格子

    目录 「题解」「美团 CodeM 资格赛」跳格子 题目描述 考场思路 思路分析及正解代码 「题解」「美团 CodeM 资格赛」跳格子 今天真的考自闭了... \(T1\) 花了 \(2h\) 都没有搞 ...

  3. 【翻译】西川善司的「实验做出的游戏图形」「GUILTY GEAR Xrd -SIGN-」中实现的「纯卡通动画的实时3D图形」的秘密,后篇

    http://www.4gamer.net/games/216/G021678/20140714079/     连载第2回的本回,  Arc System Works开发的格斗游戏「GUILTY G ...

  4. Android内存管理(4)*官方教程 含「高效内存的16条策略」 Managing Your App's Memory

    Managing Your App's Memory In this document How Android Manages Memory Sharing Memory Allocating and ...

  5. SSH连接时出现「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」解决办法

    用ssh來操控github,沒想到連線時,出現「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」,後面還有一大串英文,這時當然要向Google大神求助 ...

  6. 「Windows MFC 」「Edit Control」 控件

    「Windows MFC 」「Edit Control」 控件

  7. 「ZJOI2019」&「十二省联考 2019」题解索引

    「ZJOI2019」&「十二省联考 2019」题解索引 「ZJOI2019」 「ZJOI2019」线段树 「ZJOI2019」Minimax 搜索 「十二省联考 2019」 「十二省联考 20 ...

  8. Loj #6069. 「2017 山东一轮集训 Day4」塔

    Loj #6069. 「2017 山东一轮集训 Day4」塔 题目描述 现在有一条 $ [1, l] $ 的数轴,要在上面造 $ n $ 座塔,每座塔的坐标要两两不同,且为整点. 塔有编号,且每座塔都 ...

  9. Loj #6073.「2017 山东一轮集训 Day5」距离

    Loj #6073.「2017 山东一轮集训 Day5」距离 Description 给定一棵 \(n\) 个点的边带权的树,以及一个排列$ p\(,有\)q $个询问,给定点 \(u, v, k\) ...

随机推荐

  1. 深入谈谈 Java IOC 和 DI

    1.前言 不得不说, IOC和DI 在写代码时经常用到.还有个就是在面试时 ,面试官老喜欢问 IOC 和DI是什么的问题,都快被问吐了, 可是,仍然会让许多人说的支支吾吾. 为什么? 第一,因为这个知 ...

  2. windows更改pip源(可用)

    中国科学技术大学 : https://pypi.mirrors.ustc.edu.cn/simple 清华:https://pypi.tuna.tsinghua.edu.cn/simple 豆瓣:ht ...

  3. Spark案例练习-PV的统计

    关注公众号:分享电脑学习回复"百度云盘" 可以免费获取所有学习文档的代码(不定期更新) 云盘目录说明: tools目录是安装包res   目录是每一个课件对应的代码和资源等doc  ...

  4. 利用Word2010制作流程图

    利用Word2010制作流程图 原文链接:https://www.toutiao.com/i6483034968225235469/ 一.页面和段落的设置 启动Word2010,打开一个空白文档,并切 ...

  5. [源码解析] PyTorch 分布式之 ZeroRedundancyOptimizer

    [源码解析] PyTorch 分布式之 ZeroRedundancyOptimizer 目录 [源码解析] PyTorch 分布式之 ZeroRedundancyOptimizer 0x00 摘要 0 ...

  6. RHCSA 第八天

    1.查询ip的几种方式: ip, ifconfig, nmcli,nmtui 2.nmcli命令使用: a.在ens160网卡上新建连接static_con,并配置静态ip b.在ens160网卡上新 ...

  7. MRCTF2020 你传你🐎呢

    MRCTF2020 你传你 .htaccess mime检测 1.先尝试上传了一个文件,发现.jpg后缀的可以上传成功,但是用蚁剑连接时返回空数据 2.重新先上传一个.htaccess文件,让需要被上 ...

  8. 【涨姿势】原来golang的case <-time.After(xxx)还有这样的坑

    偶然看到这样一篇文章:<使用 pprof 排查 Golang 内存泄露>https://www.toutiao.com/i6881796351139676680/ 最后一段让我很疑惑: 修 ...

  9. Cesium1.70-介绍CesiumOSM建筑新特性

    Cesium中文网:http://cesiumcn.org/ | 国内快速访问:http://cesium.coinidea.com/ 我们很高兴宣布Cesium OSM建筑,一个覆盖整个世界的3D建 ...

  10. 【Vulnhub靶场】JANGOW: 1.0.1

    时隔这么久,终于开始做题了 环境准备 下载靶机,导入到virtualBox里面,这应该不用教了吧 开机可以看到,他已经给出了靶机的IP地址,就不用我们自己去探测了 攻击机IP地址为:192.168.2 ...