垃圾ACMer的暑假训练220720

最新推荐文章于 2024-07-22 10:02:43 发布

Hytidel

最新推荐文章于 2024-07-22 10:02:43 发布

阅读量117

点赞数

分类专栏： 2022ACM暑假训练文章标签： c++ 算法

本文链接：https://blog.csdn.net/Hytidel/article/details/125895193

版权

2022ACM暑假训练专栏收录该内容

23 篇文章 1 订阅

订阅专栏

垃圾ACMer的暑假训练220720

群论

35.1 Burnside引理与Polya定理

对集合 $G$ 和作用于 $G$ 的二元运算 $\cdot$ ,若满足如下四个性质,则称 $G$ 与 $\cdot$ 构成一个群,记作 $(G,\cdot)$ :

①封闭性:对 $\forall a,b\in G$ ,有 $a\cdot b\in G$ .

②结合律:对 $\forall a,b,c\in G$ ,有 $(a\cdot b)\cdot c=a\cdot(b\cdot c)$ .

③单位元:$\exists e\in G\ s.t.\ $对$ \forall a\in G $, 有$ a\cdot e=e\cdot a=a $, 则称$ e$为单位元.可以证明单位元是唯一的.

④逆元:对 $\forall a\in G,\exists b\in G\ s.t.\ a\cdot b=b\cdot a=e$ ,则称 $b$ 为 $a$ 的逆元,记作 $a^{-1}$ .可以证明逆元是唯一的.

设 $H$ 是 $G$ 的一个子集.若 $(H,\cdot)$ 为群,则称 $(H,\cdot)$ 为 $(G,\cdot)$ 的子群,记作 $H\leq G$ .

对 $\forall g\in G,h\in H$ :①称 $gH=g\cdot h$ 为 $H$ 在 $G$ 内关于 $g$ 的左陪集;②称 $Hg=h\cdot g$ 为 $H$ 在 $G$ 内关于 $g$ 的右陪集.

陪集的性质(以右陪集为例):

①对 $\forall g\in G$ ,有 $∣ H ∣ = ∣ H g ∣$ .

[证] 对$\forall $相异的$ h_1,h_2\in H $, 若$ h_1\cdot g=h_2\cdot g $, 两边右乘$ g^{-1} $得$ h_1=h_2 $, 矛盾 . 故$ h_1\cdot g\neq h_2\cdot g$.

②对 $\forall g\in G$ ,有 $g\in Hg$ ,进而 $Hg=H\Leftrightarrow g\in H$ .

[证] 注意到 $e\in H$ ,则 $e\cdot g\in Hg$ .

③ $Ha=Hb\Leftrightarrow a\cdot b^{-1}\in H$ .

[证] $Ha\cdot b^{-1}=H$ ,由②: $a\cdot b^{-1}\in H$ .

④若 $Ha\bigcap Hb\neq \varnothing$ ,则 $H a = H b$ ,这表明:子群 $H$ 的陪集要么为空集,要么相等.

[证] 取 $c=Ha\bigcap Hb$ ,则 $\exist h_1,h_2\in H\ s.t.\ h_1\cdot a=c=h_2\cdot b$ ,则 $a\cdot b^{-1}=h_2\cdot h_1^{-1}\in H$ .由③: $H a = H b$ .

⑤ $H$ 的全体右陪集的并为 $G$ .

[证] 对 $e\in H$ ,取遍 $g\in G$ 即证.

若 $H\leq G$ ,记 $G / H$ 为 $G$ 中所有 $H$ 的左陪集,即 $G/H=\{gH\ |\ g\in G\}$ .

若 $H\leq G$ ,记 $[G : H]$ 为 $G$ 中 $H$ 的不同陪集的数量.

设群 $(G,\cdot)$ .对 $\forall x\in G$ ,若 $\exists$ 最小的 $d\in\mathbb{Z}\ s.t.\ x^d=e$ ,则称 $d$ 为元素 $x$ 的阶.对有限群, $x$ 的阶都存在.

有限群 $G$ 的阶定义为群 $G$ 的元素个数,无限群的阶规定为 $0$ .

[Lagrange定理] 设有限群 $H$ 是有限群 $G$ 的子群,则 $∣ H ∣$ 整除 $∣ G ∣$ ,即 $H$ 的阶整除 $G$ 的阶,亦即 $|H|\cdot |G:H|=|G|$ .

[证] 由陪集的性质①④⑤:所有本质不同的陪集互不相交且阶都为 $∣ H ∣$ ,并为 $G$ ,并的阶为 $∣ G ∣$ .

[置换的双行表示法] $\sigma=\begin{pmatrix}1&2&3&4&5 \\ 2&5&4&3&1\end{pmatrix}$ 是上行到下行的一个置换,表示将排列 $1\ 2\ 3\ 4\ 5$ 映射为 $2\ 5\ 4\ 3 \ 1$ .

设集合 $N=\{1,2,\cdots,n\},M$ 为 $N$ 的若干个排列构成的集合.令群 $G(M,\cdot)$ ,其中 $\sigma$ 为置换.若 $G(M,\cdot)$ 满足群的性质,则称 $G$ 是一个置换群.

对集合 $M$ 和群 $G$ ,若给定二元函数 $\varphi(v,k)$ ,其中 $v$ 为群中的元素, $k$ 为集合中的元素,且 $\varphi(e,k)=k,\varphi(g,\varphi(s,k))=\varphi(g\cdot s,k)$ ,则称群 $G$ 作用于集合 $M$ .

对作用于集合 $X$ 的群 $G$ ,称 $x\in X$ 通过 $G$ 的元素能转移到的元素的集合为 $x$ 的轨道,记作 $G (x)$ .对 $g\in G$ ,记 $g (x)$ 为 $g$ 作用于 $x$ 的结果,即 $g(x)=\varphi(g,x)$ .

对 $g\in G$ ,称群 $G$ 中满足 $g (x) = x$ 的元素 $g$ 构成的集合为 $G$ 的稳定子,即 $G^x=\{g\ |\ g\in G,g(x)=x\}$ .

[定理35.1.1] $G^x\leq G$ .

[证] 显然 $e\in G^x$ ,且 $G^x$ 满足结合律.

封闭性:设 $f,g\in G^x$ ,则 $f (x) = x = g (x)$ ,进而 $(f\cdot g)(x)=x$ ,即 $f\cdot g\in G^x$ .

逆元:设 $g\in G^x$ ,则 $g (x) = x$ .因 $(g\cdot g^{-1})(x)=e(x)=x$ ,则 $g^{-1}(x)=x$ ,即 $g^{-1}\in G^x$ .

[例35.1.1] 给定一个 $2\times 2$ 的网格,每格可染成黑色或白色,设染色能得到的矩形的集合为 $M$ .

给定群 $G$ ,其成员:①顺时针旋转 $90^\circ$ ;②顺时针旋转 $180^\circ$ ;③顺时针旋转 $270^\circ$ ;④顺时针旋转 $0^\circ$ (模 $360^\circ$ ).

设 $0$ 表示白色, $1$ 表示黑色.对 $M$ 的一个元素 $\begin{bmatrix}1&1 \\ 0&0\end{bmatrix}$ ,其稳定子 $G^x$ 为 $\{顺时针旋转0^\circ\}$ ,

其轨道为 $\begin{bmatrix}1&1 \\ 0 &0\end{bmatrix},\begin{bmatrix}0&1 \\ 0 &1\end{bmatrix},\begin{bmatrix}0&0 \\ 1 &1\end{bmatrix},\begin{bmatrix}1&0 \\ 1 &0\end{bmatrix}$ .

[轨道-稳定子定理] $|G^x|\cdot |G(x)|=|G|$ ,即稳定子大小与轨道大小之积为群的大小.

[证] 因 $G^x\leq G$ ,由Lagrange定理: $|G^x|\cdot |G:G^x|=|G|$ .

只需证 $G:G^x|=|G(x)|$ ,即证每个 $g (x)$ 都能对应 $G:G^x]$ 中的一个左陪集或右陪集.

构造如下的对应关系:

若 $f (x) = g (x)$ ,则 $(f\cdot g^{-1})(x)=x=e(x)\in G^x$ .

由陪集的性质: $fG^x=gG^x$ ,这表明:相同的 $f (x)$ 都能对应相同的陪集,反之亦然.

对每个 $g (x)$ ,取 $gG^x$ 为其对应的陪集,则相同的 $g (x)$ 都能对应相同的陪集,故证.

设置换群 $G$ 作用于集合 $X$ .若 $\exists x,y\in X$ 在 $G$ 的作用下相等,即 $\exists f\in G\ s.t.\ f(x)=y$ ,则称 $x, y$ 属于同一等价类.

记 $X^g$ 为 $X$ 在 $g\in G$ 的作用下的不动点的的集合,即 $X^g=\{x|x\in X,g(x)=x\}$ .

[Burnside引理] 设 $A$ 和 $B$ 为有限集, $X$ 为一些从 $A$ 到 $B$ 的映射组成的集合, $G$ 是 $A$ 上的置换群,且 $X$ 中的映射在 $G$ 中的置换作用下封闭.设 $X / G$ 表示 $G$ 作用在 $X$ 上产生的所有等价类的集合,则 $\displaystyle |X/G|=\dfrac{1}{|G|}\sum_{g\in G}|X^g|$ ,且 $X^g=\{x|x\in X,g(x)=x\}$ .

[证] 因每个元素只属于一个轨道,轨道内部在群 $G$ 的作用下互达,则 $\displaystyle |X/G|=\sum_{x\in X}\dfrac{1}{|G:G^x|}$ .

由轨道-稳定子定理: $|G:G^x|=\dfrac{G}{|G^x|}$ ,则 $\displaystyle |X/G|=\sum_{x\in X}\dfrac{1}{|G:G^x|}=\dfrac{1}{|G|}\sum_{x\in X}G^x=\dfrac{1}{|G|}\sum_{g\in G}|X^g|$ .

[Polya定理] 在与Burnside引理相同的前提下,若 $X$ 为所有从 $A$ 到 $B$ 的映射,则 $\displaystyle|X/G|=\dfrac{1}{|G|}\sum_{g\in G}|B|^{c(g)}$ ,其中 $c (g)$ 表示置换 $g$ 能拆分成的不相交的循环置换的数量.

[证] Burnside引理中 $g (x) = x$ 的充要条件是: $x$ 将 $g$ 中的每个循环置换的元素都映射到 $B$ 中的同一元素.

故 $X|^g=|B|^{c(g)}$ ,即Poyla定理.

[例35.1.2] 用 $3$ 种颜色给一个正方体染色,求本质不同的方案数(经翻转后相同的两种方案视为同种).

[解1] 设 $A$ 为正方体 $6$ 个面的集合, $B$ 为 $3$ 种颜色的集合, $X$ 为直接给每个面染色,不考虑本质不同的方案的集合(共 $3^6$ 种方案), $G$ 为各种翻转操作构成的置换群, $X / G$ 为本质不同的染色方案的集合, $X^g$ 为某一翻转操作 $g$ 的所有直接染色方案中,经翻转 $g$ 后保持不变的染色方案的集合.

称正方体的 $6$ 个面分别为前、后、上、下、左、右.

$G$ 的所有置换可分为如下几类:

①不动:即恒等变换,所有直接染色方案经恒等变换都不变,对应的 $X^g|=3^6$ .

②以两相对面的中心连线为轴的 $90^\circ$ 旋转:相对面有 $3$ 种选择,旋转方向有 $2$ 种选择,则该类有 $6$ 个置换.

设选择前、后两面的中心连线为轴,则需上、下、左、右四个面颜色相同才能使旋转后不变,对应的 $X^g|=3^3$ .

③以两相对面的中心连线为轴的 $180^\circ$ 旋转:相对面有 $3$ 种选择,旋转方向的选择对置换无影响,则该类有 $3$ 个置换.

设选择前、后两面的中心连线为轴,则需上、下两面和左、右两面颜色相同才能使旋转后不变,对应的 $X^g|=3^4$ .

④以两条相对棱的终点连线为轴的 $180^\circ$ 旋转:相对棱有 $6$ 种选择,旋转方向对置换无影响,则该类有 $6$ 个置换.

设选择前、上量面的边界和下、后两面的边界为相对棱,则需前、上两面和下、后两面和左、右两面颜色相同

才能使旋转后不变,对应的 $X^g|=3^3$ .

⑤以两相对顶点连线为轴的 $120^\circ$ 旋转.相对顶点有 $4$ 种选择,旋转方向有 $2$ 种选择,则该类有 $8$ 个置换.

设选择前面的右上角和后面的左下角为相对顶点,则需前、上、有三个面和后、下、左三个面颜色相同

才能使旋转后不变,对应的 $X^g|=3^2$ .

综上,本质不同的方案数 $\dfrac{1}{1+6+3+6+8}(3^6+6\times 3^3+3\times 3^4+6\times 3^3+8\times 3^2)=57$ .

[解2] 分析解I中的④:设前、后、上、下、左、右 $6$ 个面依次编号为 $1\sim 6$ ,则置换可表示为 $\begin{pmatrix}1,3,2,4,5,6 \\ 3,1,4,2,6,5\end{pmatrix}=(1,3)\circ (2,4)\circ(5,6)$ ,其中翻转后 $1$ 面换到了 $3$ 面, $3$ 面换到了 $1$ 面,以此类推.

故 $c(g)=3,|B|^{c(g)}=3^3$ .

实际应用中,若循环置换间相互独立,则可用Polya定理;否则只能用Burnside引理.

35.1.1 Pólya 定理

题意

有 $t\ \ (1\leq t\leq 1000)$ 组测试数据.每组测试数据输入一个整数 $n\ \ (1\leq n\leq 1\mathrm{e}9)$ ,表示给定一个 $n$ 个点、 $n$ 条边的环.现用 $n$ 种颜色给每个点染色,求本质不同的染色方案数,答案对 $1\mathrm{e}9+7$ 取模.本质不同的染色方案定义为:不能通过旋转与别的染色方案相同.

思路

本质不同的 $n$ 个点的环可视为在置换群 $G=\{旋转0个,旋转1个,\cdots,旋转(n-1)个\}$ 作用下得到的等价类的数量.

定义集合 $M$ 为 $1\sim n$ 的全排列表示的环构成的集合.由Burnside引理: $\displaystyle ans=\dfrac{1}{|G|}\sum_{g\in G}|M^g|$ .

考虑每个置换对答案的贡献:

①旋转 $0$ 个的置换的不动点数为 $n^n$ ,即所有集合都合法.

②旋转 $k$ 个的置换,注意到一个元素是不动点等价于其存在一个长度为 $a$ 的循环节满足 $a\mid k$ ,而 $a\mid n$ ,则可改写判定条件为存在一个长度为 $\gcd(k,n)$ 的循环节.因该种置换中每个子串的前 $\gcd(k,n)$ 任取,则不动点数为 $n^{\gcd(k,n)}$ .

故 $\displaystyle ans=\dfrac{1}{n}\sum_{k=1}^n n^{\gcd(k,n)}\xlongequal{d=\gcd(k,n)}\dfrac{1}{n}\sum_{d\mid n}x^d\cdot \sum_{k=1}^\frac{n}{d}\left[\gcd\left(k,\dfrac{n}{d}\right)=1\right]=\dfrac{1}{n}\sum_{d\mid n}n^d\varphi\left(\dfrac{n}{d}\right)$

$\displaystyle\xlongequal{q=\frac{n}{d}}\dfrac{1}{n}\sum_{p\mid n}\varphi(p)\cdot n^\frac{n}{p}=\sum_{p\mid n}\varphi(p)\cdot n^{\frac{n}{p}-1}$ .

代码

const int MOD = 1e9 + 7;

int phi(int x) {
	int res = x;
	for (int i = 2; i <= x / i; i++) {
		if (x % i == 0) {
			res = res / i * (i - 1);
			while (x % i == 0) x /= i;
		}
	}
	if (x > 1) res = res / x * (x - 1);
	return res;
}

int main() {
	CaseT{
		int n; cin >> n;

		int ans = 0;
		for (int d = 1; d <= n / d; d++) {  // 枚举n的约数
			if (n % d)continue;

			ans = ((ll)ans + (ll)phi(d) * qpow(n, n / d - 1, MOD)) % MOD;
			if (d * d != n) ans = ((ll)ans + (ll)phi(n / d) * qpow(n, d - 1, MOD)) % MOD;
		}
		cout << ans << endl;
	}
}

35.1.2 串珠子

题意

给定 $m$ 种不同颜色的珠子,每种颜色的珠子个数足够多.现需从中选出 $n$ 个珠子,串成一个环形手链,问能制成多少种本质不同的手链.本质不同定义为:不能通过旋转或翻转与其他方案对应位置上的珠子颜色相同.

有多组测试数据.每组测试数据输入两整数 $m,n\ \ (n,m>0,nm\leq 32)$ ,最后一行输入 $0\ 0$ 表示输入结束.

对每组测试数据,输出能制成多少种本质不同的手链.

思路

考虑有多少种不同的置换及其循环数:

(1)旋转:转 $k\ \ (k=0,1,\cdots,n-1)$ 次,共 $n$ 种置换.

环上序号 $x$ 的珠子转 $k$ 次得到的序号为 $(x + k)$ .为出现循环,应使 $x+kt\equiv x\ \ (\mathrm{mod}\ n)$ ,则 $kt\equiv0\ \ (\mathrm{mod}\ n)$ .

$t$ 是不定方程 $k t + n r = 0$ 的最小正整数解.设 $d=\gcd(n,k)$ ,则 $t=\dfrac{n}{d}$ ,即转 $\dfrac{n}{d}$ 步出现循环,亦即每个循环上有 $\dfrac{n}{d}$ 个点.

因每个循环独立,则 $n$ 个点有 $\dfrac{n}{\dfrac{n}{d}}=d$ 个循环.有 $m$ 种颜色,则不动点数为 $m^d$ .

故旋转置换对答案的贡献为 $\displaystyle \sum_{k=0}^m m^{\gcd(n,k)}$ .

(2)翻转:

①奇数个珠子:对每个珠子,自己映射到自己,其余珠子左右对称配对,有 $m^{\frac{n+1}{2}}$ 个循环,则 $n$ 个珠子共 $n\cdot m^\frac{n+1}{2}$ 个循环.

②偶数个珠子:

i)若对称轴穿过相对的两珠子,则有 $m^{\frac{n}{2}+1}$ 个循环, $\dfrac{n}{2}$ 条对称轴共 $\dfrac{n}{2}\cdot m^{\frac{n}{2}+1}$ 个循环.

ii)若对称轴在相对的四珠子间,则有 $m^\frac{n}{2}$ 个循环, $\dfrac{n}{2}$ 条对称轴共 $\dfrac{n}{2}\cdot m^{\frac{n}{2}}$ 个循环.

代码

ll power(int a, int b) {  // a的b次方
	ll res = 1;
	while (b--) res *= a;
	return res;
}

int main() {
	int m, n;
	while (cin >> m >> n, m || n) {
		ll ans = 0;
		for (int i = 0; i < n; i++) ans += power(m, gcd(n, i));  // 旋转置换

		// 翻转置换
		if (n & 1) ans += n * power(m, (n + 1) / 2);
		else ans += n / 2 * (power(m, n / 2 + 1) + power(m, n / 2));
		ans /= 2 * n;
		cout << ans << endl;
	}
}

35.1.3 魔法手链

题意 ( $3\ \mathrm{s}$ )

给定编号 $1\sim m$ 的 $m$ 种不同颜色的珠子,每种颜色的珠子个数足够多.现需从中选出 $n$ 个珠子,串成一个环形手链.珠子间存在 $k$ 对排斥关系,互相排斥的两种颜色的珠子不能相邻(同种颜色的珠子也可能排斥).问能制成多少种本质不同的手链,答案对 $9973$ 取模.本质不同定义为:不能通过旋转(不考虑翻转)与其他方案对应位置上的珠子颜色相同.

有 $t\ \ (1\leq t\leq 10)$ 组测试数据.每组测试数据第一行输入整数 $n,m,k\ \ \left(1\leq n\leq 1\mathrm{e}9,\gcd(n,9973)=1,1\leq m\leq 10,0\leq k\leq\dfrac{m(m+1)}{2}\right)$ .接下来 $k$ 行每行输入两个整数 $a,b\ \ (1\leq a,b\leq m)$ ,表示颜色 $a$ 的珠子与颜色 $b$ 的珠子不能相邻.

对每组测试数据,输出能制成多少种本质不同的手链,答案对 $9973$ 取模.

思路

旋转置换:转 $k\ \ (k=0,1,\cdots,n-1)$ 次,共 $n$ 种置换.

环上序号 $x$ 的珠子转 $k$ 步得到的序号为 $(x + k)$ .为出现循环,应使 $x+kt\equiv x\ \ (\mathrm{mod}\ n)$ ,则 $kt\equiv0\ \ (\mathrm{mod}\ n)$ .

$t$ 是不定方程 $k t + n r = 0$ 的最小正整数解.设 $d=\gcd(n,k)$ ,则 $t=\dfrac{n}{d}$ ,即转 $\dfrac{n}{d}$ 步出现循环,亦即每个循环上有 $\dfrac{n}{d}$ 个点.

因每个循环独立,则 $n$ 个点有 $\dfrac{n}{\dfrac{n}{d}}=d$ 个循环.

每个循环的 $\dfrac{n}{d}$ 个点按编号顺序放在圆上,则它们两两间隔为 $d$ .

[证] 环上序号 $x$ 的珠子能转到的序号为 $(x+k),(x+2k),(x+3k)\cdots$ .

因 $d=\gcd(n,k)$ ,则它们两两间的距离是 $d$ 的倍数,它们与 $x$ 号珠子相距 $0,d,2d,3d\cdots$ .

令 $n'=\dfrac{n}{d},k'=\dfrac{k}{d}$ ,考虑在 $n^{'}$ 个点的环上转 $k^{'}$ 次,则能转到的珠子与 $x$ 号珠子相距 $0,1,2,3\cdots$ ,

这表明:每次转 $k^{'}$ 补能遍历圆上的所有点,即 $0k',1k',2k',\cdots,(n-1)k'$ 是 $0\sim n$ 的一个简化剩余系.

第一个循环的起点为 $x$ ,第二个循环的起点为 $x + 1$ , $\cdots$ ,第 $d$ 个循环的起点为 $x + d - 1$ .

若是不动点,则 $x$ 与 $x+d,x+2d,\cdots$ 的颜色相同, $x + 1$ 与 $x+1+d,x+1+2d,\cdots$ 的颜色相同, $\cdots$ ,

这表明: $[x, x + d - 1]$ 的颜色在环上循环,即环可分为 $d$ 段,每一段上的珠子颜色对应相同.

显然只需考虑 $[x, x + d - 1]$ 这一段的颜色,且与 $x + d - 1$ 相邻的点即 $x$ ,则可将 $[x, x + d - 1]$ 视为一个小环,只需求该小环的染色方案数.

设小环上珠子的编号为 $1\sim d$ ,可将 $d$ 视为 $0$ 号.枚举 $d$ 染的颜色, $d p [i] [j]$ 表示已染完前 $i$ 个珠子,且第 $(i + 1)$ 个珠子染 $j$ 色的方案数.若 $d$ 染 $1$ 色,则初始化 $dp[0][1]=1,dp[0][j]=0\ \ (j\neq 1)$ ;若 $d$ 染 $2$ 色,则初始化 $dp[0][2]=1,dp[0][j]=0\ \ (j\neq 2),\cdots$ .转移时若当前珠子染 $j$ ,上一个珠子染 $k$ ,则 $d p [i] [j] + = d p [i - 1] [k]$ .最终对 $d p [d] [1]$ 统计答案.

$n$ 最大 $1\mathrm{e}9$ ,直接暴力转移会超时.考虑矩阵乘法优化.令 $F(i)=[f(i,1)\ f(i,2)\ \cdots\ f(i,m)]$ ,则 $F(i)=F(i-1)\cdot M$ ,其中转移矩阵 $M$ 中若 $k$ 和 $j$ 相邻则为 $1$ ,否则为 $0$ .最终 $F(d)=F(0)\cdot M^d$ ,用快速幂优化,时间复杂度 $O(m^3\log n)$ . $m$ 种颜色可同时转移,即 $M$ 左乘一个 $m\times m$ 的单位矩阵,最后对矩阵的对角线求和即可.

$k$ 不能直接枚举 $0\sim (n-1)$ ,注意到方案数只取决于 $d=\gcd(n,k)$ ,则可将 $k$ 按 $d$ 分类,只需枚举 $n$ 的不超过 $1600$ 个约数即可.此时每个 $\gcd(n,k)=d$ 的 $k$ 的方案数都可用 $F(d)=F(0)\cdot M^d$ 求出.考察有多少个 $k\ s.t.\ \gcd(n,k)=d$ ,即有多少个 $\dfrac{k}{d}\ s.t.\ \gcd\left(\dfrac{n}{d},\dfrac{k}{d}\right)=1$ ,显然有 $\varphi\left(\dfrac{n}{d}\right)$ 个,暴力求时间复杂度 $O\left(\sqrt{\dfrac{n}{d}}\right)$ .

$9973$ 是素数,可直接快速幂求 $n$ 的逆元.

代码

const int MAXN = 11;  // 最大颜色数
const int MOD = 9973;
int n, m, k;  // 珠数、颜色数、限制数
struct Matrix {
	int arr[MAXN][MAXN];

	Matrix() { memset(arr, 0, so(arr)); }

	friend Matrix operator*(const Matrix& A, const Matrix& B) {
		Matrix res;
		for (int i = 1; i <= m; i++)
			for (int j = 1; j <= m; j++)
				for (int k = 1; k <= m; k++)
					res.arr[i][j] = (res.arr[i][j] + A.arr[i][k] * B.arr[k][j]) % MOD;
		return res;
	}
};

int qpow(Matrix A, int b) {
	Matrix C;
	for (int i = 1; i <= m; i++) C.arr[i][i] = 1;  // 初始化为单位矩阵

	while (b) {
		if (b & 1) C = C * A;
		A = A * A;
		b >>= 1;
	}

	int res = 0;
	for (int i = 1; i <= m; i++) res = (res + C.arr[i][i]) % MOD;  // 对矩阵的对角线统计答案
	return res % MOD;
}

int phi(int x) {
	int res = x;
	for (int i = 2; i <= x / i; i++) {
		if (x % i == 0) {
			res = res / i * (i - 1);
			while (x % i == 0) x /= i;
		}
	}
	if (x > 1) res = res / x * (x - 1);
	return res % MOD;
}

int main() {
	CaseT{
		cin >> n >> m >> k;
		
		Matrix tr;  // 转移矩阵
		for (int i = 1; i <= m; i++)
			for (int j = 1; j <= m; j++) tr.arr[i][j] = 1;  // 初始时都能转移

		while (k--) {
			int a, b; cin >> a >> b;
			tr.arr[a][b] = tr.arr[b][a] = 0;  // 有限制的不能转移
		}
				 
		int ans = 0;
		for (int d = 1; d <= n / d; d++) {  // 枚举n的约数
			if (n % d == 0) {
				ans = (ans + qpow(tr, d) * phi(n / d)) % MOD;
				if (d != n / d) ans = (ans + qpow(tr, n / d) * phi(d)) % MOD;
			}
		}
		cout << (ll)ans * qpow(n, MOD - 2, MOD) % MOD << endl;
	}
}